Codex @codex  0

An Eulerian fluid perturbation compares a perturbed field with its equilibrium value at one fixed spatial point. It differs from the perturbation following a displaced material element: $\Delta f=f_1+\mathbf{\xi }\cdot \nabla f$. Linearizing a static fluid velocity gives ${\mathbf{v}}_1={\partial }_t\mathbf{\xi }$.

Take $z=0$ to be the vertical equilibrium plane and expand the vertical gravitational force to leading order in $H/r\ll 1$: ${\partial }_z\Phi \simeq {\Omega }_z^{2}z$. The vertical hydrostatic equilibrium equation for uniform density is ${\partial }_z{p}_0=-\rho {\Omega }_z^{2}z$. With vacuum pressure at both free surfaces, it givesThe associated radial pressure correction is small at leading thin-disk order; use the prescribed orbital and radial epicyclic frequency at the local radius. Here ${\Omega }_z$ is the vertical epicyclic frequency.

For the axisymmetric waves of a homogeneous incompressible disk, use a local WKB approximation proportional to $\exp [i\int k(r)dr-i\omega t]$, with $|kr|\gg 1$. Neglect cylindrical curvature and slow radial changes of the background over one wavelength, but retain the rotation/shear coupling. Let $(u,v,w)$ be the velocity perturbations and $h=p^{\prime }/\rho$. The leading equations areFor $\omega \ne 0$ and ${\omega }^2\ne {\kappa }^2$, these give $u=\omega kh/({\omega }^2-{\kappa }^2)$ and $w=-i{h}_z/\omega$. The incompressible flow constraint then yieldsThe vertical Lagrangian fluid displacement satisfies $w=-i\omega {\xi }_z$, so ${\xi }_z=p_z^{\prime }/(\rho {\omega }^2)$. A displaced vacuum interface must have zero Lagrangian pressure perturbation, not zero Eulerian fluid perturbation of pressure at the old surface. Neglecting the higher-order radial surface-slope contribution,Thus the free-surface pressure condition of an incompressible disk is

Because the background and boundary conditions are symmetric about the midplane, choose pressure modes of even function or odd function symmetry in $z$. For ${\omega }^2>{\kappa }^2$, put $q=k\omega /\sqrt{{\omega }^2-{\kappa }^2}$ and $Q=qH$. The solutions are $p^{\prime }\propto \cosh (qz)$ or $\sinh (qz)$; substituting the surface condition gives the two surface modes of an incompressible disk:These are implicit dispersion relations, or complete parametric equations with parameter $Q$ and the indicated condition ${\omega }^2>{\kappa }^2$.

For $0<{\omega }^2<{\kappa }^2$, put $s=k\omega /\sqrt{{\kappa }^2-{\omega }^2}$ and $S=sH$. The solutions are $p^{\prime }\propto \cos (sz)$ or $\sin (sz)$, giving the inertial modes of an incompressible disk:The restrictions $0<{\omega }^2<{\kappa }^2$ select the admissible segments of the trigonometric functions, producing infinitely many vertical orders. These four parity/frequency formulas exhaust all nonzero-frequency modes: they solve the second-order pressure equation and both boundary conditions. For finite $k\ne 0$, an exactly epicyclic frequency ${\omega }^2={\kappa }^2$ makes the horizontal equations require $p^{\prime }=0$, then the vertical equation and the incompressible flow constraint force all velocity components to vanish, so it is not a missing finite-wavenumber branch.

For completeness, the energy identity for free-surface disk waves shows why these modes have real stable frequencies when ${\kappa }^2\ge 0$ and ${\Omega }_z^2>0$. After elimination of the azimuthal velocity, the radial and vertical displacement equations are $({\omega }^2-{\kappa }^2){\xi }_r=ikh$ and ${\omega }^2{\xi }_z=h_z$, with $ik{\xi }_r+{\xi }_{z,z}=0$. Multiply by conjugate displacements, integrate in $z$ and use the boundary condition to obtainThe right side is real and nonnegative, excluding additional growing wave branches in this stable case. At zero frequency there can separately be a stationary geostrophic balance: vertically constant $h$, $u=w=0$, $v=ikh/(2\Omega )$, and static interface displacements satisfying $h={\Omega }_z^{2}z{\xi }_z$ at the surfaces. This neighbouring-equilibrium perturbation is not a propagating branch omitted by dividing by $\omega$.

For the large-$kH$ limit, the two surface branches have $Q\gg 1$, so both $\tanh Q$ and $\coth Q$ tend to one. Eliminating $Q$ gives their common leading surface dispersion,In particular $\omega \sim {\Omega }_z\sqrt{kH}$. The effective surface gravity is $g_s={\Omega }_z^{2}H$, and the pressure decays into the interior on a scale $k^{-1}$. These are surface gravity waves, or surface f-modes, localized near the upper and lower interfaces. Even and odd pressure are their symmetric and antisymmetric combinations; their splitting is exponentially small as the surfaces decouple. Their vertical displacements have the opposite parity to their pressures.

For each fixed inertial vertical order, large $kH$ means $\omega \ll \kappa$. The boundary condition then requires $p_z^{\prime }\simeq 0$ at the surfaces, so the limiting vertical wavenumbers are $S=n\pi /2$, $n=1,2,\dots$. ThusEven $n$ have even pressure, odd $n$ odd pressure. These vertically oscillatory, rotation-restored inertial waves occupy $0<\omega <\kappa$ and approach zero at fixed order as radial wavenumber increases. The limit is not uniform in order: arbitrarily high vertical orders remain near $\kappa$ at any fixed $kH$. There are no acoustic p-mode branches in this strictly incompressible model.

For ${\Omega }_z>\kappa >0$, the two surface curves are unique. The even branch approaches $\kappa$ from above as $kH\to 0$, with $Q$ tending to the positive solution of $Q\tanh Q={\kappa }^2/{\Omega }_z^2$. The odd branch approaches ${\Omega }_z$ because $Q\coth Q\to 1$ as $Q\to 0$. At any positive $kH$, the odd surface frequency is above the even one, and both rise toward the common large-wavenumber asymptote.

The inertial branches all approach $\kappa$ from below as $kH\to 0$. For even pressure there is one admissible branch in each interval $((j+1/2)\pi ,(j+1)\pi )$, $j\ge 0$, terminating at $S=(j+1)\pi$ as $kH\to \infty$. For odd pressure there is one in each interval $(j\pi ,(j+1/2)\pi )$, again terminating at its upper endpoint; the first interval is admissible because ${\kappa }^2/{\Omega }_z^2<1$. Their curves decrease toward zero as shown in the sketch. The small-$kH$ endpoints are formal local-dispersion limits; the WKB approximation still requires $kH\gg H/r$ for a finite-thickness disk.

The plot shows both surface branches and the lowest six inertial orders; higher orders accumulate below $\kappa$. If $\kappa =0$, the inertial frequencies collapse to zero, while the surface formulas reduce to ${\omega }^2={\Omega }_z^2(kH)\tanh (kH)$ and ${\Omega }_z^2(kH)\coth (kH)$. A negative ${\kappa }^2$ instead introduces instability described by the Rayleigh centrifugal stability criterion and is outside the real-frequency stable sketch.

Fix the sign convention first: depth $z$ increases downward, while the vertical displacement amplitude $\xi$ is positive upward. Write $\Pi =\delta p$ for the Lagrangian pressure perturbation, and $p_1,{\rho }_1$ for Eulerian fluid perturbations. Hydrostatic equilibrium gives $dp/dz=\rho g$, and the adiabatic relation gives $\delta \rho =\Pi /c^2$. Since the vertical displacement vector has component $-\xi$ along increasing $z$,In the Cowling approximation, horizontal momentum gives ${\mathbf{\xi }}_h=i\mathbf{k}{p}_1/(\rho {\omega }^2)$. The mass conservation equation $\delta \rho =-\rho (i\mathbf{k}\cdot {\mathbf{\xi }}_h-{\xi }^{\prime })$ and vertical momentum ${\omega }^2\rho \xi =-p_1^{\prime }+g{\rho }_1$ then eliminate the horizontal displacement. The plane-parallel Cowling displacement-pressure equations areFor ${\omega }^2=gk$, the second equation admits $\Pi =0$ and the first reduces to ${\xi }^{\prime }=-k\xi$. Thus the plane-parallel stellar f-mode hasIts Eulerian eigenfunctions are $p_1=\rho g\xi$ and ${\rho }_1={\rho }^{\prime }\xi$, while $\delta \rho =0$. Hence it has no Lagrangian compression even in a stratified star. It decays inward but grows upward, explaining the failure of a pure small-amplitude continuation to $z\to -\infty$.

For the lower part of the matched chromospheric interfacial mode, put ${\omega }^2=gk(1-\epsilon )$, $\xi =e^{-k(z-z_0+\epsilon \varphi )}$ and $\Pi =\epsilon f\xi$. The pressure amplitude is of order $\epsilon$ because the limiting f-mode has identically zero Lagrangian pressure amplitude; the frequency correction perturbs that zero solution. Expand $g{k}^2/{\omega }^2=k(1+\epsilon )+O({\epsilon }^2)$ and ${\omega }^2-g^2{k}^2/{\omega }^2=-2gk\epsilon +O({\epsilon }^2)$. The leading equations areExclude the lower homogeneous component $f\propto e^{2kz}$, which would make $\Pi$ grow into the star. Integrating the remaining equation with its integrating factor givesThe accompanying correction is $\varphi (z)={\int }_{z_0}^z[1+(f/\rho )(1/g-1/(k{c}^2))]ds$, with the integrand evaluated at $s$. This fixes $\varphi (z_0)=0$ and determines the lower eigenfunctions to the stated order.

For the upper part, distinguish the rescaled sound speed by writing $\rho =\epsilon \hat{\rho }$ and $c^2={\epsilon }^{-1}{\hat{c}}^2$. Set $\Pi =\epsilon \psi$, $\psi =A{e}^{k(z-z_0+\epsilon \varphi )}$ and $\xi =\mu \psi$. To leading order the first displacement equation becomes ${\xi }^{\prime }=-k\xi -k\psi /(g\hat{\rho })$, while the pressure equation gives ${\psi }^{\prime }=k\psi +\epsilon [k\psi +2gk\hat{\rho }\xi ]+O({\epsilon }^2)$. ConsequentlyThe component $\mu \propto e^{-2kz}$ would give $\xi \propto e^{-kz}$, which diverges upward. Remove it by the upper boundary condition. The leading upper solution isIntegrating ${\varphi }^{\prime }$ with $\varphi (z_0)=0$ gives the first correction to the upper exponential. The displayed $\xi =\mu \psi$ satisfies the required decay only for compatible upper density profiles; integral convergence alone is not the full boundary condition.

To match, define the positive weighted integralsThe lower normalization is $\xi (z_0)=1$ and $f(z_0)=-2gk{I}_b$. Continuity of Lagrangian pressure perturbation makes the upper $A=f(z_0)$; continuity of fluid displacement then requires $\mu (z_0)A=1$. Since $1/\hat{\rho }=\epsilon /\rho$, the upper integral gives $\mu (z_0)=-(k/g)\epsilon I_a$. The two matching conditions therefore giveAs a sign and normalization check, two incompressible constant-density layers give $I_a=(2k{\rho }_a{)}^{-1}$ and $I_b={\rho }_b/(2k)$, so $\epsilon \simeq 2{\rho }_a/{\rho }_b$. This agrees with expansion of the exact interfacial gravity-wave dispersion relation ${\omega }^2=gk({\rho }_b-{\rho }_a)/({\rho }_b+{\rho }_a)$.

For validity, require $\epsilon \ll 1$, convergent weighted integrals, a matching level separating the low-density high-sound-speed atmosphere from the denser envelope, and small accumulated corrections $\epsilon k|\varphi |$ over the regions carrying the matching integrals. The upper expansion also needs $g/(k{c}^2)\ll 1$. The resulting upper displacement must tend to zero: for example if $\hat{\rho }\propto e^{a(z-z_0)}$ at great height, its leading amplitude is proportional to $e^{(k-a)(z-z_0)}$, requiring $a<k$; the weaker condition $a<2k$ merely makes the upper integral converge. Plane geometry, constant $g$ and the Cowling approximation must be appropriate on the wavelength scale. The physical mode amplitude must remain small compared with both its wavelength and local background scale heights, so the adiabatic fluid perturbation remains linear. If these requirements fail, neither the matching formula nor continuation of the leading expansions is controlled.

Assume a nonrotating, nonmagnetic equilibrium, inviscid adiabatic perturbations, and no dissipation. Let $g={\varphi }^{\prime }(R)>0$, so hydrostatic equilibrium is $p^{\prime }=-\rho g$. Use $\mathbf{\xi }{e}^{i\omega t}$ for the fluid displacement. To avoid a source ambiguity, write Eulerian fluid perturbations as $p_1,{\rho }_1,{\varphi }_1$, and use $\Delta$ for quantities following the displaced material element, as in Eulerian and Lagrangian fluid perturbations. They are related by $\Delta p=p_1+{\xi }_R{p}^{\prime }$ and similarly for mass density and potential.

The printed identity needs two corrections: its compression denominator must be $\gamma p$, not $\gamma \rho$, and the pressure and potential appearing in this form of the identity must be Eulerian fluid perturbations. The printed denominator has the wrong dimensions for an energy density. The claimed Lagrangian interpretation does not produce the displayed buoyancy-compression split; the following derivation specifies the consistent version.

Linearized mass conservation, the adiabatic relation, and the Poisson equation giveConsequently $p_1=-\gamma p\nabla \cdot \mathbf{\xi }-{\xi }_R{p}^{\prime }$. Linearizing the Euler equations for an inviscid fluid about the static star gives the linear adiabatic stellar oscillation equationsRequire regularity at the center, a vanishing Lagrangian surface pressure, and the usual matching to a decaying exterior gravitational potential. Because $p$ and $\rho$ vanish at the surface, the pressure boundary term vanishes for regular displacements; compactly supported trials also suffice below. The perturbed Newtonian gravitational potential solves Laplace's equation outside the star and must be included there when its field energy is computed.

Multiply the momentum equation by $-{\mathbf{\xi }}^{\ast }$ and integrate. The pressure term is $-{\int }_V{p}_1\nabla \cdot {\mathbf{\xi }}^{\ast }dV$. Let $d=\nabla \cdot \mathbf{\xi }$. The adiabatic relation and hydrostatic equilibrium implyThe pressure and equilibrium-gravity terms then combine pointwise aswhere $N^2=g[(1/\gamma )p^{\prime }/p-{\rho }^{\prime }/\rho ]$ is the stellar buoyancy frequency. The potentially troublesome cross terms cancel exactly.

For the perturbed self-gravity term, integration by parts givesThus the corrected stellar displacement energy identity isThe right side is a real quadratic form. With these conservative boundary conditions, the force operator is a self-adjoint operator. Its Rayleigh quotient characterizes squared mode frequencies, and a negative admissible trial quotient proves an unstable spectral component; it does not require the trial displacement itself to be an exact eigenfunction.

For a generic compressive disturbance of scale $\ell$ in a region with slowly varying equilibrium coefficients, estimate ${\rho }_1\sim \rho \xi /\ell$, $p_1\sim \gamma p\xi /\ell$ and ${\varphi }_1\sim 4\pi G{\rho }_1{\ell }^2$. The self-gravity and compression energy densities have orderswith adiabatic sound speed $c_s^2=\gamma p/\rho$. Under the plausible stellar estimate $c_s^2\sim G\bar{\rho }{a}^2$, global modes can have comparable gravity and pressure contributions, while their ratio for small-scale compressive modes is of order $(\rho /\bar{\rho })(\ell /a{)}^2$. Thus perturbations with $\ell \ll c_s/\sqrt{4\pi G\rho }$ admit the Cowling approximation. This estimate compares two nonzero compressive terms; it cannot alone justify neglecting gravity for a trial deliberately chosen to have $p_1=0$.

For an explicit convective trial, choose $c<d<e<b$ and a nonzero smooth bumpLet $Y$ be a normalized spherical harmonic of degree $\ell \ge 1$, let $L=\ell (\ell +1)$, and let ${\nabla }_{\Omega }$ be the unit-sphere gradient. The pressure-free localized stellar convective trial isIt is supported away from the center and surface, so it satisfies every needed mechanical boundary condition. Since ${\Delta }_{\Omega }Y=-LY$,Therefore $p_1=-\gamma p\nabla \cdot \mathbf{\xi }+\rho g{\xi }_R=0$ exactly. Its buoyancy energy is ${\int }_d^e\rho N^2{R}^2|f{|}^{2}dR<0$, whereas its kinetic norm ${\int }_d^e\rho R^2(|f{|}^2+L|h{|}^2)dR$ is positive. In the Cowling approximation this gives a negative Rayleigh quotient, and hence $N^2<0$ in an interior shell implies convective instability. Retaining perturbed self-gravity only adds a nonpositive term, so the same trial proves instability in the full conservative problem too.

There is also a quantitative small-scale check for this very trial. Its mass density amplitude is ${\rho }_1=(\rho N^2/g)fY$, independent of $\ell$. Write ${\rho }_1=s(R)Y$ and ${\varphi }_1=F(R)Y$. The radial Poisson equation givesBy the Cauchy-Schwarz inequality, $A\le 4\pi G(\int R^4|s{|}^{2}dR{)}^{1/2}(A/L{)}^{1/2}$, and so the angular-degree bound on perturbed stellar self-gravity isIt tends to zero with horizontal scale $R/\ell$, while the strictly negative buoyancy integral is fixed. This justifies the small-scale neglect of gravity even though the trial's compression term vanishes.

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.