中文

A generalized linear matrix method for normal modes in collisionless stellar disks

星系天体物理 2026-06-30 v1

摘要

We generalize the linear matrix method for computing normal modes in collisionless stellar disks to distribution functions with sharp edges at zero angular momentum (L=0L=0). The generalization adds boundary-integral terms to the matrix equation without increasing its size. We validate the method by computing m=2m=2 modes for two Kuzmin--Toomre disk models (Miyamoto nM=3n_{\rm M}=3 and Kalnajs mK=6m_{\rm K}=6 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 0.003{\sim}\,0.003 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