A generalized linear matrix method for normal modes in collisionless stellar disks
摘要
We generalize the linear matrix method for computing normal modes in collisionless stellar disks to distribution functions with sharp edges at zero angular momentum (). The generalization adds boundary-integral terms to the matrix equation without increasing its size. We validate the method by computing modes for two Kuzmin--Toomre disk models (Miyamoto and Kalnajs families) and comparing the eigenvalues with those obtained from an independent nonlinear matrix method based on logarithmic-spiral expansions. A systematic convergence study over grid resolution and harmonic truncation yields eigenvalues accurate to in both pattern speed and growth rate. Unlike the nonlinear method, the linear method naturally incorporates gravitational softening, enabling the computation of eigenmodes for softened disk models. The implementation in Julia with GPU acceleration is openly available.
引用
@article{arxiv.2606.31523,
title = {A generalized linear matrix method for normal modes in collisionless stellar disks},
author = {Evgeny V. Polyachenko and Ilia G. Shukhman},
journal= {arXiv preprint arXiv:2606.31523},
year = {2026}
}
备注
13 pages, 3 figures, 3 tables. Accepted to ApJ