Asymptotic analysis of high-frequency acoustic modes in rapidly rotating stars

The asteroseismology of rapidly rotating pulsating stars is hindered by our poor knowledge of the effect of the rotation on the oscillation properties. Here we present an asymptotic analysis of high-frequency acoustic modes in rapidly rotating stars. We study the Hamiltonian dynamics of acoustic rays in uniformly rotating polytropic stars and show that the phase space structure has a mixed character, regions of chaotic trajectories coexisting with stable structures like island chains or invariant tori. In order to interpret the ray dynamics in terms of acoustic mode properties, we then use tools and concepts developed in the context of quantum physics. Accordingly, the high-frequency acoustic spectrum is a superposition of frequency subsets associated with dynamically independent phase space regions. The sub-spectra associated with stable structures are regular and can be modelled through EBK quantization methods while those associated with chaotic regions are irregular but with generic statistical properties. The results of this asymptotic analysis are successfully confronted with the properties of numerically computed high-frequency acoustic modes. The implications for the asteroseismology of rapidly rotating stars are discussed.

💡 Research Summary

The paper tackles the long‑standing problem of interpreting the acoustic oscillation spectra of rapidly rotating stars, where centrifugal distortion and Coriolis forces break the spherical symmetry that underlies most asteroseismic analyses. The authors adopt a high‑frequency (short‑wavelength) approximation and replace the full wave equation by a Hamiltonian ray dynamics description. The stellar model is a uniformly rotating polytrope, which captures the essential structural changes induced by rapid rotation while remaining analytically tractable.

First, the Hamiltonian governing acoustic rays includes the kinetic term, the effective acoustic potential, and rotation‑induced terms (Coriolis and centrifugal). By integrating the ray equations for many initial conditions and constructing Poincaré sections, the authors map the four‑dimensional phase space onto two‑dimensional surfaces of section. The resulting portrait is not uniformly chaotic; instead, it displays a mixed phase‑space structure: stable periodic orbits surrounded by island chains, invariant tori that correspond to quasi‑integrable motion, and extensive chaotic seas where trajectories wander ergodically. The relative volumes of these regions evolve with the rotation rate – higher rotation expands the chaotic sea at the expense of regular islands.

Second, the authors translate this classical phase‑space picture into spectral predictions using semiclassical (quantum‑chaos) techniques. For rays confined to regular islands or tori, the Einstein‑Brillouin‑Keller (EBK) quantization conditions are applied. Each invariant torus carries two action variables; imposing (\oint p_i,dq_i = (n_i+\mu_i/4)h) yields a set of quantized frequencies that form nearly equally spaced ladders. These “regular” modes retain a clear pattern in both spatial eigenfunctions and frequency spacings, despite the underlying stellar deformation.

In contrast, rays that explore the chaotic sea cannot be quantized individually. Instead, the authors invoke statistical descriptions from random‑matrix theory. The spacing distribution of frequencies associated with chaotic trajectories follows the Wigner‑Dyson (GOE) law, reflecting level repulsion and universal fluctuations characteristic of quantum systems whose classical limit is chaotic. This provides a statistical, rather than deterministic, description of the “chaotic” subset of the spectrum.

Third, the theoretical framework is confronted with full numerical solutions of the three‑dimensional acoustic wave equation in the same rotating polytropic model. High‑frequency eigenmodes are computed using spectral methods, and each mode is classified by comparing its spatial structure and frequency to the predictions from ray dynamics. Modes that belong to regular islands exhibit nodal patterns that match EBK‑derived quantum numbers and have frequencies that lie on the predicted ladders with deviations well below the typical mode spacing. Chaotic modes, on the other hand, display irregular nodal geometries and frequency spacings that, when analyzed statistically, conform closely to the GOE predictions. The agreement persists across a range of rotation rates, confirming that the mixed phase‑space picture remains valid as the star approaches breakup velocity.

Finally, the authors discuss the implications for asteroseismology. Traditional mode identification schemes, which rely on spherical harmonic degree (\ell) and azimuthal order (m), break down in the rapid‑rotation regime because the eigenfunctions are no longer separable in latitude and longitude. The mixed‑phase‑space approach offers a new classification: regular modes can be labeled by their EBK quantum numbers (essentially deformed (\ell, m) equivalents), while chaotic modes are described statistically. This dual description enables the extraction of structural information—such as the internal rotation profile, the degree of centrifugal flattening, and the presence of internal barriers—from observed frequency spectra that would otherwise appear too complex. Moreover, the methodology is readily extensible to include differential rotation, magnetic fields, and non‑adiabatic effects, opening a pathway toward a comprehensive theory of oscillations in the most extreme stellar rotators.

In summary, the paper provides a rigorous asymptotic analysis that bridges classical ray dynamics, semiclassical quantization, and random‑matrix statistics to explain the high‑frequency acoustic spectra of rapidly rotating stars. It demonstrates that the spectrum is a superposition of regular, EBK‑quantized subsets and chaotic, statistically described subsets, and validates these predictions against full numerical simulations. This work constitutes a significant step forward in the theoretical foundation required for interpreting observations of rapidly rotating pulsators, such as δ Scuti and β Cephei stars, and paves the way for more accurate asteroseismic inferences in the era of high‑precision space photometry.