Self-gravitating isothermal slab ($\rho (z)={\rho }_0{\mathrm{sech}}^2(z/z_0)$)

A plane-symmetric, self-gravitating stellar population with constant vertical velocity dispersion satisfies the vertical Jeans equation and the Poisson equation for Newtonian gravity. With ${\Phi }^{\prime }(0)=0$, their solution isIts surface density is $\Sigma =2{\rho }_0{z}_0$, so $z_0={\sigma }_z^2/(\pi G\Sigma )$. Its outer density tail is exponential with scale height $z_0/2$; its central profile is not exponential.

For a self-gravitating isothermal slab with $\rho ={\rho }_0{\mathrm{sech}}^2(z/H)$ and $H^2={\sigma }^2/(2\pi G{\rho }_0)$, hydrostatic equilibrium gives the displayed Newtonian gravitational potential. A test particle has constant in-plane velocity and conserved vertical energy $E_z={\dot{z}}^2/2+\Phi (z)-\Phi (0)$. Its vertical motion is periodic for every finite $E_z>0$, with turning height $H\mathrm{arcosh}\exp [E_z/(2{\sigma }^2)]$. Near the plane it is a harmonic oscillator of squared frequency $4\pi G{\rho }_0$; at large height the force tends to the constant $2\pi G\Sigma$, where $\Sigma =2{\rho }_{0}H$. The ideal infinite slab has no finite vertical escape energy.

For an isothermal gas slab, a horizontal Fourier mode has $\delta \rho =-ik\rho {\xi }_x-(\rho {\xi }_z{)}_z$, $W=c_s^2\delta \rho /\rho +\varphi$, ${\omega }^2{\xi }_x=ikW$, ${\omega }^2{\xi }_z=W_z$ and $({\partial }_z^2-k^2)\varphi =4\pi G\delta \rho$. Multiply momentum by $\rho {\mathbf{\xi }}^{\ast }$ and integrate by parts, then use the Poisson equation for Newtonian gravity. Vanishing pressure-work and gravitational boundary terms give the displayed real quadratic identity, with $I=\int \rho |\mathbf{\xi }{|}^{2}dz>0$. Thus every admissible squared eigenfrequency is real; a negative value gives gravitational growth rather than an oscillatory instability.

For mass density $\rho ={\rho }_0{\mathrm{sech}}^2(z/H)$ with $H^2=c_s^2/(2\pi G{\rho }_0)$, a zero-frequency gravitational mode with $k\ne 0$ obeys $W=0$. Its potential satisfies ${\varphi }_{zz}+[2{\mathrm{sech}}^2(z/H)/H^2-k^2]\varphi =0$. The substitution $\tau =\tanh (z/H)$ gives the associated Legendre differential equation of degree one and order $|k|H$. Solutions decaying at both infinities exist only at $|k|H=1$, when $\varphi \propto \mathrm{sech}(z/H)$. Longer horizontal wavelengths are unstable. Neutral mass-preserving relabeling displacements and the uniform vertical translation at $k=0$ are distinct from this fragmentation threshold.