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

We study analytically the radial Schrödinger equation with long-range attractive potentials whose asymptotic behaviors are dominated by inverse power-law tails of the form V(r)=-βnr-n\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V(r)=-\beta _n r^{-n}$$\end{document} with n>2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n>2$$\end{document}. In particular, assuming that the effective radial potential is characterized by a short-range infinitely repulsive core of radius R, we derive a compact analytical formula for the threshold energy Elmax=Elmax(n,βn,R)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E^{\text {max}}_l=E^{\text {max}}_l(n,\beta _n,R)$$\end{document}, which characterizes the most weakly bound-state resonance (the most excited energy level) of the quantum system.


Introduction
The Schrödinger differential equation with inverse powerlaw 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][10] and the 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 − [10]. 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 [10]. 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 potena e-mail: shaharhod@gmail.com tial 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][11].
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 ), 1 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). 2

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 3 −h 2 2μ d 2 dr 2 +h 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 4 In addition, the repulsive core of the effective radial potential (3) dictates the inner boundary condition 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 5,6 [12] κr n 1 of small binding energies, where the characteristic length scale r n is defined by the relation 4 We shall assume, without loss of generality, that κ > 0. 5 Note that the strong inequality (7) corresponds to small resonant energies in the regime (2μβ 2/n n /h 2 ) n/(n−2) E 1 [see Eqs. (5) and (8)]. 6 It is important to point out that the special case of s-waves (l = 0) has been studied in the highly important work of Gribakin and Flambaum [12], where it was explicitly proved that the energy of the most excited energy level can be expressed in the simple form κ = 1/(a −ā), where a andā are, respectively, the s-wave scattering length and the so-called average scattering length (see, in particular, Eqs. (3) and (29) of [12]). It is worth emphasizing, however, that the interesting work presented in [12] is based on the semi-classical WKB analysis, whereas in the present paper we shall present a full quantum-mechanical treatment of the physical system. In addition, in the present paper we shall extend the interesting results of [12] to the more generic physical regime of non-spherically symmetric higher partial wave functions (that is, our analytical results, to be derived below, are valid for generic values of the dimensionless physical parameter l).

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) (see footnote 5)], 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 [13]): 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 [13]), 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 [13]): where {B 1 , B 2 } are normalization constants. 7 Using the small-argument (κr 1) asymptotic behaviors of the modified Bessel functions (see Eqs. 9.1.7 and 9.1.9 of [13]), 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 wavenumbers 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 ) 7 As we shall explicitly show below, the normalization constants 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 [13], one finds the asymptotic spatial behavior for the radial eigenfunction (17). Taking account of the boundary condition (4), which characterizes the boundstate resonances of the radial Schrödinger equation (2), one deduces from (22) the simple relation 9 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).

The resonant binding energy of the most excited energy level
The dimensionless ratio A 2 /A 1 , which 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 10 for the dimensionless resonant wavenumber which characterizes the most excited energy level (the most weakly boundstate 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)] 11 κr n 1 (27) 10 As interestingly pointed out by the anonymous referee, following [12], which determined the most excited energy level of the system in the particular case of spherically symmetric s-waves (l = 0), it may be interesting to explore the possibility to express the result (26) for the characteristic wavenumber κ of the most excited energy level (a result which is valid for generic values of the dimensionless angular parameter l) in terms of the characteristic s-wave (l = 0) scattering length a, or possibly in terms of higher partial wave scattering volumes. The exploration of this interesting possibility is beyond the scope of the present paper. 11 Since each inequality in (14) and (19) roughly corresponds to an order-of-magnitude difference between two physical quantities [that is, r n /r 10 −1 and r/(1/κ) 10 −1 in (14) and (19)], the analytically derived resonance condition (24) for the weakly bound-state resonances (the most excited energy levels) of the effective radial potential (3) is expected to be valid in the low wavenumber regime κr n 10 −2 . of small binding energies. Taking account of Eq. (26), one realizes that the small wavenumber requirement (27) is satisfied for 2 r n R where { j ν,k } k=∞ k=1 are the positive zeros of the Bessel function J ν (x) [13,14]. Defining the dimensionless small quantity one finds from (26) the expression 12,13 for the smallest resonant wavenumber which characterizes the effective binding potential (3).

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][10].
for the dimensionless threshold energy 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 noting that in the regime n 1 14 of fast decaying inverse power-law potentials, one finds from (31) the compact formula 15,16 for the characteristic threshold energy of the most excited bound-state resonance.
is worth emphasizing that the interesting work presented in [12] 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.