arxiv: v2 [quant-ph] 30 Nov 2010

Transcription

1 A new quasi-exactly solvable problem and its connection with an anharmonic oscillator Da-Bao Yang, Fu-Lin Zhang, and Jing-Ling Chen Chern Institute of Mathematics, Nankai University, Tianjin 30007, China arxiv: v [quant-ph] 30 Nov 00 (Dated: September, 08) Abstract The two-dimensional hydrogen with a linear potential in a magnetic field is solved by two different methods. Furthermore the connection between the model and an anharmonic oscillator is investigated by methods of KS transformation. PACS numbers: Electronic address: chenjl@nankai.edu.cn

2 I. INTRODUCTION The discovery of quasi-exactly solvable quantal problems was a remarkable contribution in last century[]. This kind of problem could be solved by lie algebraic methods, which was pioneered by Turbiner[]. The method was also be applied to an anharmonic oscillators by Chen et. al.[3]. Moreover a generalization of the basic lie algebraic construction was considered by Shifman[4]. You can also confer to a reviewed paper by Shifman[]. However, the quasi-exactly solvable problem could also be calculated by an analytical way[5]. This paper is concentrated on a model that is a two-dimensional hydrogen with a linear potential in a homogeneous magnetic field. This model is very meaningful among the quasiexactly problems. First, its potential is novel and complex, moreover when the coefficient of the linear potential becomes zero, its all result could go back to the simpler circumstance which have been investigated by Taut[5]. Second, the sl() algebra structure hidden in the Hamiltonian is not apparent, however it will be revealed in this paper. Third, this model is also solved by an analytical method whose result coincide with the one by the sl() method. Last but not least, a new sextic oscillator with a centrifugal potential is obtained from this model. Furthermore, the connection between them is constructed. This paper is organised as follows: in section two, the separation of variables in the corresponding schrödinger equation is displayed, so we get the radial part equation; in section three, we solve it by the lie-algebraic methods, i.e., sl(); in section four, it is solved by an analytic way, therefore, the two results coincide with each other; in section five we build its connection with the sextic oscillator with a centrifugal barrier; in section six, a conclusion is drawed. II. THE SEPARATION OF VARIABLES Now, we consider a two-dimensional hydrogen with a linear potential in a homogeneous magnetic field, where the magnetic field B is perpendicular to the plan in which the electron is located. The schrödinger equation reads [ ( p+ c A) Z r +kr]ψ = Eψ, wherecis the velocity of light and the vector potential in the symmetric gauge is A = B r. Under some operation, the above equation can be transformed to be

3 [ p ψ + c ( p A)ψ + A pψ c + A c ψ] Z ψ +krψ = Eψ, () r In polar coordinates, letting r and θ represent the radius and azimuthal angle respectively, we can get p = i( e r r + e θ r ) θ A = B r = Br e θ p A = 0. A p = i B θ A = 4 B r Substitution the above equations into the Eq.(), we get [ ψ i θ ψ +r ψ] Z ψ +krψ = Eψ, () r where = r + and ω r r r r θ L = B. Having represented the equation in a concrete c coordinate, we will separate the variables in the next paragraph. Substituting ψ(r,θ) = π e imθ R(r) (3) into the above equation, where m = 0,±,±,..., we get the radial part equation, which is ( d dr d rdr + r +kr + m Z r + m r)r = ER. (4) The above differential equation will be solved by the method of sl() algebra in the next secion. III. SOLVING THE RADIAL EQUATION WITH THE SL() ALGEBRA First of all, let us recall the essential idea that the schrödinger equation is solved by sl(). The schrödinger equation determining the stationary states of the system is Hψ(r) = Eψ(r), d H = dr +V(r). According to [], we can perform a quasi-gauge transformation 3

4 ˆ ψ(r) = ψ(r)exp( A(r)dr). (5) Then the stationary schrödinger equation can be transformed to H G ψ = E ψ, H G = ( d dx A(x)) +V(x). Furthermore H G can be represented by the combination of the partial algebraization of the generators of sl(), which is where H G = C ab T a T b + C a T a +const, (6) a,b=±,0 a=±,0 T + = jr r d dr, T 0 = j +r d dr, (7) T = d dr. It is very easy to check that the commutation relations for T ± and T 0 are indeed those inherent to the sl() algebra: [T +,T ] = T 0, [T 0,T + ] = +T +, [T 0,T ] = T. Finally, with the aid of the above commutation relations, the transformed Hamiltonian H G can be represented by a (j + ) (j + ) matrix, whose basis are { ψ} = {,r,r,,r j,r j }. Then the standard algorithm of diagonalization of the finite matrix gives us (j +) eigenvalues and the corresponding eigenfunctions of the transformed Hamiltonian. Blindly following the program outlined above, we do not find the above model that we pay attention to in the section can be solved by the methods of sl(), in another word, the Hamiltonian H isn t reducible to the quadratic combination of the sl() generators. However, we can consider a corresponding eigenvalue problem, f(r)hψ = Eψ, where f 4

5 is an appropriately chosen function of the variable r, so that fh can be reducible to the standard form (6). The above explanation maybe a litter obscure, moreover this method will be showed in more transparent terms in the following paragraphs. Let us introduce an auxiliary Hamiltonian H = γ( d dr d rdr + r + m+ m Z r ) r +γkr γx, (8) where γ and X are constant coefficients.next, we calculate the gauge -transformed Hamiltonian H G, simply stemming form H after the substitution Then d dr d A. (9) dr H G = γ( d d dr rdr +A d dr ) γa +γ r +γkr+ γa +γ A r +γm r Z r +γ(m X). (0) After a short analysis, we can get a conclusion that: when r, the term γ r +γkr will become the dominant one; whenr 0, the termγ m r will become the dominant counterpart. In order to guarantee the normalizability of the wave function, we must choose a proper form of A to cancel out the above terms. A short reflection shows that Substituting the above Eq.() into Eq.(0), one gets A = µr+δ + ν r. () H G = γ( d dr r d dr +A d dr )+ γ( µ )r +γ(k µδ)r+( γδ γδν Z) r + γ(m ν ) r +γ(+ m+µ µν δ X) Moreover, from above analysis, one knows the coefficients of r, r and are equal to 0 r respectively, which are ωl µ = 0 k µδ = 0 m ν = 0 The above equations can be solved out: µ = ± δ = k µ ν = ± m 5. (). (3)

6 Substituting Eq.() and Eq.(3) into Eq.(4), one finds that only R(r) = R(r)exp( r k r + m lnr) (4) can guarantee that the wave function converge at both sides of r 0 and r (Note that we replace ψ(r) and ψ(r) by R(r) and R(r) respectively.), that is to say µ = + δ = k ν = m Substituting Eq.(5) into Eq.(), one can gets where H G = γ[. (5) r γ H E(γ) G R = R, (6) γ d dr r d dr +(r + k m r ) d dr ] +[γ k ( m + ) Z]+γ[(+m+ m ) ( k ) X]. The reason that I divided both sides of Eq.(6) byγ is for a convenient calculation.the reader maybe puzzled by the reason why we use r to multiple H G, because r H G can reducible to a quadratic combination of the generator of sl().under the inverse transformation Eq.(9), the Eq.(6) can be transformed to be where γ HR(r) = E(γ) R(r), (7) γr R = Rexp( r k r + m lnr). Substituting Eq.(8) into Eq.(7), we can get ( d dr d rdr + r +kr + m Z +E γ r + m r)r = XR. (8) By comparison with the Eq.(4), we observe that we have actually solved the original problem provided that the parameter γ is chosen as follows Z +E γ = Z. (9) And we represented the eigenvalue by another symbol X instead. 6

7 In order to solve Eq.(8) by sl(), let s return to Eq.(6). Following what I said before, the r γ H G can be reducible to the combination of generators of sl(), which is r γ H G = C T 0 T +C T + +C 3 T 0 +C 4 T +C 0. (0) Substituting Eq.(7) into the above Eq.(0), we can get r γ H G = C r d dr C r d dr +C 3r d dr +(C 4 jc ) d dr +C jr jc 3+ C 0. By comparison with Eq.(6), the coefficients can be solved from equations C = C = C 3 = k. C 4 jc = m A further simplified calculation show that C j = (+m+ m ) X jc +C 0 = k ( + m ) Z r C = C = C 3 = k C 4 = (+j + m ). () X = (j ++m+ m ) ( k ) C 0 = k ( m + +j) Z γ Substituting the above Eq.() and Eq.(0) into Eq.(6), one gets [ T0 T T + + k T 0 (+j + m )T + k ( m + +j)] R = Z +E γ R () By use of the constraint which is Eq.(9), the above Eq.() can be transformed to [ T0 T T + + k T 0 (+j + m )T + k ( m + +j)] R = Z R (3) )j=0 7

8 Hence, we know R = a 0. (4) Substituting Eq.(4) into Eq.(4), one get R = a 0 exp( r k r + m lnr). (5) Substituting Eq.(7) and Eq.(4) into Eq.(3), we can get k ( m + ) = Z, (6) which is precisely the constraint of eigenvalues. Substituting the above Eq.(6) into the fifth Equation of Eq.(), we can get the eigenvalue X 0 = (+m+ m ) ( Z m + ) )j = It is very easy to calculate that T + = 0 T 0 = 0 T = 0 0. (7) Substitution the above Eq.(7) into Eq.(3), we can get k ( m + 3 ) ω L a ( m + ) k ( m + ) = Z a. (8) a 0 a 0 Such a equation has a nontrivial solution if the determinant of the corresponding matrix vanishes, i.e. k ( m + 3) Z ( m + ) The equation can be further simplified to be k ( m + ) Z = 0. ωlz k ( m +)Z +k ( m + 3 )( m + ) ω3 L( m + )] = 0. (9) Form it, we can draw a conclusion that the parameters and k can t be chosen arbitrarily, they must satisfy the above Eq.(9). The above Eq.(9) can be regarded as a quadratic Equation about Z. As a result, the solution of it can be expressed as Z = k( m +)± M, (30) 8

9 where M k + 4 ω3 L ( + m ). In the following paragraph, we will calculate the different eigenfunctions corresponding to different eigenvalues. a)z = k( m +)+ M, Substituting the above equation into Eq.(8), we get a 0 k + ω M L + m Substituting the above Eq.(3) into Eq.(4) a 0 ( k + M + m r). (3) R(r) m j= j= = a 0 ( k + M + m r)exp( r k r+ m lnr). (3) ω L b)z = k( m +) M, Following the above procedure, we get the radial part of the wave function, i.e., R(r) m j= j= = a 0 ( k M + m r)exp( r k r + m lnr). (33) ω L With the constraint of Eq.(9), the eigenvalue of energy becomes X = (+m+ m ) ( k ) 3)Arbitrary j According to the choice of the phase, we know the following formulation, i.e., T ± jm = (j m)(j ±m+) j,m± T 0 jm = m jm. The j ± and j 0 can be represented in the following (j +) (j +) matrix, 0 j (j ) 0 0. T = (j ) j (34) 9

10 j T 0 (j ) = (j ) j 0 j j T 0 j =... j + j Substitution the above Eq.(34), Eq.(35) and Eq.(36) into Eq.(3), so the problem is deduced to an eigenvalue and eigenfunction problem. In order to get the nontrivial eigenfunction, the determinant of the corresponding matrix must vanish, which will give a onstraint about and k.similarly, we will get the radial part of the wave function. Furthermore, the fifth equation of Eq.() tell us the eigenvalue of energy, X j = (j ++m+ m ) ( k ). To sum up, in this section, we have solved Eq.(4) by method of sl(). Moreover, we will solved it again by another method, namely, an analytical methods. (35) (36) IV. SOLVING THE SAME EQUATION WITH ANALYTICAL METHOD In order to solve Eq.(4), using the following transformation, i.e. we can get R = u r, (37) d u dr + (m 4 ) Z r u r u+mu+kru+ r u = Eu. (38) Due to the condition of square-integrability of the wave function, it is useful to consider of the limits r 0 and r of the above Eq.(38). For r 0, the term containing r is 0

11 dominant and we get the equation d u dr + (m 4 ) ru = 0. (39) Because the above Eq.(39) contains a normal singularity, so we substitute a power series u = r s i=0 b ir i into the equation and one gets (s+i)(s+i )b i r s+i +(m 4 ) b i r s+i = 0. i=0 Due to the uniqueness of the series, the coefficient of it is equal to zero, we can get i=0 (s+i)(s+i )b i +(m 4 )b i = 0. By use of the limits r 0, only the lowest-order term become dominant, so we omits the higher order terms and get s(s )b 0 +(m 4 )b 0 = 0. In order to get the nontrivial solution, we must set b 0 0, so one gets s = m + or s = m+. However, if m become negative, the solution diverges when r 0. Hence, we get s = m +. So the solution of Eq.(39) is u = r m +. (40) Having solved the solutionr 0, we consider about another situationr. Under this circumstance, the terms containing r and r in Eq.(38) become dominant. So the equation becomes d u dr + r u+kru = 0. We can verify the solution of the above equation is u = exp( r k r). (4) Considering both Eq.(40) and Eq.(4), we can write the solution of Eq.(38)as Substituting the above Eq.(4)into (38), we get u = r m + e r k r ϕ. (4) ϕ +ω rd dr L r dϕ + k dr r dϕ ( m + dr )dϕ +( m ω dr L + +m +[( m + ) k Z]ϕ = 0. k ω L E)rϕ (43)

12 The above differential equation has two nonessential singularities at r = 0,. It s index equation takes this form s(s+ m ) = 0. It s solutions are 0 and m respectively. However, if s = m, the solution diverges, in another words, it violate the demand of square-integrability of the wave function. Hence, s must be 0. Therefore, we can substitute the power series ϕ = a n r n (44) into Eq.(43)and find the recurrence relation n=0 n= a n+(n+)nr n + n= a n (n )r n + k n= a nnr n ( m + ) n=0 a n+(n+)r n +( m + +m +[( m + ) k Z] n=0 a nr n = 0 Let us discuss in detail in the following paragraphs. )n = 0 k ωl E) n= a n r n a = ( m + ) k Z m + a 0. (45) We put forward a hypothesis that the ground state exists, so we get a 0 0, a = 0, a = 0,. However, conferring to Eq.(46), if a 0 0, a may be not equal to 0, which is contradict the hypothesis. So we can draw a conclusion that a 0 = 0, that is to say, there isn t a ground state. )n = k a = m + +m E ωl a 0 + ( m + 3) k Z a. (46) ( m +) ( m +) If the first excited state exists, a 0 and a can t be equal to 0 simultaneously, a = 0, a 3 = 0,. Conferring to Eq.(49), a will have a effect on a 3, so a = 0, that is to say, conferring to Eq.(45), So the radial part of wave function is Z = ( m + ) k. R n= = r m e r k r a 0, (47)

13 which is gained by use of Eq.(37) and Eq.(4).By comparison with Eq.(5), it is found easily that the two results coincide with each other. From above analysis, we can draw a conclusion that a 0 0, a = 0, a = 0,, conferring to Eq.(46), so we get 3)n = E = (+ m +m) k k ωl. (48) a 3 = m + +m E ωl 3( m + 3) a + ( m + 5) k Z 3( m + 3) a. (49) Following the above analysis, if the second excited state exists, we can get a 3 = 0, a 4 = 0 and so on. It is necessary to demand the first term and the second term are equal to 0 respectively, conferring to Eq.(49). Moreover if a 0, it may render a 4 0, which is contradict our hypothesis. So we must let a = 0, in another word, the second term become zero already. Now let us concentrate on the first term. If a = 0, this case will be the same with the first excited state. So we can get E = (+ m +m) Substituting Eq.(45) and Eq.(50) into Eq.(46), we can get k ωl a = [( m + 3 ) k Z][( m + ) k Z] ( m + ) ( m +)( m + ) a 0. From above analysis, we know a = 0, that is to say, [( m + 3 ) k Z][( m + ) k Z] ( m + ) = 0,. (50) which coincide with Eq.(9). Therefore, regarding Z as unknown root, we can certainly get the same solution with Eq.(9), i.e., Eq.(30). Substituting Eq.(4) and Eq.(44) into Eq.(37), we can also get the radial part of the wave function, R n= = r m e r k ω r L [ ( m + ) k Z m + Substituting Eq.(30) into the above Eq.(5), we get R ± n= = r m e r k r [ which is coincide with Eq.(3) and Eq.(33). 3 k ± M m + r]a 0, r +]a 0. (5)

14 4)n k a n+ = (n+ m +m) E ωl (n+)( m + +n) a n + ( m + +n) k Z (n+)( m + +n) a n. (5) Following the similar analysis above, also conferring to [5], it is a sufficient condition to guarantee normalizability of the eigenfunctions, i.e., E n = (n+ m +m) k ωl, (53) and For a 0 0, we can deduce a n = F( m,n,e,,k,z)a 0 = 0. F( m,n,e,,k,z) = 0. (54) Substituting Eq.(53)into above Eq.(54), we obtain F( m,n, (n+ m +m) k ωl,,k,z) = 0. (55) From it, we obtain a constraint for and k, in another word, and k can t be chosen arbitrarily. Furthermore, from Eq.(45), Eq.(46) and Eq.(5), we can get the polynomials of ϕ, i.e. Eq.(44). By use of Eq.(37) and Eq.(4), we can get the corresponding R n. In conclusion, we have solved the Eq.(4) by an analytical method. In the following section, we will build the connection between the radial part equation and the sextic oscillator with a centrifugal potential. V. KS TRANSFORMATION Via KS transformation, i.e., the Eq.() becomes r = ρ θ = ϕ, (56) [ 4ρ ρ + ρ ρ + ρ ϕ )ψ iω L ϕ ψ + ρ4 ψ] Z ρ ψ +kρ ψ = Eψ. ( Substituting ψ(ρ,ϕ) = π e i mϕ χ(ρ) (57) 4

15 into the above equation, where m = 0,±,±,..., we get ( d dρ + d ρdρ m ρ )χ+( mω L 4E)ρ χ+4kρ 4 ψ +ωl ρ6 χ = 4Zχ. In order to cancel the item which contains ρ, substituting into the above equation, one obtains d ζ dρ + 4 m 8 χ = ρ ζ (58) ρ ζ +( m 4E)ρ ζ +4kρ 4 ψ +ωl ρ6 ζ = 4Zζ, (59) which is the so-called sextic oscillator with a centrifugal barrier [6]. By comparison between Eq.(3) and Eq.(57), we can get m = m, (60) and χ(ρ) = R(r). (6) Substituting Eq.(58) into the above Eq.(6), we obtain ζ = ρr(r), (6) which is precisely the wave function of the sextic oscillator with a centrifugal barrier. Moreover, using Eq.(56) and Eq.(60), we can replace the parameters in Eq.(6) by the counterparts used in Eq.(59). In section two and three, we have calculate Z by two different methods. Therefore, it is very easy to get the eigenvalue 4Z of Eq.(59). VI. CONCLUSION AND ACKNOWLEDGMENTS In this paper, we have researched the two-dimensional hydrogen with a linear potential in a homogeneous magnetic field by two different methods, i.e., lie algebraic and analytic methods respectively. Furthermore we have also build the connection between the radial part of the schrödinger equation and the sextic oscillator with centrifugal barrier. This work was supported in part by NSF of China(Grants No and No ). [] M. A. Shifman, Int. J. Mod. A 4 (989)

16 [] A. V. Turbiner, Commun. Math. Phys. 8 (988) 467. [3] J. L. Chen, L. C. Kwek, C. H. Oh, and Y. Liu, J. Phys. A: Math. Gen. 34 (00) [4] M. A. Shifman, Int. J. Mod. A 4 (989) 33. [5] M. Taut, J. Phys. A Math. Gen. 8 (995) 08. [6] G. Levai, Phys. Rev. C 69 (004)