🔭 astrophysics

Near-degeneracy effects in Quadrupolar Mixed Modes. From an Asymptotic Description to Data Fitting

This paper presents a new asymptotic formulation and Bayesian fitting framework for quadrupolar (l=2) mixed modes that accounts for near-degeneracy effects, enabling precise, model-independent measurements of internal rotation and asymptotic parameters in evolved solar-like stars.

Original authors: B. Liagre, A. Desai, L. Einramhof, L. Bugnet

Published 2026-03-18

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: B. Liagre, A. Desai, L. Einramhof, L. Bugnet

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Listening to a Star's Heartbeat

Imagine a star is like a giant, glowing drum. When it vibrates, it creates sound waves that travel through its interior. By listening to these "songs" (a field called asteroseismology), astronomers can figure out what's happening inside the star, just like a doctor uses an ultrasound to see inside a human body.

For a long time, astronomers have been very good at listening to the "low notes" of these stars (called dipolar modes or $\ell=1$ ). These notes have told us that the cores of old stars spin much faster than their outer layers, but not as fast as our computer models predicted. This suggests there is some invisible "brake" or "mixer" inside the star that we don't fully understand yet.

The Problem: The Quiet, Messy High Notes

The star also sings "high notes" (called quadrupolar modes or $\ell=2$ ). These high notes are actually even more valuable because they can tell us about different parts of the star's interior. However, they are very hard to hear for two reasons:

They are whisper-quiet: They are much fainter than the low notes, so they get lost in the cosmic static.
They are "glitchy": When the star spins, these high notes usually split into a neat group of frequencies. But because the star's core and envelope spin at different speeds, these groups of notes start to overlap and interfere with each other. It's like two singers trying to harmonize, but one is slightly out of tune, creating a messy, distorted sound that is hard to decipher.

This distortion is called a "near-degeneracy effect." Until now, trying to analyze these messy high notes required building a massive, complex 3D computer model of the specific star, which is slow and difficult.

The Solution: A New "Decoder Ring"

The authors of this paper (Liagre et al.) have invented a new mathematical "decoder ring." Instead of building a whole new computer model for every star, they developed a set of analytical formulas (like a recipe) that can predict exactly how these messy, overlapping high notes should look.

The Analogy:
Imagine you are trying to identify a specific car engine sound in a noisy parking lot.

The Old Way: You build a full-scale simulation of the car, the weather, and the pavement to figure out what the engine should sound like, then compare it to the recording.
The New Way: The authors wrote a formula that says, "If the engine is spinning at speed X and the gears are interacting like Y, the sound will look like this specific pattern." You can plug the numbers in, and it tells you exactly what to expect, without needing the full simulation.

What They Did

Derived the Math: They took the complex physics of how stars spin and derived a new formula that accounts for the "messy overlap" (near-degeneracy) of the high notes.
Built a Fitting Tool: They put this formula into a powerful computer program that looks at the star's entire song at once (a "global fit"). It doesn't just pick out one note; it tries to fit the whole melody to find the best match.
Tested on Two Stars:
- KIC 7341231: They tested their new method on a star that had already been studied using the old, heavy computer-modeling method. Their new, faster method gave the exact same results for the star's rotation speed, but with much higher precision (smaller error bars).
- KIC 8179973: They applied the method to a brand-new star that had never been analyzed this way before. They successfully measured its rotation speed and internal structure for the first time.

Why This Matters

Speed and Precision: This new method is much faster and gives more precise measurements of how fast the star's core and surface are spinning.
Unlocking New Data: Because they can now handle the "messy" high notes, astronomers can use all the data from space telescopes (like Kepler and TESS), not just the easy-to-read low notes.
Solving the Mystery: By getting better measurements of how stars spin, we can finally figure out what that invisible "brake" is inside them. Is it magnetic fields? Waves? Turbulence? This new tool brings us one step closer to solving that mystery.

In a Nutshell

The authors found a way to listen to the "whispering, glitchy" high notes of dying stars. They created a new mathematical tool that untangles these messy sounds without needing heavy computer simulations. This allows them to measure the spin of stars' cores with incredible accuracy, helping us understand the hidden physics that governs how stars evolve.

1. Problem Statement

The internal dynamics of evolved solar-like stars (subgiants and red giants) are primarily probed using asteroseismology, specifically through mixed modes (coupled acoustic $p$ -modes and gravity $g$ -modes). While dipolar ( $\ell=1$ ) mixed modes have been extensively used to measure core and envelope rotation rates, they reveal rotation gradients that are weaker than predicted by hydrodynamical models, suggesting missing angular momentum transport mechanisms.

Extending this analysis to quadrupolar ( $\ell=2$ ) mixed modes is crucial because they sample different regions of the stellar interior (particularly the convective zone) and provide additional constraints. However, $\ell=2$ modes are rarely characterized due to two main challenges:

Low Amplitude: The coupling between $p$ and $g$ cavities is weaker for $\ell=2$ than $\ell=1$ , resulting in lower signal-to-noise ratios in the power spectral density (PSD).
Near-Degeneracy Effects: In rotating stars, adjacent rotational multiplets can become close in frequency. When the frequency separation is comparable to the rotation rate, near-degeneracy occurs. This breaks the symmetry of rotational splitting predicted by first-order perturbation theory, creating complex, asymmetric patterns that are difficult to identify and interpret using standard methods.

Existing frameworks either treat these effects implicitly (requiring numerical stellar models) or calculate them in a decoupled basis ( $\pi$ - $\gamma$ ), making direct application to observational data fitting difficult.

2. Methodology

The authors developed a comprehensive framework to model and fit $\ell=2$ mixed modes explicitly accounting for near-degeneracy effects without relying on numerical stellar models.

A. Theoretical Formalism

Asymptotic Derivation: The authors derived a new asymptotic formulation for the off-diagonal terms of the rotation matrix in the mixed-mode basis.
Near-Degeneracy Terms: They explicitly calculated the interaction terms ( $\gamma_{c,ij}$ and $\gamma_{e,ij}$ ) between modes of adjacent radial orders using the JWKB approximation for horizontal eigenfunctions in the $g$ -mode cavity.
Closed-Form Solution: The resulting expression (Eq. 27) provides an analytical estimate of the near-degeneracy coupling using standard asymptotic parameters ( $\Delta\Pi$ , $\zeta$ , frequencies), allowing for the solution of the full quadratic eigenvalue problem directly in the mixed-mode basis.
Two-Zone Rotation Model: The formalism assumes a simplified rotation profile with a core ( $\Omega_{core}$ ) and an envelope ( $\Omega_{env}$ ), calculating diagonal (symmetric splitting) and off-diagonal (asymmetric distortion) terms.

B. Data Analysis Framework

Global Bayesian Fitting: Instead of "peakbagging" (identifying individual peaks first), the authors implemented a global fit of the entire Power Spectral Density (PSD) using a Monte Carlo Markov Chain (MCMC) approach.
Synthetic PSD Generation: They used the sloscillations code to generate synthetic spectra containing:
- Pure $p$ -modes and $g$ -modes based on asymptotic relations.
- Mixed modes coupled via the factor $q_\ell$ .
- Rotational splitting (symmetric for $\ell=1$ , full quadratic solution for $\ell=2$ ).
- Realistic mode heights, linewidths, and granulation backgrounds.
Priors and Constraints: The framework uses priors derived from $\ell=1$ fits to constrain $\ell=2$ parameters (e.g., relating $\Delta\Pi_2$ to $\Delta\Pi_1$ and deriving bounds for the coupling factor $q_2$ based on inclination).

3. Key Contributions

First Asymptotic Description of $\ell=2$ Near-Degeneracy: The paper presents the first analytical, asymptotic formulation of off-diagonal near-degeneracy terms specifically for mixed quadrupolar modes in the mixed-mode basis.
Model-Independent Fitting: The method allows for the direct fitting of $\ell=2$ modes without requiring a pre-computed numerical stellar model to compute coupling kernels.
Global Bayesian Pipeline: A robust pipeline capable of simultaneously fitting $\ell=0, 1, 2$ modes, distinguishing between signal and noise, and resolving overlapping multiplets.
New Constraints on Coupling: The authors provide methods to constrain the $\ell=2$ coupling factor ( $q_2$ ) and phase ( $\epsilon_{g,2}$ ) using observational data.

4. Results

The framework was applied to two Kepler targets:

KIC 7341231 (Validation Target)

Comparison: The results were compared against a previous study by Deheuvels et al. (2017) that used numerical modeling.
Agreement: The asymptotic fit yielded core and envelope rotation rates ( $\Omega_{core} \approx 781$ nHz, $\Omega_{env} \approx 53$ nHz) consistent with the numerical model but with significantly higher precision (uncertainties reduced by factors of ~2–3).
Statistical Evidence: Including mixed $\ell=2$ modes reduced the Bayesian Information Criterion (BIC) by 3,726 points compared to a model using only $\ell=1$ mixed modes and $\ell=2$ pure $p$ -modes, providing overwhelming statistical evidence for the extended model.
Identification: The method successfully identified $\ell=2$ modes across four radial orders and resolved overlapping modes of different degrees ( $\ell=0, 1, 2$ ) that standard methods struggle to separate.

KIC 8179973 (New Discovery)

First Characterization: This star was analyzed for the first time, revealing strong near-degeneracy effects in its quadrupolar modes.
Parameters: The authors measured the first set of asymptotic parameters for this star, including $\Delta\Pi_2$ , $q_2$ , and $\epsilon_{g,2}$ .
Rotation: Derived rotation rates of $\Omega_{core} \approx 672$ nHz and $\Omega_{env} \approx 49$ nHz, consistent with young red giant expectations.
Robustness: The fit remained robust even with a less informative prior on $q_2$ , demonstrating the method's reliability.

5. Significance and Future Outlook

Precision Asteroseismology: The global fitting approach significantly reduces uncertainties in rotation rate estimates compared to standard peakbagging methods, offering better constraints for angular momentum transport studies.
Beyond $\ell=1$ : This work paves the way for using $\ell=2$ modes as standard seismic probes, effectively doubling the number of available oscillation probes and providing sensitivity to different dynamical kernels.
Magnetic Field Detection: The framework sets the stage for disentangling rotation-induced asymmetries from magnetic signatures in $\ell=2$ modes. Since magnetic fields can suppress mode amplitudes and alter period spacings, having a precise asymptotic model for $\ell=2$ is essential for future magnetic field topology studies in red giants.
Mission Readiness: The methodology is designed to handle data from current (Kepler, TESS) and future (PLATO) missions, enabling large-scale, homogeneous seismic analyses of subgiants and early red giants.

In summary, this paper bridges the gap between theoretical asymptotic descriptions and observational data analysis, providing a powerful, model-independent tool to unlock the internal dynamics of stars using quadrupolar mixed modes.