Highly excited bound-state resonances of short-range inverse power-law potentials

We study analytically the radial Schr\"odinger equation with long-range attractive potentials whose asymptotic behaviors are dominated by inverse power-law tails of the form $V(r)=-\beta_n r^{-n}$ with $n>2$. In particular, assuming that the effective radial potential is characterized by a short-range infinitely repulsive core of radius $R$, we derive a compact {\it analytical} formula for the threshold energy $E^{\text{max}}_l=E^{\text{max}}_l(n,\beta_n,R)$ which characterizes the most weakly bound-state resonance (the most excited energy level) of the quantum system.


I. INTRODUCTION
The Schrödinger differential equation with inverse power-law attractive potentials has attracted the attention of physicists and mathematicians since the early days of quantum mechanics. In particular, long-range power-law potentials play a key role in theoretical models describing the physical interactions of atoms and molecules (see [1][2][3][4][5][6][7][8][9] and references therein).
It is well known that the attractive Coulombic potential is characterized by an infinite spectrum {E k } k=∞ k=0 of stationary bound-state resonances with the asymptotic property E k→∞ → 0 − [9]. On the other hand, attractive radial potentials whose asymptotic spatial behaviors are dominated by inverse power-law decaying tails of the form can only support a finite number of bound-state resonances [9]. In particular, it is interesting to note that, for generic values of the physical parameters n and β n , the discrete energy spectrum of an attractive inverse power-law potential of the form (1) terminates at some finite non-zero energy E max (n, β n ) [1][2][3][4][5][6][7][8][9][10].
The main goal of the present paper is to present a simple and elegant mathematical technique for the calculation of the most excited energy levels E max (n, β n ) [11] which characterize the family (1) of attractive inverse power-law potentials. In particular, below we shall derive a compact analytical formula for the threshold (maximal) energies E max (n, β n ) which characterize the most weakly bound-state resonances (the most excited energy levels) of the radial Schrödinger equation with the inverse power-law attractive potentials (1) [12].

II. DESCRIPTION OF THE SYSTEM
We shall analyze the physical properties of a quantum system whose stationary resonances are determined by the radial Schrödinger equation [13] − 2 2µ d 2 dr 2 + 2 l(l + 1) where the effective radial potential V (r) in (2) is characterized by a long-range inverse power-law attractive part and a short-range infinitely repulsive core. Specifically, we shall consider a composed radial potential of the form The bound-state (E < 0) resonances of the Schrödinger differential equation (2) that we shall analyze in the present paper are characterized by exponentially decaying radial eigenfunctions at spatial infinity: where [14] In addition, the repulsive core of the effective radial potential (3) dictates the inner boundary condition ψ l (r = R) = 0 (6) for the characteristic radial eigenfunctions. The Schrödinger equation (2), supplemented by the radial boundary conditions (4) and (6), determine the discrete spectrum of bound-state eigen-wavenumbers {κ(n, β n , R)} [or equivalently, the discrete spectrum of binding energies E(n, β n , R)] which characterize the effective radial potential (3). As we shall explicitly show in the next section, the most weakly bound-state resonance (that is, the most excited energy level) which characterizes the quantum system (3) can be determined analytically in the regime [15][16][17] κr n ≪ 1 (7) of small binding energies, where the characteristic length-scale r n is defined by the relation . (8)

III. THE RESONANCE EQUATION AND ITS REGIME OF VALIDITY
In the present section we shall analyze the radial Schrödinger equation which determines the spatial behavior of the bound-state eigenfunctions ψ l (r) in the regime r > R. As we shall explicitly show below, the characteristic radial equation (9) can be solved analytically in the two asymptotic radial regions r ≪ 1/κ and r ≫ r n . We shall then show that, for small resonant energies in the regime κr n ≪ 1 [see (7) [15]], one can use a functional matching procedure in the overlapping region r n ≪ r ≪ 1/κ in order to determine the binding energies {E(n, r n , R, l)} [or equivalently, the eigen-wavenumbers {κ(n, r n , R, l)}] which characterize the marginally bound-state resonances of the radial Schrödinger equation (2) with the effective binding potential (3). We shall first solve the Schrödinger equation (9) in the radial region in which case one may approximate (9) by The general solution of the radial differential equation (11) can be expressed in terms of the Bessel functions of the first and second kinds (see Eq. 9.1.53 of [18]): where {A 1 , A 2 } are normalization constants to be determined below. Using the small-argument (r n /r ≪ 1) asymptotic behaviors of the Bessel functions (see Eqs. 9.1.7 and 9.1.9 of [18]), one finds from (12) the expression for the radial eigenfunction which characterizes the weakly-bound (highly-excited) states of the Schrödinger differential equation (9) in the intermediate radial region We shall next solve the Schrödinger equation (9) in the radial region in which case one may approximate (9) by The general solution of the radial differential equation (16) can be expressed in terms of the Bessel functions of the first and second kinds (see Eq. 9.1.49 of [18]): where {B 1 , B 2 } are normalization constants [19]. Using the small-argument (κr ≪ 1) asymptotic behaviors of the modified Bessel functions (see Eqs. 9.1.7 and 9.1.9 of [18]), one finds from (17) the expression for the radial eigenfunction which characterizes the weakly bound-state resonances (the highly-excited states) of the Schrödinger differential equation (9) in the intermediate radial region Interestingly, for weakly bound-state resonances (that is, for small resonant wave-numbers in the regime κr n ≪ 1), the two expressions (13) and (18) for the characteristic eigenfunction ψ l (r) of the radial Schrödinger equation (9) are both valid in the intermediate radial region r n ≪ r ≪ 1/κ [see Eqs. (14) and (19)]. Note, in particular, that these two analytical expressions for the radial eigenfunction ψ l (r) are characterized by the same functional (radial) behavior. One can therefore express the coefficients {B 1 , B 2 } of the radial solution (17) in terms of the coefficients {A 1 , A 2 } of the radial solution (12) by matching the two mathematical expressions (13) and (18) for the characteristic radial eigenfunction ψ l (r) in the intermediate radial region r n ≪ r ≪ 1/κ. This functional matching procedure yields the relations [20] and We are now in a position to derive the resonance equation which determines the binding energies {E(n, r n , l)} [or equivalently, the eigen-wavenumbers {κ(n, r n , l)}] of the weakly-bound (highly-excited) states which characterize the radial Schrödinger equation (2) with the effective radial potential (3). Using Eqs. 9.2.1 and 9.2.2 of [18], one finds the asymptotic spatial behavior for the radial eigenfunction (17). Taking cognizance of the boundary condition (4), which characterizes the bound-state resonances of the radial Schrödinger equation (2), one deduces from (22) Substituting Eqs. (20) and (21) into (23), one obtains the characteristic resonance equation for the highly-excited bound-state resonances which characterize the Schrödinger equation (2) with the effective radial potential (3).

IV. THE RESONANT BINDING ENERGY OF THE MOST EXCITED ENERGY LEVEL
The dimensionless ratio A 2 /A 1 that appears in the resonance equation (24) can be determined by the inner boundary condition (6) which is dictated by the short-range repulsive part of the effective radial potential (3). In particular, substituting (12) into (6), one finds Substituting the dimensionless ratio (25) into the resonance equation (24), one finally finds the expression [22] κr n = (−1) l+1 n − 2 for the dimensionless resonant wave-number which characterizes the most excited energy level (the most weakly bound-state resonance) of the radial Schrödinger equation (2) with the effective binding potential (3). It is worth emphasizing again that the analytically derived resonance equation (24) is valid in the regime [see (14) and (19) where {j ν,k } k=∞ k=1 are the positive zeros of the Bessel function J ν (x) [18,24]. Defining the dimensionless small quantity for the smallest resonant wave-number which characterizes the effective binding potential (3).

V. SUMMARY
We have studied analytically the Schrödinger differential equation with attractive radial potentials whose asymptotic behaviors are dominated by inverse power-law tails of the form V (r) = −β n r −n with n > 2. These long-range radial potentials are of great importance in physics and chemistry. In particular, they provide a quantitative description for the physical interactions of atoms and molecules [1][2][3][4][5][6][7][8][9].
Using a low-energy matching procedure, we have derived the analytical expression [see Eqs.
As a consistency check, it is worth mentioning that, in the special case of spherically symmetric (l = 0) wave functions, Eq. (24) reduces to the semi-classical result κ = 1/(a −ā) of [17], where a is the s-wave scattering length andā = πr n (n − 2) cot[π/(n − 2)]/Γ 2 [1/(n − 2)] [30]. In addition, it is worth emphasizing that the interesting work presented in [17] for the l = 0 case is based on the semi-classical WKB analysis, whereas in the present paper we have presented a full quantum-mechanical treatment of the physical system which is valid for generic values of the dimensionless angular momentum parameter l.

ACKNOWLEDGMENTS
This research is supported by the Carmel Science Foundation. I thank Oded Hod for helpful discussions. I would also like to thank Yael Oren, Arbel M. Ongo, Ayelet B. Lata, and Alona B. Tea for stimulating discussions.