The Contour Method: a New Approach to Finding Modes of Nonadiabatic Stellar Pulsations
Abstract
The contour method is a new approach to calculating the nonadiabatic pulsation frequencies of stars. These frequencies can be found by solving for the complex roots of a characteristic equation constructed from the linear nonadiabatic stellar pulsation equations. A complexroot solver requires an initial trial frequency for each nonadiabatic root. A standard method for obtaining initial trial frequencies is to use a star's adiabatic pulsation frequencies, but this method can fail to converge to nonadiabatic roots, especially as the growth and/or damping rate of the pulsations becomes large. The contour method provides an alternative way to obtain initial trial frequencies that robustly converges to nonadiabatic roots, even for stellar models with extremely nonadiabatic pulsations and thus high growth/damping rates. We describe the contour method implemented in the GYRE stellar pulsation code and use it to calculate the nonadiabatic pulsation frequencies of $10\,{M}_{\odot }$ and $20\,{M}_{\odot }$ β Cephei star models, and of a $0.9\,{M}_{\odot }$ extreme helium star model.
 Publication:

The Astrophysical Journal
 Pub Date:
 August 2020
 DOI:
 10.3847/15384357/aba748
 arXiv:
 arXiv:2006.13223
 Bibcode:
 2020ApJ...899..116G
 Keywords:

 Asteroseismology;
 Stellar oscillations;
 Computational methods;
 Astronomy software;
 73;
 1617;
 1965;
 1855;
 Astrophysics  Solar and Stellar Astrophysics
 EPrint:
 12 pages, 14 figures, Accepted to ApJ