Separability of the Planar $1/\rho^{2}$ Potential In Multiple Coordinate Systems

With a number of special Hamiltonians, solutions of the Schr\"{o}dinger equation may be found by separation of variables in more than one coordinate system. The class of potentials involved includes a number of important examples, including the isotropic harmonic oscillator and the Coulomb potential. Multiply separable Hamiltonians exhibit a number of interesting features, including"accidental"degeneracies in their bound state spectra and often classical bound state orbits that always close. We examine another potential, for which the Schr\"{o}dinger equation is separable in both cylindrical and parabolic coordinates: a $z$-independent $V\propto 1/\rho^{2}=1/(x^{2}+y^{2})$ in three dimensions. All the persistent, bound classical orbits in this potential close, because all other orbits with negative energies fall to the center at $\rho=0$. When separated in parabolic coordinates, the Schr\"{o}dinger equation splits into three individual equations, two of which are equivalent to the radial equation in a Coulomb potential---one equation with an attractive potential, the other with an equally strong repulsive potential.


Introduction
There are certain special Hamiltonians for which the Schrödinger equation is separable in more than one coordinate system. The spectra of these Hamiltonians exhibit what are known as "accidental" degeneracies, and a number of the Hamiltonians are extremely important. Well known examples of systems that are separable in more than one set of coordinates include the free particle (separable in any coordinates for which the Robertson condition is satisfied [1,2]), the three-dimensional isotropic harmonic oscillator (separable in spherical, ellipsoidal, cylindrical, and rectangular coordinates), the Coulomb potential (separable in spherical, prolate ellipsoidal [3], and parabolic coordinates), and the constant magnetic field (separable in rectangular or cylindrical coordinates, depending on the gauge, but in either case with a free choice of the location of the origin). All these examples are also well known for the degeneracies present in their spectra. A slightly less well known system is a particular anisotropic harmonic oscillator, with potential That this energy spectrum possesses accidental degeneracy is obvious; however, it is not as well known that, just like the Coulomb problem, this potential problem is separable in parabolic coordinates. Separability in parabolic coordinates is guaranteed if the Hamiltonian involves just a standard nonrelativistic kinetic term and a potential V of the form for some functions f and g. The quantities η and ξ are two of the parameters of the parabolic coordinate system (η, ξ, φ), where (Surfaces of constant η or ξ are orthogonal paraboloids of revolution, each with its focus at the origin.) For the anisotropic harmonic oscillator (1), the two functions are f (s) = g(s) = 1 2 Mω 2 0 s 3 . The use of parabolic coordinates in the Coulomb problem is also especially convenient for dealing with the Stark effect. While the Hamiltonian remains separable in spherical coordinates when an external magnetic field is applied, it remains separable in the parabolic coordinates with a external electric field present.
Classically, all the potentials we have mentioned with accidental degeneracy are well known for another feature. All their bound orbits close. In this paper, we shall look at another potential that, in a sense, shares this classical feature. Bertrand's theorem is normally taken to hold that the only two central potentials for which all the bound orbits close are the attractive Coulomb potential and the isotropic harmonic oscillator potential. However, this is not quite accurate; there are other examples of potentials with the stated property, but they are typically discounted because they do not have a full spectrum of bound states. For example, all the bounds orbits in a constant potential, which exerts no force, close-precisely because there are no bound states. A system with a charged particle in a constant magnetic field also evades the strong restriction imposed by Bertrand's Theorem [4], since the force in this instance is not derived from a central potential. The orbits in the magnetic field are right circular helices. The velocity parallel to the magnetic field is a constant of the motion; only when the velocity in that direction vanishes are the orbits truly bound-in which case they are closed circles. The potential considered in this paper is in a similar category to the two examples just mentioned. We shall see that the attractive V ∝ 1/ρ 2 potential possesses few bound orbits (for an appropriate interpretation of the meaning of "bound"), but those that it does possess are circular, with exactly zero energy (as expected from the virial theorem).
Separability of the Schrödinger equation in a given set of coordinates means that the energy of a given state may be expressed as a function of three quantum numbers, one corresponding to each coordinate. If separations in multiple coordinate systems are possible, there must be multiple formulas for the energy, based on different sets on quantum numbers. The eigenstates with fixed quantum numbers are generally different in different coordinates; an eigenfunction in one coordinates must be a linear combination of eigenfunctions in the other coordinates. The only way this is possible for an energy eigenstate is if there are multiple eigenstates with exactly the same energy. This shows why separability is multiple coordinate systems requires the presence of accidental degeneracy. Moreover, it is worth noting that the same argument can be applied even to the degeneracy of system with a generic central potential V (r), which is separable in spherical coordinates only. A system with angular momentum ℓ possesses a (2ℓ + 1)-fold degeneracy, which is actually related to the fact that the spherical coordinates may be chosen with their polar axis pointing in any direction. The Schrödinger equation is thus separable in an infinite number of different spherical coordinate systems.
However, the existence of accidental degeneracy does not absolutely require that a system be separable in multiple coordinate systems. Any anisotropic three-dimensional harmonic oscillator for which the frequencies of the motions along the three coordinate direction are rational multiples will have classical orbits that eventually close and degeneracies in its quantum mechanical spectrum. So there are systems for which the degeneracy is seemingly too "sporadic" to be indicative of any deeper underlying symmetry principle at work.
The existence of alternative bases of quantum numbers is also related, of course, the existence of additional observables that commute with the Hamiltonian. For the spherical harmonic oscillator and the Coulomb potential, these extra conserved quantities are well known. The harmonic oscillator has separate, commuting Hamiltonians governing the motion along the three orthogonal axes. For the Coulomb problem, there is the Runge-Lenz vector, which points out the direction of the major axis of a bound state elliptical orbit. For the free particle, with its extensive degeneracy, the additional conserved quantity is the momentum itself. Using the algebras generated with the inclusion of any of these conserved quantities, it is possible to determine the bound state spectra of these problems using operator algebra alone.
In the Coulomb problem, part of the accidental degeneracy associated with the additional operators that commute with the Hamiltonian persists even in the relativistic Dirac theory, although the separability in parabolic coordinates actually does not carry over. Another potential, albeit a potential in only one dimension, that is also amenable to similar operator methods is the V = V 0 sech 2 ax potential; the eigenstates for potentials of different depths V 0 are related by operators [5], and as with the other potentials solvable by operator methods, the operators involved may be interpreted as elements of a (0 + 1)dimensional supersymmetry algebra [6,7,8,9]. (An excellent introductory treatment of quantum-mechanical supersymmetry, especially as applied to the Coulomb problem, may be found in the lecture notes, Ref. [10].) The harmonic oscillator, Coulomb, and sech 2 ax potentials just mentioned are all among the shape-invariant potentials, and the main object of study in this paper will be yet another such one. There are a wide variety of mathematical tools that may be useful for solving and addressing questions about these shape-invariant potentials [11,12].
The harmonic oscillator and Coulomb systems are also known for the fact that certain classical and semiclassical approximations yield exact results when applied to these systems. The Bohr-Sommerfeld quantization rule derived from the Wentzel-Kramers-Brillouin (WKB) approximation gives the exact energies for a harmonic oscillator system in one dimension-and hence also for an isotropic or anisotropic harmonic oscillator in any number of dimensions. Harmonic oscillators also have coherent states, with the zeropoint uncertainties in position and momentum added onto a classically orbiting wave function centroid. For the Coulomb potential, there is the fact that the full nonperturbative scattering cross section is the same as the cross sections derived classically or from the first-order Born approximation. Moreover, for the V = V 0 sech 2 ax potential also, a certain approximation is exact; the potential is reflectionless, so the classical reflection coefficient is always precisely correct. The exact successes of these various approximations are highly appealing features of these special potential types, although it is not clear whether we should expect anything similar when dealing with more esoteric multiply separable potentials. This paper will examine the 1/ρ 2 = 1/(x 2 + y 2 ) potential in three spatial dimensions. The attractive version of the potential has been observed physically, in the interaction of a long charged wire with a polarizable atom [13]. Section 2 introduces the features that make this potential special classically. Finding analytical solutions of the equations of motion is easy, and there are relatively few bound, stable orbits, but those that do exist close. The simplicity of the classical problem leads us to suspect that the corresponding quantum mechanical theory may also exhibit special properties.
In section 3, we demonstrate that the Hamiltonian with a 1/ρ 2 potential belongs to the elite family that are separable in multiple coordinate systems-cylindrical and parabolic, in this case. The scattering state wave functions of this potential in two space dimensions have already been studied [14], but it is only with the inclusion of the third dimension that the dual separations become possible. So new symmetry phenomena are expected to be found in the three-dimensional axisymmetric potential, although it is not really clear in advance what special behavior could be expected. If the potential possessed bound states, we would obviously expect them to have additional degeneracies. Unfortunately, however, the attractive 1/ρ 2 potential does not support a stable spectrum of bound negative-energy states; the attractive singularity at ρ = 0 is too strong. For the scattering states of a repulsive potential, we might hope, based on the behavior of other multiply separable Hamiltonian systems, to uncover some interesting behavior. The special features of the 1/ρ 2 potential that we have identified are discussed in 4, and our conclusions are summarized in section 5.

Classical Features of the 1/ρ Potential
Central potential problems in three (or more) dimensions may be reduced to two-dimensional problems using the conservation of angular momentum. The differential equation for the orbital curve r(φ) = [u(φ)] −1 then has the simple, well-known form where F (r) = −dV /dr is the radial force. The u term on the right-hand side of (5) corresponds to the centrifugal term in the effective potential governing the radial motion. For potentials V ∝ r n , the orbital shape can be expressed in terms of elementary trigonometric functions for n = 2, 0, −1, or −2 [4]. These are, respectively, the harmonic oscillator, the free particle, the Coulomb potential, and the inverse square potential of interest here. For certain other integral and rational values of the exponent n, the solutions may be expressed in terms of elliptic, hypergeometric, or other progressively more general functions. Because it is straightforward to reduce the central force problem from three to two dimensions, there is very little practical difference at the classical level between computations with a three-dimensional central potential V (r) and the two-dimensional analogue V (ρ). (The same cannot be said in quantum mechanics though; for example, a twodimensional attractive potential always possesses at least one bound state, but in three dimensions the bound state need not be present.) The classical radial equation of motion in two-dimensional space, with a potential V (ρ) = κ/ρ 2 , using the effective potential (depending on the angular momentum L = L z ), is This is clearly just as solvable with κ = 0 as it is for the free particle (κ = 0) case. The equation for the orbital curve is correspondingly The full solutions in three dimensions are simply the two-dimensional (ρ, φ) motion superimposed upon uniform motion in the z-direction, z(t) = z(0) +ż(0)t. The classical approach we are using corresponds to the solution of the Schrödinger problem in cylindrical coordinates. The nature of the solutions for ρ(φ) = [u(φ)] −1 depends on the sign of 1 + 2κM/L 2 , and thus on the sign of κ. If κ > 0, the coefficient in parentheses on the right-hand side of (8) is automatically positive. The only possible trajectories in this repulsive potential are scattering orbits In a coordinate system with the phase angle δ = 0, the radial coordinate diverges at φ = ±π/2 1 + 2κM/L 2 , corresponding to a classical scattering angle which depends on the angular momentum, but not separately on the energy-a consequence of the scale invariance of the problem. This expression can also be cast in terms of the impact parameter b in the xy-plane, via L = √ 2ME ′ b, where E ′ > 0 is the energy of the in-plane motion (so that the total energy is E = E ′ + 1 2 Mż 2 ). For attractive potentials with κ < 0, there are similar scattering orbits when the energy is positive (which means L 2 > 2|κ|M). The scattering angle is again given by (10); the absolute value present in that formula, which was superfluous for the repulsive potential, is needed in the attractive case to give a nonnegative ϕ scat .
The angle ϕ scat represents the scattering angle in the xy-plane. When the uniform motion in the third dimension is included, it is also possible to describe the total scattering angle ϑ scat . Since the potential in three dimensions is not spherically symmetric, the scattering behavior does not depend solely on an impact parameter (or equivalently, for fixed energy, on an angular momentum). Instead, we shall describe the incoming trajectory of a particle by a directionΘ, together with the angular momentum component L = L z . Choosing an appropriate orientation for the x and y coordinates,Θ iŝ The z-axis around which the potential V (ρ) is symmetric and the incoming trajectory (along which the particle would travel if it were not deflected) are generally skew lines. Their distance of closest approach to one-another is given by the in-plane impact param- After the scattering, the in-plane component of the velocity has been rotated through an angle ±ϕ scat , making the outgoing direction vector Therefore, the three-dimensional scattering angle ϑ scat is where ϕ scat is still a function of L or b, according to (10). Naturally, when the motion is planar (ż = 0), this gives ϑ scat = ϕ scat . Conversely, when the velocity in the z-direction (which does not change) predominates, ϑ scat → 0. Apart from the sign difference inside the absolute value in (10), there is another important difference between the attractive and repulsive regimes. When κ > 0, the scattering angle is limited to the range 0 ≤ ϕ scat < π; the trajectory never crosses itself. In contrast, when κ < 0, the angle ϕ scat may be arbitrarily large. When the potential is attractive, the particle may orbit around the center any number of times before it escapes again to infinity. The resulting two-dimensional trajectories cross over themselves repeatedly. This is quite different than the behavior seen in classical Rutherford scattering, in which the trajectories for attractive and repulsive potentials are represented by the two disjoint branches of the same hyperbola. However, this behavior, with the number of times the orbital curve intersects itself increasing as the total energy approaches zero, is by no means unique to the attractive 1/ρ 2 potential, but is in fact fairly generic.
The κ < 0 scattering orbits, with more and more revolutions around the origin as the energy decreases, are approaching the limit of perfectly circular orbits, which occur when the energy vanishes at L 2 = −2κM. Any attractive potential will have classical circular orbits. However, for the potential we are interested in, it turns out that these circular orbits are, in a certain meaningful sense, the only bound orbits. If the total energy is negative, then the quantity in parentheses in (8) is negative, and the orbital solution u(φ) becomes a linear combination of equiangular spirals, so that or depending on whether the both endpoints lie at ρ = 0 or one at ρ = 0 and the other at ρ = ∞. Note that all the bound orbits in two dimensions do therefore (in a certain sense) close, because the circular orbits are the only persistent bound orbits. Alternatively, taking the negative-energy solutions of (8) as linear combinations ρ(φ) has the form (14) In this regime, the time t f required for the particle to reach ρ = 0 is This is behavior is known as "falling to the center." Note that makingρ(0) more negative while keeping ρ(0) fixed decreases t f , which is clearly correct; the effects of varying ρ(0) independently are less intuitively obvious.

Two Separations of the Schrödinger Equation
We now turn our attention to the quantum theory. The classical theory with an attractive potential was dominated by the falling to the center. The quantum mechanical wave functions in the presence of the κ < 0 potential exhibit their own manifestation of this phenomenon. If ψ( ρ ) is an eigenfunction of the two-dimensional time-independent Schrödinger equation (with energy eigenvalue E ′ 0 ), then for any real number α, ψ(α ρ ) is also an eigenfunction, with energy α 2 E ′ 0 . Thus, if there is a normalizable eigenstate with energy E ′ 0 < 0, then there must be eigenstates with arbitrarily negative energies. By compressing the wave function closer to the attractive singularity at ρ = 0, the energy may be made as negative as we wish, meaning that there cannot be a stable Hilbert space of quantum states; the energy is not bounded below.
In contrast, when κ > 0, the energy eigenvalues are always positive. We may still dilate the wave function to decrease its energy, but the energies remain bounded from below by zero; and it is, of course, no surprise that all positive energies are allowed in this scattering system.
It is possible, via one of several renormalization procedures, to let the strength of the attractive 1/r 2 potential go to zero, in such a way that there is a reasonable physical spectrum (with exactly one bound state) [15,16,17]. However, the resulting Hamiltonian is not self-adjoint, in spite of it having a naively Hermitian appearance [18]; moreover, the scale invariance of the solutions is broken by an anomaly. Most analyses of the regularized Hamiltonian have focused on the three-dimensional 1/r 2 potential, although the general character of the solutions appears to the independent of the dimensionality [15,19,20,21]. Our results in this paper might be extended to this renormalized regime; in fact, it would be very interesting to see how the renormalization would interact with the separation of the Schrödinger equation in parabolic coordinates. However, this regime lies beyond the scope of the present work.

Cylindrical Coordinates
In order to have a quantum theory with well-defined wave functions, without additional regularization of the potential in the vicinity of ρ = 0, the potential we shall consider in the remainder of our analysis is with repulsive K > 0. (The potential strength κ has been rescaled to avoid unnecessary factors of and M.) With this potential, the Schrödinger equation in cylindrical coordinates is where E = 2ME/ 2 , with E being the total energy. With a separable ansatz, this reduces to a single-variable Schrödinger equation for P (ρ), This is just the usual Bessel's equation that arises for a free particle in two dimensions, except with the indices of the Bessel function solutions changed to √ m 2 + K. (Note that m 2 + K just corresponds to the classical quantity L 2 z + 2Mκ, measured in units of .) The general solution is thus If the Bessel function index √ m 2 + K happens to be an integer, the usual replacement of the linearly dependent J − √ m 2 +K by the Neumann function N √ m 2 +K is required. However, only the Bessel function with positive index is regular at ρ = 0 (and thus permitted).
The scattering theory of these solutions is straightforward. The scattering by a threedimensional 1/r 2 potential is worked out in [22]. One surprising result is that the classical limit only exists for strong potentials, as the classical scattering cross section is linear (never quadratic) in the strength of the potential. The two-dimensional 1/ρ 2 version is completely analogous, merely using the formalism for partial wave scattering in two dimensions [23]. The partial wave expansion for an incoming plane wave is in terms of the free wave radial functions. The Bessel functions have the limiting behavior J ν (s) = 2 πs cos s − νπ 2 − π 4 , and the scattering state wave function may be written The non-free wave functions with J √ m 2 +K (qρ) are simply phase shifted by The fact that δ m is independent of the energy for each partial wave is another consequence of the scale invariance. Moreover, as noted, the classical limit corresponds to K ≫ |m|.

Parabolic Coordinates
Unlike the scattering solution in cylindrical coordinates, the solution of of the 1/ρ 2 potential in parabolic coordinates is not a standard problem. A 1/ρ 2 potential has, however, previously been considered algebraically, as a perturbation added to the Coulomb Hamiltonian (which, as noted above, is also parabolic separable) [24]. In parabolic coordinates, the Laplacian is Whether an eigenfunction is separable in cylindrical coordinates, with ψ = P (ρ)Φ(φ)Z(z) or in parabolic coordinates ψ = H(η)Ξ(ξ)Φ(φ), we may take it to be an eigenfunction of L z , Φ(φ) = e imφ . Noting that ηξ = ρ 2 , taking this azimuthal dependence reduces the Laplacian plus potential in the parabolic coordinates to Once again, and not coincidentally, the inclusion of the potential corresponds to the change m 2 → m 2 + K.
With Φ(φ) factored out, the remaining Schrödinger equation can be written in the separation form Letting the first bracketed term in (30) be equal to a constant C, the ordinary differential equation for H is and with C → −C in the equation for Ξ.
Since the wave function is complex, it may not be automatically clear whether C should be real or complex. Note that a purely imaginary C gives the real and imaginary parts of the solutions to the ordinary differential equations definite behavior under inversions of the variables, η → −η or ξ → −ξ. However, this behavior is not actually physically mandated by the theory, because the physical space is limited to the parameter region where both η and ξ are nonnegative. It will, however, necessarily be the case that only a one-parameter family of C values will correspond to physically meaningful states. Any separable energy eigenfunction in three dimensions is determined (up to phase and normalization) by the values of the three real quantum numbers. In this system, we have the physical observables represented by m and E, so the choice of C must provide exactly one additional real degree of freedom. Since with a real-valued C, the separate differential equations for H and Ξ can be cast as eigenvalue equations for Hermitian operators, a real C is a sufficient condition for having equations that yield bases of wave functions with asymptotic forms that are continuum normalizable. Thus, a real C gives the correct one-parameter family of solutions.
The linearly independent solutions of (31) are expressible in terms of the confluent hypergeometric fuctions 1 F 1 (a; b; s), In order to have regularity at the origin (where η = ξ = 0), we must have the solution H = h + (η), and for Ξ, (33) Unfortunately, the overlap integrals giving the weights needed to write the wave functions H(η)Ξ(ξ) as superpositions of the P (ρ)Z(z) in cylindrical coordinates are intractable in the general case.
The asymptotic behavior of 1 F 1 (a; b; s) for |s| → ∞ and − 3π 2 < arg s < π 2 is For the solutions h + and h ′ + , the relevant values of a and b − a(= a * ) always have real parts 1 2 + 1 2 √ m 2 + K, which means that the first and second terms in (34) are of the same magnitude when when η or ξ is large. The large η behavior of the confluent hypergeometric function appearing in h + (η) is accordingly [using the phase convention that s = e iπ 2

√
Eη, corresponding to that in (34)], In the intermediate formula (36), arg Γ(a) has been used to abbreviate the complex argu- It follows that the asymptotic behavior of the wave function (r → ∞, but away from from the z-axis, where η or ξ will vanish) is using ηξ = ρ 2 and η ξ = cot 2 θ 2 . The limiting form (39) away from the z-axis is clearly normalizable as a continuum state. On the other hand, in the vicinity of the z-axis, either h + (η) or h ′ + (ξ) is close to 1, while the other function-and the wave function ψ as a whole-scale as ∼ 1/ |z|, which is again normalizable behavior. This confirms that our earlier choice of a real separation constant C was the correct one for the physical wave function solutions.

Special Features
Remarkably, the separated equation (31) for H can actually be cast in nearly the same form as the radial Schrödinger equation for a Coulomb potential. Letting U 1 (η) = √ ηH(η), (31) becomes The ordinary differential equation of U 2 (ξ) = √ ξΞ(ξ) is identical, except for the switch C → −C, equivalent to interchanging an attractive Coulomb potential with a repulsive one. Moreover, the normalization condition for the wave function, sets the same kinds of of constraints on how quickly the functions U 1 and U 2 must decay at spatial infinity as in the Coulomb problem. The equivalence also immediately explains the presence of the log( √ Eη) and arg Γ(a) terms in the argument of the cosine in (37), since these same kinds of terms appear in the phases of Coulomb waves.
The transformation of the separated parts of the Schrödinger equation into Coulomblike forms opens up a number of tools that can be used to further analyze the wave function solutions. However, those tools may play different roles in the analysis of the 1/ρ 2 potential than in the study of the 1/r potential. For example, there are (0 + 1)dimensional supersymmetry transformations that carry solutions of the radial Schrödinger equation in the Coulomb problem to other radial solutions with the same energies but different values of the angular momentum (changing l ↔ l + 1) [9,10]. Applied to (40), these transformations would still leave the energy affected (and also the separation constant C unchanged), but the strength of the potential would be modified through a change to the quantity m 2 + K, which combines the z-component of angular momentum with the strength of the repulsive potential. This is analogous to the situation with the one-dimensional sech 2 ax potential, where the supersymmetry transformations connect potentials with the same functional form, but of different depths.
Another interesting feature of the solution in parabolic coordinates stems from the fact that the choice of coordinate system breaks the translation invariance along the zdirection. This symmetry is manifestly present in the equations of motion in a cylindrical coordinate system, where z is a cyclic coordinate. Since z never enters the dynamics explicitly, there is nothing special about the location of z = 0. The translation symmetry is obscured somewhat in the parabolic coordinates, but it must still exist. If ψ 1 = H(η)Ξ(ξ)e imφ is a solution of the Schrödinger equation, then must also be a solution, since it is simply a translate of ψ 1 along the z-direction. The degeneracy of these states is analogous to the energy degeneracy of the Landau levels for a charged particle moving in the plane perpendicular to a constant magnetic field. The magnitude of the degeneracy is proportional to the area of the plane, since the origin of the coordinates may be located anywhere in the plane.
The presence of C in the eigenfunction equation (40) also appears to break the scale invariance of the problem, since the Hermitian operator on the left-hand side contains C, which has units of (length) −1 . However, since C is merely a separation constant, which can take any real value, a rescaling ρ → α ρ (and thus η → αη, ξ → αξ) may be accompanied by C → α −1 C. Since the same C appears in the equations for H and Ξ, this restores the physical scaling invariance.
Finally, the separation of the quantum-mechanical problem in parabolic coordinates can give some insight about the classical behavior in that coordinate system. Because the η and ξ portions of the Schrödinger equation are the same as those for two Coulomb problems, one attractive and the other equally repulsive, we can apply the normal methods for solving the Kepler problem to the classical time evolution of a particle's (η, ξ, φ) coordinates. Recalling the that the classical limit applies when K is large, we may neglect the −1 in the m 2 + K − 1 appearing in (40). Then, restoring the factors of 2 /2M, (40) corresponds to a classical limit of

Conclusions
Hamiltonians that are amenable to separation of variables methods in more that one coordinate system have a number of important properties. These include accidental degeneracies in their bound state spectra, and classical behavior that typically involves bound orbits that always close. The nonrelativistic Hamiltonians for number of important physical systems, such as the hydrogen atom and the charged particle in a constant magnetic field, are multiply separable in this way. These features are also tied to the usefulness of operator methods in solving these Hamiltonians. The z-independent 1/ρ 2 potential in three dimensions is obviously separable in cylindrical coordinates (ρ, φ, z), and we have shown that it is also separable in parabolic coordinates (η, ξ, φ). Although parabolic coordinates are not used nearly as frequently as rectangular, spherical, and cylindrical coordinate systems, they were already known to be useful for addressing certain aspects of the Coulomb problem. The attractive 1/ρ 2 potential is too strong to support a stable set of bound states, but the repulsive version is well behaved. When separated in parabolic coordinates, the one-dimensional Schrödinger equations for the component functions H(η) and Ξ(ξ) have the same forms as the radial Schrödinger equation in the Coulomb problem, although the strength of the Coulomblike term is set by the separation constant C, so that one equation features the effective equivalent of an attractive potential, which the other has a repulsive potential of equal magnitude.
The asymptotic behavior of the of full wave function ψ(η, ξ, φ) in parabolic coordinates is thus determined by the limiting behavior of a product of Coulomb waves, albeit ones that typically have nonintegral values for what would normally be the angular momentum parameter ℓ. Moreover, although we have focused on the scattering states of a repulsive 1/ρ 2 potential, it is clear that for sufficiently large values of L z = m , the scattering wave function in an attractive 1/ρ 2 potential will have essentially the same structures. They will still be products of Coulomb radial functions in η and ξ, with unconventional values of ℓ and equal and opposite effective potential strengths.
While the separation of variables in cylindrical coordinates keeps the translation symmetry along the z-direction and the scale invariance of the problem manifest, both of these invariances are obscured in the parabolic coordinate system, which picks a particular z = 0 location about with the η and ξ coordinate surfaces are focused. The fact that these important features are hidden in the parabolic coordinates formalism suggests that there may be yet other interesting features of these potentials still to be uncovered. In any case, the planar 1/ρ 2 potential, as well as being a system of real physical significance [13], appears to be a fruitful theoretical laboratory for understanding the structure of mechanics in parabolic coordinates and the behavior of multiply separable quantum systems.