Non-Adiabatic Pulsation Analyses in Stars

• Non-adiabatic pulsation analyses are a framework that incorporates entropy perturbations to model energy exchanges and predict mode stability in stars.
• They elucidate the roles of opacity (κ–γ) and nuclear (ε) mechanisms in driving and damping oscillations, accurately defining instability strips.
• Computational methods using linearized stellar models validate the theory against observed pulsations in low-mass white dwarfs and classical variable stars.

Non-Adiabatic Pulsation Analyses

Non-adiabatic pulsation analyses are a cornerstone of modern stellar astrophysics, providing the theoretical and computational framework for understanding the driving, damping, stability, amplitude growth, and observable characteristics of stellar oscillations in contexts where the adiabatic approximation is insufficient. By explicitly including the energy exchanges between oscillations and the stellar medium—incorporating radiative, convective, and nuclear energy flows—these analyses enable the prediction of instability strips, mode amplitudes, and the interpretation of rich oscillation spectra observed in variable stars, especially low-mass white dwarfs, classical pulsators, and convective envelope objects. The scope of non-adiabatic theory includes both self-excited pulsators (classical variables) and stochastically excited, solar-like oscillations, with the underlying physics varying according to region in the HR diagram and the local properties of the stellar model.

1. Physical Foundation and Governing Equations

Non-adiabatic theory departs from the adiabatic limit by explicitly retaining all terms associated with variations of entropy, allowing for net energy exchange between a pulsating mode and the stellar environment. The linearized set of equations encompasses mass conservation, momentum, Poisson’s equation (gravity), the perturbed equation of state, and critically, the full linearized energy equation with entropy perturbations ($\delta s \neq 0$): $$T \frac{\partial \delta s}{\partial t} = \delta \left( \epsilon - \frac{\partial L}{\partial m} \right)$$ where $\epsilon$ is nuclear energy generation per unit mass and $L$ is the luminosity at mass coordinate $m$. The solutions yield complex eigenfrequencies, $\omega = \omega_r + i\omega_i$, with $\omega_i$ quantifying mode growth ($\omega_i > 0$) or damping ($\omega_i < 0$). Perturbative terms also incorporate opacity derivatives with respect to temperature and density: $$\frac{\delta \kappa}{\kappa} = \kappa_T \frac{\delta T}{T} + \kappa_\rho \frac{\delta \rho}{\rho}$$ This formalism is necessary to analyze both the local processes (e.g., κ- and γ-mechanisms) and global stability properties (e.g., work integral balance) of oscillatory modes.

2. Mode Excitation and Damping Mechanisms

Two principal classes of mechanisms emerge in non-adiabatic analyses:

• Opacity (κ–γ) Mechanism: Operating most efficiently in partial ionization zones (e.g., hydrogen and helium in low-mass white dwarfs, iron-group “Z-bump” in hot pulsators), this mechanism relies on localized enhancements of opacity during compression, resulting in heat trapping and periodic driving when the local thermal time matches the modal period. The mathematical criterion for driving involves the quantity:

$\kappa_T + \frac{\kappa_\rho}{\Gamma_3 - 1}$

Here, mode driving requires this sum to increase outward, as confirmed in He-core white dwarfs and BLAPs.

• Epsilon (ε) Mechanism: In models with ongoing nuclear burning (e.g., ELM white dwarfs with thick hydrogen envelopes), the temperature sensitivity of energy generation can excite low-order g-modes if significant temperature perturbations overlap with the burning region.

Damping is primarily associated with radiative diffusion in deep layers, the time-dependent modulation of convective flux (requiring time-dependent convection models), and turbulent viscosity, with the net effect quantified by the imaginary component of mode eigenfrequencies or by the sign of the work integral at the stellar surface.

3. Computational Methods and Observational Validation

Non-adiabatic mode stability calculations require detailed input models, typically derived from evolutionary tracks including binary evolution, mass loss, and diffusion, as in the construction of He-core WD sequences via binary interaction. Pulsation properties (radial and nonradial, low- and high-order modes) are computed using established codes with the frozen convection approximation (neglecting perturbation of convective flux), which, for white dwarfs, reliably predicts observed instability strip boundaries.

Key diagnostic quantities include the growth rate parameter ($\eta = -\frac{\Im(\sigma)}{\Re(\sigma)}$), the $e$-folding time ($\tau_e = 1/|\Im(\sigma)|$), differential and integrated work functions ($dW/dr$, $W$), and ranges of unstable mode periods, each sensitive to stellar mass, $T_{\rm eff}$, convective efficiency (as parameterized by MLT flavor: ML1, ML2, ML3), and mode degree $\ell$. Contemporary analyses yield excellent agreement between theoretical instability domains (in both period and temperature) and observed periods for ELMVs, including the role of the ε-mechanism in stars observed to have short-period modes indicative of residual nuclear burning.

4. Parameter Dependencies and Mode Instability Domains

Non-adiabatic instability domains are acutely sensitive to several key stellar parameters:

• Stellar Mass: Lower-mass (especially $M_*\lesssim 0.18\,M_\odot$) He-core white dwarfs present more extended unstable $g$-mode period ranges and are susceptible to ε-driven low-order mode excitation.
• Effective Temperature ($T_{\rm eff}$) and Convection Efficiency: The blue edge of the instability strip shifts to higher temperature for higher mass and more efficient convection. The period range of unstable modes broadens toward lower $T_{\rm eff}$. Different MLT parameterizations yield blue edge differences up to 1300 K.
• Mode Degree ($\ell$): The blue edge shifts slightly for $g$-modes with increasing $\ell$, but instability domains for $p$- and radial modes are largely insensitive to $\ell$.
• H-Envelope Mass and Evolutionary History: The mass of the hydrogen envelope, determined by evolutionary pathway (binary case vs. single evolution), governs the presence and strength of ongoing nuclear burning and the excitation of short-period modes via the ε-mechanism.

Observed ELMVs lie within the predicted instability domains, affirming the theoretical models. The location of observed periods relative to predicted blue edges further supports specific convective parameterizations (favoring ML2/ML3 and 3D-corrected atmospheric parameters in some cases).

5. Theoretical and Observational Implications

Non-adiabatic pulsation analyses provide critical insight into the physics of both the driving and selection of observable modes, the structure of the instability strips, and the temporal evolution of mode amplitudes with respect to stellar cooling timescales. In the context of ELM white dwarfs, these analyses establish:

• Predominant κ–γ excitation of a dense spectrum of radial and non-radial modes due to partial H-ionization.
• Supplemental excitation of low-order $g$-modes by the ε-mechanism in models with thick H envelopes, offering indirect evidence for ongoing H burning in the coolest observed WDs.
• Consistency of observable $g$-mode periods and their temperature dependence with theoretical predictions, and the ability to explain short-period $g$-modes via nuclear driving in certain targets.

The sensitivity to convective efficiency underscores the importance of accurate convective parameterization in evolutionary and atmospheric models.

6. Methodological Advances and Limitations

While standard non-adiabatic analyses for white dwarfs employ the frozen convection approximation, capturing the blue edge with high accuracy, parameterization of time-dependent convection remains an open area, especially for stars with substantial convective envelopes or in contexts where amplitude saturation or non-linear effects are significant. Instability domain predictions are robust to first order but detailed amplitude and modal selection (including red edge definition and mode visibility) are still affected by the specifics of convective physics and limitations of linear theory. For variable ELMVs, all analytically excited modes have $e$-folding timescales much shorter than evolutionary cooling, implying rapid growth to observable amplitudes.

---

Summary Table: Mode Instability Domains in Low-Mass He-Core White Dwarfs

Parameter	Influence on Instability Domain
Mass ($M_*$)	Lower mass: wider $g$-mode range, stronger ε-driving
$T_{\rm eff}$	Higher mass/MLT = hotter blue edge; range broadens at low $T_{\rm eff}$
Convective Efficiency	ML3 > ML2 > ML1: blue edge shifts hotter and period limits altered
Mode degree ($\ell$)	Slight effect for $g$-mode blue edge; not for $p$-/radial modes
H Envelope Mass	Thicker → more efficient nuclear-driven modes

---

Non-adiabatic pulsation analyses thus form the theoretical basis for asteroseismic inferences about the internal structure, evolution, and atmospheric properties of low-mass white dwarfs and related variable stars. They act as a stringent test of evolutionary theory and of the adequacy of mixing-length convective prescriptions in the context of the observed ELMV instability strip.