Quantum Exactly Solvable Radial Systems on Curved Space
A Family of Exactly Solvable Radial
Quantum Systems on Space of NonConstant
Curvature with Accidental Degeneracy
in the Spectrum
Orlando RAGNISCO and Danilo RIGLIONI
O. Ragnisco and D. Riglioni
Dipartimento di Fisica Universitá Roma Tre and Istituto Nazionale di Fisica Nucleare,
Sezione di Roma, I00146 Roma, Italy
\Emailragnisco@fis.uniroma3.it, riglioni@fis.uniroma3.it
\URLaddresshttp://webusers.fis.uniroma3/~ragnisco/
Received October 05, 2010, in final form December 07, 2010; Published online December 15, 2010
A novel family of exactly solvable quantum systems on curved space is presented. The family is the quantum version of the classical Perlick family, which comprises all maximally superintegrable 3dimensional Hamiltonian systems with spherical symmetry. The high number of symmetries (both geometrical and dynamical) exhibited by the classical systems has a counterpart in the accidental degeneracy in the spectrum of the quantum systems. This family of quantum problem is completely solved with the techniques of the SUSYQM (supersymmetric quantum mechanics). We also analyze in detail the ordering problem arising in the quantization of the kinetic term of the classical Hamiltonian, stressing the link existing between two physically meaningful quantizations: the geometrical quantization and the position dependent mass quantization.
superintegrable quantum systems; curved spaces; PDM and LB quantisation
(81S10; 81R12; 31C12)
1 Introduction
A Hamiltonian system of degrees of freedom is said to be completely integrable, in the Arnold–Liouville sense, if it possesses functionally independent globally defined and single valued integrals of motion in involution [2]. In analogy to the classical mechanics, a quantum mechanical system with degrees of freedom is called integrable if it is equipped with a set of algebraically independent observables (selfadjoint operators densely defined on a suitable Hilbert space) all of them commuting with the Hamiltonian itself. The maximum number of independent observables commuting with the Hamiltonian is ; in that case out of this set one can extract at least two (complete) sets of commuting observables, both containing the Hamiltonian. Such a system is called maximally superintegrable (M.S.).
From a classical point of view these systems are characterized by the so called “closed orbit property” for bounded motions, and from a quantum point of view these systems present a maximal degeneracy in the discrete spectrum. Over the years the M.S. systems have been extremely useful models in the description of physical phenomena, as for instance the well known oscillator and Kepler radial problems. Recently, exactly solvable 1dimensional systems have been proposed, above all in quantum mechanics, as models describing electronic properties in condensed matter physics, for instance in semiconductors and quantum dots [3, 4, 5, 6, 7, 8]; the main feature of these systems is the presence of an effective mass in the Schrödinger equation, that can be considered as a positiondependent particle mass. Since, from a classical point of view, there exists a complete classification of all radial M.S. systems (flat case Bertrand theorem [9], curved space Perlick theorem [10]), the main goal of this paper is to “translate” these classical systems in the quantum context. Speaking about quantum M.S. systems many authors have recently stressed the existing relations of such quantum systems with systems studied in supersymmetric quantum mechanics (SUSYQM) [11] and in fact in this paper we will achieve our goal by using the powerful tools of the (SUSYQM) in particular the intertwining techniques that in the end allow us to obtain the spectrum and the eigenfunctions in a pure algebraic way.
2 Perlick’s family of classical M.S. systems
Let us consider the standard radial Hamiltonian describing a particle subject to a central force:
A theorem proved by Bertrand in the nineteenth century [9] establishes that the only radial potentials such that all bounded orbits are closed turn out to be:
We can obtain a wider class of such systems if we take in exam the more general Hamiltonian:
(2.1) 
The kinetic term of this new class of Hamiltonian systems can be regarded as that of a particle moving according to a geodesic motion in a complete Riemannian manifold with the metric:
Considering this new scenario, Perlick managed to find a classification of the pair of functions admitting stable closed orbit property [10]. This classification is indeed the most general [12] and it contains as particular cases other systems with the closed orbit property as for example the so called multifold Kepler systems and Darboux III system [13]. The Perlick’s systems can be classified in two parametric families that can be viewed as a deformation of the two original Bertrand systems:
Family I (Kepler type)
(2.2) 
Family II (oscillator type)
In the above formulas, is just an additive constant, and indeed it could be harmlessly removed. However we prefer to keep it in order to be consistent with the formalism used in our previous papers [14, 12]. These systems were proven to be M.S. in a rigorous way in 2009 [14] and define a complete classification of all radial problems in a Riemannian manifold with radial symmetry.
3 The generalized Kepler problem
As said before, the Hamiltonians with radial symmetry turn out to be of physical interest above all in the quantum context. So, this section will be devoted to give a short introduction of the quantization problem of systems defined in a non flat space and to provide the tools needed to solve them exactly. The main problem is how to define the quantum analog of the kinetic term in (2.1), since there is an order ambiguity in its quantization . This problem always arises when in the kinetic energy term there is a dependence on the coordinates [15, 4, 7], which comes from the nonzero curvature of the space. So, our Hamiltonian could be quantized for example as:
(3.1) 
Since our goal is to preserve the M.S. property of the classical Perlick’s systems at the quantum level, we will set aside, for a moment, the choice of the parameters , and start directly from a M.S. quantum system. To this aim we choose the Perlick I space with the parameter ; this system is universally known as the generalized Kepler problem, namely the Kepler problem on a space of constant curvature [16, 17, 18, 19, 20, 21]. As many other M.S. quantum systems [11], this system is indeed linked to a socalled shapeinvariant system.
3.1 Shapeinvariant systems and intertwining technique
Let us introduce the 1dimensional hyperbolic Kepler Hamiltonian, namely the radial part of the 3dimensional Schrödinger equation on a space of constant negative curvature (“hyperbolic space” [17]):
It should be noted that the parameter is linked to the scalar curvature through the relation and in fact in the limit we recover the radial part of the usual flat Kepler problem. This Hamiltonian system is a member of a small, though extremely relevant, family of systems called shapeinvariant systems: the reason for this name will be clear soon.
These systems can be solved easily by the “intertwining technique” [22], which can be seen as a generalization of the Dirac ladderoperator technique. To clarify this point, let us define the “prepotential” function:
and the following two ladder operators:
The Hamiltonian turns out to be factorized by the ladder operators:
Since we have:
therefore we can read as an eigenfunctions of the Hamiltonian .
The crucial property turns up by intertwining the two ladder operators:
which shows that we get the same Hamiltonian system, apart from the value of the parameter that turns out to be changed to . This is indeed the exceptional property of the shape invariant systems. This property allows us to find in a purely algebraic way the spectrum and the eigenvectors. In fact, let us define the following relations:
Multiplying by from the left:
and finally:
therefore can be considered the first excited state of .
By iterating this procedure it is possible to gain the whole set of bound states eigenfunctions as well as their eigenvalues. In particular the eigenvalues read:
while, up to a normalization factor, the finite set of the eigenfunctions will be:
As stated at the beginning of this section this system can be seen as the deformation of the radial part of the hydrogen atom Hamiltonian: in fact, in the limit , its eigenvalues (and eigenfunctions) turn out to be exactly those of the hydrogen atom:
3.2 Hyperbolic Kepler and Perlick systems
The radial part of the hyperbolic Kepler system is an exactly solvable 1dimensional quantum system. In this subsection we are going to see how to link this 1dimensional problem to one of the 3dimensional problems associated to the metric (2.2). This link can be found easily by making the following point transformation:
yielding
namely the quantization (3.1) with parameters , .
The scalar product will change to the new one:
Finally, in order to get a geometrically coherent Hamiltonian system, we perform a similarity transformation that (of course) doesn’t change the spectrum, but modifies the expressions of the eigenfunctions
We define this new operator as the Laplace–Beltrami (LB) Hamiltonian:
with eigenfunctions:
orthogonal to each other with the new weight function in the semiline :
The problem solved above is still a 1dimensional problem on the semiline. Then in order to establish an exact correspondence between this 1dimensional problem and the generalized Kepler problem we have to “upgrade” this 1dimensional problem to a 3dimensional one. We define the standard 3dimensional angular momentum operator in angular variables:
and their spherical harmonic eigenfunctions:
It is possible now write down the 3dimensional Hamiltonian operator:
We can recognize this Hamiltonian to be exactly the quantum version of the 3dimensional generalized Kepler system in which the kinetic part turns out to be the Laplace–Beltrami operator of the Perlick metric of type I with . This solves in a natural way the ordering problem of the quantization
(3.2) 
We underline also the presence of a (maximal) degeneracy in the spectrum, a property which is characteristic of a M.S. quantum system
4 The conformal Perlick problem
Let us rewrite the general Perlick I metric in a conformal way performing a point transformation:
(4.1)  
The classical Hamiltonian related to this new system of coordinates turns out to be:
Seen from this perspective, the physical meaning of this classical system appears to be twofold: in fact we could regard it as describing either a particle with a constant mass in a curved space or a particle with a nonconstant mass in a flat space.
In Section 3 we have investigated the hyperbolic Kepler system, namely the particular case of the Perlick I systems with . Therefore in this case the transformation (4.1) is:
In this new coordinate system the Hamiltonian (3.2) turns out to be:
The above Hamiltonian is of course the LB quantization of the generalized Kepler in the new coordinate system. Now the question is: it is possible to get a “physical”, though “nongeometrical” quantization of this Perlick I system without loosing the M.S. property? The answer to this question is affirmative, and indeed in the following we will introduce another M.S. quantization of this problem.
To this end, let us first rewrite the scalar product in this new variable:
with two possible domains:
Following the strategy of the previous section we get:
we define now the new Hamiltonian:
(4.2)  
We recognize easily the differential part of this operator as the radial part of the Laplacian operator; this is the so called Schrödinger quantization [23] which amounts to solve the order ambiguity in the quantization of the classical system putting the conformal factor in front of the differential part, in other words the quantizations (3.1) with parameters , . Therefore in this case the Laplace–Beltrami quantization and the Schrödinger quantization turn out to be linked through a similarity transformation, up to a constant shift proportional to the scalar curvature of this system. The set of eigenfunctions and eigenvectors are linked to each other, through easy algebraic manipulations. In Section 3 we showed the way to get by an algebraic computation the whole set of eigenfunctions of this Hamiltonian operator.
Before closing this section we wish to summarize all the algebraic manipulations we have carried out, and to display the general eigenfunctions written in a smarter way in term of Jacobi orthogonal polynomials
Laplace–Beltrami vs Schrödinger quantization
The above link between the Laplace–Beltrami and the Schrödinger quantization is far from being a lucky occurrence. Let us introduce a general conformal metric:
The scalar curvature of this conformal metric turns out to be:
The classical Hamiltonian of a free particle in such space is:
Computing the Laplace–Beltrami quantization we get:
while the Schrödinger quantization leads to:
Now by a direct computation it is possible to verify that the two quantum operators are related by a similarity transformation, apart for an additive function which bears a deep geometrical meaning, inasmuch as it represents the scalar curvature of the space:
Therefore the spectrum of the two different quantizations coincide up to a constant shift only in the case of constant curvature: otherwise the scalar curvature term becomes a radial potential term and cannot any more be considered as a shift in the eigenvalues.
5 Nonconstant curvature quantum Perlick systems
The machinery showed up to now plays a fundamental role in the quantization of the general Perlick problem I, namely the metric (2.2) with . The general case is, geometrically speaking, a highly nontrivial space because the scalar curvature is no more constant as it was in the generalized Kepler:
In the previous sections we performed the quantization following two different schemes, but since the scalar curvature turns out to be positiondependent the quantizations are no more similarityequivalent and we need to make a choice. The most suitable choice, which keeps the solvability property, turns out to be the quantization (3.1) with parameters , :
Performing the standard separation of variables:
the Schrödinger problem is reduced to:
Now let us apply some algebraic manipulations already used in the past sections:
This turns out to be exactly the Hamiltonian operator (4.2) if we set: . Performing this substitution we found the new spectrum:
Let us consider as a rational number, i.e. , . The two different wave functions , with , are associated with the same energy eigenvalue, and therefore the degeneracy is still present in the general case (but rational).
Now it is crucial to underline that, since , are natural numbers, the degeneracy of the spectrum can be preserved only if is a rational number too; but this is exactly the restriction requested also in the classical case for the closed orbit condition.
Of course it is still possible to translate these results to the Laplace–Beltrami Hamiltonian by a suitable similarity transformation. Clearly, in the LB case we would find a slightly different Hamiltonian with a potential correction term proportional to the scalar curvature:
6 Conclusions
In the present paper, the first of a series on the radial quantum M.S., we have completely solved the first family of the quantum Perlick’s systems with parameter , giving contextually a few insights about the ordering problem. The case with is also doable and pretty much similar to the one investigated in the present paper; the main difference consists just in the set of bounded states that turns out to be infinite.
All the problems related to this family turn out to be solvable by the factorization technique (SUSY quantum mechanics). We have also provided the general solutions of the eigenfunctions in terms of the Jacobi orthogonal polynomials.
The property of solvability and degeneracy in the spectrum holds indeed for the second family as well, providing therefore a complete classifications of all the radial quantum M.S. Hamiltonians.
The results about the second family will be presented in a forthcoming paper.
Acknowledgments
We wish to thank our colleagues and friends A. Ballesteros, A. Enciso and F.J. Herranz for illuminating discussions and crucial suggestions about the content of this paper. The results reported here have been obtained in the framework of the INFNMICINN collaboration 2010, and the related research activity has been partially supported by the Italian MIUR, through the PRIN 2008 research project n.20082K9KXZ/005.
[1]Referencesref
Footnotes
 This paper is a contribution to the Proceedings of the Workshop “Supersymmetric Quantum Mechanics and Spectral Design” (July 18–30, 2010, Benasque, Spain). The full collection is available at http://www.emis.de/journals/SIGMA/SUSYQM2010.html
References
 Arnold V.I., Mathematical methods of classical mechanics, 2nd ed., Graduate Texts in Mathematics, Vol. 60, SpringerVerlag, New York, 1989.
 Gritsev V.V., Kurochkin Yu.A., Model of excitations in quantum dots based on quantum mechanics in space of constant curvature, Phys. Rev. B 64 (2001), 035308, 9 pages.
 Quesne C., Tkachuk V.M., Deformed algebras,positiondependent effective masses and curved spaces: an exactly solvable Coulomb problem, J. Phys. A: Math. Gen. 37 (2004), 4267–4281, mathph/0403047.
 Bagchi B., Quesne C., Roychoudury R., PseudoHermitian versus Hermitian positiondependentmass Hamiltonians in a perturbative framework, J. Phys. A: Math. Gen. 39 (2006), L127–L134, quantph/0511182.
 Cruz y Cruz S., Negro J., Nieto L.M., Classical and quantum positiondependent mass harmonic oscillators, Phys. Lett. A 369 (2007), 400–406.
 Mazharimousavi S., Mustafa O., Flatland positiondependent mass: polar coordinates, separability and exact solvability, SIGMA 6 (2010), 088, 8 pages, arXiv:1003.3003.
 Tanaka T., fold supersimmetry in quantum systems with positiondependent mass, J. Phys. A: Math. Gen. 39 (2006), 219–234, quantph/0509132.
 Bertrand J., Théorème relatif au mouvement d’un point attiré vers un centre fixe, C. R. Math. Acad. Sci. Paris 77 (1873), 849–853.
 Perlick V., Bertrand spacetimes, Classical Quantum Gravity 9 (1992), 1009–1021.
 Marquette I., Superintegrability with third order integrals of motion, cubic algebras, and supersymmetric quantum mechanics. II. Painlevé trascendent potentials, J. Math. Phys. 50 (2009), 095202, 18 pages, arXiv:0811.1568.
 Ballesteros A., Enciso A., Herranz F.J., Ragnisco O., Bertrand spacetimes as Kepler/oscillator potentials, Classical Quantum Gravity 25 (2008), 165005, 13 pages, arXiv:0803.3430.
 Iwai T., Katayama N., Multifold Kepler systems – dynamical systems all of whose bounded trajectories are closed, J. Math. Phys. 36 (1995), 1790–1811.
 Ballesteros A., Enciso A., Herranz F.J., Ragnisco O., Hamiltonian systems admitting a Runge–Lenz vector and an optimal extension of Bertrand’s theorem to curved manifolds, Comm. Math. Phys. 290 (2009), 1033–1049, arXiv:0810.0999.
 von Roos O., Positiondependent effective masses in semiconductor theory, Phys. Rev. B 27 (1983), 7547–7752.

Schrödinger E.,
A method of determining quantum mechanical eigenvalues and eigenfunctions,
Proc. Roy. Irish Acad. Sect. A 46 (1940), 9–16.
Schrödinger E., Further studies on solving eigenvalue problems by factorization, Proc. Roy. Irish Acad. Sect. A 46 (1941), 183–206.  Infield L., Schild A., A note on the Kepler problem in a space of constant negative curvature, Phys. Rev. 67 (1945), 121–122.
 Higgs P.W., Dynamical symmetries in a spherical geometry. I, J. Phys. A: Math. Gen. 12 (1979), 309–323.
 Leemon H.I., Dynamical symmetries in a spherical geometry. II, J. Phys. A: Math. Gen. 12 (1979), 489–501.
 Kalnins E.G., Miller W. Jr., Pogosyan G.S., Coulomboscillator duality in spaces of constant curvature, J. Math. Phys. 41 (2000), 2629–2657, quantph/9906055.
 Cariñena J.F., Rañada M.F., Santander M., Central potentials on spaces of constant curvature: the Kepler problem on the twodimensional sphere and hyperbolic plane , J. Math. Phys. 46 (2005), 052702, 25 pages, mathph/0504016.
 Bagrov V.G., Samsonov B.F., Darboux transformation, factorization, and supersymmetry in onedimensional quantum mechanics, Theoret. and Math. Phys. 104 (1995), 1051–1060.
 Iwai T., Uwano Y., Katayama N., Quantization of the multifold Kepler system, J. Math. Phys. 37 (1996), 608–624.