Past exam of the mathematics course of the University of Cambridge / 2015 / iii / Paper 58 / 3 / c

Let $g=Gm/r^2>0$, $H_P=-dr/d\log P>0$, $\nabla =d\log T/d\log P$, and ${\nabla }_{\mu }=d\log \mu /d\log P$. Consider a small radial fluid displacement $\xi$ whose sound-crossing time is short enough to keep its pressure equal to its surroundings. Its composition is frozen and its heat exchange negligible. This is the local stellar convective stability test.

For gas-pressure-dominated ideal gas matter, $\rho \propto \mu P/T$. The ambient gradient is thereforeThe displaced element conserves specific entropy and mean molecular weight, giving $(d\log \rho /dr{)}_{\mathrm{parcel}}=-(1-{\nabla }_{\mathrm{ad}})/H_P$, where ${\nabla }_{\mathrm{ad}}=({\Gamma }_2-1)/{\Gamma }_2$. The parcel-minus-environment mass density difference isIts buoyancy acceleration is minus $g$ times this difference. ThusPositive stellar buoyancy frequency squared gives a restoring force. The gas-pressure-dominated Ledoux criterion is thereforeEquality is marginal in this ideal adiabatic test. For monatomic gas ${\nabla }_{\mathrm{ad}}=2/5$. An inward increase of mean molecular weight has ${\nabla }_{\mu }>0$ and stabilizes the layer; uniform composition recovers the Schwarzschild criterion. For a radiative layer one substitutes the stellar radiative temperature gradient, ${\nabla }_{\mathrm{rad}}=3\kappa LP/(16\pi acGm{T}^4)$.

More generally define $\delta =-(\partial \log \rho /\partial \log T{)}_{P,\mu }$ and $\phi =(\partial \log \rho /\partial \log \mu {)}_{P,T}$. The same displacement argument gives $N^2=(g/H_P)[\delta ({\nabla }_{\mathrm{ad}}-\nabla )+\phi {\nabla }_{\mu }]$ and stability for $\nabla <{\nabla }_{\mathrm{ad}}+(\phi /\delta ){\nabla }_{\mu }$. Both coefficients equal one in the gas-dominated ideal-gas limit used above. This is local dynamical stability against convection; it is distinct from global radial collapse and from instabilities requiring heat or composition diffusion.