Representation of solutions to the one-dimensional Schrödinger equation in terms of Neumann series of Bessel functions
This provides a novel numerical method for solving spectral problems in quantum mechanics and related fields, offering improved accuracy for large eigenvalue sets.
The paper presents a new representation of solutions to the one-dimensional Schrödinger equation as Neumann series of Bessel functions, enabling uniform approximation in the spectral parameter and accurate computation of large sets of eigendata with nondeteriorating accuracy.
A new representation of solutions to the equation $-y"+q(x)y=ω^2 y$ is obtained. For every $x$ the solution is represented as a Neumann series of Bessel functions depending on the spectral parameter $ω$. Due to the fact that the representation is obtained using the corresponding transmutation operator, a partial sum of the series approximates the solution uniformly with respect to $ω$ which makes it especially convenient for the approximate solution of spectral problems. The numerical method based on the proposed approach allows one to compute large sets of eigendata with a nondeteriorating accuracy.