Abstract:In this paper we consider the Sturm-Liouville equation -y"+qy = lambda*y on the half line (0,infinity) under the assumptions that x=0 is a regular singular point and nonoscillatory for all real lambda, and that either (i) q is L_1 near x=infinity, or (ii) q' is L_1 near infinity with q(x) --> 0 as x --> infinity, so that there is absolutely continuous spectrum in (0,infinity). Characterizations of the spectral density function for this doubly singular problem, similar to those obtained in [12] and [13] (when the left endpoint is regular) are established; corresponding approximants from the two algorithms in [12] and [13] are then utilized, along with the Frobenius recurrence relations and piecewise trigonometric - hyperbolic splines, to generate numerical approximations to the spectral density function associated with the doubly singular problem on (0,infinity). In the case of the radial part of the separated hydrogen atom problem, the new algorithms are capable of achieving near machine precision accuracy over the range of lambda from 0.1 to 10000, accuracies which could not be achieved using the SLEDGE software package.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to estimate the spectral density function for the Sturm - Liouville problem with two singular endpoints. Specifically, the authors consider the Sturm - Liouville equation on the half - line $(0,\infty)$:
\[
-y'' + q(x)y=\lambda y
\]
where $x = 0$ is a regular singular point and is non - oscillatory for all real $\lambda$, and the assumptions at $x=\infty$ include two cases: (i) $q$ belongs to the $L^1$ space near $x = \infty$; (ii) $q'$ belongs to the $L^1$ space near $x=\infty$ and $q(x)\to0$ as $x\to\infty$, thus ensuring the existence of an absolutely continuous spectrum in the interval $(0,\infty)$.
The main objectives of the paper are two - fold:
1. Under appropriate assumptions, extend the analysis in [12] and [13] to show that the spectral density function $f(\lambda)=\rho'(\lambda)$ associated with equation (1.2) can be expressed as:
\[
f(\lambda)=\frac{1}{\pi}\left[P(x,\lambda)\phi(x,\lambda)^2+Q(x,\lambda)\phi(x,\lambda)\phi'(x,\lambda)+R(x,\lambda)\phi'(x,\lambda)^2\right]
\]
where $(P(\cdot,\lambda),Q(\cdot,\lambda),R(\cdot,\lambda))^T$ is the unique solution of the initial - value problem at $x = \infty$.
2. Extend the numerical algorithms in [12] and [13] to calculate the spectral density function $f(\lambda)$. Approximate the solution $\phi(x,\lambda)$ and its derivative by using the Frobenius power series solution, and perform numerical calculations using piecewise trigonometric - hyperbolic splines.
In addition, the paper also presents the application of the new numerical algorithm on several test examples, including the problem of separating the radial part of the hydrogen atom, demonstrating the superiority of the new algorithm in terms of accuracy.