 Solution

Take a static, nonrotating, nonmagnetic spherical equilibrium with $\nabla P=-\rho \nabla \Phi$ and ${\nabla }^2\Phi =4\pi G\rho$. Use a fluid displacement $\mathbf{\xi }(\mathbf{r})e^{-i\omega t}$. Write $p_1,{\rho }_1,{\varphi }_1$ for Eulerian fluid perturbations. The Lagrangian pressure perturbation and corresponding mass density change obey $\Delta P=p_1+\mathbf{\xi }\cdot \nabla P$, $\Delta \rho ={\rho }_1+\mathbf{\xi }\cdot \nabla \rho$. Linearizing mass conservation, the Euler equations for an inviscid fluid and Poisson equation for Newtonian gravity givesThere is no equilibrium acceleration to multiply a perturbed mass density. The adiabatic process condition, with composition carried by the parcel, closes the system:These are the complete linear adiabatic stellar oscillation equations, including the perturbation of self-gravity. Neglect of heat exchange is appropriate when oscillation periods are short compared with relevant thermal relaxation times; it does not determine nonadiabatic excitation or damping.

For the pressure and buoyancy modes, separate angular dependence using spherical harmonics:The horizontal amplitude ${\xi }_h$ has dimensions of length. Angular differentiation gives $\nabla \cdot \mathbf{\xi }=[r^{-2}(r^2{\xi }_r{)}^{\prime }-\ell (\ell +1){\xi }_h/r]Y_{\ell }^m$, and tangential momentum gives ${\omega }^{2}r{\xi }_h=p_{\ell }/\rho +{\varphi }_{\ell }$. Define the adiabatic sound speed, stellar buoyancy frequency and Lamb frequency byThe adiabatic mass density relation becomes ${\rho }_{\ell }/\rho =p_{\ell }/(\rho c_s^2)+(N^2/g){\xi }_r$. Substitution produces a radial form of the full oscillation equations for nonzero $\omega$:The apparent $N^2/g$ factor is evaluated through its defining gradient at the centre rather than by dividing two zeros.

Regularity at the centre excludes singular solutions. At a free surface the Lagrangian pressure perturbation vanishes, $\Delta P=0$; outside the star the gravitational perturbation decays as $r^{-\ell -1}$. For a model with mass density tending to zero at its surface, continuity of ${\varphi }_{\ell }$ and its radial derivative gives ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R$. If the equilibrium mass density jumps to vacuum, include the displaced-surface mass sheet: the outward-minus-inward derivative jump is $4\pi G\rho (R){\xi }_r(R)$, so the interior condition is ${\varphi }_{\ell }^{\prime }(R)=-(\ell +1){\varphi }_{\ell }(R)/R-4\pi G\rho (R){\xi }_r(R)$. An atmospheric boundary condition can replace the ideal free surface.

These conditions make ${\omega }^2$ an eigenvalue, not an arbitrary local sound frequency. With conservative boundary conditions the adiabatic operator is self-adjoint, giving real ${\omega }^2$; negative values describe instability. The frequencies depend on $P(r),\rho (r),{\Gamma }_1(r)$, self-gravity, stratification and boundaries, as well as angular degree $\ell$ and radial order. In a spherical nonrotating star they are degenerate in $m$. For homologous equilibrium structures, their scale isthus the typical oscillation time measures inverse square root of mean mass density, while individual frequencies probe the interior adiabatic sound speed and stellar buoyancy frequency profiles.

For radial modes, write ${\xi }_r=r\eta$. Eliminating the pressure and gravitational perturbations gives the radial stellar pulsation equationThis is a Sturm-Liouville problem. Multiplication by $\eta$ and integration, with vanishing boundary terms, gives its Rayleigh quotientFor constant ${\Gamma }_1>4/3$, $P^{\prime }<0$ makes both numerator contributions nonnegative. At ${\Gamma }_1=4/3$, a homologous displacement is neutral; for constant ${\Gamma }_1<4/3$ the same trial displacement makes the quotient negative. For the uniform-density stellar model, $\eta =$ constant is an exact mode and ${\omega }^2=(3{\Gamma }_1-4)GM/R^3$. With varying ${\Gamma }_1$, the integral criterion, rather than a universal pointwise threshold, controls radial stability.

The Cowling approximation neglects ${\varphi }_{\ell }$ while retaining the equilibrium gravitational field. It is useful for short-wavelength modes, but is not needed for the full derivation above. In a locally slowly varying region, take both remaining amplitudes proportional to $e^{i\int k_{r}dr}$ and retain the leading derivative terms. ThenEliminating either amplitude gives the acoustic-gravity propagation relationPositive $k_r^2$ is oscillatory propagation; negative $k_r^2$ means an evanescent wave. The high-frequency branch, ${\omega }^2>S_{\ell }^2,N^2$, describes stellar acoustic modes, restored chiefly by compressibility and pressure. At frequencies well above $N$, this gives ${\omega }^2\simeq c_s^2[k_r^2+\ell (\ell +1)/r^2]$. The low-frequency propagating branch in stable stratification, ${\omega }^2<N^2,S_{\ell }^2$, instead describes stellar gravity modes, restored by buoyancy. At $\ell =0$ there is no such nonradial gravity-wave cavity.

A stellar acoustic mode is trapped between an inner turning point near $\omega =S_{\ell }$ and an outer reflecting region. Low-$\ell$ modes penetrate deeply; radial modes reach the centre. Higher-degree modes turn farther out. Standing waves require the WKB quantization conditionwhere the phase $\alpha$ depends on the central or turning-point behaviour and surface reflection. For high radial order and small degree, the leading acoustic travel-time result is the large frequency separationThe $\ell /2$ term is the leading central angular phase shift; smaller frequency separations depend on detailed interior gradients. The near-surface mass density stratification sets the acoustic cutoff frequency. In a plane-parallel isothermal atmosphere with mass density scale height $H_{\rho }$, ${\omega }_{\mathrm{ac}}\simeq c_s/(2H_{\rho })$. Modes below this cutoff can reflect and form a cavity; waves sufficiently above it escape and need an outgoing-wave boundary condition. Nonadiabatic damping, driving and rotation alter real-star mode properties, but the adiabatic frequency problem isolates their dependence on the equilibrium structure.