Appendix/Ramblings/51BiPolytropeStability/NoAnalytic

Analog of Bonnor-Ebert Limiting Pressure

As has been demonstrated in an accompanying discussion, the mass of a pressure-truncated, n = 5 polytrope is,

$\frac{M}{M_{\mathrm{S}\mathrm{W}\mathrm{S}}}$

$=$

$\left[\right(\frac{3\cdot 5^3}{2^2\pi })\frac{({\xi }_e^2/3{)}^3}{(1+{\xi }_e^2/3{)}^4}{]}^{1/2}$

where,

$M_{\mathrm{S}\mathrm{W}\mathrm{S}}$	$\equiv$	$\left[\right(\frac{n+1}{nG}{)}^{3/2}{K}_n^{2n/(n+1)}{P}_{\mathrm{e}}^{(3-n)/[2(n+1)]}{]}_{n=5}$
	$\equiv$	$(\frac{2\cdot 3}{5G}{)}^{3/2}{K}^{5/3}{P}_{\mathrm{e}}^{-1/6}$

Do Not Confine Search to Analytic Eigenvector

Overview

STEP01:
Develop an algorithm (for Excel) that numerically integrates the LAWEs from the center to the surface of a $(n_c,n_e)=(5,1)$ bipolytrope, for an arbitrary specification of the three parameters:   ${\mu }_e/{\mu }_c,{\xi }_i$, and ${\sigma }_c^2$.

• Enforce the proper interface matching condition(s) at the interface location, ${\xi }_i$.
• Note that in general, for an arbitrarily chosen set of the three parameter values, the resulting surface displacement function will not match the desired boundary condition.

STEP02:
Fix your chosen value of the parameter pair, $({\mu }_e/{\mu }_c,{\xi }_i)$, and vary ${\sigma }_c^2$ until the proper surface boundary condition is realized.

• In an accompanying discussion, we claim to have identified at what point along various ${\mu }_e/{\mu }_c$ sequences the fundamental mode of radial oscillation becomes unstable — that is, when ${\sigma }_c^2=0$. For a given choice of ${\mu }_e/{\mu }_c$, it would be wise to begin our eigenvector search at a value of ${\xi }_i<[{\xi }_i{]}_{\mathrm{F}\mathrm{M}}$, as specified in the following table:

  Marginally Unstable Fundamental Modes
  $\frac{{\mu }_e}{{\mu }_c}$	$[{\xi }_i{]}_{\mathrm{F}\mathrm{M}}$
  1	1.6686460157
  $\frac{1}{2}$	2.27925811317
  0.345	2.560146865247
  $\frac{1}{3}$	2.582007485476
  0.309	2.6274239687695
  $\frac{1}{4}$	2.7357711469398
  See orange-colored triangular markers in the associated Figure 4

• Keep steadily raising the value of the interface location until you find the 1^st overtone mode; a related discussion (with animation) shows the results of this type of search in the context of isolated n = 1 polytropes. Our expectation is that, if this mode is unstable, the model will coincide with the turning point along the equilibrium sequence and its eigenvector will essentially overlap with the eigenvector found using the B-KB74 conjecture. At the same time, the square-of-the-eigenfrequency for the fundamental mode will be very negative.

STEP03:
Regarding analytically specified eigenvectors that satisfy the governing LAWES …

• If we force ${\sigma }_c^2=0$ in the core, we have shown that a parabolic-shaped eigenfunction satisfies the LAWE of the core. We expect this eigenfunction to precisely overlay the numerically determined, marginally unstable displacement function in both the case of the unstable fundamental mode and the case of the unstable 1^st overtone.
• If we force ${\sigma }_c^2=0$ in the envelope, we have derived a different — $\cos (\eta -B)$ dependent — eigenfunction that satisfies the LAWE of the envelope. However, this proves to be irrelevant in the context of our bipolytrope because the derived eigenfunction does not match the physically relevant surface boundary condition.

Renormalized LAWE

As presented, for example, in a parallel discussion, in terms of our original ( * ) parameter normalizations, the polytropic LAWE takes the form,

$0$	$=$	$\frac{d^{2}x}{dr{*}^2}+\{\frac{4}{r^{*}}-(\frac{{\rho }^{*}}{P^{*}}\left)\frac{M_r^{*}}{(r^{*}{)}^2}\right\}\frac{dx}{dr*}+\left(\frac{{\rho }^{*}}{P^{*}}\right)\{\frac{{\omega }^2}{{\gamma }_{\mathrm{g}}G{\rho }_c}+(\frac{4}{{\gamma }_{\mathrm{g}}}-3\left)\frac{M_r^{*}}{(r^{*}{)}^3}\right\}x$
	$=$	$\frac{d^{2}x}{dr{*}^2}+\{4-(\frac{{\rho }^{*}}{P^{*}}\left)\frac{M_r^{*}}{(r^{*})}\right\}\frac{1}{r^{*}}\frac{dx}{dr*}+\left(\frac{{\rho }^{*}}{P^{*}}\right)\{\frac{2\pi {\sigma }_c^2}{3{\gamma }_{\mathrm{g}}}-\frac{{\alpha }_{\mathrm{g}}{M}_r^{*}}{(r^{*}{)}^3}\}x,$

where,

${\rho }^{*}$	$\equiv$	$\frac{\rho }{{\rho }_0}$	;	$r^{*}$	$\equiv$	$\frac{r}{[K_c^{1/2}/(G^{1/2}{\rho }_0^{2/5})]}$
$P^{*}$	$\equiv$	$\frac{P}{K_c{\rho }_0^{6/5}}$	;	$M_r^{*}$	$\equiv$	$\frac{M_r}{[K_c^{3/2}/(G^{3/2}{\rho }_0^{1/5})]}$

New Normalization $\tilde{\rho }$ $\equiv$ $\rho \left[\right(\frac{K_c}{G}{)}^{3/2}\frac{1}{M_{\mathrm{t}\mathrm{o}\mathrm{t}}}{]}^{-5};$ $\tilde{P}$ $\equiv$ $P\left[K_c^{-10}{G}^9{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^6\right];$ $\tilde{r}$ $\equiv$ $r\left[\right(\frac{K_c}{G}{)}^{5/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{-2}],$ ${\tilde{M}}_r$ $\equiv$ $\frac{M_r}{M_{\mathrm{t}\mathrm{o}\mathrm{t}}};$ $\tilde{H}$ $\equiv$ $H\left[K_c^{-5/2}{G}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}\right].$

Switching to the new normalization, where it is understood that,

$(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}$

$=$

$(\frac{{\mu }_e}{{\mu }_c}{)}^{-2}{{𝓂}}_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}},$

where,

${{𝓂}}_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}$

$\equiv$

$\left(\frac{2}{\pi }{)}^{1/2}{\theta }_i^{-1}\right(-{\eta }^2\frac{d\varphi }{d\eta }{)}_s=(\frac{2}{\pi }{)}^{1/2}\frac{A{\eta }_s}{{\theta }_i},$

we find the following relevant relations:

$(r^{*}{)}^2$	$=$	$r^2\left[\frac{G{\rho }_0^{4/5}}{K_c}\right]={\tilde{r}}^2\left[\right(\frac{K_c}{G}{)}^{5/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{-2}{]}^{-2}\left[\frac{G{\rho }_0^{4/5}}{K_c}\right]$
	$=$	${\tilde{r}}^2\left[\right(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}{]}^4;$
$\frac{{\rho }^{*}}{P^{*}}$	$=$	$\left(\frac{\rho }{{\rho }_0}\right)\left[\frac{K_c{\rho }_0^{6/5}}{P}\right]=\left(\frac{\rho }{P}\right)\left[K_c{\rho }_0^{1/5}\right]$
	$=$	$\frac{\tilde{\rho }}{\tilde{P}}\left[\right(\frac{K_c}{G}{)}^{3/2}\frac{1}{M_{\mathrm{t}\mathrm{o}\mathrm{t}}}{]}^5\left[K_c^{-10}{G}^9{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^6\right]\left[K_c{\rho }_0^{1/5}\right]$
	$=$	$\frac{\tilde{\rho }}{\tilde{P}}\left[\right(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}];$
$\frac{M_r^{*}}{r^{*}}$	$=$	$\frac{M_r}{r}\left\{\frac{[K_c^{1/2}/(G^{1/2}{\rho }_0^{2/5})]}{[K_c^{3/2}/(G^{3/2}{\rho }_0^{1/5})]}\right\}=\frac{M_r}{r}\left[\right(\frac{G}{K_c}\left){\rho }_0^{-1/5}\right]$
	$=$	$\frac{{\tilde{M}}_r}{\tilde{r}}\left[\right(\frac{K_c}{G}{)}^{5/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{-2}\left]\right[\left(\frac{G}{K_c}\right){\rho }_0^{-1/5}]M_{\mathrm{t}\mathrm{o}\mathrm{t}}$
	$=$	$\frac{{\tilde{M}}_r}{\tilde{r}}\left[\right(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}{]}^{-1}.$

Therefore, in terms of the renormalized variables the LAWE becomes,

$0$	$=$	$\frac{d^{2}x}{dr{*}^2}+\{4-(\frac{{\rho }^{*}}{P^{*}}\left)\frac{M_r^{*}}{(r^{*})}\right\}\frac{1}{r^{*}}\frac{dx}{dr*}+\left(\frac{{\rho }^{*}}{P^{*}}\right)\{\frac{2\pi {\sigma }_c^2}{3{\gamma }_{\mathrm{g}}}-\frac{{\alpha }_{\mathrm{g}}{M}_r^{*}}{(r^{*}{)}^3}\}x$
	$=$	$\left[\right(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}{]}^{-4}\frac{d^{2}x}{d{\tilde{r}}^2}+\{4-(\frac{\tilde{\rho }}{\tilde{P}}\left)\frac{{\tilde{M}}_r}{\tilde{r}}\right\}\left[\right(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}{]}^{-4}\frac{1}{\tilde{r}}\frac{dx}{d\tilde{r}}+\frac{\tilde{\rho }}{\tilde{P}}\left[\right(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}\left]\right\{\frac{2\pi {\sigma }_c^2}{3{\gamma }_{\mathrm{g}}}-\frac{{\alpha }_{\mathrm{g}}{\tilde{M}}_r}{(\tilde{r}{)}^3}\left[\right(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}{]}^{-5}\}x.$

After multiplying through by

$\left[\right(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}{]}^4$

$=$

$\left[\right(\frac{{\mu }_e}{{\mu }_c}{)}^{-2}{{𝓂}}_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}{]}^4$

we have,

$0$	$=$	$\frac{d^{2}x}{d{\tilde{r}}^2}+\{4-(\frac{\tilde{\rho }}{\tilde{P}}\left)\frac{{\tilde{M}}_r}{\tilde{r}}\right\}\frac{1}{\tilde{r}}\frac{dx}{d\tilde{r}}+\frac{\tilde{\rho }}{\tilde{P}}\left\{\frac{2\pi {\sigma }_c^2}{3{\gamma }_{\mathrm{g}}}\right[(\frac{G}{K_c}{)}^{3/2}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}{\rho }_0^{1/5}{]}^5-\frac{{\alpha }_{\mathrm{g}}{\tilde{M}}_r}{(\tilde{r}{)}^3}\}x$
	$=$	$\frac{d^{2}x}{d{\tilde{r}}^2}+\{4-(\frac{\tilde{\rho }}{\tilde{P}}\left)\frac{{\tilde{M}}_r}{\tilde{r}}\right\}\frac{1}{\tilde{r}}\frac{dx}{d\tilde{r}}+\frac{\tilde{\rho }}{\tilde{P}}\left\{\frac{2\pi {\sigma }_c^2}{3{\gamma }_{\mathrm{g}}}\right[\left(\frac{{\mu }_e}{{\mu }_c}{)}^{-10}{{𝓂}}_{\mathrm{s}\mathrm{u}\mathrm{r}\mathrm{f}}^5\right]-\frac{{\alpha }_{\mathrm{g}}{\tilde{M}}_r}{(\tilde{r}{)}^3}\}x.$

See Also