A radial basis of massless potential-density pairs
In the authors' words
In the matrix method for linear perturbations of spherical stellar systems, the perturbed potential and density are expanded over a set of potential-density pairs. For radial perturbations, mass conservation makes the coefficient of in the perturbed potential vanish, so that the potential decays faster than . Starting from the Hernquist--Ostriker family, we construct in closed form a set of potential-density pairs whose potentials decay as with a prescribed , the tail containing all subsequent integer powers. For each pair is individually massless, and the potential of the -th pair is expressed through a single Jacobi polynomial. We derive the combination weights, the leading tail coefficients and the Gram matrix analytically; the Gram matrix is banded with half-width (tridiagonal at ). The -th potential element has exactly nodes, spread from to , so that the set resolves both the centre and the far periphery. As a test, the expansion of the dilation-mode potential of the isochrone model converges exponentially and carries no parasitic mass at any truncation, in contrast to the standard set.
Appeared: Wednesday, September 23. arXiv. Preprint, not yet peer-reviewed.
Authors' comment: 6 pages, 3 figures. Submitted to Astron.Rep