A Novel Approach to Solve Quasiexactly Solvable Pauli Equation

The spectra for some specific forms of external magnetic fields in Pauli equation are obtained in the framework of the asymptotic iteration method (AIM). AIM is applied to find the solution of Pauli equation.When the method is applied to quasiexactly solvable systems, it not only easily gives the corresponding spectrum, but also produces accurate results for the eigenvalues of the system having sl(2) symmetry.


Introduction
It is well known that motion of the charged particles in electromagnetic fields described by Pauli equation is one of the fascinating phenomena in modern physics, providing a general approach to understanding the properties of electrons.The Hamiltonian for this system when ℏ =  = 2  =  = 1 is where   ,   are the momentum operators,  = 2 is the gyromagnetic ratio, and A and   are the vector potential of the electromagnetic field and the Pauli matrix, respectively.The Pauli equation has been solved exactly for uniform magnetic field   =   = 0,   =  0 , leading the Landau levels.Moreover, it was proved in [1] that for any general magnetic field   = (, ) perpendicular to the  plane, the ground state is solved exactly as determined in [2][3][4].The general result of [1] can be viewed as a very special case of the recently discovered QES systems that only a finite spectrum to be solved with algebraic methods [5][6][7][8][9].The Pauli equation has been solved in a quasiexactly solvable model [10] for possible (nonuniform) magnetic fields [11].If not impossible, it seems difficult to find other forms of the magnetic field such that the Pauli equation could be exactly solved.Therefore, the forms of magnetic fields, for which Pauli equation could be exactly solved, are indeed rather rare.Therefore, a natural question arising at this point is can the asymptotic iteration method be applied to these systems?The answer to this question is the main topic of this paper.An asymptotic iteration technique which improves both analytical and numerical solutions of the problems has been introduced by C ¸iftc ¸i et al. [12].Afterwards, it has been developed for many physical systems in order to obtain the whole spectra [13,14].It has been also suggested to solve the Schrödinger equation and developed for some quantum optical systems [15].Later, a considerable attention has been paid to AIM [12,[16][17][18][19].The results obtained by this method fit with the solutions of many exactly solvable differential equations in the form of the Schrödinger equation.It also gives accurate result for the nonexactly solvable Schrödinger equation such as sextic oscillator, cubic oscillator, and deformed coulomb potential, which are important in applications to many problems in physics.Although AIM has been applied to solve the Schrödinger equation, its application to the solution of the Pauli equation is still missing.It is well known that the AIM is very powerful in describing many physical problems.It improves both analytical and numerical solutions and also it is useful for understanding the nature of physical structures.
One can show that the termination relation Δ  () = 0 in the corresponding method reflects effectively the peculiar properties of the functions  0 and  0 .After appropriate change of the variable, their hidden polynomial nature can be revealed.On these polynomials of finite order, the finitedimensional representations of the hidden (2, ) algebraic structure underlying the quasiexact solvability are realized.Also with this hidden polynomial structure, the nonlinear supersymmetry discussed in [10] can be related.

Formalism of the Asymptotic Iteration Method
The AIM is proposed to solve the second-order differential equations and details can be found in [12].In this section we systematically apply the method to obtain eigenvalues of the Pauli equation for various nonuniform magnetic fields.Consider the following second-order differential equation: where  0 ̸ = 0 and  0 ,  0 are functions defined in  ∞ (, ) and they have sufficiently many continuous derivatives.  and   denote derivatives of  with respect to .It is easy to show that ( + 2)th derivative of the function  can be written as where   and   are given by the recurrence relations then the solution of ( 2) can be written as follows: where  1 and  2 are the integration constants.In calculating the parameters in (4), when  = 0, we take the initial conditions as  −1 = 1 and  −1 = 0 [14] and the termination condition Before going further we note that, in this paper, we will determine magnetic fields that lead to QES solvable potential forms in Pauli equation.We consider the magnetic fields only in asymmetric gauge.Solution of the Pauli equation including nonuniform magnetic fields has recently been studied in the framework of the QES systems [10,11].Ho and Roy presented a general procedure for determining possible (nonuniform) magnetic fields such that the Pauli equation becomes QES with an underlying (2) symmetry.But our paper includes a part of the motivation those QES systems which can be solved by using AIM and it gives the whole spectrum with accurate results.
The following section includes the solution of Pauli Equation within the framework of the AIM.
In this case components of the magnetic field B with asymmetric gauge are given by where () is an arbitrary function of  and   () = ()/.The Pauli Hamiltonian then takes the form The variables  and  can be separated by introducing the wavefunction Here  (−∞ <  < ∞) are the eigenvalues of   , and () = ( ) is a two-component function of .Then the Hamiltonian, , can be written as It is obvious that the Hamiltonian has been written in the one-dimensional form of the Schrödinger equation with the potential.This final form is the SUSY form of the Pauli equation with the superpotential () = () + .Now, our aim is to determine the magnetic fields that lead to the QES potential forms.By defining the magnetic fields in appropriate forms that construct QES potentials, we will try to apply AIM to the resultant Hamiltonian.
If we apply AIM quantization condition, that is, from (4) by substituting  0 and  0 , in order to obtain polynomial solutions we see that the parameters must be equal to  = −2, V = −2 and  values are given in Table 1; we will show the results for type 1 and  =  =  = 1 for simplicity.Note that there are two types of solutions that give us the polynomials, for example, for  = 2, V = −2 (as given above), and V = (6/)(( + 2)(2 + 3) − ).If we use the first value of V we get the results given in Table 1, but one can use the second value of  and get different results which give us the polynomials.We will not give this second situation here because of complication on the calculations.In this paper, we use only V = −2 type solution.If we use  = 2(+1)−, V = 2 + , ( + 1) = ( + 1) + 2,  = −2, and V = −2, we have  = 2(++1) and  = −2(+1).If we go back to definitions of , , , and , we find the following expression for the energy as follows: If the parameters in (18) are mapped to the parameters in energy values in [10,11], the energy values can be in accord.In Table 1, we will show the results for type 1 and  =  =  = 1 for simplicity.Similarly, in Table 2, we will show the results of different values of  for  =  = 1 for type 2 solution.

Class IX.
Let us assume that the magnetic field takes the form of With this magnetic field, we can write the following eigenvalue problem for  − after replacing  by : where  = ,  = √  2 − ,  = 2(1 + ),  = 2, and  = ( + 1) + 2.This potential form is the Class IX potential form.By applying the same procedure performed in the previous potential form, when  → 0 and  → ∞, () ∼  +1 exp(−/ − ), and after substituting () =  +1 exp(−/ − )() into (20), we have Here  = 2(+1)−, V = 2+, and (+1) = +2.In this form of differential equation, if we apply AIM to this final Advances in Mathematical Physics equation with  0 = 2(−(+1)/−/ 2 ) and  0 = /−(2+ )/ 3 , we will obtain the exact solutions of (20) under some constraints on the parameters , , , and .It can be easily seen that the solutions of this form of potential are equal to the solution of the first example except for definitions of , , , and .In this type of potential, the definitions of   and  for the type 1 solution can be obtained as below:

Results and Discussion
The energy obtained for Class IX potential in the form of parameters , , , , and  agrees with those in [20,21].
Similar to the Class I potential, the values obtained for the second type solution in Table 2 are valid for this potential.The energy in (22) depends only on , , and .For simplicity, the energy values for  = 1 are also tabulated in Table 1.
Similar to the procedure presented in this study, it is possible to obtain the same solutions for other class potentials described by Turbiner [6] except for Class IV and Class V potentials.For these two potentials, an appropriate superpotential to construct the potentials has not been written because of the constant  in the ()+.The other six classes satisfy the requirements of using AIM which is the differential equations conditions being transformable to a specific form in (2), with quite restrictive asymptotic conditions on the coefficient functions of () in (5).

Conclusions
In this work, we have shown that the AIM method can also be used to obtain the QES energies obtained in [10,11] by using the known forms of the magnetic field B with asymmetric gauge obtained in the same paper.We follow these corresponding studies and consider the magnetic field in that form.The procedure presented here makes full use of the close connection between QES systems and SUSY-like form.Additionally, an attempt to give the QES solution of the Pauli equation by using asymptotic iteration technique is presented.The application of AIM is considered for Class I and Class IX QES systems.By the same procedure, AIM can be applied to the other QES potentials in one dimension except for Class IV and Class V.The whole spectrum of the QES systems is presented.
Finally, we will try to apply this method to other QES potentials in all branches of physics.