Two sets of potential-density basis pairs for the study of radial perturbations in collisionless spherical stellar systems

Clicks: 9
ID: 328767
2026
Article Quality & Performance Metrics
Overall Quality
Not rated
Combines reader engagement with the AI quality analysis. This article has not been analysed, so there is no overall score — reader engagement is measured and shown alongside.
AI Quality Assessment
Not analyzed
Readership in this journal
Emerging

Ranked #1,275 of 1,339 articles by views in monthly notices of the royal astronomical society

Most read Least read

Bar heights use a square-root scale. Only the 120 most-read articles are drawn; the journal has 1,339 in total.

Mint this article as an NFT
Not yet minted

Create a permanent, verifiable on-chain record of this article on the Scimatic Network. The NFT is held in your Journament account, and you can withdraw it to your own wallet at any time.

5 SUSD one-off · no wallet required
Abstract
Abstract The Kalnajs matrix method is a widely used framework for studying the global linear stability and possible Landau damping of collisionless stellar systems. However, for radial perturbations (l = 0) in three-dimensional spherical models with infinite boundaries, standard biorthogonal sets such as the Clutton-Brock basis often converge slowly. This stems from a physical constraint: mass-conserving radial modes force the perturbed potential to decay faster than a point mass at infinity, whereas individual Clutton-Brock elements carry a fictitious net mass and decay only as $\mathcal {O}(1/r)$, producing unphysical asymptotic tails. We construct two new families of potential-density basis pairs that are free of these tails by construction. The first modifies the Clutton-Brock set: a specific linear combination of adjacent elements analytically cancels the leading $\mathcal {O}(1/r)$ term, giving potentials that decay as $\mathcal {O}(1/r^3)$ and densities as $\mathcal {O}(1/r^5)$. This breaks strict biorthogonality but yields a compact tridiagonal Gram matrix that enters the response equation at negligible additional cost. The second family is built from Jacobi polynomials that embed the required $\mathcal {O}(1/r^2)$ potential decay directly into their construction while retaining strict diagonal biorthogonality. Numerical tests demonstrate convergence that is uniform in radius for both expansions: the truncation error decreases exponentially when the asymptotic tail of the expanded potential matches the parity of the basis elements, and algebraically otherwise. Both bases reduce the dimension of the response matrix required for a given accuracy and are suited to studying radial perturbations in open stellar systems.
Reference Key
openalex_W7208730293 Use this key to autocite in the manuscript while using SciMatic Manuscript Manager or Thesis Manager
Authors
Journal monthly notices of the royal astronomical society
Year 2026
DOI
10.1093/mnras/stag1717
URL
Keywords Keywords not found

Citations

No citations found. To add a citation, contact the admin at info@scimatic.org

No comments yet. Be the first to comment on this article.