ABSTRACT
We discuss the application of the discrete variable representation (DVR) to Schrödinger problems which involve singular Hamiltonians. Unlike recent authors who invoke transformations to rid the eigenvalue equation of singularities at the cost of added complexity, we show that an approach based solely on an orthogonal polynomial basis is adequate, provided the Gauss-Lobatto or Gauss-Radau quadrature rule is used. This ensures that the mesh contains the singular points and by simply discarding the DVR functions corresponding to those points, all matrix elements become well behaved, the boundary conditions are satisfied, and the calculation is rapidly convergent. The accuracy of the method is demonstrated by applying it to the hydrogen atom. We emphasize that the method is equally capable of describing bound states and continuum solutions.
