User:Tohline/SSC/Stability BoundedCompositePolytropes

Instabilities in Bounded and Composite Polytropes

Unbounded, Complete Polytropes

Free-Energy Function and Its Derivatives

The free-energy function that is relevant to a discussion of the structure and stability of unbounded configurations having polytropic index, <math>~n</math>, has the form,

<math>~\mathcal{G}(x)</math>

<math>~=</math>

<math> -ax^{-1} +b x^{-3/n} + \mathcal{G}_0 \, , </math>

where <math>~x \equiv R/R_\mathrm{SWS}</math> identifies the radius of the configuration and <math>\mathcal{G}_0</math> is an arbitrary constant. If the coefficients, <math>~a, b</math>, and <math>~c</math>, are held constant while varying the configuration's size, we see that,

<math>~\frac{d\mathcal{G}}{dx}</math>

<math>~=</math>

<math> x^{-2} \biggl[ a - \frac{3b}{n}\cdot x^{(n-3)/n} \biggr] \, , </math>

and,

<math>~\frac{d^2\mathcal{G}}{dx^2}</math>

<math>~=</math>

<math> x^{-3} \biggl[ -2a + \frac{3(3+n)b}{n^2}\cdot x^{(n-3)/n} \biggr] \, . </math>

In terms of the system's mass and its structural form factors, <math>~\mathfrak{f}_M</math>, <math>~\mathfrak{f}_A</math>, and <math>~\mathfrak{f}_W</math>, the two relevant coefficients are,

<math>~a</math>	<math>~=</math>	<math>~\frac{3}{5}\biggl(\frac{n+1}{n}\biggr) \biggl[ \biggl( \frac{M}{M_\mathrm{SWS}}\biggr) \cdot \frac{1}{\mathfrak{f}_M} \biggr]^2 \mathfrak{f}_W \, ,</math>
<math>~b</math>	<math>~=</math>	<math>~n \biggl( \frac{3}{4\pi}\biggr)^{1/n} \biggl[ \biggl( \frac{M}{M_\mathrm{SWS}}\biggr) \cdot \frac{1}{\mathfrak{f}_M} \biggr]^{(n+1)/n} \mathfrak{f}_A \, . </math>

Equilibrium Radius

A configuration's equilibrium radius, <math>~x_\mathrm{eq}</math>, can be determined in one of two ways:

Extrema in the Free Energy

Equilibria are identified by extrema in the free-energy function. Setting <math>d\mathcal{G}/dx = 0</math>, we find,

<math>~x_\mathrm{eq}</math>	<math>~=</math>	<math> \biggl(\frac{an}{3b} \biggr)^{n/(n-3)} </math>
<math>~\Rightarrow ~~~~~ R_\mathrm{eq}^{n-3}</math>	<math>~=</math>	<math>\frac{R_\mathrm{SWS}^{n-3}}{M_\mathrm{SWS}^{n-1}} \biggl( \frac{4\pi}{3\cdot 5^n}\biggr) \biggl(\frac{n+1}{n}\biggr)^n \biggl[ \frac{\mathfrak{f}_W^n \mathfrak{f}_M^{1-n}}{\mathfrak{f}_A^n} \biggr] M^{n-1} </math>
	<math>~=</math>	<math>\biggl( \frac{4\pi}{3\cdot 5^n}\biggr) \biggl[ \frac{\mathfrak{f}_W^n \mathfrak{f}_M^{1-n}}{\mathfrak{f}_A^n} \biggr] G^n K^{-n} M^{n-1} \, . </math>

If one assumes that the equilibrium system has no internal structure — that is, that the interior density and temperature are uniform throughout — then <math>~\mathfrak{f}_M = \mathfrak{f}_A = \mathfrak{f}_W = 1</math> and this derived expression gives a good estimate of the equilibrium radius, given any choice of the pair of parameters, <math>~M</math> and <math>~K</math>.

Detailed Force Balance

Alternatively, a solution of the,

Lane-Emden Equation

<math>~\frac{1}{\xi^2} \frac{d}{d\xi}\biggl( \xi^2 \frac{d\Theta_H}{d\xi} \biggr) = - \Theta_H^n</math>

gives information regarding the detailed interior structure of the equilibrium system via knowledge of the properties of the Lane-Emden function, <math>~\Theta_H(\xi)</math>, as well as an exact expression for the equilibrium radius, namely,

<math>~R_\mathrm{eq}^{n-3} = \frac{3^n}{(n+1)^n} \biggl(\frac{4\pi}{3}\biggr) G^n K^{-n} M^{n-1} \cdot \frac{\mathfrak{f}_M^{1-n}}{\xi_1^{2n}} \, , </math>

where,

<math>~\mathfrak{f}_M \equiv \frac{\bar\rho}{\rho_c} = \biggl(- \frac{3\Theta_H^'}{\xi} \biggr)_{\xi_1} \, .</math>

Together

Now, once the Lane-Emden function, <math>~\Theta_H</math>, is known from a detailed force-balance model, the other two structural form factors also can be straightforwardly determined. For unbounded, complete polytropes, the relevant expressions are,

<math>~\mathfrak{f}_W</math>	<math>~=</math>	<math>~\frac{3^2\cdot 5}{5-n} \biggl[ \frac{\Theta_H^'}{\xi} \biggr]^2_{\xi_1} \, ,</math>
<math>~\mathfrak{f}_A</math>	<math>~=</math>	<math>~\frac{3(n+1) }{(5-n)} ~\biggl[ \Theta_H^' \biggr]^2_{\xi_1} \, .</math>

Plugging the appropriate ratio of these two functions, namely,

<math>~\biggl[\frac{\mathfrak{f}_W}{\mathfrak{f}_A} \biggr]^n</math>

<math>~=</math>

<math>~\biggl( \frac{3\cdot 5}{n+1} \biggr)^n \biggl[ \frac{1}{\xi} \biggr]^{2n}_{\xi_1} \, ,</math>

into the expression for the equilibrium radius obtained from the free-energy analysis gives precisely the same answer as was obtained from the detailed force-balance analysis. Using either method of determination we conclude, therefore, that,

<math>~\frac{K^n R_\mathrm{eq}^{n-3}}{G^n M^{n-1}} = \frac{4\pi}{(n+1)^n} \cdot \xi_1^{-(n+1)} (-\Theta_H^')^{1-n}_{\xi_1} \, . </math>

Stability

The procedure that has been used to obtain a detailed force-balanced model of unbounded polytropes cannot readily be extended to provide a stability analysis of such systems. However, the free-energy analysis can be readily extended. If the first derivative of the free-energy function is zero — that is, if you have identified an equilibrium configuration — and the second derivative of the free-energy function for that configuration is positive, then the equilibrium system is dynamically stable. If, however, the second derivative is negative, then the equilibrium system is dynamically unstable.

Using the expression for the second derivative of the free-energy function derived above, we deduce that equilibrium configurations are dynamically unstable when,

<math>~\biggl[\frac{3(3+n)b}{n^2}\cdot x^{(n-3)/n} - 2a \biggr]_{x_\mathrm{eq}} </math>	<math>~<</math>	<math> ~0 </math>
<math>~\Rightarrow ~~~~~ x_\mathrm{eq}^{(n-3)/n} </math>	<math>~<</math>	<math> ~\frac{2an^2}{3(3+n)b} </math>
<math>~\Rightarrow ~~~~~ \frac{an}{3b} </math>	<math>~<</math>	<math> ~\frac{2an^2}{3(3+n)b} </math>
<math>~\Rightarrow ~~~~~ n </math>	<math>~></math>	<math> ~3 \, . </math>

Bounded (Pressure-Truncated) Polytropes

Throughout this section we will rely on the following definitions of two normalization constants:

<math>M_\mathrm{SWS} \equiv \biggl( \frac{n+1}{nG} \biggr)^{3/2} K_n^{2n/(n+1)} P_\mathrm{e}^{(3-n)/[2(n+1)]} \, ,</math>

<math> R_\mathrm{SWS} \equiv \biggl( \frac{n+1}{nG} \biggr)^{1/2} K_n^{n/(n+1)} P_\mathrm{e}^{(1-n)/[2(n+1)]} \, . </math>

Free-Energy Function and Its Derivatives

The free-energy function that is relevant to a discussion of the structure and stability of a pressure-truncated configuration having polytropic index, <math>~n</math>, has the form,

<math>~\mathcal{G}(x)</math>

<math>~=</math>

<math> -ax^{-1} +b x^{-3/n} + c x^3 + \mathcal{G}_0 \, , </math>

where <math>~x \equiv R/R_\mathrm{SWS}</math> identifies the radius of the configuration and <math>\mathcal{G}_0</math> is an arbitrary constant. If the coefficients, <math>~a, b</math>, and <math>~c</math>, are held constant while varying the configuration's size, we see that,

<math>~\frac{d\mathcal{G}}{dx}</math>	<math>~=</math>	<math> ax^{-2} - \frac{3b}{n}\cdot x^{-(3+n)/n} + 3c x^2 </math>
	<math>~=</math>	<math> x^{-2} \biggl[ a - \frac{3b}{n}\cdot x^{(n-3)/n} + 3c x^4 \biggr] \, , </math>

and,

<math>~\frac{d^2\mathcal{G}}{dx^2}</math>

<math>~=</math>

<math> x^{-3} \biggl[ -2a + \frac{3(3+n)b}{n^2}\cdot x^{(n-3)/n} + 6c x^4 \biggr] \, . </math>

In terms of the system's mass and its structural form factors, <math>~\tilde\mathfrak{f}_M</math>, <math>~\tilde\mathfrak{f}_A</math>, and <math>~\tilde\mathfrak{f}_W</math>, the three relevant coefficients are,

<math>~a</math>	<math>~=</math>	<math>~\frac{3}{5}\biggl(\frac{n+1}{n}\biggr) \biggl[ \biggl( \frac{M}{M_\mathrm{SWS}}\biggr) \cdot \frac{1}{\tilde\mathfrak{f}_M} \biggr]^2 \tilde\mathfrak{f}_W \, ,</math>
<math>~b</math>	<math>~=</math>	<math>~n \biggl( \frac{3}{4\pi}\biggr)^{1/n} \biggl[ \biggl( \frac{M}{M_\mathrm{SWS}}\biggr) \cdot \frac{1}{\tilde\mathfrak{f}_M} \biggr]^{(n+1)/n} \tilde\mathfrak{f}_A \, , </math>
<math>~c</math>	<math>~=</math>	<math>~\frac{4\pi}{3} \, . </math>

Equilibrium Radius

A configuration's equilibrium radius, <math>~x_\mathrm{eq}</math>, can be determined in one of two ways:

Extrema in the Free Energy

Equilibria are identified by extrema in the free-energy function, that is, by setting <math>d\mathcal{G}/dx = 0</math>. Hence, <math>~x_\mathrm{eq}</math> is given by the root(s) of the polynomial expression that is often referred to as the,

Scalar Virial Theorem

<math>~x^{(n-3)/n}_\mathrm{eq} </math>

<math>~=</math>

<math>\frac{n}{b} \biggl[\frac{a}{3} + c\cdot x^4_\mathrm{eq} \biggr] \, . </math>

Plugging the definitions of the three free-energy coefficients into this virial expression gives the mass-radius relationship for pressure-truncated, polytropic equilibrium configurations, namely,

<math>~\biggl( \frac{3}{4\pi} \biggr)^{1/n} \tilde\mathfrak{f}_A \biggl[ \biggl( \frac{M}{M_\mathrm{SWS}} \biggr) \frac{1}{\tilde\mathfrak{f}_M} \biggr]^{(n+1)/n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} </math>

<math>~=</math>

<math>\frac{1}{5} \cdot \tilde\mathfrak{f}_W \biggl( \frac{n+1}{n}\biggr) \biggl[\biggl( \frac{M}{M_\mathrm{SWS}} \biggr) \frac{1}{\tilde\mathfrak{f}_M} \biggr]^{2} + \frac{4\pi}{3} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^4 \, . </math>

If one assumes that the equilibrium system has no internal structure — that is, that the interior density and temperature are uniform throughout — then <math>~\mathfrak{f}_M = \mathfrak{f}_A = \mathfrak{f}_W = 1</math> and this derived expression provides a good, approximate mass-radius relationship for configurations having a given <math>~K</math> that are truncated by a surrounding, massless medium exerting a pressure, <math>~P_e</math>, on the system.

Detailed Force Balance

From the detailed force-balance analysis presented in Appendix B of Steven W. Stahler (1983), we see that the mass, <math>~M</math>, associated with the equilibrium radius, <math>~R_\mathrm{eq}</math>, of bounded (pressure-truncated) polytropic spheres is given through the following pair of parametric relations:

<math> ~\frac{M}{M_\mathrm{SWS}} </math>	<math>~=~</math>	<math> \biggl( \frac{n^3}{4\pi} \biggr)^{1/2} \tilde\Theta_H^{(n-3)/2} \tilde\xi^2 (-\tilde\Theta_H^') \, , </math>
<math> ~\frac{R_\mathrm{eq}}{R_\mathrm{SWS} } </math>	<math>~=~</math>	<math> \biggl( \frac{n}{4\pi} \biggr)^{1/2} \tilde\xi \tilde\Theta_H^{(n-1)/2} \, . </math>

Together

As was realized in the case of unbounded polytropes, once the Lane-Emden function, <math>~\tilde\Theta_H(\xi)</math>, and its radial derivative, <math>~\tilde\Theta_H^'(\xi)</math>, are known from a detailed force-balance analysis, all three structural form factors can be straightforwardly determined. For bounded (pressure-truncated) polytropes, the relevant expressions are,

<math>~\tilde\mathfrak{f}_M</math>	<math>~=</math>	<math>~ - \frac{3\tilde\Theta_H^'}{\tilde\xi} \, ,</math>
<math>\tilde\mathfrak{f}_W</math>	<math>~=</math>	<math>\frac{3\cdot 5}{(5-n)\tilde\xi^2} \biggl[\tilde\Theta_H^{n+1} + 3 (\tilde\Theta_H^')^2 - \tilde\mathfrak{f}_M \tilde\Theta_H \biggr] \, , </math>
<math>~ \tilde\mathfrak{f}_A </math>	<math>~=</math>	<math>~\frac{1}{(5-n)} \biggl\{ 6\tilde\Theta_H^{n+1} + (n+1) \biggl[3 (\tilde\Theta_H^')^2 - \tilde\mathfrak{f}_M \tilde\Theta_H \biggr] \biggr\} \, . </math>

We have noticed, as well, that the relationship between <math>~\tilde\mathfrak{f}_A</math> and <math>~\tilde\mathfrak{f}_W</math> is relatively simple, namely,

<math>~ \tilde\mathfrak{f}_A </math>

<math>~=</math>

<math>~ \tilde\Theta_H^{n+1} +\biggl[ \frac{(n+1)}{3\cdot 5} \biggr] \tilde\xi^2\cdot \tilde\mathfrak{f}_W \, . </math>

As it turns out, it is sufficient to use this relationship — rather than the separate, explicit definitions of <math>~\tilde\mathfrak{f}_A</math> and <math>~\tilde\mathfrak{f}_W</math> — in order to obtain identical expressions for the equilibrium radius from both the detailed force-balance and free-energy analyses. By way of demonstration, let's plug this simpler expression for <math>~\tilde\mathfrak{f}_A</math> into the virial equilibrium relation.

<math>~\biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^4</math>	<math>~=</math>	<math> \biggl( \frac{3}{4\pi} \biggr)^{(n+1)/n} \biggl[ \tilde\Theta_H^{n+1} +\frac{(n+1)}{3\cdot 5} \cdot \tilde\xi^2 \tilde\mathfrak{f}_W \biggr] \biggl[ \biggl( \frac{M}{M_\mathrm{SWS}} \biggr) \frac{1}{\tilde\mathfrak{f}_M} \biggr]^{(n+1)/n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} - ~\frac{3}{5 \cdot 4\pi} \cdot \tilde\mathfrak{f}_W \biggl( \frac{n+1}{n}\biggr) \biggl[\biggl( \frac{M}{M_\mathrm{SWS}} \biggr) \frac{1}{\tilde\mathfrak{f}_M} \biggr]^{2} </math>
	<math>~=</math>	<math> \biggl[ \tilde\Theta_H^{n+1} +\frac{(n+1)}{3\cdot 5} \cdot \tilde\xi^2 \tilde\mathfrak{f}_W \biggr] \mathfrak{M}^{(n+1)/n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} - ~\frac{4\pi}{3\cdot 5} \cdot \tilde\mathfrak{f}_W \biggl( \frac{n+1}{n}\biggr) \mathfrak{M}^{2} </math>
	<math>~=</math>	<math> \biggl[ \tilde\Theta_H^{n} \mathfrak{M}\biggr]^{(n+1)/n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} + ~\biggl[ \mathfrak{M}^{(n+1)/n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} - ~\frac{4\pi}{n} \biggl( \frac{\mathfrak{M}}{\tilde\xi}\biggr)^{2} \biggr] \frac{(n+1)}{3\cdot 5} \cdot \tilde\xi^2 \tilde\mathfrak{f}_W \, , </math>

where, we have temporarily adopted the shorthand notation,

<math>~\mathfrak{M}</math>

<math>~\equiv</math>

<math>~\frac{3}{4\pi} \biggl( \frac{M}{M_\mathrm{SWS}} \biggr) \frac{1}{\tilde\mathfrak{f}_M} = \biggl( \frac{M}{M_\mathrm{SWS}} \biggr) \frac{\tilde\xi}{4\pi (- \tilde\Theta_H^')} \, .</math>

Regrouping terms from this last virial expression, we also can write,

<math>~\biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^4 \biggl\{1 - \biggl[ \tilde\Theta_H^{n} \mathfrak{M} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{-3} \biggr]^{(n+1)/n} \biggr\}</math>

<math>~=</math>

<math> ~\mathfrak{M}^2 \biggl\{ \frac{n\tilde\xi^2}{4\pi} \cdot \mathfrak{M}^{(1-n)/n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} - ~1 \biggr\} \frac{4\pi(n+1)}{3\cdot 5n} \cdot \tilde\mathfrak{f}_W \, . </math>

Now, this relation will be satisfied without having to plug in the fully explicit definition of <math>~\tilde\mathfrak{f}_W</math> if we demand that the expressions inside of the curly braces on both sides of the equation are zero, independently. This demand implies that (from the LHS),

<math>~\tilde\Theta_H^{n} \mathfrak{M} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{-3} </math>

<math>~=</math>

<math> ~1 \, ,</math>

and, independently (from the RHS),

<math>~ \frac{n\tilde\xi^2}{4\pi} \cdot \mathfrak{M}^{(1-n)/n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} </math>

<math>~=</math>

<math> ~1 \, . </math>

This pair of conditions, in turn, imply that,

<math>~ \biggl[ \tilde\Theta_H^{n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{-3} \biggr]^{(n-1)/n} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} </math>	<math>~=</math>	<math> ~\frac{4\pi}{n\tilde\xi^2} </math>
<math>~\Rightarrow ~~~~~ \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^{2} </math>	<math>~=</math>	<math> ~\frac{n}{4\pi} \cdot \tilde\xi^{2} \tilde\Theta_H^{n-1} \, ; </math>

and,

<math>~\mathfrak{M} </math>	<math>~=</math>	<math> ~\biggl(\frac{n}{4\pi} \cdot \tilde\xi^{2} \tilde\Theta_H^{n-1}\biggr)^{3/2} \tilde\Theta_H^{-n} = \biggl(\frac{n}{4\pi}\biggr)^{3/2} \cdot \tilde\xi^{3} \tilde\Theta_H^{(n-3)/2} </math>
<math>~\Rightarrow~~~~~\biggl( \frac{M}{M_\mathrm{SWS}} \biggr) </math>	<math>~=</math>	<math> ~ \biggl(\frac{n^3}{4\pi}\biggr)^{1/2} (- \tilde\xi^{2}\tilde\Theta_H^')\tilde\Theta_H^{(n-3)/2} \, , </math>

which match exactly the pair of parametric equations relating <math>~M</math> to <math>~R_\mathrm{eq}</math> that Stahler derived via his detailed force-balance analysis. We note as well that, after eliminating the parameter, <math>~P_e</math> — which appears in the definition of both normalizations, <math>~M_\mathrm{SWS}</math> and <math>~R_\mathrm{SWS}</math> — using either method of determination we have,

<math>~\frac{K^n R_\mathrm{eq}^{n-3}}{G^n M^{n-1}} = \frac{4\pi}{(n+1)^n} \cdot \tilde\xi^{-(n+1)} (-\tilde\Theta_H^')^{1-n} \, . </math>

Stability

As before, if the first derivative of the free-energy function is zero — that is, if you have identified an equilibrium configuration — and the second derivative of the free-energy function for that configuration is positive, then the equilibrium system is dynamically stable. If, however, the second derivative is negative, then the equilibrium system is dynamically unstable.

Using the expression for the second derivative of the free-energy function derived above, we deduce that equilibrium configurations are dynamically unstable when,

<math>~\biggl[\frac{3(3+n)b}{n^2}\cdot x^{(n-3)/n} - 2a + 6cx^4\biggr]_{x_\mathrm{eq}} </math>	<math>~<</math>	<math> ~0 </math>
<math>~\Rightarrow ~~~~~ x_\mathrm{eq}^{(n-3)/n} </math>	<math>~<</math>	<math> ~\frac{n^2}{3(3+n)b}\biggl[ 2a - 6cx^4_\mathrm{eq} \biggr] </math>
<math>~\Rightarrow ~~~~~ \frac{n}{b} \biggl[\frac{a}{3} + c\cdot x^4_\mathrm{eq} \biggr] </math>	<math>~<</math>	<math> ~\frac{n^2}{3(3+n)b}\biggl[ 2a - 6cx^4_\mathrm{eq} \biggr] </math>
<math>~\Rightarrow ~~~~~ (3+n) [a + 3c x^4_\mathrm{eq} ] </math>	<math>~<</math>	<math> ~n[ 2a - 6cx^4_\mathrm{eq} ] </math>
<math>~\Rightarrow ~~~~~ [3c(3+n) + 6nc] x^4_\mathrm{eq} </math>	<math>~<</math>	<math> ~[2n - (3+n)]a </math>
<math>~\Rightarrow ~~~~~ x^4_\mathrm{eq} </math>	<math>~<</math>	<math> ~\frac{(n-3)a}{9c(n+1)} \, . </math>

Plugging in the expressions for the free-energy coefficients, <math>~a</math> and <math>~c</math>, this in turn implies that an equilibrium system will be dynamically unstable if,

<math>~\biggl(\frac{R_\mathrm{eq}}{R_\mathrm{SWS}}\biggr)^4 </math>	<math>~<</math>	<math> ~\frac{(n-3)}{20\pi n} \biggl[ \biggl( \frac{M}{M_\mathrm{SWS}}\biggr) \cdot \frac{1}{\tilde\mathfrak{f}_M} \biggr]^2 \tilde\mathfrak{f}_W </math>
<math>~\Rightarrow ~~~~~ \biggl(\frac{R_\mathrm{eq}}{R_\mathrm{SWS}}\biggr)^4 \biggl( \frac{M}{M_\mathrm{SWS}}\biggr)^{-2}</math>	<math>~<</math>	<math> ~\biggl[ \frac{(n-3)}{12\pi n (5-n)} \biggr] \biggl[\frac{1}{\tilde\xi(-\tilde\Theta_H^')^2} \biggr] \biggl[\tilde\xi \tilde\Theta_H^{n+1} + 3\tilde\xi (\tilde\Theta_H^')^2 - 3(-\tilde\Theta_H^')\tilde\Theta_H \biggr] </math>
<math>~\Rightarrow ~~~~~ \frac{P_e R_\mathrm{eq}^4}{GM^2}</math>	<math>~<</math>	<math> ~\biggl[ \frac{(n-3)}{12\pi (5-n)(n+1)} \biggr] \biggl[\frac{1}{\tilde\xi(-\tilde\Theta_H^')^2} \biggr] \biggl[\tilde\xi \tilde\Theta_H^{n+1} + 3\tilde\xi (\tilde\Theta_H^')^2 - 3(-\tilde\Theta_H^')\tilde\Theta_H \biggr] \, . </math>

Composite Polytropes (Bipolytropes)

Throughout this section we will rely on the following definitions of three normalization constants:

<math>~R_\mathrm{norm}</math>	<math>~=</math>	<math>~\biggl[ \biggl( \frac{G}{K} \biggr)^n M_\mathrm{tot}^{n-1} \biggr]^{1/(n-3)} \, ,</math>
<math>~P_\mathrm{norm}</math>	<math>~=</math>	<math>~\biggl[ \frac{K^{4n}}{G^{3(n+1)} M_\mathrm{tot}^{2(n+1)}} \biggr]^{1/(n-3)} \, , </math>
<math>~E_\mathrm{norm}</math>	<math>~=</math>	<math>~ \biggl[ K^n G^{-3}M_\mathrm{tot}^{n-5} \biggr]^{1/(n-3)} \, .</math>

Free-Energy Function and Its Derivatives

The free-energy function that is relevant to a discussion of the structure and stability of a composite polytropic configurations of index, <math>~n_c</math>, for the core and, <math>~n_e</math>, for the envelope, has the form,

<math>~\mathcal{G}(x)</math>

<math>~=</math>

<math> -(a_c + a_e) x^{-1} +b_c x^{-3/n_c} +b_e x^{-3/n_e} + \mathcal{G}_0 \, , </math>

where <math>~x \equiv R/R_\mathrm{norm}</math> identifies the radius of the configuration and <math>\mathcal{G}_0</math> is an arbitrary constant. If the coefficients, <math>~a_c, a_e, b_c</math>, and <math>~b_e</math>, are held constant while varying the configuration's size, we see that,

<math>~\frac{d\mathcal{G}}{dx}</math>	<math>~=</math>	<math> (a_c + a_e)x^{-2} - \frac{3b_c}{n_c}\cdot x^{-(3+n_c)/n_c} - \frac{3b_e}{n_e}\cdot x^{-(3+n_e)/n_e} </math>
	<math>~=</math>	<math> x^{-2} \biggl[ (a_c + a_e) - \frac{3b_c}{n_c}\cdot x^{(n_c-3)/n_c} - \frac{3b_e}{n_e}\cdot x^{(n_e-3)/n_e} \biggr] \, , </math>

and,

<math>~\frac{d^2\mathcal{G}}{dx^2}</math>

<math>~=</math>

<math> x^{-3} \biggl[ -2(a_c + a_e) + \frac{3(3+n_c)b_c}{n_c^2}\cdot x^{(n_c-3)/n_c} + \frac{3(3+n_e)b_e}{n_e^2}\cdot x^{(n_e-3)/n_e} \biggr] \, . </math>

The four relevant coefficients are,

<math>~a_c</math>	<math>~\equiv</math>	<math> \biggl( \frac{\nu}{q^3} \biggr) \int_0^{q} 3\biggl[\frac{M_r(y)}{M_\mathrm{tot}} \biggr]_\mathrm{core} \biggl[ \frac{\rho(y)}{\bar\rho} \biggr]_\mathrm{core} y dy \, , </math>
<math>~a_e</math>	<math>~\equiv</math>	<math> \biggl( \frac{1-\nu}{1-q^3} \biggr) \int_{q}^{1} 3\biggl[\frac{M_r(y)}{M_\mathrm{tot}} \biggr]_\mathrm{env} \biggl[ \frac{\rho(y)}{\bar\rho} \biggr]_\mathrm{env} y dy \, , </math>
<math>~b_c</math>	<math>~\equiv</math>	<math> \frac{4\pi n_c}{3} \biggl[ \frac{P_{ic} x^{3(n_c+1)/n_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \int_0^q 3\biggl[\frac{1 - p_c(y)}{1-p_c(q)} \biggr] y^2 dy = \frac{4\pi n_c}{3} \biggl[ \frac{P_{ic} x^{3(n_c+1)/n_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} q^3 s_\mathrm{core} \, , </math>
<math>~b_e</math>	<math>~\equiv</math>	<math> \frac{4\pi n_e}{3} \biggl[ \frac{P_{ie} x^{3(n_e+1)/n_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \int_q^1 3\biggl[1 - p_e(y) \biggr] y^2 dy = \frac{4\pi n_e}{3} \biggl[ \frac{P_{ie} x^{3(n_e+1)/n_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} (1-q^3)s_\mathrm{env} \, . </math>

Specifically for analytically definable composite polytropes having <math>~(n_c, n_e) = (5, 1)</math>,

<math>~a_c</math>	<math>~\equiv</math>	<math> \frac{3}{5} \biggl[\biggl(\frac{\nu}{q^3} \biggr)^2 \biggl( 1 + \ell_i^2 \biggr)^{3} \biggr]_\mathrm{eq} \biggl(\frac{5}{2^4}\biggr) \biggl( \frac{q}{\ell_i}\biggr)^{5} \biggl[ \ell_i \biggl(\ell_i^4 - \frac{8}{3} \ell_i^2 -1 \biggr) (1 + \ell_i^2)^{-3} + \tan^{-1}(\ell_i) \biggr] \, , </math>
<math>~a_e</math>	<math>~\equiv</math>	<math> \biggl( \frac{\mu_e}{\mu_c} \biggr)^2 \biggl[ \frac{3^2 A^2 }{b_\eta^4 } \biggr] \biggl( \frac{\nu^2 \theta_i^4}{q^6} \biggr) \frac{1}{8b_\eta} \biggl[6 b_\eta y - 3\sin[2(b_\eta y-B)] - 4b_\eta y \sin^2(b_\eta y - B) \biggr]_q^1 </math>
	<math>~=</math>	<math> \biggl( \frac{\nu }{q^2 \ell_i} \biggr)^2 \frac{3A^2 }{8b_\eta^3} \biggl[2 b_\eta y - \sin[2(b_\eta y-B)] - \frac{4}{3}\cdot b_\eta y \sin^2(b_\eta y - B) \biggr]_q^1 \, , </math>
<math>~q^3 s_\mathrm{core}</math>	<math>~\equiv</math>	<math> ~\frac{3}{2^3}\biggl( \frac{q}{\ell_i} \biggr)^3 (1 + \ell_i^2)^3 \biggl[ \ell_i (\ell_i^4 - 1 )(1+\ell_i^2)^{-3} + \tan^{-1}(\ell_i) \biggr] \, , </math>
<math>~(1-q^3)s_\mathrm{env}</math>	<math>~\equiv</math>	<math> ~\biggl( \frac{P_{ic}}{P_{ie}} \biggr) \frac{3A^2}{4b_\eta^3} \biggl[ 2b_\eta y -\sin[2(b_\eta y - B)] \biggr]_q^1 \ , </math>

where,

<math>~\ell</math>	<math>~\equiv~</math>	<math>~\frac{\xi}{\sqrt{3}} \, ,</math>
<math>~A</math>	<math>~=~</math>	<math>~\frac{\eta_i}{\sin(\eta_i - B)} \, ,</math>
<math>~(\eta_s - B)</math>	<math>~=~</math>	<math>~\pi </math>       …	<math>~[\mathrm{hence}, ~\sin^2(\eta_s - B) = 0 ]</math>
<math>~\eta_i</math>	<math>~=~</math>	<math>~3\biggl( \frac{\mu_e}{\mu_c} \biggr) \ell_i ( 1 + \ell_i^2 )^{-1}\, ,</math>
<math>~b_\eta</math>	<math>~=~</math>	<math>~\eta_s = \frac{\eta_i}{q} \, .</math>

Let's try putting these terms all together in the context of the virial theorem.

<math>a_c - \frac{3b_c}{n_c}\cdot x^{(n_c-3)/n_c} </math>	<math>~=</math>	<math> \frac{3}{2^4}\biggl( \frac{q}{\ell_i} \biggr)^3 \biggl\{ \biggl(\frac{\nu}{q^2 \ell_i} \biggr)^2 \biggl[ - \frac{8}{3} \ell_i^3 + ( 1 + \ell_i^2 )^{3}\cdot \mathfrak{L}_1 \biggr] - 8\pi \cdot x^4 \biggl[ \frac{P_{ic} }{P_\mathrm{norm}} \biggr]_\mathrm{eq} \biggl[ (1 + \ell_i^2)^3 \cdot \mathfrak{L}_1 \biggr]\biggr\} </math>
	<math>~=</math>	<math> \frac{3}{2^4}\biggl( \frac{q}{\ell_i} \biggr)^3 \biggl\{ - \frac{2^3 \nu^2 \ell_i}{3 q^4 } +\biggl[ \biggl(\frac{\nu}{q^2 \ell_i} \biggr)^2 - 8\pi \cdot x^4 \biggl( \frac{P_{ic} }{P_\mathrm{norm}} \biggr)_\mathrm{eq} \biggr] (1 + \ell_i^2)^3 \cdot \mathfrak{L}_1 \biggr\} </math>
	<math>~=</math>	<math> - \frac{ \nu^2 }{2 q \ell_i^2} +\biggl[ \biggl(\frac{\nu}{q^2 \ell_i} \biggr)^2 - 8\pi \cdot x^4 \biggl( \frac{P_{ic} }{P_\mathrm{norm}} \biggr)_\mathrm{eq} \biggr] \frac{3}{2^4}\biggl( \frac{q}{\ell_i} \biggr)^3 (1 + \ell_i^2)^3 \cdot \mathfrak{L}_1 \, , </math>

<math>a_e - \frac{3b_e}{n_e}\cdot x^{(n_e-3)/n_e} </math>	<math>~=</math>	<math> \frac{3A^2 }{8b_\eta^3} \biggl\{ \biggl( \frac{\nu }{q^2 \ell_i} \biggr)^2 \biggl[\mathfrak{L}_2- \frac{4}{3}\cdot b_\eta y \sin^2(b_\eta y - B) \biggr]_q^1 - 8\pi \cdot x^4 \biggl[ \frac{P_{ic} }{P_\mathrm{norm}} \biggr]_\mathrm{eq} \biggl[ \mathfrak{L}_2\biggr]_q^1 \biggr\} </math>
	<math>~=</math>	<math> \frac{3A^2 }{8b_\eta^3} \biggl\{ - \frac{2^2 \nu^2 b_\eta }{3 q^4 \ell_i^2} \biggl[ y \sin^2(b_\eta y - B) \biggr]_q^1 +\biggl[ \biggl( \frac{\nu }{q^2 \ell_i} \biggr)^2 - 8\pi \cdot x^4 \biggl( \frac{P_{ic} }{P_\mathrm{norm}} \biggr)_\mathrm{eq}\biggr] \biggl[ \mathfrak{L}_2\biggr]_q^1 \biggr\} </math>
	<math>~=</math>	<math> - \frac{\nu^2 }{ 2q^4 \ell_i^2} \cdot \frac{A^2 }{b_\eta^2} \biggl[ y \sin^2(b_\eta y - B) \biggr]_q^1 +\biggl[ \biggl( \frac{\nu }{q^2 \ell_i} \biggr)^2 - 8\pi \cdot x^4 \biggl( \frac{P_{ic} }{P_\mathrm{norm}} \biggr)_\mathrm{eq}\biggr] \frac{3A^2 }{8b_\eta^3} \biggl[ \mathfrak{L}_2\biggr]_q^1 </math>
	<math>~=</math>	<math> + \frac{\nu^2 }{ 2q \ell_i^2} +\biggl[ \biggl( \frac{\nu }{q^2 \ell_i} \biggr)^2 - 8\pi \cdot x^4 \biggl( \frac{P_{ic} }{P_\mathrm{norm}} \biggr)_\mathrm{eq}\biggr] \frac{3A^2 }{8b_\eta^3} \biggl[ \mathfrak{L}_2\biggr]_q^1 \, , </math>

where the last step has been made by appreciating that,

<math>~\frac{A^2 }{b_\eta^2} \biggl[ y \sin^2(b_\eta y - B) \biggr]_q^1</math>	<math>~=</math>	<math> ~\biggl(\frac{q}{\eta_i}\biggr)^2\biggl[ \frac{\eta_i}{\sin(\eta_i-B)} \biggr]^2 \biggl[ \sin^2(b_\eta - B) - q \sin^2(b_\eta q - B) \biggr] </math>
	<math>~=</math>	<math> ~\biggl(\frac{q}{\eta_i}\biggr)^2\biggl[ \frac{\eta_i}{\sin(\eta_i-B)} \biggr]^2 \biggl[ \cancelto{0}{\sin^2(\eta_s - B)} - q \sin^2(\eta_i - B) \biggr] </math>
	<math>~=</math>	<math> ~-q^3 \, , </math>

and we have adopted the shorthand notation,

<math>~\mathfrak{L}_1</math>	<math>~\equiv</math>	<math>~ \ell_i (\ell_i^4 - 1 )(1+\ell_i^2)^{-3} + \tan^{-1}(\ell_i) \, ,</math>
<math>~\mathfrak{L}_2</math>	<math>~\equiv</math>	<math>~2b_\eta y -\sin[2(b_\eta y - B)] \, .</math>

Adding these two expressions together and setting the result equal to zero (as prescribed by the virial theorem) gives,

<math>~0 </math>	<math>~=</math>	<math>~ \biggl[ \biggl(\frac{\nu}{q^2 \ell_i} \biggr)^2 - 8\pi \biggl(\frac{R_\mathrm{eq}}{R_\mathrm{norm}} \biggr)^4 \biggl( \frac{P_{i} }{P_\mathrm{norm}} \biggr)_\mathrm{eq} \biggr] \biggl\{ \frac{3}{2^4}\biggl( \frac{q}{\ell_i} \biggr)^3 (1 + \ell_i^2)^3 \cdot \mathfrak{L}_1 +\frac{3A^2 }{8b_\eta^3} \biggl[ \mathfrak{L}_2\biggr]_q^1 \biggr\} </math>
<math>~\Rightarrow ~~~~ \biggl(\frac{R_\mathrm{eq}}{R_\mathrm{norm}} \biggr)^4 \biggl( \frac{P_{i} }{P_\mathrm{norm}} \biggr)_\mathrm{eq}</math>	<math>~=</math>	<math>~ \frac{1}{8\pi} \biggl(\frac{\nu}{q^2 \ell_i} \biggr)^2 </math>
<math>~\Rightarrow ~~~~ \frac{R_\mathrm{eq}^4 P_{i} }{GM_\mathrm{tot}^2} </math>	<math>~=</math>	<math>~ \frac{3}{8\pi} \biggl(\frac{\nu}{q^2 \xi_i} \biggr)^2 \, . </math>

Related Discussions

• Constructing BiPolytropes
• Analytic description of BiPolytrope with <math>(n_c, n_e) = (5,1)</math>
• Bonnor-Ebert spheres
• Schönberg-Chandrasekhar limiting mass
• Relationship between Bonnor-Ebert and Schönberg-Chandrasekhar limiting masses

• Wikipedia introduction to the Lane-Emden equation
• Wikipedia introduction to Polytropes

© 2014 - 2021 by Joel E. Tohline |   H_Book Home   |   YouTube   | Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS | Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation