format: text

text:
# USACH-FM-02/08 hep-th/0212117 Dolan-Grady Relations and Noncommutative Quasi-Exactly Solvable Systems

## 1 Introduction

The Dolan-Grady relations (DGR) arose 20 years ago in a very interesting paper , where Dolan and Grady considered a class of systems described by a Hamiltonian given as a linear combination of the two operators obeying some nonlinear relations. The relations, afterward called by their names, guarantee the existence of an infinite set of mutually commuting charges which includes the Hamiltonian. This naturally connects the construction to the integrable systems . A rich algebraic structure associated with the Dolan-Grady relations was revealed by Perk  and Davies  who showed that the two operators satisfying DGR generate recursively the infinite-dimensional Onsager algebra. The latter appeared in a seminal work  on exact solution of the planar Ising model. The set of commuting charges constructed by Dolan and Grady is linear in generators of the Onsager algebra.

Recently it was observed  that the Dolan-Grady relations arise as necessary conditions for the anomaly-free quantization of pseudo-classical systems with nonlinear holomorphic supersymmetry ($`n`$-HSUSY). In accordance with general algebraic construction of the $`n`$-HSUSY , the Hamiltonian depending on a natural parameter $`n`$ is a quadratic form in two mutually conjugate operators subjected to DGR, whereas the supercharges are (anti) holomorphic polynomials of the order $`n`$ in the same operators. It was shown that the $`n`$-HSUSY Hamiltonian belongs to an infinite set of mutually commuting even operators being quadratic in the Onsager algebra generators, and that these even integrals together with the odd supercharges form a polynomial superalgebra of the order $`n`$. Realization of this construction in various physical systems  revealed an intimate relation of the $`n`$-HSUSY with the quasi-exactly solvable (QES) systems .

Though the quadratic set of commuting integrals admits any real and even complex values for the parameter $`n`$ playing the role of a coupling constant, such an extension breaks the nonlinear supersymmetry. Nevertheless, in Ref.  it was observed that the quasi exact solvability of the systems associated with the quadratic set of integrals of motion roots not in the $`n`$-HSUSY construction itself, but rather in the Dolan-Grady relations. In particular, it was shown that the corresponding systems can be quasi-exactly solvable not only for integer values of the coupling constant $`n`$.

In the present paper, exploiting the Dolan-Grady relations we consider a class of $`U(1)`$ gauge invariant quantum mechanical systems on a two-dimensional noncommutative manifold, and demonstrate the quasi exact solvability of their spectral problem for any value of the coupling constant.

The layout of the paper is as follows. In Section 2 the basic information on the Dolan-Grady relations and Onsager algebra necessary for further analysis is presented in the context of integrable systems, $`n`$-HSUSY and low-dimensional matrix models. Realization of the Dolan-Grady relations in terms of the generalized deformed oscillator algebra is given in Section 3. In Section 4 we construct a $`U(1)`$ gauge invariant quantum mechanical system with noncommutative coordinates generating the generalized deformed oscillator algebra. In Section 5 we demonstrate that by a reduction to finite-difference equations, the spectral problem of the model on the fuzzy sphere and noncommutative hyperbolic plane can be partially solved for some values of the coupling constant. In Section 6 we develop an algebraic covariant approach to deal with the spectral problem, and show that really our noncommutative system is quasi-exactly solvable for any value of the coupling constant. The obtained results are summarized in Section 7. In Appendix we discuss a deformation of the $`𝔰𝔩(2,)`$ scheme associated with the finite-difference quasi-exactly solvable equations.

## 2 Dolan-Grady relations

### 2.1 Dolan-Grady relations and Onsager algebra

The Dolan-Grady relations  can be represented in the form

$`[Z,[Z,[Z,\overline{Z}]]]`$ $`=\omega ^2[Z,\overline{Z}],`$ $`[\overline{Z},[\overline{Z},[Z,\overline{Z}]]]`$ $`=\overline{\omega }^2[Z,\overline{Z}].`$ (2.1)
Here $`Z`$ and $`\overline{Z}`$ are some operators which we shall call generating elements while $`\omega `$ and $`\overline{\omega }`$ are some constants. The nonlinear relations (2.1) can be used to construct the infinite set of mutually commuting operators for some integrable systems. In the most elegant way they can be represented in terms of the Onsager algebra  spanned by the generators $`A_n`$ and $`G_n`$ having the commutation relations

$`[A_m,A_n]`$ $`=4G_{mn},`$ $`[G_m,A_n]`$ $`=2A_{n+m}2A_{nm},`$ $`[G_m,G_n]`$ $`=0,`$ (2.2)
where $`m,n`$. In fact, the operators $`A_0\overline{Z}`$ and $`A_1Z`$ are the generating elements of the Onsager algebra obeying the Dolan-Grady relations and all the other operators are defined recursively by the relations (2.2) . The values of the constants $`\omega `$ and $`\overline{\omega }`$ can be fixed by rescaling the generators. Then, the infinite set of the mutually commuting operators constructed by Dolan and Grady  is represented as

$$2J_m=A_m+A_m+\lambda \left(A_{m+1}+A_{1m}\right),$$
(2.3)
where $`\lambda `$. The commutativity of the operators $`J_m`$ follows directly from the algebra (2.2). Now, treating the operator $`J_0=A_0+\lambda A_1`$ as a Hamiltonian, one obtains, generally, an integrable system with infinite number of conserved commuting charges. The charges of the form (2.3) appear in various models, such as the Baxter eight-vertex , the two-dimensional Ising  and $`Z_N`$ spin models . The power of the algebraic formulation (2.2), (2.3) is that it does not refer to the number of dimensions or to the nature of the space-time manifold which can be lattice, continuum or loop space.

### 2.2 Structure of $`n`$-HSUSY

In the systems associated with the linear set of integrals (2.3) the operators $`Z`$ ($`A_1`$) and $`\overline{Z}`$ ($`A_0`$) are supposed to be Hermitian. In Ref.  it was shown that the set (2.3) is not a unique set of mutually commuting operators which can be constructed in terms of generators of the Onsager algebra. There exists a set of commuting operators which are quadratic in the generators of the algebra (6.2). Such commuting operators naturally appear in the scheme of the $`n`$-HSUSY , which we discuss shortly below. Note, that for the systems associated with the quadratic set of commuting operators, the $`Z`$, $`\overline{Z}`$ and the $`\omega `$, $`\overline{\omega }`$ are usually supposed to be mutually conjugate operators and complex constants, respectively.

The Hamiltonian of the system associated with the quadratic set of commuting charges (see Eq. (6.2) for their explicit form) reads as

$$H_\lambda =\frac{1}{4}\{\overline{Z},Z\}\frac{\lambda }{4}[Z,\overline{Z}],$$
(2.4)
where a real parameter $`\lambda `$ serves as a coupling constant, while mutually conjugate operators $`Z`$ and $`\overline{Z}`$ are supposed to obey the nonlinear Dolan-Grady relations (2.1)111In principle, one can also deal with Hermitian operators $`Z`$ and $`\overline{Z}`$. Then the Hamiltonian (and the whole set of the integrals (6.2)) is Hermitian if the coupling constant is pure imaginary . However, here we will not discuss such systems..

When the coupling constant is integer, $`\lambda =\pm n`$, $`n_+`$, the system (2.4) can be extended to the supersymmetric system with the Hamiltonian $`_n=H_n\sigma _{}\sigma _++H_n\sigma _+\sigma _{}`$, $`\sigma _\pm =\frac{1}{2}(\sigma _1\pm i\sigma _2)`$. In this case the Dolan-Grady relations (2.1) guarantee the existence of the odd integrals of motion (supercharges)

$`Q`$ $`=\sigma _+{\displaystyle \underset{k=0}{\overset{n1}{}}}\left(Z{\displaystyle \frac{\omega }{2}}(2kn+1)\right),`$ $`\overline{Q}`$ $`=\sigma _{}{\displaystyle \underset{k=0}{\overset{n1}{}}}\left(\overline{Z}{\displaystyle \frac{\overline{\omega }}{2}}(2kn+1)\right).`$ (2.5)
The supercharges generate a nonlinear superalgebra of the order $`n`$. Its structure and further details of the algebraic construction underlying the nonlinear holomorphic supersymmetry can be found in Ref. .

The advantage of the supersymmetric formulation is that the zero modes of the supercharge $`Q`$ (or $`\overline{Q}`$) form an invariant subspace with respect to the action of the Hamiltonian $`H_n`$ (or $`H_n`$). Furthermore, as a rule the zero modes are easier calculable than the eigenfunctions of the Hamiltonian due to the factorized structure of the supercharges. When the number of the zero modes is finite, then the Hamiltonian is quasi-exactly solvable operator. The information on the functional form of the zero modes is very helpful for solving the corresponding spectral problem.

Although the system (2.4) formally admits an infinite number of integrals, nevertheless, for any given realization of the generating elements in a system with finite number of degrees of freedom this infinite set is reduced to a finite set of independent integrals. Such systems being not integrable in general case, however, reveal a profound connection with quasi-exactly solvable systems . We shall return to this aspect below.

### 2.3 Matrix models

Here we show that the Dolan-Grady relations can be revealed in some low dimensional matrix models.

So, let us consider a matrix model given by the action

$$S=𝖳𝗋\left(\frac{1}{2}[X_\mu ,X_\nu ][X^\mu ,X^\nu ]+R_{\mu \nu }X^\mu X^\nu \right),$$
(2.6)
where $`X_\mu `$ are Hermitian operators on a Hilbert space and $`R_{\mu \nu }`$ are $`c`$-numbers. The last term can be treated as describing the interaction of the vector field $`X_\mu `$ with an external symmetric field $`R_{\mu \nu }`$. One can also think that its diagonal part contains the mass term of the vector field. The matrix models of the form (2.6) are widely used in the context of string theories. In the case of positively definite metric the action (2.6) corresponds to the potential part of the BFSS model  and defines static (vacuum) configurations of the theory. On the other hand, for the Lorentz metric the action can be interpreted as an effective theory of the reduced model of Yang-Mills theory  with the specific interaction. Actually, action (2.6) represents a slightly modified massive matrix model .

The action (2.6) leads to the equations of motion

$$[X^\nu ,[X_\mu ,X_\nu ]]+R_{\mu \nu }X^\nu =0.$$
(2.7)
In the two-dimensional case Eq. (2.7) can be represented in the form

$`[X_{},[X_{},X_+]]+R_+X_{}+R_{}X_+`$ $`=0,`$
$`[X_+,[X_+,X_{}]]+R_+X_++R_{++}X_{}`$ $`=0,`$
where $`X_\pm =\frac{1}{\sqrt{2}}(X_1\pm X_2)`$ or $`X_\pm =\frac{1}{\sqrt{2}}(X_1\pm iX_2)`$ for the Lorentz or Euclidean metric, respectively. In the first case $`R_{++}`$ and $`R_{}`$ are real constants, while in the latter $`R_{}^{}=R_{++}`$. It is easy to see that any solution of Eq. (LABEL:para) also obeys the Dolan-Grady relations (2.1) with the identification:

$`X_{}`$ $`Z,`$ $`X_+`$ $`\overline{Z},`$ $`R_{}`$ $`\omega ^2,`$ $`R_{++}`$ $`\overline{\omega }^2.`$
Moreover, note that in the case of positively definite metric and when $`R_{++}=R_{}=0`$, the equations (LABEL:para) coincide with the basic commutations relations for the single-mode parafermion .

Now let us consider the following modification of the three-dimensional matrix model

$$S=𝖳𝗋\left(\frac{1}{2}[X_j,X_k]^2\frac{4}{3}i\alpha ϵ_{jkl}X_jX_kX_l+R_{jk}X_jX_k\right),$$
(2.9)
where $`\alpha `$ and $`R_{jk}`$ are real numbers and we imply positively definite metric. This model can be treated as the matrix model corresponding to the usual Yang-Mills theory with Chern-Simon term. As before, the quadratic term in action (2.9) represents a specific interaction with an external symmetric field. The equations of motion are

$$[[X_k,X_j],X_k]i\alpha ϵ_{jkl}[X_k,X_l]+R_{jk}X_k=0.$$
(2.10)
Imposing the condition

$$[X_0,X_\pm ]=\pm \beta X_\pm $$
(2.11)
with $`\beta `$ on the components $`X_0`$, $`X_\pm `$, the equation for $`j=0`$ is satisfied identically for $`\beta =\alpha `$, $`R_{0i}=0`$. The equations of motion for $`X_\pm `$ have the form (LABEL:para) with the substitution $`R_+R_++\alpha `$. Therefore, the Dolan-Grady relations arise in this model as well. Again, like in model (2.6), for $`R_{\pm \pm }=0`$ the equations of motion for the components $`X_\pm `$ take the form of basic commutation relations of the single-mode parafermion.

## 3 Generalized deformed oscillator

### 3.1 Main definitions

Another basic ingredient we shall use in what follows is the generalized deformed oscillator (GDO). Therefore, in this subsection we present a brief description of the construction (see also Refs. ).

The generalized deformed oscillator is defined by the algebra generated by the set of operators $`\{\text{11I},z,\overline{z},N\}`$ satisfying the relations

$`[N,z]`$ $`=\theta z,`$ $`[N,\overline{z}]`$ $`=\theta \overline{z},`$ $`\overline{z}z`$ $`=F(N),`$ $`z\overline{z}`$ $`=F(N+\theta ),`$ (3.1)
where $`N`$ is the number operator and the structure function $`F(.)`$ is an analytic function with the properties $`F(0)=0`$ and $`F(n\theta )>0`$, $`n=1,\mathrm{}`$. The structure function $`F(.)`$ is characteristic to the deformation scheme. Here we have introduced a positive real parameter $`\theta `$, which will serve later on as a measure of noncommutativity.

The operators $`z`$, $`\overline{z}`$ obey the commutation relation

$$[z,\overline{z}]=\theta \mathrm{\Delta }_+F(N).$$
(3.2)
By symbols $`\mathrm{\Delta }_\pm `$ we denote the forward and backward difference derivatives

$`\mathrm{\Delta }_+f(N)`$ $`={\displaystyle \frac{f(N+\theta )f(N)}{\theta }},`$ $`\mathrm{\Delta }_{}f(N)`$ $`={\displaystyle \frac{f(N)f(N\theta )}{\theta }}.`$ (3.3)
The generalized deformed oscillator algebra can be naturally represented on the Fock space of the eigenstates of the number operator $`N`$, $`N|n=n\theta |n`$, $`n|m=\delta _{nm}`$, $`n=1,\mathrm{}`$, with the vacuum state defined as $`z|0=0`$. The eigenstates are given by

$`|n`$ $`={\displaystyle \frac{1}{\sqrt{F(n\theta )!}}}\overline{z}^n|0,`$ $`F(n\theta )!`$ $`{\displaystyle \underset{k=1}{\overset{n}{}}}F(k\theta ).`$ (3.4)
The operators $`z`$ and $`\overline{z}`$ are the annihilation and creation operators of the oscillator algebra,

$`z|n`$ $`=\sqrt{F(n\theta )}|n1,`$ $`\overline{z}|n`$ $`=\sqrt{F((n+1)\theta )}|n+1.`$ (3.5)
In what follows we shall denote the Fock space defined by (3.4), (3.5) as $`_F`$. On this space the operators $`z`$ and $`\overline{z}`$ are mutually adjoint, $`\overline{z}=z^{}`$.

Note that when the characteristic function obeys the condition

$$F((p+1)\theta )=0$$
(3.6)
for some $`p`$, the creation-annihilation operators satisfy the nilpotency relations $`z^{p+1}=\overline{z}^{p+1}=0`$. This means that the corresponding representation is $`(p+1)`$-dimensional and we have a parafermionic type system of the order $`p`$.

### 3.2 The Dolan-Grady relations in terms of GDO

The simplest realization of the Dolan-Grady relations (2.1) can be obtained in terms of the generalized deformed oscillator:

$$Z=z,\overline{Z}=\overline{z}.$$
(3.7)
As a consequence of the natural grading of the GDO algebra defined by the operator $`N`$, for this representation of the generating elements $`Z`$ and $`\overline{Z}`$ it is possible to realize the Dolan-Grady relations only in its contracted form, $`\omega =0`$ ,

$`[Z,[Z,[Z,\overline{Z}]]]`$ $`=0,`$ $`[\overline{Z},[\overline{Z},[Z,\overline{Z}]]]`$ $`=0.`$ (3.8)
On the other hand, these relations lead to the restriction on the structure function in the form of the finite-difference equation:

$$\mathrm{\Delta }_+^3F(N)=0.$$
(3.9)
Since the structure function has to vanish at zero, the general solution of the equation (3.9) is

$$F(N)=c_1N^2+c_0\theta N,$$
(3.10)
where ad interim $`c_0,c_1`$. Three different cases corresponding to the solution (3.10) should be distinguished:

1. The case $`c_1=0`$, $`c_00`$ corresponds to the usual oscillator given by the Heisenberg algebra. The nonlinear supersymmetry for such a system was discussed in Ref. .
2. For $`c_1>0`$ the polynomial (3.10) is a structure function when the parameters obey the inequality $`c_1+c_0>0`$. Then the generalized oscillator realizes the half-bounded infinite-dimensional discrete series of the unitary representations of the $`𝔰𝔩(2,R)`$ algebra. The generators are, correspondingly, $`L_0=\theta \mathrm{\Delta }_+F(N)`$, $`L_{}=z`$, $`L_+=\overline{z}`$, or $`L_0=\theta \mathrm{\Delta }_+F(N)`$, $`L_{}=\overline{z}`$, $`L_+=z`$ for the series $`𝒟_\alpha ^+`$ or $`𝒟_\alpha ^{}`$  with $`\alpha =(c_0+c_1)/(2c_1)`$.
3. In the case $`c_1<0`$, for the solution (3.10) to be a structure function one has to impose the additional condition (3.6). This leads to quantization of one of the parameters, $`c_0=(p+1)c_1`$, and results in the $`(p+1)`$-dimensional representation of the $`𝔰𝔲(2)`$ algebra with the generators $`L_{}=z`$, $`L_+=\overline{z}`$ and $`L_3=\theta \mathrm{\Delta }_+F(N)`$. This is the single-mode parafermion case.

## 4 Noncommutative quantum mechanics and DGR

Let us consider a 2D quantum mechanical system with noncommutative complex-like coordinates $`z`$, $`\overline{z}`$ obeying the GDO algebra (3.1) realized on the Hilbert space $`_F`$. The operator-valued functions $`\mathrm{\Psi }(\overline{z},z)`$ on $`_F`$ with a finite norm with respect to the scalar product

$$(\mathrm{\Phi },\mathrm{\Psi })\mathrm{𝖳𝗋}__F\left(\mathrm{\Phi }^{}(\overline{z},z)\mathrm{\Psi }(\overline{z},z)\right)$$
(4.1)
will be treated by us as the states of the noncommutative system. Here the symbols $`\mathrm{𝖳𝗋}__F`$ and  denote the usual trace and Hermitian conjugation operations on the space $`_F`$ (3.4). The operator-valued functions with the finite norm themselves form a Hilbert space, which we shall denote as $`\widehat{}_F`$.

The corresponding derivative operators on the space $`\widehat{}_F`$ can be defined as

$$\theta ^1𝖺𝖽\overline{z},\overline{}\theta ^1𝖺𝖽z.$$
(4.2)
These operators possess the usual conjugation property $`^{}=\overline{}`$ with respect to the scalar product (4.1). It is worth noting that the derivative operators do not commute if the structure function $`F(N)`$ is different from a linear one. Indeed, $`[,\overline{}]=\theta ^1𝖺𝖽\mathrm{\Delta }_+F(N)`$.

Let us introduce a $`U(1)`$ gauge interaction associated with an external magnetic field. The noncommutative system admits the following two types of gauge transformations:

$`\mathrm{\Psi }^{}(\overline{z},z)`$ $`=U(\overline{z},z)\mathrm{\Psi }(\overline{z},z),`$ $`\mathrm{\Psi }^{}{}_{}{}^{}(\overline{z},z)`$ $`=\mathrm{\Psi }^{}(\overline{z},z)U^{}(\overline{z},z),`$ (4.3)
or

$$\mathrm{\Psi }^{}(\overline{z},z)=U(\overline{z},z)\mathrm{\Psi }(\overline{z},z)U^{}(\overline{z},z)$$
(4.4)
with $`U(\overline{z},z)`$ being a unitary operator. The scalar product (4.1) is invariant with respect to the both transformations. In the commutative limit, $`\theta 0`$, the transformation (4.3) is reduced to the usual $`U(1)`$ gauge transformation while (4.4) becomes trivial.

We realize the generating elements from (2.4) as the gauge covariant derivatives

$`Z`$ $`=+\rho \left[A(\overline{z},z)\right],`$ $`\overline{Z}`$ $`=\overline{}+\rho \left[\overline{A}(\overline{z},z)\right],`$ (4.5)
where the operator $`\rho [𝒪]`$, $`𝒪\{A(\overline{z},z),\overline{A}(\overline{z},z)\}`$, acts on $`\mathrm{\Psi },\mathrm{\Psi }^{}`$ as

$$\rho [𝒪]\mathrm{\Psi }=𝒪\mathrm{\Psi },\rho [𝒪]\mathrm{\Psi }^{}=\mathrm{\Psi }^{}𝒪$$

manifesting the fact that the states $`\mathrm{\Psi }`$ and $`\mathrm{\Psi }^{}`$ carry the charges of the opposite sign. The operator-valued functions $`A`$ and $`\overline{A}=A^{}`$ are the components of the noncommutative vector gauge potential corresponding to the case of stationary magnetic field, and we assume that the gauge coupling constant (“electric charge”) and the mass parameter are equal to unit, $`e=m=1`$. Under the both gauge transformations (4.3) and (4.4) the components of the noncommutative vector potential are transformed in the usual way

$`A^{}`$ $`=UAU^{}+UU^{},`$ $`\overline{A}^{}`$ $`=U\overline{A}U^{}U\overline{}U^{}.`$
The commutator of the generating elements has the form

$``$ $`[Z,\overline{Z}]=\theta ^1𝖺𝖽\mathrm{\Delta }_+F(N)+B(\overline{z},z).`$ (4.6)
We will treat

$$B(\overline{z},z)\overline{A}+\overline{}A+[A,\overline{A}]$$
(4.7)
as an external quasi-magnetic field since in the commutative limit the commutator in (4.7) disappears and it turns into a true magnetic field. Nevertheless, in the noncommutative case the quantity (4.7) transforms as a connection under the gauge transformations. Actually, if one considers the operator $`N`$ as the third dependent coordinate, then the quantity (4.7) is the $`N`$-th component of the gauge connection.

Let us suppose that the magnetic field depends on $`N`$ only. This corresponds to an axially symmetric field in the commutative limit. The following choice of the gauge,

$$A(\overline{z},z)=\overline{z}f(N),\overline{A}(\overline{z},z)=f(N)z,$$
(4.8)
where $`f(.)`$ is a real function, guarantees such a dependence of the magnetic field. Indeed, in the gauge (4.8) the quasi-magnetic field (4.7) acquires the form

$$B(N)=\mathrm{\Delta }_+\left(F(N)f(N\theta )\left(2\theta f(N\theta )\right)\right).$$
(4.9)
Now let us analyze the restrictions which follow from the Dolan-Grady relations. In the representation (4.5) the l.h.s. of Eq. (2.1) can be rewritten as follows

$$l.h.s.(\text{2.1})=\theta ^1𝖺𝖽\left(\overline{z}^2\mathrm{\Delta }_+^3F(N)\right)+\overline{z}^2\left(1\theta f(N+\theta )\right)\left(1\theta f(N)\right)\mathrm{\Delta }_+^2B(N)$$

$$+\overline{z}^2\left(f(N)+f(N+\theta )(1\theta f(N))\right)\mathrm{\Delta }_+^3F(N).$$

The comparison of the first terms in (4.6) and (LABEL:cdgrZ) leads to the conclusion that the Dolan-Grady relations (2.1) with $`\omega 0`$ can be satisfied only when the structure function of the GDO algebra is linear. On the other hand, for $`\omega =0`$ the expression (LABEL:cdgrZ) vanishes only when the structure function obeys the equation (3.9). This leads to the following constraint on the quasi-magnetic field:

$$\mathrm{\Delta }_+^2B(N)=0.$$
(4.11)
In other words, the quasi-magnetic field has to be at most linear in $`N`$. The function $`f(N)`$ appeared in the gauge (4.8) can be found from (4.9) taking into account the constraint (4.11).

The classification of the GDO algebra with the quadratic structure function (3.10) was presented in the previous section. In accordance with it, the noncommutative coordinates $`z`$, $`\overline{z}`$ and the operator $`N`$ form the $`𝔰𝔲(2)`$, $`𝔰𝔩(2,)`$, or the Heisenberg algebras. Therefore, the noncommutative system given by the Hamiltonian (2.4) in the representation (4.5) can be treated, correspondingly, as the system on the fuzzy sphere $`S_\theta ^2`$, on the noncommutative hyperbolic plane $`H_\theta ^2`$, or on the noncommutative plane $`_\theta ^2`$.

Note that since the structure function is quadratic, the simplest solution for $`f(N)`$ is

$$f(N)=f=const.$$
(4.12)
Indeed, one can verify that in this case the quasi-magnetic field (4.7) is a linear function,

$$B(N)=f(2\theta f)(2c_1N+(c_0+c_1)\theta ),$$
(4.13)
and hence, it obeys the equation (4.11). Below we shall consider only the solution (4.12), (4.13). Later on we shall show that for the given solution the expression (4.13) coincides with a covariant definition of the magnetic field in the gauge (4.8).

From the explicit form of the magnetic field (4.13) it follows that the parameter $`f`$ has two special values $`f=0`$ and $`f=2\theta ^1`$. The first case is rather trivial. It corresponds to the model without gauge interaction (a free model (2.4) on the corresponding noncommutative space). In the second case the magnetic field (4.13) vanishes. Nevertheless, the gauge interaction still is present in the model since the corresponding connection (4.8) with $`f0`$ is not trivial. This resembles the famous Aharonov-Bohm effect but its origin roots here in the noncommutativity of the configuration space. As we shall see, there is also another critical value of the parameter, $`f=\theta ^1`$. For this value the quasi-magnetic field is nontrivial while the system is degenerate. Such a case, in turn, resembles the situation with existence of the critical value of the constant external magnetic field being coherent with the value of the parameter of noncommutativity (see, e.g. ).

In the next section we shall discuss the spectral problem for the noncommutative system (2.4) on the fuzzy sphere and noncommutative hyperbolic plane.

## 5 Reduction to finite-difference QES systems

Here we show that the noncommutative system given by the Hamiltonian (2.4) is quasi-exactly solvable for some values of the coupling constant. It is demonstrated by a reduction of the Hamiltonian to one-dimensional finite-difference QES operators.

First of all, one can notice that in the gauge (4.8) the operator $`𝖺𝖽N`$ commutes with the Hamiltonian. One can refer this property to the existence in the system of an “axial” symmetry. Therefore, the eigenstates of the Hamiltonian can be represented in the form

$$\mathrm{\Psi }_L(\overline{z},z)=\overline{z}^m\psi (N),\text{or}\text{ }\mathrm{\Psi }_R(\overline{z},z)=\psi (N)z^m,$$
(5.1)
where $`m_+`$. Evidently, these operator functions are eigenstates of the operator $`𝖺𝖽N`$. We shall call the functions $`\mathrm{\Psi }_L`$ and $`\mathrm{\Psi }_R`$ the left and right modes.

Now it is clear enough that the corresponding 2D spectral equation with the eigenstates of the form (5.1) can be reduced to a one-dimensional problem. But unlike the commutative case, here the corresponding 1D system is represented by a finite-difference equation. One notes that such a situation is typical for any noncommutative system with the “axial” symmetry.

The Hamiltonian (2.4) is associated with the nonlinear holomorphic supersymmetry for integer values of the coupling constant. The zero modes of the supercharge $`Q`$ form an invariant subspace of the Hamiltonian. We will use their functional form as an anzatz for investigating the spectral problem with arbitrary coupling constant.

Let us temporarily set $`\lambda =n`$. For the left mode wave functions represented in the factorized form

$$\mathrm{\Psi }_L(\overline{z},z)=\overline{z}^m\phi (N)\varphi (N),$$

with $`\varphi (N)`$ obeying the equation

$$(\mathrm{\Delta }_++f)\varphi (N)=0,$$
(5.2)
the zero mode equation $`Z^n\mathrm{\Psi }_L=0`$ is reduced to

$$(1\theta f)^n\mathrm{\Delta }_+^n\phi (N)=0.$$

For $`\theta f1`$ this equation gives the polynomial solution $`\phi (N)=P_{n1}(N)`$, where $`P_k(.)`$ denotes a polynomial of the $`k`$-th degree. For $`\theta f=1`$ there are no zero modes at all, i.e. there are no bounded (normalized) states in the system.

The existence of the zero modes means that the Hamiltonian is quasi-exactly solvable when the coupling constant $`\lambda `$ is a nonnegative integer. Now we investigate the question on existence of the zero modes for an arbitrary value of the coupling constant. For the left modes the Hamiltonian can be reduced, $`H_\lambda \mathrm{\Psi }_L(\overline{z},z)H_{\lambda ,m}^L\psi (N),`$ to the one-dimensional finite-difference operator

$`H_{\lambda ,m}^L`$ $`=(\theta f1)\left(F(N)\mathrm{\Delta }_+\mathrm{\Delta }_{}+(m+1)(N)\mathrm{\Delta }_+\right)+f^2F(N)`$
$`+{\displaystyle \frac{1}{2}}f\left((2m+\lambda +1)\theta f2(m+\lambda )\right)\mathrm{\Delta }_+F(N)+const`$
with the function $`(N)`$ defined by

$`(N)`$ $`=\mathrm{\Delta }_+F(N)+{\displaystyle \frac{m\theta }{2}}\mathrm{\Delta }_+^2F(N).`$ (5.4)
For $`\theta f=1`$ the reduced Hamiltonian (5) is a multiplicative operator. The generating elements of the form (4.5) are multiplicative operators as well, $`Z\mathrm{\Psi }\mathrm{\Psi }\overline{z}`$, $`\overline{Z}\mathrm{\Psi }\mathrm{\Psi }z`$. Actually, this representation of the Dolan-Grady relations becomes equivalent to the representation (3.7). However, the Hamiltonian (2.4) acting on the Hilbert space $`\widehat{}_F`$ has a trivial dynamics in this case. Therefore, later on we shall always assume $`\theta f1`$.

The leading term of the “potential” in (5) is proportional to the structure function $`F(N)`$. We represent the function $`\psi (N)`$ in the factorized form

$$\psi (N)=\phi (N)\varphi (N),$$
(5.5)
where the function $`\varphi (N)`$ is supposed to obey the equation

$$(\mathrm{\Delta }_++\gamma )\varphi (N)=0$$
(5.6)
with a real parameter $`\gamma `$. This equation is a slight modification of Eq. (5.2). The factorization (5.5) allows us to discard the leading term in the equation for $`\phi (N)`$ when the constant $`\gamma `$ obeys the algebraic equation

$$\frac{\gamma ^2}{1\theta \gamma }=\frac{f^2}{1\theta f},$$
(5.7)
that gives a possibility to look for the solutions $`\phi (N)`$ in a polynomial form.

The first evident solution of the equation (5.7)

$$\gamma =f$$
(5.8)
corresponds to the zero modes found above. Indeed, the reduction

$$H_{\lambda ,m}^L\psi (N)H_{\lambda ,m}^{}\phi (N)$$

results in the one-dimensional finite-difference operator

$$H_{\lambda ,m}^{}=F(N)\mathrm{\Delta }_+\mathrm{\Delta }_{}(m+1)(1\theta f)^2(N)\mathrm{\Delta }_++f(2\theta f)𝒯_\lambda +const.$$
(5.9)
Here, the function $`(N)`$ is defined by (5.4), while the operator $`𝒯_\lambda `$ is given by

$$𝒯_\lambda =F(N)\mathrm{\Delta }_+\lambda \mathrm{\Delta }_+F(N).$$
(5.10)
Let $`𝒫_k(x)`$, $`k_+`$, denotes the space of all the polynomials in $`x`$ up to the order $`k`$. On such spaces the finite-difference derivatives act similarly to the usual derivative, $`\mathrm{\Delta }_\pm :𝒫_k(N)𝒫_{k1}(N)`$. Using this fact and that the function $`F(N)`$ is quadratic while $`(N)`$ is linear, one can notice that the spaces $`𝒫_k(N)`$ are invariant subspaces of the operator (5.9) except for the third term. On the other hand, the space $`𝒫_n(N)`$ is an invariant subspace of the operator $`𝒯_\lambda `$ for any $`\lambda =n_+`$. As a result, for nonnegative integer values of the coupling constant the operator (5.9) is quasi-exactly solvable .

The second solution of the equation (5.7) is

$$\gamma =\frac{f}{\theta f1}.$$
(5.11)
As a consequence, the reduced Hamiltonian acting on $`\phi (N)`$ reads as

$`H_{\lambda ,m}^{\prime \prime }`$ $`=(1\theta f)^2(F(N)\mathrm{\Delta }_+\mathrm{\Delta }_{}+(m+1)(N)\mathrm{\Delta }_+)f(2\theta f)𝒯_{\lambda 2m}+const.`$ (5.12)
Therefore, this operator is quasi-exactly solvable if one imposes the condition $`\lambda +2m_{}`$ on the coupling constant.

Now we pass over to the discussion of the right modes. Unfortunately, in this case we failed to find the zero modes of the supercharge $`Q`$ in a closed form. Therefore, we shall investigate the conjectural partial algebraization of this part of the spectrum from viewpoint of the reduction of the Hamiltonian.

For the right modes (5.1) the 2D Hamiltonian is reduced to the finite-difference operator

$`H_{\lambda ,m}^R`$ $`=(\theta f1)\left(F(N)\mathrm{\Delta }_+\mathrm{\Delta }_{}+(m+1)(N)\mathrm{\Delta }_+\right)+f^2F(N)`$
$`+{\displaystyle \frac{1}{2}}f\left((\lambda +1)\theta f+2(m\lambda )\right)\mathrm{\Delta }_+F(N)+const.`$
This Hamiltonian can be obtained from (5) by means of the formal change $`\lambda \lambda 2m`$. Using this property, one can reproduce all the corresponding results for the right modes. The factorization (5.5) with the conditions (5.6) and (5.7) for the solution (5.8) leads to the finite-difference operator containing the term proportional to the operator $`𝒯_{\lambda 2m}`$. Hence, for $`\lambda 2m_+`$ the polynomial space $`𝒫_{\lambda 2m}(N)`$ is an invariant subspace of the obtained operator, i.e. it is QES operator. For the second solution (5.11) the reduced Hamiltonian contains $`𝒯_\lambda `$, and so, it is QES operator for non-positive integer values of the coupling constant.

All the reduced finite-difference operators, (5.9), (5.12) and those for the right modes, can be represented as QES operators of the form (LABEL:gQES) with the following coefficients:

1. The solution (5.8):

$$a_1=0,a_2=0,b_+^+=\frac{c_1}{\theta }(1\theta f)^2,b_0^+=c_0,b_{}^+=0,b_+^{}=\frac{c_1}{\theta },$$

$$b_0^{}=\left(c_0+2(m+1)c_1\right)(1\theta f)^2,b_{}^{}=\theta (m+1)\left(c_0+(m+1)c_1\right)(1\theta f)^2,$$

where $`𝖭=\lambda `$ or $`𝖭=\lambda 2m`$ for the left or right modes, respectively;

2. The solution (5.11):

$$a_1=0,a_2=0,b_+^+=\frac{c_1}{\theta },b_0^+=c_0(1\theta f)^2,b_{}^+=0,$$

$$b_+^{}=\frac{c_1}{\theta }(1\theta f)^2,b_0^{}=c_02(m+1)c_1,b_{}^{}=\theta (m+1)\left(c_0+(m+1)c_1\right),$$

where $`𝖭=\lambda 2m`$ or $`𝖭=\lambda `$ for the left or right modes, respectively.

It is worth noting that for $`\theta f=2`$ the both operators (5.9) and (5.12) are exactly solvable since the term proportional to $`𝒯_\lambda `$ vanishes.

Let us briefly discuss a normalizability condition for the found solutions. Since the function $`\phi (N)`$ is a polynomial, the normalizability has to be provided by the factor $`\varphi (N)`$. The solution of the equation (5.6) has the form $`\varphi (N)(1\theta \gamma )^N`$. Therefore, the norm of the corresponding state $`\mathrm{\Psi }\widehat{}_F`$ can be represented as

$`\mathrm{\Psi }_{\widehat{}_F}`$ $`=𝖳𝗋__F(\mathrm{\Psi }^{}\mathrm{\Psi })=C_𝒩{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}n^K(1\theta \gamma )^{2n}+\mathrm{}.`$ (5.14)
Here the dots denote less divergent terms in the sense that if the first term in (5.14) is convergent, then the others are convergent as well. The factor $`C_𝒩`$ is related to the normalization constant and the integer parameter $`K`$ is defined by the type of solution. The requirement of convergence of the series leads to the condition

$`|1\theta \gamma |<1.`$ (5.15)
Therefore, we have to impose the following restrictions on the values of the parameter $`f`$: $`\theta f(0,2)`$ for the solution (5.8), and $`\theta f(\mathrm{},0)(2,\mathrm{})`$ for the solution (5.11).

Of course, in the case of the fuzzy sphere the norm (5.14) is finite because of the finitness of the corresponding Hilbert space, and no additional condition appears for it.

## 6 Intrinsic algebra and quasi exact solvability

In this section we show that the noncommutative system under consideration really is quasi-exactly solvable for any value of the coupling constant. The demonstration is realized in terms of an intrinsic algebra.

### 6.1 Intrinsic algebra

In general, the operators $`Z_0Z`$ and $`\overline{Z}_0\overline{Z}`$ together with the contracted Dolan-Grady relations (3.8) recursively generate the infinite-dimensional contracted Onsager algebra :

$`[Z_k,\overline{Z}_l]`$ $`=B_{k+l+1},`$ $`[Z_k,B_l]`$ $`=Z_{k+l},`$ $`[B_k,\overline{Z}_l]`$ $`=\overline{Z}_{k+l},`$
$`[Z_k,Z_l]`$ $`=0,`$ $`[\overline{Z}_k,\overline{Z}_l]`$ $`=0,`$ $`[B_k,B_l]`$ $`=0,`$
where $`k,l_+`$ and $`B_00`$ is implied. This algebra can be extended discarding the last condition, i. e. by treating $`B_0`$ as a nontrivial generator. Then the extended algebra will be generated by the set $`Z`$, $`\overline{Z}`$ and $`B_0`$. The element $`B_0`$ splits up the contracted Onsager algebra into the three subalgebras of the grade $`+1`$, $`0`$ and $`1`$, given, correspondingly, by the span of $`\overline{Z}_m`$, $`B_m`$ and $`Z_m`$. In other words, the generator $`B_0`$ is a grading element which provides the extended algebra with the triangular decomposition. This grading is different from that considered in Ref. .

The algebra (LABEL:cOA) admits the infinite set of the commuting quadratic charges

$$𝖩_\lambda ^k=\frac{1}{2}\underset{p=1}{\overset{k}{}}\{\overline{Z}_{p1},Z_{kp}\}\frac{1}{2}\underset{p=1}{\overset{k1}{}}B_pB_{kp}\frac{\lambda }{2}B_k$$
(6.2)
with $`k_+`$, which contains the Hamiltonian (2.4), $`𝖩_\lambda ^1=2H_\lambda `$, and the grading operator, $`𝖩_\lambda ^0B_0`$.

In the gauge (4.8) the generating elements of the form (4.5) read as

$`Z`$ $`=\theta ^1𝖺𝖽\overline{z}+f\rho [\overline{z}],`$ $`\overline{Z}=`$ $`\theta ^1𝖺𝖽z+f\rho [z].`$
These operators generate the following algebra intrinsic to the noncommutative system:

$`[Z,\overline{Z}]`$ $`=2c_1\left(S+f\left(2\theta f\right)G\right),`$ $`[𝒟,\overline{𝒟}]`$ $`=2c_1\theta G,`$
$`[Z,\overline{𝒟}]`$ $`=2c_1\left(1\theta f\right)G,`$ $`[G,𝒟]`$ $`=\theta 𝒟,`$
$`[𝒟,\overline{Z}]`$ $`=2c_1\left(1\theta f\right)G,`$ $`[\overline{𝒟},G]`$ $`=\theta \overline{𝒟},`$
$`[Z,G]`$ $`=\left(1\theta f\right)𝒟,`$ $`[S,\overline{𝒟}]`$ $`=\overline{𝒟},`$ (6.3)
$`[G,\overline{Z}]`$ $`=\left(1\theta f\right)\overline{𝒟},`$ $`[𝒟,S]`$ $`=𝒟,`$
$`[Z,S]`$ $`=Z,`$ $`[S,\overline{Z}]`$ $`=\overline{Z},`$
$$[S,G]=[Z,𝒟]=[\overline{Z},\overline{𝒟}]=0,$$

where the operators $`\overline{𝒟}`$, $`𝒟`$, $`G`$ and $`S`$ in the chosen gauge are

$`𝒟`$ $`=\overline{z},`$ $`\overline{𝒟}`$ $`=z,`$ $`G`$ $`=N+\theta {\displaystyle \frac{c_0+c_1}{2c_1}},`$ $`S`$ $`=\theta ^1𝖺𝖽N.`$ (6.4)
The algebra has the following two Casimir operators:

$`C_1`$ $`=\{\overline{𝒟},𝒟\}2c_1G^2,`$
$`C_2`$ $`=\{\theta \overline{Z}+(1\theta f)\overline{𝒟},\theta Z+(1\theta f)𝒟\}2c_1\left(\theta S+G\right)^2.`$
The algebra (6.1) is $`𝔰𝔩(2,)𝔰𝔩(2,)`$ ($`c_1>0`$) or $`𝔰𝔲(2)𝔰𝔲(2)`$ $`(c_1<0)`$. In the chosen gauge it can be shown by representing the generators of the algebra in the form $`L_z,L_{\overline{z}},L_N`$ and $`R_z,R_{\overline{z}},R_N`$, where $`L_a`$ and $`R_a`$ denote the left and right multiplications. In the given representation the Casimir operators (6.1) are equal to the same number,

$`C_1`$ $`=C_2={\displaystyle \frac{\theta ^2}{2c_1}}\left(c_1^2c_0^2\right).`$ (6.6)
It is necessary to note that the values $`\theta f=0,\mathrm{\hspace{0.17em}1},\mathrm{\hspace{0.17em}2}`$ are special from the viewpoint of the algebra (6.1). Although the algebraic content of it remains the same, the algebra is not generated by the operators $`Z`$ and $`\overline{Z}`$, i. e. it is not intrinsic any more. In these cases the generating elements produce the $`𝔰𝔩(2,)`$ or $`𝔰𝔲(2)`$ only.

The operators $`𝒟`$, $`\overline{𝒟}`$ and $`G`$ can be treated as dependent covariant coordinates of a noncommutative membrane. Indeed, the equivalent set of coordinates

$`X_1`$ $`={\displaystyle \frac{𝒟+\overline{𝒟}}{\sqrt{8|c_1|}}},`$ $`X_2`$ $`=i{\displaystyle \frac{𝒟\overline{𝒟}}{\sqrt{8|c_1|}}},`$ $`X_3`$ $`=G`$
obeys the condition

$$X_1^2+X_2^2\epsilon X_3^2=\epsilon \theta ^2\alpha \left(1\alpha \right)=\epsilon r^2,$$
(6.7)
where $`\epsilon =𝗌𝗂𝗀𝗇c_1`$ and $`\alpha =(c_1+c_0)/(2c_1)(c_1>0)`$ or $`\alpha =p/2`$ ($`c_1<0`$), $`p_+`$. So, the operators $`X_i`$ form the set of noncommutative covariant coordinates of the fuzzy sphere $`S_\theta ^2`$ ($`c_1<0`$), or of the noncommutative hyperbolic plane $`H_\theta ^2`$ ($`c_1>0`$). The operators $`Z`$, $`\overline{Z}`$ and $`S`$, in their turn, can be interpreted as the covariant derivatives. Note that in the noncommutative system under consideration, unlike the constant magnetic field case , the covariant coordinates are not proportional to the covariant derivatives.

The intrinsic algebra gives us the possibility to define the magnetic field in a covariant way. In general, if the derivatives form a Lie algebra, $`[_i,_j]=f_{ij}{}_{}{}^{k}_{k}^{}`$, then the corresponding gauge covariant field strength can be defined as $`iF_{ij}=[𝒟_i,𝒟_j]f_{ij}{}_{}{}^{k}𝒟_{k}^{}`$ (e. g., see Ref. ), where $`𝒟_i`$ are gauge covariant derivatives. From this definition and the intrinsic algebra it follows that in the present case only the “transversal” component ($`i=3`$) of the magnetic field, $`B_i=\frac{1}{2}f_i{}_{}{}^{jk}F_{jk}^{}`$, does not vanish. In the gauge (4.8) this component exactly coincides with the definition (4.13). It is interesting to note that like in the commutative case , the magnetic field is a linear function of the “transversal” coordinate.

In terms of the generators of the algebra (6.1) the higher operators of the contracted Onsager algebra are

$`Z_m`$ $`=(2c_1)^m\left(Z+\theta ^1(1\theta f)(1(1\theta f)^{2m})𝒟\right),`$
$`B_m`$ $`=(2c_1)^m\left(S+\theta ^1(1(1\theta f)^{2m})G\right),`$
where $`\theta f0,1,2`$. Applying this form of the generators of the Onsager algebra to the infinite set of the commuting charges (6.2), one can obtain

$`𝖩_\lambda ^k={\displaystyle \frac{(2c_1)^{k1}}{\theta f(2\theta f)}}\left(2\left(1(1\theta f)^{2k}\right)H_\lambda c_1\left((1\theta f)^2(1\theta f)^{2k}\right)S(S\lambda )\right)+\mathrm{},`$
where dots denote a linear combination of the Casimir operators (6.1). This representation leads to the conclusion that the Hamiltonian (2.4) and the operator $`S`$ form a complete set of nontrivial integrals of the system under consideration.

### 6.2 The commutative limit

Our aim is to apply the found intrinsic algebra of the noncommutative system with the Hamiltonian ( 2.4) for investigation of its spectrum. The algebraic scheme we shall discuss is an analogue of the Fock space construction. It is based on the introduction of the “ground” vectors in the Hilbert space of the system (2.4). These vectors have special properties with respect to the action of the generators of the intrinsic algebra.

First of all, let us note that in the commutative limit $`\theta 0`$ the algebra (6.1) is converted into the intrinsic algebra of the system (2.4) living on the (pseudo)sphere . Indeed, with the rescaling

$`G`$ $`{\displaystyle \frac{G}{4c_1f}},`$ $`𝒟`$ $`{\displaystyle \frac{𝒟}{2f}},`$ $`\overline{𝒟}`$ $`{\displaystyle \frac{\overline{𝒟}}{2f}},`$ $`c_12\beta ,`$
where the parameter $`\beta `$ is related to the radius of the (pseudo)sphere, the algebra (6.1) changes into the intrinsic algebra appeared in Ref.  in the limit $`\theta 0`$.

Since the commutative case is much simpler than the noncommutative one, we shall start the discussion of the covariant algebraic construction from it.

We introduce a vector (“ground” state) of the Hilbert space which obeys the conditions

$`Z\mathrm{\Psi }^{(0)}`$ $`=𝒟\alpha (G)\mathrm{\Psi }^{(0)},`$ $`\overline{Z}\mathrm{\Psi }^{(0)}`$ $`=\overline{𝒟}\beta (G)\mathrm{\Psi }^{(0)},`$ $`S\mathrm{\Psi }^{(0)}`$ $`=s\mathrm{\Psi }^{(0)}.`$ (6.8)
Here $`\alpha (.)`$ and $`\beta (.)`$ are some functions, which will be defined later, and $`s`$. We start with the real eigenvalues of the operator $`S`$ since no quantization condition follows from the intrinsic algebra. Besides, we know  that the eigenvalues do are real in general, since the magnetic flux gives a contribution to them and for a noncompact surface it is not quantized.

From the commutation relations for the generating elements $`Z`$ and $`\overline{Z}`$ the following constraint on the functions $`\alpha (G)`$ and $`\beta (G)`$ appears:

$`{\displaystyle \frac{1}{2}}(G^2{\displaystyle \frac{c_1^2}{4}})\varrho (G)+G\varrho (G)=G4\beta ^2s,`$ (6.9)
where $`\varrho (G)=\alpha (G)+\beta (G)`$. The solution of the differential equation is

$`\varrho (G)=1{\displaystyle \frac{8\beta ^2s}{G+\frac{c_1}{2}}}{\displaystyle \frac{4\beta ^2c_1m}{G^2\frac{c_1^2}{4}}}.`$ (6.10)
The parameter $`m`$ is the integration constant. One should regard this constant as an integer since it corresponds to the quantum number of the planar angular momentum of the system .

We will look for the eigenstates of the Hamiltonian in the form

$`\mathrm{\Psi }=\sigma (G)\mathrm{\Psi }^{(0)},`$ (6.11)
where we suppose that $`\sigma (G)`$ is a polynomial. Then the Hamiltonian (2.4) is reduced to the second order differential operator

$`H_\lambda ^{}`$ $`=4\beta ^2\left(G^2{\displaystyle \frac{c_1^2}{4}}\right){\displaystyle \frac{d^2}{dG^2}}+\left(\left(G^2{\displaystyle \frac{c_1^2}{4}}\right)\left(\alpha (G)\beta (G)\right)8\beta ^2G\right){\displaystyle \frac{d}{dG}}`$
$`+G(2f+\lambda 1)+(G^2{\displaystyle \frac{c_1^2}{4}})(\alpha (G)+{\displaystyle \frac{\alpha (G)\beta (G)}{4\beta ^2}})+const`$
acting on the function $`\sigma (G)`$.

Let us fix the residuary arbitrariness in the definition of the functions $`\alpha (G)`$ and $`\beta (G)`$ requiring that the operator (LABEL:comH) is quasi-exactly solvable. Therefore, it should be represented in terms of the generators of the $`𝔰𝔩(2,)`$ algebra

$`J_+`$ $`=G^2{\displaystyle \frac{d}{dG}}(𝖭1)G,`$ $`J_0`$ $`=G{\displaystyle \frac{d}{dG}}{\displaystyle \frac{𝖭1}{2}},`$ $`J_{}`$ $`={\displaystyle \frac{d}{dG}},`$ (6.13)
where if the real parameter $`𝖭`$ is a nonnegative integer, it defines dimensionality of the corresponding finite-dimensional representation. Moreover, the operator (LABEL:comH) is a quadratic polynomial in the generators (6.13) since it is a second order differential operator .

Eventually one can conclude that the operator (LABEL:comH) can be represented in the form

$`T`$ $`=c_1^2\beta ^2J_{}^24\beta ^2J_0^2+a_+J_++a_0J_0+a_{}J_{}+const,`$ (6.14)
where the parameters $`a_\pm `$, $`a_0`$ and $`𝖭`$ have the following eight solutions:

$`𝖭`$ $`=ϵ\lambda ,`$ $`a_+`$ $`=ϵ,`$ $`a_0`$ $`=4ϵ\beta ^2\left(2s\lambda \right),`$ $`a_{}`$ $`={\displaystyle \frac{ϵ}{4}}c_1\left(c_1+16\beta ^2(ms)\right),`$
$`𝖭`$ $`=ϵ(\lambda m),`$ $`a_+`$ $`=ϵ,`$ $`a_0`$ $`=4ϵ\beta ^2(2sm\lambda ),`$ $`a_{}`$ $`={\displaystyle \frac{ϵ}{4}}c_1(c_116\beta ^2s),`$
$`𝖭`$ $`=ϵ(\lambda 2s),`$ $`a_+`$ $`=ϵ,`$ $`a_0`$ $`=4ϵ\lambda \beta ^2,`$ $`a_{}`$ $`={\displaystyle \frac{ϵ}{4}}c_1(c_116\beta ^2(ms)),`$
$`𝖭`$ $`=ϵ(\lambda +m2s),`$ $`a_+`$ $`=ϵ,`$ $`a_0`$ $`=4ϵ(m\lambda )\beta ^2,`$ $`a_{}`$ $`={\displaystyle \frac{ϵ}{4}}c_1(c_1+16\beta ^2s),`$ (6.15)
where $`ϵ=\pm 1`$. So, if the parameters of the system take such values that $`𝖭_+`$ for any of the solutions, then the system is quasi-exactly solvable. Any QES part of the spectrum is a superstructure over the corresponding “ground” state. Therefore, one can interpret these states as some QES “vacua” of the system. Here and in what follows we will understand under the QES “vacuum” a state over which one can construct a certain number of states of the QES part of the spectrum following an oscillator-like procedure based on an intrinsic algebra.

The number of the $`𝔰𝔩(2,)`$ representations (6.14) coincides with the number of the solutions found in Ref. . At the same time, it is worth noting that the method based on the intrinsic algebra is completely covariant since all the generators of the algebra are covariant quantities.

We would like also to remind that this commutative system possesses a discrete symmetry $`_2_2_2`$, which intertwines the solutions (6.15) between themselves. The detailed discussion of the symmetry can be found in Ref. .

### 6.3 The noncommutative case

To investigate the spectral problem for the noncommutative system with the Hamiltonian (2.4) in terms of the intrinsic algebra (6.1), we introduce a “ground” state obeying the conditions (6.8) but with the only change

$$S\mathrm{\Psi }^{(0)}=m\mathrm{\Psi }^{(0)},m.$$
(6.16)
One can treat this number as an “azimuthal” quantum number since it corresponds to the existence of “axial” symmetry in the system with operator $`S`$ as its generator. We assume this in accordance with the representation (5.1) of the eigenstates. Therefore, in the noncommutative case the conserved quantity $`S`$ is an analogue of the angular momentum operator. We would like to note that although we use formally the same conditions (6.8), one should remember that the operators $`𝒟`$, $`\overline{𝒟}`$ and $`G`$ do not commute and the ordering of the operators is important.

From the commutation relations of the generating elements (6.1) and the definition of the “ground” state (6.8) it follows that the functions $`\alpha (G)`$ and $`\beta (G)`$ are not independent but obey the constraint

$$\begin{array}{c}2c_1G\left(\left(1\theta f\right)\left(\alpha (G)+\beta (G)\right)+\theta \alpha (G)\beta (G)\right)+\overline{𝒟}𝒟\left(1\theta f+\theta \alpha (G)\right)\mathrm{\Delta }_+\beta (G)\hfill \\ \hfill +𝒟\overline{𝒟}\left(1\theta f+\theta \beta (G)\right)\mathrm{\Delta }_{}\alpha (G)=2c_1\left(m+f\left(2\theta f\right)G\right).\end{array}$$
(6.17)
This is the noncommutative analogue of the restriction (6.9).

For investigation of the spectral problem it is convenient to represent the Hamiltonian as the linear in $`Z`$, $`\overline{Z}`$ form:

$`H_\lambda `$ $`=2\theta ^1\left(1\theta f\right)\left(𝒟\overline{Z}+\overline{𝒟}Z\right)+2c_1\theta ^1f\left(2\theta f\right)G(G\lambda \theta )`$
$`+4c_1\theta ^1GS+2c_1S(S\lambda )\theta ^2C_1\left(1\theta f\right)^2+\theta ^2C_2.`$
Using the definition of the “ground” state (6.8), the spectral problem for the Hamiltonian in this form is reduced to that of the finite-difference operator

$`𝖧_\lambda `$ $`=`$ $`2(\theta f1)(𝒟\overline{𝒟}(1\theta f+\theta \beta (G))\mathrm{\Delta }_+\mathrm{\Delta }_{}+(\overline{𝒟}𝒟\alpha (G)𝒟\overline{𝒟}\beta (G)`$ (6.19)
$`+2c_1(1\theta f)G)\mathrm{\Delta }_++\theta ^1(\overline{𝒟}𝒟\alpha (G)+𝒟\overline{𝒟}\beta (G)))`$
$`+2c_1f\theta ^1\left(2\theta f\right)G(G\lambda \theta )+4c_1m\theta ^1G+const`$
acting on the function $`\sigma (G)`$. Here the operators $`\mathrm{\Delta }_\pm `$ are defined as the difference derivatives (3.3) but with $`N`$ changed for $`G`$. The operator (6.19) is a finite-difference operator in the variable $`G`$ since $`\overline{𝒟}𝒟`$ and $`𝒟\overline{𝒟}`$ are the quadratic polynomials in $`G`$,

$`\overline{𝒟}𝒟`$ $`=c_1\left(G+{\displaystyle \frac{\theta }{2}}\right)^2+{\displaystyle \frac{\theta ^2c_0}{4c_1}},`$ $`𝒟\overline{𝒟}`$ $`=c_1\left(G{\displaystyle \frac{\theta }{2}}\right)^2+{\displaystyle \frac{\theta ^2c_0}{4c_1}},`$
as this follows from the intrinsic algebra (6.1) and the form of the Casimir operators (6.1), (6.6). The evident advantage of the representation (LABEL:linH) is that the finite-difference operator (6.19) is linear in the functions $`\alpha (G)`$ and $`\beta (G)`$.

Representing the functions $`\alpha (G)`$ and $`\beta (G)`$ in the form

$$\alpha (G)=\frac{\stackrel{~}{\alpha }(G)}{\overline{𝒟}𝒟},\beta (G)=\frac{\stackrel{~}{\beta }(G)}{𝒟\overline{𝒟}},$$

one can see that if the functions $`\stackrel{~}{\alpha }(G)`$ and $`\stackrel{~}{\beta }(G)`$ are polynomials, then the operator (6.19) is a finite-difference operator with polynomial coefficient functions and vice versa.

Let us fix the form of the functions $`\stackrel{~}{\alpha }(G)`$ and $`\stackrel{~}{\beta }(G)`$ by demanding the Hamiltonian (LABEL:linH) to be a QES operator of the form (LABEL:gQES) for some values of the parameters. Moreover, we require that the functions to obey the constraint (LABEL:linH). In comparison with the commutative case the noncommutative one is drastically different. There are four solutions for the functions $`\stackrel{~}{\alpha }(G)`$ and $`\stackrel{~}{\beta }(G)`$ providing the quasi exact solvability of the Hamiltonian (LABEL:linH). The explicit form of the function is not important. The corresponding solutions for the coefficients in the representation (LABEL:gQES) are

$$a_{1,2}=0,b_+^\pm =\frac{2c_1}{\theta }(1\theta f)^{1\pm \mu },b_0^\pm =2c_1\left(2j+\mu \left((1\pm 1)m\lambda \right)\right)(1\theta f)^{1\mu },$$

$$b_{}^\pm =\pm \frac{\theta (c_1ϵc_0)}{2c_1}\left(ϵc_0\left(4j+2\mu \left((1\pm \mu )m\lambda \right)+1\right)c_1\right)(1\theta f)^{1\mu },$$
(6.20)
where $`ϵ=\pm 1`$ and $`\mu =\pm 1`$ parametrize the set of solutions. The main difference of these solutions from those of the commutative case (6.15) is the absence of restrictions on the parameter $`𝖭=2j+1`$. Therefore, the operator (6.19) is quasi-exactly solvable since we can always choose the parameter $`𝖭`$ to be natural.

The functions $`\stackrel{~}{\alpha }(G)`$ and $`\stackrel{~}{\beta }(G)`$ corresponding to the solution (6.20) depend on the parameter $`j`$. Therefore, one can consider that it labels the “ground” states as well,

$$\mathrm{\Psi }^{(0)}=\mathrm{\Psi }_{m,j}^{(0)}.$$

Here we recovered the dependence of the state on the quantum number $`m`$ (6.16). Hence, unlike the commutative case, the noncommutative system admits the infinite number of the “ground” states for every value of the quantum number $`m`$. These states serve as QES “vacua” of the model.

Like the commutative case, the noncommutative system has a discrete symmetry. Indeed, evidently that the substitution

$$c_0c_0$$
(6.21)
converts the solutions (6.20) with $`ϵ=1`$ to those with $`ϵ=1`$ and vice versa. This transformation is similar to the transformations of the commutative case since it converts one QES operator (6.19) into another of the same type. Besides, in this case there is another discrete transformation

$`\theta `$ $`\theta ,`$ $`f`$ $`f,`$ $`\lambda `$ $`\lambda ,`$ $`m`$ $`m.`$ (6.22)
This mapping induces the following change of the parameters (6.20) and the operators (A.5):

$$b_\alpha ^\pm (\mu )b_\alpha ^{}(\mu ),J_\alpha ^\pm J_\alpha ^{},$$

where $`\alpha =0,\pm `$. The transformation of the generators is provoked by the change of the sign of the parameter $`\theta `$. This follows from the definitions of the forward and backward derivatives (3.3) and of the generators (A.5) themselves. In this way, the transformation (6.22) intertwines the QES operators corresponding to the solutions (6.20) with $`\mu =1`$ and $`\mu =1`$, respectively. The whole group of the discrete symmetry consists of four elements including the identity mapping. The transformations (6.21) and (6.22) are associated with the generators of a finite generated Abelian group . This discrete group is $`_2_2`$ since all its elements are involutive. So, like in the commutative case, this discrete symmetry is a kind of duality.

## 7 Discussion and outlook

Let us summarize briefly the obtained results and discuss some open problems that deserve further attention.

The Dolan-Grady relations provide a rich algebraic structure, which arises in numerous branches of mathematical physics. These relations originally appeared in the framework of various integrable models . Later it was observed that the Dolan-Grady relations are necessary conditions of the anomaly-free quantization  of a pseudo-classical 1D systems with the nonlinear holomorphic supersymmetry . This observation revealed the intimate connection of the relations with the QES systems. The algebraic origin of the construction allowed us to apply it to the two-dimensional systems . In the present paper we demonstrated that

* The Dolan-Grady relations arise in low-dimensional matrix models.

The matrix models naturally lead to appearance of the noncommutative geometry . Using this observation as a hint, we applied the algebraic scheme described in Ref.  to the two-dimensional systems on a noncommutative background subjected to an external $`U(1)`$ gauge interaction. We started by postulating that the coordinates form the generalized deformed oscillator algebra . The Hamiltonian of the system was constructed in terms of generating elements realized as gauge covariant derivatives on the noncommutative manifold. We showed that the generating elements obey the Dolan-Grady relations only when

* The noncommutative coordinates generate $`𝔰𝔲(2)`$, $`𝔰𝔩(2,)`$ or the Heisenberg algebras.

This allowed us to identify the resulting noncommutative manifold, correspondingly, with the fuzzy sphere $`S_\theta ^2`$, the noncommutative hyperbolic plane $`H_\theta ^2`$ or the noncommutative plane $`_\theta ^2`$. The cases of $`𝔰𝔲(2)`$ and $`𝔰𝔩(2,)`$ were discussed in detail.

We showed that the parameter $`f`$ defining the intensity of the quasi-magnetic field has the three critical values: $`f=0`$, $`f=\theta ^1`$ and $`f=2\theta ^1`$. The first case is trivial since it corresponds to the absence of the gauge interaction ($`B=0`$). The second case is analogous to the situation for the particle on the fuzzy sphere and the noncommutative hyperbolic plane in the presence of the constant magnetic field $`B=\theta ^1`$  (see also Ref.  for the case of the noncommutative plane). For this special value the systems are degenerate. In the third case, $`f=2\theta ^1`$, the quasi-magnetic field vanishes while the gauge interaction still is present in the model since the corresponding connection is not trivial, $`f0`$ (see Eqs. (4.8), (4.12)). This resembles the famous Aharonov-Bohm effect, but here its origin roots in the noncommutativity of the configuration space. Note that there is no analogue of such a phenomenon in the constant magnetic field case.

By the reduction of the Hamiltonian to one-dimensional finite-difference operators, it was showed that the constructed noncommutative systems are quasi-exactly solvable for some values of the coupling constant. In the commutative limit the systems on the fuzzy sphere and noncommutative hyperbolic plane lead to QES systems on the corresponding Riemann surfaces .

The algebraic content of the given model was revealed by the construction of the intrinsic algebra. Using the notion of the intrinsic algebra, in Ref.  we proposed a general algebraic method of constructing the covariant form of the integrals of motion. Developing this idea, here we demonstrated how

* The spectral problem can be treated in terms of the intrinsic algebra.

The advantage of this approach is its complete covariance. For example, one does not need to impose any gauge conditions since the intrinsic algebra operates with covariant quantities. The application of this method allowed us to show that the noncommutative systems under consideration are quasi-exactly solvable for any value of the coupling constant.

In the noncommutative case the reduction of the Hamiltonian and the application of intrinsic algebra to the spectral problem lead to some finite-difference equations. To investigate their quasi exact solvability,

* A deformation of the $`𝔰𝔩(2,)`$ scheme for finite-difference QES systems was proposed.

In this approach the $`𝔰𝔩(2,)`$ algebra is changed for its nonlinear analogue (A.6). We showed that the deformed scheme is not equivalent to the approach based on the usual $`𝔰𝔩(2,)`$ algebra  if one requires for finite-difference QES operators to be the polynomials in the generators of the corresponding algebras.

Surprisingly enough, but the system (2.4) in a non-constant magnetic field on the noncommutative plane (the coordinates form the Heisenberg algebra) is more complicated than the systems on the fuzzy sphere and the noncommutative hyperbolic plane investigated in this paper. Therefore, it is of interest to investigate this system and verify if it inherits such a property as the quasi exact solvability. Furthermore, it would be interesting to apply the observed analogy between the matrix equations and the parafermionic trilinear commutation relations  to construct new solutions. Such a consideration will be presented elsewhere .

Acknowledgements

We thank W. Orrick and J. Perk for useful comments. M.P. is also grateful to A. Turbiner for valuable discussions and kind hospitality during his visit to Universidad Nacional Autonoma de Mexico, where part of this work was done. The work was supported by the grants 1010073 and 3000006 from FONDECYT (Chile) and by DICYT (USACH).

## Appendix A Deformed $`𝔰𝔩(2,)`$ scheme

Here we discuss a deformation of the $`𝔰𝔩(2,)`$ scheme related to one-dimensional finite-difference operators, and compare the approaches based on the deformed algebra and on the usual one.

In Ref.  Smirnov and Turbiner proposed the approach to finite-difference QES equations based on the following representation of the Heisenberg algebra:

$$a=\mathrm{\Delta }_+,b=x(1\theta \mathrm{\Delta }_{}),$$
(A.1)
where $`x`$ is a real variable and $`\mathrm{\Delta }_\pm `$ are the forward and backward difference operators (3.3) with respect to this variable. These operators obey the standard commutation relation, $`[a,b]=1`$.

Like in the standard approach to 1D differential QES equations, an operator for the corresponding spectral problem is supposed to be a polynomial in the $`𝔰𝔩(2,)`$ generators

$`J_+`$ $`=b^2a(𝖭1)b,`$ $`J_0`$ $`=ba{\displaystyle \frac{𝖭1}{2}},`$ $`J_{}`$ $`=a.`$ (A.2)
When $`𝖭`$, the space of polynomials $`𝒫_{𝖭1}(x)`$ is invariant with respect to the action of these operators. Hence, any polynomial in the generators (A.2) is an invariant operator on this space as well, i. e. it is a QES operator.

We are interested in the spectral problem for a three-point finite-difference equation,

$$A(x)\mathrm{\Psi }(x+\theta )+B(x)\mathrm{\Psi }(x)+C(x)\mathrm{\Psi }(x\theta )=E\mathrm{\Psi }(x).$$
(A.3)
Ref.  claims that all such QES equations with finite number of polynomial eigenfunctions are classified via the cubic polynomial element of the universal enveloping algebra of $`𝔰𝔩(2,)`$ taken in the representation (A.1):

$$𝖳=A_+\left(J_++\theta J_0^2\right)+A_1J_0^2\left(\theta J_{}+1\right)+A_2J_0J_{}+A_3J_0+A_4J_{}+const,$$
(A.4)
where the coefficients are arbitrary real parameters. Below we shall show that this operator is not the most general among QES operators corresponding to the three-point problem (A.3).

Let us introduce the finite-difference operators which are linear in the finite-difference derivatives,

$`J_+^\pm `$ $`=x^2\mathrm{\Delta }_\pm 2jx,`$ $`J_0^\pm `$ $`=x\mathrm{\Delta }_{}j,`$ $`J_{}^\pm `$ $`=\mathrm{\Delta }_{}.`$ (A.5)
The half-integer parameter $`j`$ ($`2j_+`$) is an analogue of “spin” in the usual $`𝔰𝔩(2,)`$ scheme . These linearly independent operators are finite-difference analogues of the $`𝔰𝔩(2,)`$ generators, which can be reproduced in the commutative limit, $`\theta 0`$. The operators (A.5) have another useful property of the $`𝔰𝔩(2,)`$ generators. They have the invariant finite-dimensional subspace, which is the space of polynomials $`𝒫_{2j}(x)`$, $`2j_+`$, $`𝖽𝗂𝗆𝒫_{2j}(x)=2j+1=𝖭`$. This allows us to use the generators (A.5) for the construction of finite-difference QES operators in the manner analogous to the continuous case. The exhaustive investigation of this construction lies out of the scope of the paper and here we will discuss just some aspects of it.

The both sets $`J_\mu ^+`$ and $`J_\mu ^{}`$, $`\mu =0,\pm `$, form the nonlinear subalgebras of the form

$`[J_{}^\pm ,J_+^\pm ]`$ $`=2J_0^\pm \pm \theta (2j1)J_{}^\pm ,`$
$`[J_0^\pm ,J_+^\pm ]`$ $`=J_+^\pm \theta \left(2J_+^\pm J_{}^\pm +(2j+1)(J_0^\pm +j)\right),`$ (A.6)
$`[J_{}^\pm ,J_0^\pm ]`$ $`=J_{}^\pm T_{},`$
where the operators

$$T_\pm =1\pm \theta J_{}^{}$$
(A.7)
are the discrete forward and backward translation generators, $`T_\pm \psi (x)=\psi (x\pm \theta )`$. The intertwining commutation relations between the sets $`J_\mu ^+`$ and $`J_\mu ^{}`$ are

$`[J_{}^{},J_+^\pm ]`$ $`=\left(2J_0^{}\pm \theta J_{}^{}\right)T_\pm ,`$ $`[J_+^+,J_+^{}]`$ $`=\theta \left(2J_+^+J_0^+J_+^++(2j+1)J_+^{}\right),`$
$`[J_0^{},J_+^\pm ]`$ $`=\left(J_+^\pm \pm \theta (J_0^{}+j)\right)T_\pm ,`$ $`[J_0^+,J_0^{}]`$ $`=J_0^+J_0^{},`$
$`[J_{}^{},J_0^\pm ]`$ $`=J_{}^{},`$ $`[J_{}^+,J_{}^{}]`$ $`=0.`$ (A.8)
It is obvious that in the commutative limit, $`\theta 0`$, the nonlinear commutation relations (A.6), (A.8) tend to the usual $`𝔰𝔩(2,)`$ algebra. Furthermore, for the realization (A.5) of the nonlinear algebra the following quadratic identities hold:

$`J_0^\pm J_0^\pm J_+^{}J_{}^\pm J_0^\pm \pm \theta \left(J_0^\pm +j\right)J_{}^\pm j(j+1)=0,`$
$`\theta \left(J_0^++j\right)\left(J_0^{}+j1\right)+J_+^{}J_+^+=0,`$
$`J_+^+J_0^+J_+^{}J_0^{}(j+1)\left(J_+^+J_+^{}\right)=0,`$
$`J_+^{}T_+J_+^++2j\theta \left(J_0^{}+j\right)=0,`$ (A.9)
$`\left(J_0^{}+j\right)T_{}J_0^+j=0,`$
$`J_{}^{}T_{}J_{}^+=0.`$
Taking into account these identities and the commutation relations (A.6), (A.8), the most general three-point operator (A.3), quadratic in the generators (A.5), can be written as

$`H_{QES}`$ $`=`$ $`a_1J_+^+J_+^{}+a_2J_+^{}J_0^{}+b_+^+J_+^++b_0^+J_0^++b_{}^+J_{}^+`$
$`+b_+^{}J_+^{}+b_0^{}J_0^{}+b_{}^{}J_{}^{}+const.`$
By the construction this operator is quasi-exactly solvable.

The generators of the $`𝔰𝔩(2,)`$ algebra in the representation (A.2) can be represented as polynomials in the operators (A.5),

$`J_+`$ $`=\left(J_+^{}\theta (J_0^++j)\right)T_{},`$ $`J_0`$ $`=J_0^+,`$ $`J_{}`$ $`=J_{}^{}.`$ (A.11)
Hence, the QES scheme based on the $`𝔰𝔩(2,)`$ generators (A.2) is included into that based on the operators (A.5). On the other hand, the operators (A.5) can also be represented in terms of (A.2) via the two last equalities in (A.11) and

$`J_+^+`$ $`=J_++\theta \left(J_0+j\right)\left(\theta \left(J_0+j\right)J_{}+2J_01\right),`$ $`J_{}^+`$ $`=J_{}\left(1+\theta J_{}\right)^1,`$
$`J_+^{}`$ $`=J_+T_++\theta (J_0+j),`$ $`J_0^{}`$ $`=J_0+\theta \left(J_0+j\right)J_{}.`$
According to these relations the shift operator $`T_+`$ (A.7) is linear in the generator $`J_{}`$. As we see, the correspondence is not of a polynomial form for the operator $`J_{}^+`$. Therefore, the $`𝔰𝔩(2,)`$ scheme (A.2) is not equivalent to that based on the operators (A.5) if one restricts QES operators to a polynomial in the generators of the (A.2) and (A.5) forms. One can also notice that the 5-parametric (up to overall shift) QES operator (A.4) is not the most general QES operator in the framework of $`𝔰𝔩(2,)`$ scheme. Indeed, for the case $`b_3=0`$ the QES operator (LABEL:gQES) can be represented as a 5-th degree polynomial in the generators (A.2) with 7 arbitrary real parameters.

In addition, we note that if one imposes the following restrictions on the parameters of the QES operator (LABEL:gQES)

$`a_1`$ $`=0,`$ $`a_2`$ $`=0,`$ $`b_+^{}`$ $`=b_+^+,`$ (A.13)
then it becomes exactly solvable.