SSC/Stability/NeutralMode

LAWE

Most General Form

In an accompanying discussion, we derived the so-called,

Adiabatic Wave (or Radial Pulsation) Equation

$\frac{d^{2}x}{d{r}_0^2}+[\frac{4}{r_0}-(\frac{g_0{\rho }_0}{P_0}\left)\right]\frac{dx}{d{r}_0}+\left(\frac{{\rho }_0}{{\gamma }_{\mathrm{g}}{P}_0}\right)[{\omega }^2+(4-3{\gamma }_{\mathrm{g}})\frac{g_0}{r_0}]x=0$

whose solution gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, self-gravitating fluid configurations. The boundary condition conventionally used in connection with the adiabatic wave equation is,

$r_0\frac{d\ln x}{d{r}_0}$

$=$

$\frac{1}{{\gamma }_g}(4-3{\gamma }_g+\frac{{\omega }^2{R}^3}{G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}})$        at         $r_0=R.$

Polytropic Configurations

Part 1

If the initial, unperturbed equilibrium configuration is a polytropic sphere whose internal structure is defined by the function, $\theta (\xi )$, that provides a solution to the,

Lane-Emden Equation

$\frac{1}{{\xi }^2}\frac{d}{d\xi }\left({\xi }^2\frac{d{\mathrm{\Theta }}_H}{d\xi }\right)=-{\mathrm{\Theta }}_H^n$

then,

$r_0$	$=$	$a_n\xi ,$
${\rho }_0$	$=$	${\rho }_c{\theta }^n,$
$P_0$	$=$	$K{\rho }_0^{(n+1)/n}=K{\rho }_c^{(n+1)/n}{\theta }^{n+1},$
$g_0$	$=$	$\frac{GM(r_0)}{r_0^2}=\frac{G}{r_0^2}\left[4\pi a_n^3{\rho }_c\right(-{\xi }^2\frac{d\theta }{d\xi }\left)\right],$

where,

$a_n$

$=$

$[\frac{(n+1)K}{4\pi G}\cdot {\rho }_c^{(1-n)/n}{]}^{1/2}.$

Hence, after multiplying through by $a_n^2$, the above adiabatic wave equation can be rewritten in the form,

$\frac{d^{2}x}{d{\xi }^2}+[\frac{4}{\xi }-\frac{g_0}{a_n}(\frac{a_n^2{\rho }_0}{P_0}\left)\right]\frac{dx}{d\xi }+\left(\frac{a_n^2{\rho }_0}{{\gamma }_{\mathrm{g}}{P}_0}\right)[{\omega }^2+(4-3{\gamma }_{\mathrm{g}})\frac{g_0}{a_n\xi }]x$

$=$

$0.$

In addition, given that,

$\frac{g_0}{a_n}$

$=$

$4\pi G{\rho }_c(-\frac{d\theta }{d\xi }),$

and,

$\frac{a_n^2{\rho }_0}{P_0}$

$=$

$\frac{(n+1)}{(4\pi G{\rho }_c)\theta }=\frac{a_n^2{\rho }_c}{P_c}\cdot \frac{{\theta }_c}{\theta },$

we can write,

$\frac{d^{2}x}{d{\xi }^2}+\left[\frac{4-(n+1)V(\xi )}{\xi }\right]\frac{dx}{d\xi }+\left[{\omega }^2\right(\frac{a_n^2{\rho }_c}{{\gamma }_g{P}_c})\frac{{\theta }_c}{\theta }-(3-\frac{4}{{\gamma }_g})\cdot \frac{(n+1)V(x)}{{\xi }^2}]x$

$=$

$0,$

where we have adopted the function notation,

$V(\xi )$

$\equiv$

$-\frac{\xi }{\theta }\frac{d\theta }{d\xi }.$

Part 2

Drawing from an accompanying discussion, we have the following:

Polytropic LAWE (linear adiabatic wave equation)

	$0=\frac{d^{2}x}{d{\xi }^2}+[4-(n+1)Q]\frac{1}{\xi }\cdot \frac{dx}{d\xi }+(n+1)\left[\right(\frac{{\sigma }_c^2}{6{\gamma }_g})\frac{{\xi }^2}{\theta }-\alpha Q]\frac{x}{{\xi }^2}$
where:    $Q(\xi )\equiv -\frac{d\ln \theta }{d\ln \xi },$    ${\sigma }_c^2\equiv \frac{3{\omega }^2}{2\pi G{\rho }_c},$     and,     $\alpha \equiv (3-\frac{4}{{\gamma }_{\mathrm{g}}})$

All physically reasonable solutions are subject to the inner boundary condition,

$\frac{dx}{d\xi }=0$         at         $\xi =0,$

but the relevant outer boundary condition depends on whether the underlying equilibrium configuration is isolated (surface pressure is zero), or whether it is a "pressure-truncated" configuration. As is the case with the pressure-truncated isothermal spheres, discussed above, if the polytropic configuration is truncated by the pressure, $P_e$, of a hot, tenuous external medium, then the solution to the LAWE is subject to the outer boundary condition,

$-\frac{d\ln x}{d\ln \xi }=3$         at         $\xi =\tilde{\xi }.$

But, for isolated polytropes — see the supporting derivation, below — the sought-after solution is subject to the more conventional boundary condition,

$-\frac{d\ln x}{d\ln \xi }=\left(\frac{3-n}{n+1}\right)+\frac{n{\sigma }_c^2}{6(n+1)}\left[\frac{\xi }{{\theta }^{\text{'}}}\right]$         at         $\xi ={\xi }_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}.$

See Also