= Solution
Assume a nonrotating, nonmagnetic equilibrium, inviscid adiabatic perturbations, and no dissipation. Let $g=\phi'(R)>0$, so <hydrostatic equilibrium> is $p'=-\rho g$. Use $\boldsymbol\xi e^{i\omega t}$ for the <fluid displacement>. To avoid a source ambiguity, write <Eulerian fluid perturbations> as $p_1,\rho_1,\phi_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_Rp'$ and similarly for <mass density> and potential.
The printed identity needs two corrections: its compression denominator must be \b[$\gamma p$, not $\gamma\rho$], and the <pressure> and potential appearing in this form of the identity must be \b[<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> give
$$
\boxed{\rho_1=-\nabla\cdot(\rho\boldsymbol\xi),\qquad \Delta p=-\gamma p\nabla\cdot\boldsymbol\xi,\qquad \nabla^2\phi_1=4\pi G\rho_1.}
$$
Consequently $p_1=-\gamma p\nabla\cdot\boldsymbol\xi-\xi_Rp'$. Linearizing the <Euler equations for an inviscid fluid> about the static star gives the <linear adiabatic stellar oscillation equations>
$$
-\omega^2\rho\boldsymbol\xi=-\nabla p_1-\rho_1\nabla\phi-\rho\nabla\phi_1.
$$
Require 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 $-\boldsymbol\xi^*$ and integrate. The <pressure> term is $-\int_Vp_1\nabla\cdot\boldsymbol\xi^*\,dV$. Let $d=\nabla\cdot\boldsymbol\xi$. The adiabatic relation and <hydrostatic equilibrium> imply
$$
d=\frac{\rho g\xi_R-p_1}{\gamma p},\qquad \rho_1=\frac{\rho p_1}{\gamma p}-\left(\frac{\rho^2g}{\gamma p}+\rho'\right)\xi_R.
$$
The <pressure> and equilibrium-gravity terms then combine pointwise as
$$
-p_1d^*+g\rho_1\xi_R^*=\frac{|p_1|^2}{\gamma p}-\left(\frac{\rho^2g^2}{\gamma p}+g\rho'\right)|\xi_R|^2=\frac{|p_1|^2}{\gamma p}+\rho N^2|\xi_R|^2,
$$
where $N^2=g[(1/\gamma)p'/p-\rho'/\rho]$ is the <stellar buoyancy frequency>. The potentially troublesome cross terms cancel exactly.
For the perturbed <self-gravity> term, integration by parts gives
$$
\int_V\rho\boldsymbol\xi^*\cdot\nabla\phi_1\,dV=\int_V\rho_1^*\phi_1\,dV=\frac1{4\pi G}\int_{\mathbb R^3}\phi_1\nabla^2\phi_1^*\,dV=-\frac1{4\pi G}\int_{\mathbb R^3}|\nabla\phi_1|^2\,dV.
$$
Thus the corrected <stellar displacement energy identity> is
$$
\boxed{\omega^2\int_V\rho|\boldsymbol\xi|^2\,dV=-\frac1{4\pi G}\int_{\mathbb R^3}|\nabla\phi_1|^2\,dV+\int_V\left(\frac{|p_1|^2}{\gamma p}+\rho N^2|\xi_R|^2\right)dV.}
$$
The 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 $\phi_1\sim4\pi G\rho_1\ell^2$. The <self-gravity> and compression energy densities have orders
$$
|E_g|/V\sim4\pi G\rho^2\xi^2,\qquad E_p/V\sim\rho c_s^2\xi^2/\ell^2,\qquad \frac{|E_g|}{E_p}\sim\frac{4\pi G\rho\ell^2}{c_s^2},
$$
with <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 bump
$$
f(R)=\begin{cases}\varepsilon\exp\!\left[-\dfrac{(e-d)^2}{(R-d)(e-R)}\right],&d<R<e,\\0,&\text{otherwise.}\end{cases}
$$
Let $Y$ be a normalized <spherical harmonic> of degree $\ell\geq1$, let $L=\ell(\ell+1)$, and let $\nabla_\Omega$ be the unit-sphere gradient. The <pressure-free localized stellar convective trial> is
$$
\boxed{\boldsymbol\xi=fY\mathbf e_R+h\nabla_\Omega Y,\qquad h=\frac R L\left[\frac{(R^2f)'}{R^2}-\frac{\rho g}{\gamma p}f\right].}
$$
It is supported away from the center and surface, so it satisfies every needed mechanical boundary condition. Since $\Delta_\Omega Y=-LY$,
$$
\nabla\cdot\boldsymbol\xi=\left[\frac{(R^2f)'}{R^2}-\frac LRh\right]Y=\frac{\rho g}{\gamma p}fY.
$$
Therefore $p_1=-\gamma p\nabla\cdot\boldsymbol\xi+\rho g\xi_R=0$ exactly. Its <buoyancy> energy is $\int_d^e\rho N^2R^2|f|^2dR<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 \b[$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 $\phi_1=F(R)Y$. The radial <Poisson equation> gives
$$
A=\int_0^\infty\left(R^2|F'|^2+L|F|^2\right)dR=-4\pi G\int_0^\infty R^2F^*s\,dR.
$$
By the <Cauchy-Schwarz inequality>, $A\leq4\pi G(\int R^4|s|^2dR)^{1/2}(A/L)^{1/2}$, and so the <angular-degree bound on perturbed stellar self-gravity> is
$$
|E_g|=\frac A{4\pi G}\leq\frac{4\pi G}{\ell(\ell+1)}\int_d^eR^4|s|^2dR.
$$
It 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.
Back to article page
