Bipolytrope Generalization

Setup

In a more general context, we have discussed a Gibbs-like free-energy function of the generic form,

<math> \mathfrak{G} = W_\mathrm{grav} + \mathfrak{S}_\mathrm{therm} + T_\mathrm{kin} + P_e V + \cdots </math>

Here we are interested in examining the free energy of isolated, nonrotating, spherically symmetric bipolytropes, so we can drop the term that accounts for the influence of an external pressure and we can drop the kinetic energy term. But we need to consider separately the contributions to the reservoir of thermodynamic energy by the core and envelope. In particular, we will assume that compressions/expansions occur adiabatically, but that the core and the envelope evolve along separate adiabats — <math>~\gamma_c</math> and <math>~\gamma_e</math>, respectively.

Review of Isolated Polytrope

If we were configuring isolated polytropes — instead of bipolytropes — the free-energy expression would be, simply,

<math> \mathfrak{G} = W_\mathrm{grav} + \mathfrak{S}_A \, , </math>

and, following the detailed steps presented in our introductory discussion of the free energy of spherically symmetric, configurations, properly normalized expressions for the two contributing energy terms would be,

<math>~\frac{W_\mathrm{grav}}{E_\mathrm{norm}}</math>	<math>~=</math>	<math> - \int_0^{\chi = R_\mathrm{edge}^} 3\biggl[\frac{M_r(r^)}{M_\mathrm{tot}} \biggr] r^* \rho^* dr^* </math>
	<math>~=</math>	<math> - \biggl\{ \biggl( \frac{\rho_c}{\bar\rho} \biggr) \int_0^{1} 3x \biggl[\frac{M_r(x)}{M_\mathrm{tot}} \biggr] \biggl[ \frac{\rho(x)}{\rho_c} \biggr] dx \biggr\}_\mathrm{eq} \cdot \chi^{-1} \, , </math>
<math>~\frac{\mathfrak{S}_A}{E_\mathrm{norm}}</math>	<math>~=</math>	<math>~\frac{1}{({\gamma_g}-1)} \int_0^{\chi=R_\mathrm{edge}^} 4\pi (r^)^2 P^* dr^* </math>
	<math>~=</math>	<math>~\frac{1}{({\gamma_g}-1)} \biggl( \frac{4\pi}{3} \biggr)^{1-\gamma} \biggl\{ \biggl( \frac{\rho_c}{\bar\rho} \biggr)^{\gamma} \int_0^{1} 3x^2 \biggl[ \frac{P(x)}{P_c} \biggr] dx \biggr\}_\mathrm{eq} \cdot \chi^{3-3\gamma} \, ,</math>

where,

<math>~\frac{M_r(x)}{M_\mathrm{tot} } </math>

<math>~ \biggl( \frac{\rho_c}{\bar\rho} \biggr)_\mathrm{eq} \int_0^{x} 3x^2 \biggl[ \frac{\rho(x)}{\rho_c} \biggr] dx \, .</math>

We note that, because <math>~M_r(x)/M_\mathrm{tot} = 1</math> in the limit, <math>~x \rightarrow 1</math>, we can write,

<math>~\biggl( \frac{\rho_c}{\bar\rho} \biggr)_\mathrm{eq} </math>

<math>~ \biggl[ \int_0^{1} 3x^2 \biggl( \frac{\rho(x)}{\rho_c} \biggr) dx \biggr]^{-1} \, ,</math>

or, if desired, the central-to-mean density ratio in one or both energy terms could be replaced by a term involving the normalized central pressure and the dimensionless equilibrium radius, <math>~\chi_\mathrm{eq}</math>, via the relation,

<math> \biggl[ \biggl(\frac{3}{4\pi} \biggr) \frac{\rho_c}{\bar\rho} \biggr]_\mathrm{eq}^{\gamma} </math>

<math>~\biggl[ \biggl( \frac{P_c}{P_\mathrm{norm}} \biggr) \chi^{3\gamma} \biggr]_\mathrm{eq} \, .</math>

Bipolytrope

When considering an isolated, spherically symmetric bipolytropic configuration, each energy term will be made up of separate contributions coming from the core and envelope, that is,

<math> \mathfrak{G} = (W_\mathrm{grav} + \mathfrak{S}_A)_\mathrm{core} + (W_\mathrm{grav} + \mathfrak{S}_A)_\mathrm{env} \, . </math>

Partitioning the Mass

The core will be principally defined in terms of two dimensionless parameters — <math>~q</math> and <math>~\nu</math> — which are, respectively, the core's radius relative to the bipolytrope's total radius, and the core's mass relative to the total mass of the bipolytropic configuration, specifically,

<math>~q</math>	<math>~\equiv</math>	<math>~x_i = \frac{r_i}{R_\mathrm{edge}} \, ,</math>
<math>~\nu</math>	<math>~\equiv</math>	<math>~\frac{M_\mathrm{core}}{M_\mathrm{tot}} \, .</math>

Given the separate (equilibrium) density profiles of the core and the envelope while sticking to the notation used in our introductory discussion, we can write,

<math>(\mathrm{For}~0 \leq r^* \leq r_i^*)</math> <math>~M_r </math>	<math>~=</math>	<math> \biggl( \frac{4\pi}{3} \biggr) R_\mathrm{norm}^3 \rho_\mathrm{norm} \int_0^{r^} 3 (r^)^2 \rho_\mathrm{core}^* dr^* </math>
	<math>~=</math>	<math> M_\mathrm{tot} \cdot \chi^3 \int_0^{x} 3 \biggl[ \frac{\rho_\mathrm{core}(x)}{{\bar\rho}_\mathrm{core}} \biggr] \biggl[ \frac{M_\mathrm{core}/(x_i R_\mathrm{edge})^3}{M_\mathrm{tot}/R^3_\mathrm{norm}} \biggr] x^2 dx </math>
	<math>~=</math>	<math> M_\mathrm{tot} \biggl( \frac{\nu}{q^3} \biggr) \int_0^{x} 3 \biggl[ \frac{\rho(x)}{\bar\rho} \biggr]_\mathrm{core} x^2 dx \, ; </math>
<math>(\mathrm{For}~r_i^* \leq r^* \leq R_\mathrm{edge}^*)</math> <math>~M_r </math>	<math>~=</math>	<math> M_\mathrm{core} + \biggl( \frac{4\pi}{3} \biggr) R_\mathrm{norm}^3 \rho_\mathrm{norm} \int_{r_i^}^{r^} 3 (r^)^2 \rho_\mathrm{env}^ dr^* </math>
	<math>~=</math>	<math>M_\mathrm{core} + M_\mathrm{tot} \cdot \chi^3 \int_{x_i}^{x} 3 \biggl[ \frac{\rho_\mathrm{env}(x)}{{\bar\rho}_\mathrm{env}} \biggr] \biggl\{ \frac{M_\mathrm{env}/[ (1-x_i^3) R_\mathrm{edge}^3]}{M_\mathrm{tot}/R^3_\mathrm{norm}} \biggr\} x^2 dx </math>
	<math>~=</math>	<math> M_\mathrm{tot} \biggl\{\nu + \biggl( \frac{1-\nu}{1-q^3} \biggr) \int_{x_i}^{x} 3 \biggl[ \frac{\rho(x)}{\bar\rho} \biggr]_\mathrm{env} x^2 dx \biggr\} \, . </math>

Example Bipolytrope Mass Profiles

Envelope and core of <math>~(n_c, n_e) = (0, 0) </math> bipolytrope:

In this case, <math>~\rho_\mathrm{core}(x) = \rho_c = </math> constant — hence, also, <math>~[\rho(x)/\bar\rho]_\mathrm{core} = 1</math> — and <math>~\rho_\mathrm{env}(x) = \rho_e = </math> constant — hence, also, <math>~[\rho(x)/\bar\rho]_\mathrm{env} = 1</math> — but in general <math>~\rho_e \ne \rho_c</math>. Performing the separate integrals to obtain expressions for <math>~M_r(r)</math> inside the core and the envelope, we obtain:

<math>(\mathrm{For}~0 \leq x \leq q)</math> <math>~M_r </math>	<math>~=</math>	<math> M_\mathrm{tot} \biggl( \frac{\nu}{q^3} \biggr) \int_0^{x} 3x^2 dx = \nu M_\mathrm{tot} \biggl( \frac{x}{q} \biggr)^3 \, ; </math>
<math>(\mathrm{For}~q \leq x \leq 1)</math> <math>~M_r </math>	<math>~=</math>	<math> M_\mathrm{tot} \biggl\{\nu + \biggl( \frac{1-\nu}{1-q^3} \biggr) \int_{q}^{x} 3 x^2 dx \biggr\} = M_\mathrm{core} + (1-\nu) M_\mathrm{tot}\biggl( \frac{x^3 - q^3}{1-q^3} \biggr)\, . </math>

When <math>~x = q</math>, both expressions give,

<math>~M_\mathrm{core} = \nu M_\mathrm{tot} \, ,</math>

as they should. We deduce, as well, that the mass contained in the envelope is,

<math>~M_\mathrm{env}</math>

<math>~M_\mathrm{tot} - M_\mathrm{core} = (1-\nu) M_\mathrm{tot} \, ,</math>

and that the volumes occupied by the core and envelope are, respectively,

<math>~V_\mathrm{core}</math>	<math>~=</math>	<math>~q^3 R_\mathrm{edge}^3 \, ,</math>
<math>~V_\mathrm{env}</math>	<math>~=</math>	<math>~(1- q^3) R_\mathrm{edge}^3 \, .</math>

Hence, the ratio of envelope density to core density is,

<math>~\frac{\rho_e}{\rho_c} = \frac{\bar\rho_\mathrm{env}}{\bar\rho_\mathrm{core}}</math>

<math>~ \frac{M_\mathrm{env}/V_\mathrm{env}}{M_\mathrm{core}/V_\mathrm{core}} = \frac{q^3(1-\nu)}{\nu(1-q^3)} \, . </math>

These relations should be compared to — and ultimately must match — the prescriptions for <math>~M_r</math> that have been presented elsewhere in connection with detailed force-balance models of <math>~(n_c, n_e) = (0, 0)</math> bipolytropes and in our introductory discussion of the virial stability of bipolytropes.

Core of <math>~(n_c, n_e) = (5, 1)</math> bipolytrope:

The core has <math>~n_c = 5 \Rightarrow \gamma_c = 1+1/n_c = 6/5</math>. Referring to the general relation derived above, and using <math>~\rho_0</math> to represent the central density, we can write,

<math>(\mathrm{For}~0 \leq x \leq q)</math> <math>~M_r </math>

<math> M_\mathrm{tot} \biggl( \frac{\nu}{q^3} \biggr) \biggl( \frac{\rho_0} {{\bar\rho}_\mathrm{core}}\biggr)_\mathrm{eq} \int_0^{x} 3 \biggl[ \frac{\rho(x)}{\rho_0} \biggr]_\mathrm{core} x^2 dx \, . </math>

Drawing on the derivation of detailed force-balance models of <math>~(n_c, n_e) = (5, 1)</math> bipolytropes, the density profile throughout the core is,

<math>~\biggl[ \frac{\rho(\xi)}{\rho_0} \biggr]_\mathrm{core}</math>

<math>~\biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-5/2} \, ,</math>

where the dimensionless radial coordinate is,

<math>~\biggl[ \frac{G \rho_0^{4/5}}{K_c} \biggr]^{1/2} \biggl( \frac{2\pi}{3} \biggr)^{1/2} r \, .</math>

Switching to the normalizations that have been adopted in the broad context of our discussion of configurations in virial equilibrium and inserting the adiabatic index of the core <math>~(\gamma_c = 6/5)</math> into all normalization parameters, we have,

<math>~R_\mathrm{norm} = \biggl[ \biggl(\frac{G}{K_c} \biggr) M_\mathrm{tot}^{2-\gamma} \biggr]^{1/(4-3\gamma)}</math>	<math>~\Rightarrow</math>	<math>~R_\mathrm{norm} = \biggl( \frac{G^5 M_\mathrm{tot}^4}{K_c^5} \biggr)^{1/2} \, ,</math>
<math>~\rho_\mathrm{norm} = \frac{3}{4\pi} \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2} \biggr]^{1/(4-3\gamma)}</math>	<math>~\Rightarrow</math>	<math>~\rho_\mathrm{norm} = \frac{3}{4\pi} \biggl( \frac{K_c^{3}}{G^3 M_\mathrm{tot}^2} \biggr)^{5/2} \, .</math>

Hence, we can rewrite,

<math>~\xi</math>	<math>~=</math>	<math>~\biggl( \frac{r}{R_\mathrm{norm}} \biggr) \biggl( \frac{\rho_0}{\rho_\mathrm{norm}} \biggr)^{2/5} \biggl[ \frac{G }{K_c} \biggr]^{1/2} \biggl( \frac{2\pi}{3} \biggr)^{1/2} R_\mathrm{norm} \rho_\mathrm{norm}^{2/5}</math>
	<math>~=</math>	<math>~r^* (\rho_0^*)^{2/5} \biggl[ \frac{G }{K_c} \biggr]^{1/2} \biggl( \frac{2\pi}{3} \biggr)^{1/2} \biggl( \frac{G^5 M_\mathrm{tot}^4}{K_c^5} \biggr)^{1/2} \biggl( \frac{3}{4\pi} \biggr)^{2/5} \biggl( \frac{K_c^{3}}{G^3 M_\mathrm{tot}^2} \biggr)</math>
	<math>~=</math>	<math> ~r^* (\rho_0^)^{2/5} \biggl[ \biggl( \frac{2\pi}{3} \biggr)^{5} \biggl( \frac{3}{4\pi} \biggr)^{4} \biggr]^{1/10} = r^ (\rho_0^*)^{2/5} \biggl[ \frac{\pi}{2^3 \cdot 3}\biggr]^{1/10} \, . </math>

Now, following the same approach as was used in our introductory discussion and appreciating that our aim here is to redefine the coordinate, <math>~\xi</math>, in terms of normalized parameters evaluated in the equilibrium configuration, we will set,

<math>~r^*</math>	<math>~\rightarrow~</math>	<math> ~ x \chi_\mathrm{eq} \, ; </math>
<math>~\rho_0^*</math>	<math>~\rightarrow~</math>	<math> \biggl[ \frac{\rho_0}{\bar\rho} \biggr]_\mathrm{core} \biggl( \frac{{\bar\rho}_\mathrm{core}}{\rho_\mathrm{norm}} \biggr) = \biggl[ \frac{\rho_0}{\bar\rho} \biggr]_\mathrm{core} \frac{\nu M_\mathrm{tot}/(q^3 R_\mathrm{edge}^3)_\mathrm{eq}}{M_\mathrm{tot}/R_\mathrm{norm}^3} = \frac{\nu}{q^3} \biggl[ \frac{\rho_0}{\bar\rho} \biggr]_\mathrm{core} \chi_\mathrm{eq}^{-3} \, . </math>

Then we can set,

in which case,

<math>~\biggl[ \frac{\rho(x)}{\rho_0} \biggr]_\mathrm{core}</math>

<math>~\biggl( 1 + a_\xi x^2 \biggr)^{-5/2} \, ,</math>

where the coefficient,

<math>~(3a_\xi)^{1/2}</math>	<math>~\equiv</math>	<math>~ \chi_\mathrm{eq} \biggl[ \frac{\nu}{q^3} \biggl( \frac{\rho_0}{\bar\rho} \biggr)_\mathrm{core} \chi_\mathrm{eq}^{-3} \biggr]^{2/5} \biggl( \frac{\pi}{2^3 \cdot 3}\biggr)^{1/10} =\chi_\mathrm{eq}^{-1/5} \biggl[ \frac{\nu}{q^3} \biggl( \frac{\rho_0}{\bar\rho} \biggr)_\mathrm{core} \biggr]_\mathrm{eq}^{2/5} \biggl( \frac{\pi}{2^3 \cdot 3}\biggr)^{1/10} </math>
<math>\Rightarrow~~~~a_\xi</math>	<math>~\equiv</math>	<math>~ \frac{1}{3} \biggl\{ \chi_\mathrm{eq}^{-1/5} \biggl[ \frac{\nu}{q^3} \biggl( \frac{\rho_0}{\bar\rho} \biggr)_\mathrm{core} \biggr]_\mathrm{eq}^{2/5} \biggl( \frac{\pi}{2^3 \cdot 3}\biggr)^{1/10} \biggr\}^2 = \chi_\mathrm{eq}^{-2/5} \biggl[ \frac{\nu}{q^3} \biggl( \frac{\rho_0}{\bar\rho} \biggr)_\mathrm{core} \biggr]_\mathrm{eq}^{4/5} \biggl( \frac{\pi}{2^3 \cdot 3^6}\biggr)^{1/5} \, . </math>

We therefore have,

<math>(\mathrm{For}~0 \leq x \leq q)</math> <math>~M_r </math>	<math>~=</math>	<math> M_\mathrm{tot} \biggl[ \frac{\nu}{q^3} \biggl( \frac{\rho_0} Template:\bar\rho\biggr)_\mathrm{core} \biggr]_\mathrm{eq} \int_0^{x} 3 \biggl( 1 + a_\xi x^2 \biggr)^{-5/2} x^2 dx </math>
	<math>~=</math>	<math> M_\mathrm{tot} \biggl[ \frac{\nu}{q^3} \biggl( \frac{\rho_0} Template:\bar\rho\biggr)_\mathrm{core} \biggr]_\mathrm{eq} \biggl[ x^3\biggl( 1 + a_\xi x^2 \biggr)^{-3/2} \biggr] \, . </math>

Note that, when <math>~x \rightarrow q</math>, <math>~M_r \rightarrow M_\mathrm{core} = \nu M_\mathrm{tot}</math>. Hence, this last expression gives,

<math>~\nu M_\mathrm{tot}</math>	<math>~=</math>	<math> M_\mathrm{tot} \biggl[ \frac{\nu}{q^3} \biggl( \frac{\rho_0} Template:\bar\rho\biggr)_\mathrm{core} \biggr]_\mathrm{eq} \biggl[ q^3\biggl( 1 + a_\xi q^2 \biggr)^{-3/2} \biggr] </math>
<math>\Rightarrow~~~~\biggl[\biggl( \frac{\rho_0} Template:\bar\rho\biggr)_\mathrm{core} \biggr]_\mathrm{eq}</math>	<math>~=</math>	<math> \biggl( 1 + a_\xi q^2 \biggr)^{3/2} \, . </math>

Hence, finally,

<math>(\mathrm{For}~0 \leq x \leq q)</math> <math>~M_r </math>

<math> \nu M_\mathrm{tot} \biggl( \frac{x^3}{q^3} \biggr) \biggl[ \frac{ 1 + a_\xi x^2 }{ 1 + a_\xi q^2 } \biggr]^{-3/2} \, ; </math>

and the coefficient, <math>~a_\xi</math>, will be determined only after the equilibrium radius, <math>~\chi_\mathrm{eq}</math>, has been determined, via the relation,

<math>~\chi_\mathrm{eq}^{2} </math>

<math>~\biggl( \frac{\pi}{2^3 \cdot 3^6}\biggr) \biggl( \frac{\nu}{q^3} \biggr)^{4} \biggl( 1 + a_\xi q^2 \biggr)^{6} a_\xi^{-5} \, . </math>

Separate Contributions to Gravitational Potential Energy

Given the separate (equilibrium) density and <math>~M_r</math> profiles of the core and the envelope while sticking to the notation used in our introductory discussion, we can write,

<math>~W_\mathrm{grav}\biggr\|_\mathrm{core}</math>	<math>~=</math>	<math> - E_\mathrm{norm} \int_0^{r_i^} 3\biggl[\frac{M_r(r^)}{M_\mathrm{tot}} \biggr] r^* \rho^* dr^* </math>
	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{\nu}{q^3} \biggr) \int_0^{x_i} 3x \biggl[\frac{M_r(x)}{M_\mathrm{tot}} \biggr]_\mathrm{core} \biggl[ \frac{\rho(x)}{\bar\rho} \biggr]_\mathrm{core} dx </math>
<math>~W_\mathrm{grav}\biggr\|_\mathrm{env}</math>	<math>~=</math>	<math> - E_\mathrm{norm} \int_{r_i^}^{\chi = R_\mathrm{edge}^} 3\biggl[\frac{M_r(r^)}{M_\mathrm{tot}} \biggr] r^ \rho^* dr^* </math>
	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{1-\nu}{1-q^3} \biggr) \int_{x_i}^{1} 3x \biggl[\frac{M_r(x)}{M_\mathrm{tot}} \biggr]_\mathrm{env} \biggl[ \frac{\rho(x)}{\bar\rho} \biggr]_\mathrm{env} dx </math>

Example Bipolytrope Gravitational Potential Energies

Core and envelope of <math>~(n_c, n_e) = (0, 0)</math> bipolytrope:

Let's do the core first. In this case, <math>~\rho_\mathrm{core}(x) = \rho_c = </math> constant — hence, also, <math>~[\rho(x)/\bar\rho]_\mathrm{core} = 1</math>. As shown elsewhere, the corresponding <math>~M_r</math> function is,

<math>~\biggl[ \frac{M_r(x)}{M_\mathrm{tot}}\biggr]_\mathrm{core} </math>

<math> \biggl( \frac{\nu}{q^3} \biggr) x^3 \, . </math>

Hence,

<math>~W_\mathrm{grav}\biggr\|_\mathrm{core}</math>	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{\nu}{q^3} \biggr) \int_0^{q} 3\biggl( \frac{\nu}{q^3} \biggr)x^4 dx </math>
	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{\nu}{q^3} \biggr)^2 \biggl( \frac{3}{5} q^5 \biggr) </math>

Now, let's do the envelope. In this case, <math>~\rho_\mathrm{env}(x) = \rho_e = </math> constant; hence, also, <math>~[\rho(x)/\bar\rho]_\mathrm{env} = 1</math>. As shown elsewhere, the corresponding <math>~M_r</math> function is,

<math>~\biggl[ \frac{M_r(x)}{M_\mathrm{tot}} \biggr]_\mathrm{env} </math>

<math> \nu + \biggl(\frac{1-\nu}{1-q^3} \biggr) (x^3 - q^3) \, . </math>

Hence,

<math>~W_\mathrm{grav}\biggr\|_\mathrm{env}</math>	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{1-\nu}{1-q^3} \biggr) \biggl\{ \int_{q}^{1} \biggl[ \nu -q^3 \biggl(\frac{1-\nu}{1-q^3} \biggr) \biggr]3x dx + \int_{q}^{1} \biggl[ \biggl(\frac{1-\nu}{1-q^3} \biggr) \biggr] 3x^4 dx \biggr\} </math>
	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{1-\nu}{1-q^3} \biggr) \biggl\{ \frac{3}{2} \biggl[ \nu -q^3 \biggl(\frac{1-\nu}{1-q^3} \biggr) \biggr] (1-q^2) + \frac{3}{5} \biggl[ \biggl(\frac{1-\nu}{1-q^3} \biggr) \biggr] (1-q^5) \biggr\} </math>
	<math>~=</math>	<math> - \frac{3}{5} \biggl(\frac{\nu^2}{q} \biggr) E_\mathrm{norm} \cdot \chi^{-1} \biggl[ \frac{1}{\nu} \biggl( \frac{1-\nu}{1-q^3} \biggr)\biggr] \biggl\{ \frac{5}{2} \biggl[ 1 - \frac{q^3}{\nu} \biggl(\frac{1-\nu}{1-q^3} \biggr) \biggr] (q-q^3) + \biggl[ \frac{q}{\nu} \biggl(\frac{1-\nu}{1-q^3} \biggr) \biggr] (1-q^5) \biggr\} </math>
	<math>~=</math>	<math> - \frac{3}{5} \biggl(\frac{\nu^2}{q} \biggr) E_\mathrm{norm} \cdot \chi^{-1} \biggl[ \frac{q^3}{\nu} \biggl( \frac{1-\nu}{1-q^3} \biggr)\biggr] \biggl\{ \frac{5}{2} \biggl[ 1 - \frac{q^3}{\nu} \biggl(\frac{1-\nu}{1-q^3} \biggr) \biggr] \biggl(\frac{1}{q^2}-1 \biggr) + \biggl[ \frac{q^3}{\nu} \biggl(\frac{1-\nu}{1-q^3} \biggr) \biggr] \biggl( \frac{1}{q^5}-1\biggr) \biggr\} \, . </math>

Realizing from our mass segregation derivation that,

<math>~\frac{q^3}{\nu} \biggl( \frac{1-\nu}{1-q^3} \biggr) = \frac{\rho_e}{\rho_c} \, ,</math>

this last expression can be rewritten as,

<math>~W_\mathrm{grav}\biggr\|_\mathrm{env}</math>	<math>~=</math>	<math> - \frac{3}{5} \biggl(\frac{\nu^2}{q} \biggr) E_\mathrm{norm} \cdot \chi^{-1} \biggl[ \frac{\rho_e}{\rho_c} \biggr] \biggl\{ \frac{5}{2} \biggl[ 1 - \frac{\rho_e}{\rho_c} \biggr] \biggl(\frac{1}{q^2}-1 \biggr) + \biggl[ \frac{\rho_e}{\rho_c} \biggr] \biggl( \frac{1}{q^5}-1\biggr) \biggr\} </math>
	<math>~=</math>	<math> - \frac{3}{5} \biggl(\frac{\nu^2}{q} \biggr) E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl\{ \frac{5}{2}\biggl(\frac{1}{q^2}-1 \biggr) +\biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl[ \frac{1}{q^5}-1 + \frac{5}{2} \biggl( 1-\frac{1}{q^2}\biggr)\biggr] \biggr\} \, . </math>

So, when put together to obtain the total gravitational potential energy, we have,

<math>~W_\mathrm{grav} = W_\mathrm{grav}\biggr|_\mathrm{core} + W_\mathrm{grav}\biggr|_\mathrm{env}</math>

<math> - \frac{3}{5} E_\mathrm{norm} \cdot \chi^{-1} \biggl(\frac{\nu^2}{q} \biggr) f(\nu,q) \, , </math>

where,

<math>~\equiv</math>

<math> 1+ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^2}-1 \biggr) +\biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \frac{1}{q^5}-1 + \frac{5}{2} \biggl( 1-\frac{1}{q^2}\biggr)\biggr] \, . </math>

(This result agrees with Tohline's earlier derivations in other sections of this H_Book, which may now be erased to avoid repetition.)

Core of <math>~(n_c, n_e) = (5, 1)</math> bipolytrope:

Borrowing from our separate discussion of the mass distribution in this type of bipolytrope, we can rewrite the general expression for the gravitational potential energy in the core as,

<math>~W_\mathrm{grav}\biggr\|_\mathrm{core}</math>	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl[ \frac{\nu}{q^3} \biggl(\frac{\rho_0}{\bar\rho} \biggr)_\mathrm{core} \biggr]_\mathrm{eq} \int_0^{q} 3x \biggl[\frac{M_r(x)}{M_\mathrm{tot}} \biggr]_\mathrm{core} \biggl[ \frac{\rho(x)}{\rho_0} \biggr]_\mathrm{core} dx </math>
	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl[ \frac{\nu}{q^3} \biggl( 1 + a_\xi q^2 \biggr)^{3/2} \biggr]_\mathrm{eq} \int_0^{q} 3x \biggl\{ \nu \biggl( \frac{x^3}{q^3} \biggr) \biggl[ \frac{ 1 + a_\xi x^2 }{ 1 + a_\xi q^2 } \biggr]^{-3/2} \biggr\} \biggl( 1 + a_\xi x^2 \biggr)^{-5/2} dx </math>
	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl[ 3\biggl(\frac{\nu}{q^3} \biggr)^2 \biggl( 1 + a_\xi q^2 \biggr)^{3} \biggr]_\mathrm{eq} \int_0^{q} x^4 \biggl( 1 + a_\xi x^2 \biggr)^{-4} dx </math>
	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl[ 3\biggl(\frac{\nu}{q^3} \biggr)^2 \biggl( 1 + a_\xi q^2 \biggr)^{3} \biggr]_\mathrm{eq} \biggl\{ \frac{a_\xi^{1/2} q(3a_\xi^2 q^4 - 8a_\xi q^2 - 3) + 3(a_\xi q^2 +1)^3 \tan^{-1}(a_\xi^{1/2} q)}{48 a_\xi^{5/2}(a_\xi q^2 + 1)^3} \biggr\} </math>
	<math>~=</math>	<math> - E_\mathrm{norm} \cdot \chi^{-1} \biggl[\biggl(\frac{3}{2^4}\biggr) a_\xi^{-5/2}\biggl(\frac{\nu}{q^3} \biggr)^2 \biggl( 1 + a_\xi q^2 \biggr)^{3} \biggr]_\mathrm{eq} \biggl[ a_\xi^{1/2} q(a_\xi^2 q^4 - \frac{8}{3}a_\xi q^2 - 1) (a_\xi q^2 +1)^{-3} + \tan^{-1}(a_\xi^{1/2} q) \biggr] \, . </math>

But, also from our earlier discussion, we can write,

<math>~a_\xi^{-5/2} \biggl( \frac{\nu}{q^3} \biggr)^2 (1 + a_\xi q^2)^3</math>

<math>~\chi_\mathrm{eq} \biggl( \frac{2^3 \cdot 3^6}{\pi} \biggr)^{1/2} \, .</math>

Hence,

<math>~\biggl( \frac{W_\mathrm{grav}}{E_\mathrm{norm}} \biggr)_\mathrm{core} </math>

<math> - \frac{\chi_\mathrm{eq}}{\chi} \biggl( \frac{3^8}{2^5\pi} \biggr)^{1/2} \biggl[ a_\xi^{1/2} q(a_\xi^2 q^4 - \frac{8}{3}a_\xi q^2 - 1) (a_\xi q^2 +1)^{-3} + \tan^{-1}(a_\xi^{1/2} q) \biggr] \, . </math>

After making the substitution, <math>~(a_\xi^{1/2} q) \rightarrow x_i</math>, this expression agrees with a result for the dimensionless energy, <math>~W^*_\mathrm{core}</math>,derived by Tohline in the context of detailed force-balanced bipolytropes.

Separate Thermodynamic Energy Reservoirs

In our introductory discussion of the free energy function for spherically symmetric configurations, we developed expressions that define the separate contributions to the thermodynamic energy reservoir that pertain to the core and envelope of bipolytropic configurations. In that discussion we pointed out that, in general for the core, the pressure drops monotonically from a value of <math>~P_0</math> at the center of the configuration according to an expression of the form,

<math>~P_\mathrm{core}(x) = P_0 [1 - p_c(x)]</math> for <math>~0 \leq x \leq q \, ,</math>

and that, for the envelope, the pressure drops monotonically from a value of <math>~P_{ie}</math> at the interface according to an expression of the form,

<math>~P_\mathrm{env}(x) = P_{ie} [1 - p_e(x)]</math> for <math>~q \leq x \leq 1 \, ,</math>

where <math>~p_c(x)</math> and <math>~p_e(x)</math> are both dimensionless functions that will depend on the equations of state that are chosen for the core and envelope, respectively. By prescription, the pressure in the envelope must drop to zero at the surface of the bipolytropic configuration, hence, we should expect that <math>~p_e(1) = 1</math>. Furthermore, by prescription, the pressure in the core will drop to a value, <math>~P_{ic}</math>, at the interface, so we can write,

In equilibrium — that is, when <math>~R_\mathrm{edge} = R_\mathrm{eq}</math> — we will demand that the pressure at the interface be the same, whether it is referenced in the core or in the envelope, that is, we will demand that <math>~P_{ic} = P_{ie} \, .</math> It is therefore strategically advantageous to rewrite the expression for the run of pressure through the core in terms of the pressure at the interface rather than in terms of the central pressure; specifically,

<math>~P_\mathrm{core}(x) = P_{ic} \biggl[\frac{1 - p_c(x)}{1-p_c(q)} \biggr] \, .</math>

KEY QUESTION: How should piecewise pressure be normalized?

Here we will operate under two constraints when relating the central pressure, <math>~P_0</math>, to the core pressure at the interface, <math>~P_{ic}</math>. First, the pressure in the core will drop monotonically from its central value according to the following very general prescription,

<math>~P_\mathrm{core}(x) = P_0 [1 - p_c(x)] \, .</math>

Second, it is the pressure at the interface, <math>~P_{ic} = P_\mathrm{core}(q)</math>, that will be determined by the virial equilibrium condition. There are two possible ways to determine the central pressure from knowledge of <math>~P_{ic}</math>, but which is the physically correct method to embrace?

Case A: Addition …

<math>~P_0 = P_{ic} + \Delta P \, ,</math> where <math>\Delta P \equiv P_0 p_c(q) \, ,</math>

in which case, we should write that,

<math>~P_\mathrm{core}(x) = ( P_{ic} + \Delta P ) - P_0 p_c(x) \, .</math>

Case B: Multiplication …

<math>~P_{ic} = P_0 [1-p_c(q)] ~~~~\Rightarrow~~~~~ P_0 = \frac{P_{ic}}{[1-p_c(q)]} \, ,</math>

in which case, we should write that,

<math>~P_\mathrm{core}(x) = P_{ic} \biggl[ \frac{ 1-p_c(x) }{ 1-p_c(q) } \biggr] \, .</math>

Prior to August 2014, we have been naively implementing "Case A," effectively assuming that the quantity, <math>~\Delta P</math> (as well as <math>~P_{ic}</math>), is held fixed as we search for the equilibrium value of <math>~P_0</math>. See, for example, the comment dated 12 February 2014 in connection with my discussions with Kundan Kadam where we set,

<math>~\Delta P = \Pi q^2 \, .</math>

But we now suspect that "Case B" is the proper approach to embrace because, once the parameter <math>~q</math> has been specified, it allows for the function, <math>~P_\mathrm{core}(x)</math>, to scale with the system size in exactly the same way as the interface pressure scales with size.

With these generic expressions for the pressure profile in hand, the separate components of the thermodynamic energy reservoir derived in our introductory discussion are,

<math>~\biggl( \frac{\mathfrak{S}_A}{E_\mathrm{norm}} \biggr)_\mathrm{core}</math>	<math>~=</math>	<math> \frac{4\pi }{({\gamma_c}-1)} \biggl[ \frac{P_{ic} \chi^{3\gamma_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_c} \int_0^q \biggl[\frac{1 - p_c(x)}{1-p_c(q)} \biggr] x^2 dx </math>
<math>~\biggl( \frac{\mathfrak{S}_A}{E_\mathrm{norm}} \biggr)_\mathrm{env}</math>		<math> \frac{4\pi }{({\gamma_e}-1)} \biggl[ \frac{P_{ie} \chi^{3\gamma_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_e} \int_q^1 \biggl[1 - p_e(x) \biggr] x^2 dx \, . </math>

Example Bipolytrope Thermodynamic Energy Reservoirs

Envelope and core of <math>~(n_c, n_e) = (0, 0) </math> bipolytrope:

According to our derivation of the properties of detailed force-balance <math>~(n_c, n_e) = (0, 0) </math> bipolytropes, in this case the pressure throughout the core is defined by the dimensionless function,

<math>~\biggl( \frac{2\pi}{3} \biggr) \xi^2 \, ,</math>

and the pressure throughout the envelope is defined by the dimensionless function,

<math>\frac{2\pi}{3} \biggl(\frac{\rho_e}{\rho_0}\biggr) \frac{P_0}{P_{ie} } \biggl[ \frac{\rho_e}{\rho_0} (\xi^2 - \xi_i^2) - 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \xi_i^3\biggl( \frac{1}{\xi} - \frac{1}{\xi_i}\biggr) \biggr] \, , </math>

where, for both functions,

<math>~\xi</math>	<math>~\equiv</math>	<math>~\biggl[ \biggl( \frac{G\rho_0^2}{P_0} \biggr)^{1/2} R_\mathrm{edge} \biggr]_\mathrm{eq} x</math>
	<math>~=</math>	<math>~R_\mathrm{norm} \biggl[ \biggl( \frac{G \rho_0^{2- {\gamma_c} }}{K_c} \biggr)^{1/2} \chi \biggr]_\mathrm{eq} x</math>
	<math>~=</math>	<math>~\chi_\mathrm{eq} \biggl[ \frac{G R_\mathrm{norm}^2 \rho_\mathrm{norm}^{2- {\gamma_c} }}{K_c} \biggr]^{1/2} (\rho_0^*)^{(2- {\gamma_c} )/2}_\mathrm{eq} \cdot x</math>
	<math>~=</math>	<math>~\chi_\mathrm{eq} \biggl( \frac{3}{4\pi} \biggr)^{1/2} \biggl[ \biggl( \frac{\rho_0}{\bar\rho} \biggr)_\mathrm{core} \frac{M_\mathrm{core}/M_\mathrm{tot}}{R_\mathrm{core}^3/R_\mathrm{norm}^3} \biggr]^{(2- {\gamma_c} )/2}_\mathrm{eq} x</math>
	<math>~=</math>	<math>~\chi_\mathrm{eq} \biggl( \frac{3}{4\pi} \biggr)^{1/2} \biggl( \frac{\rho_0}{\bar\rho} \biggr)_\mathrm{core}^{(2- {\gamma_c} )/2} \biggl[ \frac{\nu}{q^3} \cdot \chi^{-3} \biggr]^{(2- {\gamma_c} )/2}_\mathrm{eq} x</math>
	<math>~=</math>	<math>~\chi_\mathrm{eq}^{(3 {\gamma_c} -4)/2} \biggl( \frac{3}{4\pi} \biggr)^{1/2} \biggl[ \frac{\nu}{q^3} \biggr]^{(2- {\gamma_c} )/2}_\mathrm{eq} x \, .</math>

[In this last step, we set <math>~(\rho_0/\bar\rho)_\mathrm{core} \rightarrow 1</math> because we are examining the case of a uniform density core.] So, defining the coefficient,

<math>~\equiv</math>

<math>~\frac{1}{2}~ \chi_\mathrm{eq}^{3 {\gamma_c} -4} \biggl( \frac{\nu}{q^3} \biggr)^{2- {\gamma_c} } \, ,</math>

such that,

<math>~\biggl( \frac{3}{2\pi} \cdot b_\xi \biggr)^{1/2} x \, ,</math>

and remembering that, at the interface, <math>~x \rightarrow x_i = q</math>, so <math>~\xi_i = (3b_\xi/2\pi)^{1/2} q</math>, the two dimensionless pressure functions become,

and,

<math>b_\xi\biggl(\frac{\rho_e}{\rho_0}\biggr) \frac{P_0}{P_{ie} } \biggl[ \frac{\rho_e}{\rho_0} (x^2 - q^2) - 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^3\biggl( \frac{1}{x} - \frac{1}{q}\biggr) \biggr] \, . </math>

The desired integrals over these pressure distributions therefore give,

<math>~\int_0^q \biggl[\frac{1 - p_c(x)}{1-p_c(q)} \biggr] x^2 dx</math>	<math>~=</math>	<math>~\biggl[ \frac{1}{1-b_\xi q^2} \biggr] \int_0^q (1-b_\xi x^2)x^2 dx</math>
	<math>~=</math>	<math>~\biggl[ \frac{1}{1-b_\xi q^2} \biggr] \biggl[ \frac{1}{3}\cdot q^3 - \biggl( \frac{b_\xi}{5} \biggr) q^5 \biggr] </math>
	<math>~=</math>	<math>~\frac{q^3}{3} \biggl[ \frac{1}{1-b_\xi q^2} \biggr] \biggl[ 1 - \biggl( \frac{3b_\xi}{5} \biggr) q^2 \biggr] = \frac{q^3}{3} \biggl( \frac{P_0}{P_{ic}} \biggr) \biggl[ 1 - \biggl( \frac{3b_\xi}{5} \biggr) q^2 \biggr] \, ;</math>
<math>~\int_q^1 \biggl[1 - p_e(x) \biggr] x^2 dx</math>	<math>~=</math>	<math>~\frac{1}{3}(1-q^3) - b_\xi\biggl(\frac{\rho_e}{\rho_0}\biggr) \frac{P_0}{P_{ie} } \int_q^1 \biggl[ \frac{\rho_e}{\rho_0} (x^2 - q^2) - 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^3\biggl( \frac{1}{x} - \frac{1}{q}\biggr) \biggr] x^2 dx</math>
	<math>~=</math>	<math>~\frac{1}{3}(1-q^3) - b_\xi\biggl(\frac{\rho_e}{\rho_0}\biggr) \frac{P_0}{P_{ie} } \biggl[ \frac{\rho_e}{\rho_0} \biggl( \frac{x^5}{5} - \frac{q^2 x^3}{3} \biggr) - 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^3\biggl( \frac{x^2}{2} - \frac{x^3}{3q}\biggr) \biggr]_q^1</math>
	<math>~=</math>	<math>~\frac{1}{3}(1-q^3) - \frac{b_\xi}{3} \biggl(\frac{\rho_e}{\rho_0}\biggr) \frac{P_0}{P_{ie} } \biggl\{ \biggl[ \frac{\rho_e}{\rho_0} \biggl( \frac{3}{5} - q^2 \biggr) - \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^2\biggl( 3q - 2\biggr) \biggr] </math>
		<math>~ - \biggl[ \frac{\rho_e}{\rho_0} \biggl( -\frac{2}{5} \biggr)q^5 - \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^5 \biggr] \biggr\} </math>
	<math>~=</math>	<math>~\frac{1}{3}(1-q^3) - \frac{b_\xi}{3} \biggl(\frac{\rho_e}{\rho_0}\biggr) \frac{P_0}{P_{ie} } \biggl\{ \biggl[ q^2(2-3q) +q^5\biggr] + \frac{\rho_e}{\rho_0}\biggl[ \biggl( \frac{3}{5} - q^2 \biggr) + q^2\biggl( 3q - 2\biggr) +\frac{2q^5}{5} -q^5\biggr] \biggl\} </math>
	<math>~=</math>	<math>~\frac{1}{3}(1-q^3) - \frac{b_\xi}{3} \biggl(\frac{\rho_e}{\rho_0}\biggr) \frac{P_0}{P_{ie} } \biggl[ (2q^2 - 3q^3 +q^5) + \frac{3}{5} \cdot \frac{\rho_e}{\rho_0} ( 1 - 5q^2 + 5q^3 - q^5 ) \biggr] \, . </math>

Finally, then, we have,

<math>~\biggl( \frac{\mathfrak{S}_A}{E_\mathrm{norm}} \biggr)_\mathrm{core}</math>	<math>~=</math>	<math> \frac{4\pi/3 }{({\gamma_c}-1)} \biggl[ \frac{P_{ic} \chi^{3\gamma_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_c} \biggl\{ \biggl( \frac{P_0}{P_{ic}} \biggr) \biggl[ q^3 - \biggl( \frac{3b_\xi}{5} \biggr) q^5 \biggr] \biggr\} </math>
<math>~\biggl( \frac{\mathfrak{S}_A}{E_\mathrm{norm}} \biggr)_\mathrm{env}</math>	<math>~=</math>	<math> \frac{4\pi/3 }{({\gamma_e}-1)} \biggl[ \frac{P_{ie} \chi^{3\gamma_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_e} \biggl\{ (1-q^3) - b_\xi\biggl(\frac{\rho_e}{\rho_0}\biggr) \frac{P_0}{P_{ie} } \biggl[ (2q^2 - 3q^3 +q^5) + \frac{3}{5} \cdot \frac{\rho_e}{\rho_0} ( 1 - 5q^2 + 5q^3 - q^5 ) \biggr] \biggr\} \, . </math>

Sanity Check

SANITY CHECK

In equilibrium, we will demand that <math>~P_{ie} = P_{ic}</math>. If, in addition, we set <math>~\gamma_e = \gamma_c</math>, the sum of these two terms gives,

<math>~\frac{\mathfrak{S}_A}{E_\mathrm{norm}}</math>	<math>~=</math>	<math> \frac{4\pi/3 }{({\gamma_c}-1)} \biggl[ \frac{P_{0} \chi^{3\gamma_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_c} \biggl\{ \biggl[ q^3 - \biggl( \frac{3b_\xi}{5} \biggr) q^5 \biggr] + \frac{P_{ic} }{P_0} (1-q^3) </math>
		<math> - b_\xi\biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ (2q^2 - 3q^3 +q^5) + \frac{3}{5} \cdot \frac{\rho_e}{\rho_0} ( 1 - 5q^2 + 5q^3 - q^5 ) \biggr] \biggr\} </math>
	<math>~=</math>	<math> \frac{4\pi/3 }{({\gamma_c}-1)} \biggl[ \frac{P_{0} \chi^{3\gamma_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_c} \biggl\{ \biggl[ q^3 - \biggl( \frac{3b_\xi}{5} \biggr) q^5 \biggr] + (1-b_\xi q^2)(1-q^3) </math>
		<math> - b_\xi\biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ (2q^2 - 3q^3 +q^5) + \frac{3}{5} \cdot \frac{\rho_e}{\rho_0} ( 1 - 5q^2 + 5q^3 - q^5 ) \biggr] \biggr\} </math>
	<math>~=</math>	<math> \frac{4\pi/3 }{({\gamma_c}-1)} \biggl[ \frac{P_{0} \chi^{3\gamma_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_c} \biggl\{ 1-b_\xi q^2\biggl( 1 - \frac{2}{5} q^3 \biggr) - b_\xi\biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ (2q^2 - 3q^3 +q^5) + \frac{3}{5} \cdot \frac{\rho_e}{\rho_0} ( 1 - 5q^2 + 5q^3 - q^5 ) \biggr] \biggr\} \, . </math>

In equilibrium, we also expect (from the virial theorem) that,

<math>~2S_\mathrm{therm} + W_\mathrm{grav}</math>	<math>~=</math>	<math>~0 </math>
<math>~\Rightarrow~~~~~3(\gamma_c-1) \frac{\mathfrak{S}_A}{E_\mathrm{norm}} </math>	<math>~=</math>	<math>~- \frac{W_\mathrm{grav}}{E_\mathrm{norm}} \, .</math>

This, in turn, implies,

<math> 4\pi \biggl[ \frac{P_{0} }{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{4} \biggl\{ 1-b_\xi q^2\biggl( 1 - \frac{2}{5} q^3 \biggr) - b_\xi\biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ (2q^2 - 3q^3 +q^5) + \frac{3}{5} \cdot \frac{\rho_e}{\rho_0} ( 1 - 5q^2 + 5q^3 - q^5 ) \biggr] \biggr\} </math>

<math> \frac{3}{5}\biggl(\frac{\nu^2}{q} \biggr) f(\nu,q) \, . </math>

Now, if we also set the density ratio, <math>\rho_e/\rho_c = 1</math>, the result should match the structure of an isolated, uniform density sphere. To check this, note that the term inside curly braces on the left-hand-side becomes,

<math>~\biggl\{ ~~~ \biggr\}</math>	<math>~=</math>	<math>~ 1-b_\xi \biggl[ q^2\biggl( 1 - \frac{2}{5} q^3 \biggr) + (2q^2 - 3q^3 +q^5) + \frac{3}{5} \cdot ( 1 - 5q^2 + 5q^3 - q^5 ) \biggr] </math>
	<math>~=</math>	<math>~ 1-b_\xi \biggl[\frac{3}{5} + q^2 \biggl( 1 + 2 - 3 \biggr) + q^3 \biggl( -3 + 3 \biggr) + q^5 \biggl( -\frac{2}{5} + 1 - \frac{3}{5} \biggr) \biggr] </math>
	<math>~=</math>	<math>~ 1- \frac{3}{5} b_\xi \, , </math>

while the function that appears on the right-hand-side becomes,

<math>~f(\nu,q)</math>	<math>~=</math>	<math> 1+\frac{5}{2} \biggl(\frac{1}{q^2}-1 \biggr) +\biggl[ \frac{1}{q^5}-1 + \frac{5}{2} \biggl( 1-\frac{1}{q^2}\biggr)\biggr] </math>
	<math>~=</math>	<math> \frac{1}{q^5} \, . </math>

Hence, the statement of virial equilibrium becomes,

<math> 4\pi \biggl[ \frac{P_{0} }{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{4} \biggl\{ 1- \frac{3}{5} b_\xi \biggr\} </math>

<math> \frac{3}{5}\biggl(\frac{\nu}{q^3} \biggr)^2 </math>

Old Stuff

<math>~\mathfrak{G}</math>	<math>~=</math>	<math>~W_\mathrm{grav} + \mathfrak{S}_A\biggr\|_\mathrm{core} + \mathfrak{S}_A\biggr\|_\mathrm{env} </math>
	<math>~=</math>	<math> ~W_\mathrm{grav} + \biggl[ \frac{2}{3(\gamma_c - 1)} \biggr] S_\mathrm{core} + \biggl[ \frac{2}{3(\gamma_e - 1)} \biggr] S_\mathrm{env} \, . </math>

In addition to the gravitational potential energy, which is naturally written as,

<math>~W_\mathrm{grav}</math>

<math>~- \frac{3}{5} \biggl( \frac{GM_\mathrm{tot}^2}{R} \biggr) \cdot \mathfrak{f}_{WM} \, ,</math>

it seems reasonable to write the separate thermal energy contributions as,

<math>~S_\mathrm{core}</math>	<math>~=</math>	<math> ~\frac{3}{2}\biggl[ M_\mathrm{core} \biggl( \frac{P_{ic}}{\rho_{ic}} \biggr) \biggr] s_\mathrm{core} = \frac{3}{2} \biggl[ \nu M_\mathrm{tot} P_{ic} \biggl( \frac{\rho_{ic}}{\bar\rho} \biggr)^{-1} \biggl( \frac{4\pi R^3}{3 M_\mathrm{tot}} \biggr) \biggr] s_\mathrm{core} = 2\pi R^3 P_{ic} \biggl[ q^3 s_\mathrm{core} \biggr] \, ,</math>
<math>~S_\mathrm{env}</math>	<math>~=</math>	<math> ~\frac{3}{2}\biggl[ M_\mathrm{env} \biggl( \frac{P_{ie}}{\rho_{ie}} \biggr) \biggr] s_\mathrm{env} = \frac{3}{2} \biggl[ (1-\nu) M_\mathrm{tot} P_{ie} \biggl( \frac{\rho_{ie}}{\bar\rho} \biggr)^{-1} \biggl( \frac{4\pi R^3}{3 M_\mathrm{tot}} \biggr) \biggr] s_\mathrm{env} = 2\pi R^3 P_{ie} \biggl[ (1-q^3) s_\mathrm{env} \biggr] \, ,</math>

where the subscript "<math>i</math>" means "at the interface," and <math>~\mathfrak{f}_{WM},</math> <math>~s_\mathrm{core},</math> and <math>~s_\mathrm{env}</math> are dimensionless functions of order unity (all three functions to be determined) akin to the structural form factors used in our examination of isolated polytropes.

While exploring how the free-energy function varies across parameter space, we choose to hold <math>~M_\mathrm{tot}</math> and <math>~K_c</math> fixed. By dimensional analysis, it is therefore reasonable to normalize all energies, length scales, densities and pressures by, respectively,

<math>~E_\mathrm{norm}</math>	<math>~\equiv</math>	<math>~\biggl[ \frac{G^{3(\gamma_c-1)} M_\mathrm{tot}^{5\gamma_c-6}}{K_c} \biggr]^{1/(3\gamma_c -4)} \, ,</math>
<math>~R_\mathrm{norm}</math>	<math>~\equiv</math>	<math>~\biggl[ \biggl( \frac{K_c}{G} \biggr) M_\mathrm{tot}^{\gamma_c-2} \biggr]^{1/(3\gamma_c -4)} \, ,</math>
<math>~\rho_\mathrm{norm}</math>	<math>~\equiv</math>	<math>~\frac{3}{4\pi} \biggl[ \frac{G^3 M_\mathrm{tot}^2}{K_c^3} \biggr]^{1/(3\gamma_c -4)} \, ,</math>
<math>~P_\mathrm{norm}</math>	<math>~\equiv</math>	<math>~\biggl[ \frac{G^{3\gamma_c} M_\mathrm{tot}^{2\gamma_c}}{K_c^4} \biggr]^{1/(3\gamma_c -4)} \, .</math>

As is detailed below — first, here, and via an independent derivation, here — quite generally the expression for the normalized free energy is,

<math>~\mathfrak{G}^* \equiv \frac{\mathfrak{G}}{E_\mathrm{norm}}</math>	<math>~=</math>	<math> - \frac{3}{5} \biggl( \frac{GM_\mathrm{tot}^2}{E_\mathrm{norm}} \biggr) \biggl( \frac{1}{R} \biggr) \cdot \mathfrak{f}_{WM} + \biggl[ \frac{4\pi q^3 s_\mathrm{core} }{3(\gamma_c - 1)} \biggr] \biggl[ \frac{ R^3 P_{ic} }{E_\mathrm{norm}} \biggr] + \biggl[ \frac{4\pi (1-q^3) s_\mathrm{env} }{3(\gamma_e - 1)} \biggr] \biggl[ \frac{ R^3 P_{ie} }{E_\mathrm{norm}} \biggr] </math>
	<math>~=</math>	<math> - \frac{3\cdot \mathfrak{f}_{WM}}{5} \chi^{-1} + \biggl[ \frac{4\pi q^3 s_\mathrm{core} }{3(\gamma_c - 1)} \biggr] \biggl[ \frac{ R_\mathrm{norm}^4 P_\mathrm{norm} }{E_\mathrm{norm} R_\mathrm{norm}} \biggr] \biggl[ \biggl( \frac{P_{ic}}{P_\mathrm{norm}} \biggr) \chi^3 \biggr] </math>
		<math> + \biggl[ \frac{4\pi (1-q^3) s_\mathrm{env} }{3(\gamma_e - 1)} \biggr] \biggl[ \frac{ R_\mathrm{norm}^4 P_\mathrm{norm} }{E_\mathrm{norm} R_\mathrm{norm}} \biggr] \biggl[ \biggl( \frac{P_{ie}}{P_\mathrm{norm}} \biggr) \chi^3 \biggr] </math>
	<math>~=</math>	<math> - \frac{3\cdot \mathfrak{f}_{WM}}{5} \chi^{-1} + \biggl[ \frac{4\pi q^3 s_\mathrm{core} }{3(\gamma_c - 1)} \biggr] \biggl[ \frac{P_{ic} \chi^{3\gamma_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3 -3\gamma_c} + \biggl[ \frac{4\pi (1-q^3) s_\mathrm{env} }{3(\gamma_e - 1)} \biggr] \biggl[ \frac{P_{ie} \chi^{3\gamma_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3 -3\gamma_e} </math>

where we have introduced the parameter, <math>~\nu \equiv M_\mathrm{core}/M_\mathrm{tot}</math>. After defining the normalized (and dimensionless) configurarion radius, <math>~\chi \equiv R/R_\mathrm{norm}</math>, we can write the normalized free energy of a bipolytrope in the following compact form:

<math>~\mathfrak{G}^*</math>

<math>~ -~ 3\mathcal{A} \chi^{-1} - \frac{\mathcal{B}}{(1-\gamma_c)} ~\chi^{3-3\gamma_c} - \frac{\mathcal{C}}{(1-\gamma_e)} ~\chi^{3-3\gamma_e} \, ,</math>

where,

<math>~\mathcal{A}</math>	<math>~\equiv</math>	<math>~\frac{1}{5} \cdot \mathfrak{f}_{WM} \, ,</math>
<math>~\mathcal{B}</math>	<math>~\equiv</math>	<math> \biggl( \frac{4\pi }{3} \biggr) q^3 s_\mathrm{core} \biggl[ \frac{P_{ic} \chi^{3\gamma_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \, , </math>
<math>~\mathcal{C}</math>	<math>~\equiv</math>	<math>\biggl( \frac{4\pi }{3} \biggr) (1-q^3) s_\mathrm{env} \biggl[ \frac{P_{ie} \chi^{3\gamma_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \, .</math>

As is further detailed below, the second expression for the coefficient, <math>~\mathcal{C}</math>, ensures that the pressure at the "surface" of the core matches the pressure at the "base" of the envelope; but it should only be employed after an equilibrium radius, <math>~\chi_\mathrm{eq}</math>, has been identified by locating an extremum in the free energy.

Simplest Bipolytrope

Familiar Setup

As has been shown in an accompanying presentation, for an <math>~(n_c, n_e) = (0, 0)</math> bipolytrope,

<math>~\mathfrak{f}_{WM}</math>	<math>~\equiv</math>	<math>~\frac{\nu^2}{q} \cdot f \, ,</math>
<math>~s_\mathrm{core}</math>	<math>~\equiv</math>	<math> 1 + \Lambda_\mathrm{eq} \, , </math>
<math>~(1-q^3) s_\mathrm{env}</math>	<math>~\equiv</math>	<math> (1-q^3) + \Lambda\biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] \, ,</math>

and where (see, for example, in the context of its original definition, or another, separate derivation),

<math>~\Lambda_\mathrm{eq} </math>

<math>

\frac{1}{5} \biggl( \frac{\nu}{q} \biggr) \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{1-\gamma_c}

\chi_\mathrm{eq}^{3\gamma_c - 4} </math>

<math>\frac{2}{5(g^2-1)} = \biggl\{ \frac{5}{2}\biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( 1-q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] \biggr\}^{-1} \, ,</math>

and where (see the associated discussion of relevant mass integrals),

<math> \frac{\rho_c}{\bar\rho} = \frac{\nu}{q^3} \, ; ~~~~~ \frac{\rho_e}{\bar\rho} = \frac{1-\nu}{1-q^3} \, ; ~~~~~ \frac{\rho_e}{\rho_c} = \frac{q^3(1-\nu)}{\nu (1-q^3)} ~~\Rightarrow ~~~ \frac{q^3}{\nu} = \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3 \, . </math>

Cleaner Virial Presentation

In an effort to show the similarity in structure among the several energy terms, we have also found it useful to write their expressions in the following forms:

<math>~W_\mathrm{grav}</math>	<math>~=</math>	<math>~- \frac{3}{5} \biggl( \frac{GM_\mathrm{tot}^2}{R} \biggr) \frac{\nu^2}{q} \cdot f = - 4\pi P_i R^3 \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q} \cdot \lambda_i f \, ,</math>
<math>~S_\mathrm{core}</math>	<math>~=</math>	<math>~2\pi P_{ic} R^3 \biggl[ q^3 + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q} \cdot \lambda_{ic} \biggr] \, ,</math>
<math>~S_\mathrm{env}</math>	<math>~=</math>	<math>~2\pi P_{ie} R^3 \biggl[ (1-q^3) + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q} \cdot \lambda_{ie} \mathfrak{F} \biggr] \, ,</math>

where (see an associated discussion or the original derivation),

<math> f\biggl(q, \frac{\rho_e}{\rho_c}\biggr) = 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \frac{1}{q^5} \biggl[ (q^3 - q^5 ) + \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl( \frac{2}{5} - q^3 + \frac{3}{5}q^5 \biggr) \biggr] \, , </math>

and where,

<math>~\lambda_i</math>

<math>~\equiv</math>

<math>~\frac{GM_\mathrm{tot}^2}{R^4 P_i} \, ,</math>

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

<math>~\equiv</math>

<math>~ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \frac{1}{q^5} \biggl[ (-2q^2 + 3q^3 - q^5) + \frac{3}{5} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-1 +5q^2 - 5q^3 + q^5) \biggr] </math>

<math>~ \frac{1}{\lambda_{ie}} \biggl( \frac{2^2 \cdot 5\pi}{3} \biggr) \frac{q(1-q^3)}{\nu^2} (s_\mathrm{env} -1)

\, .</math>

This also means that the three key terms used as shorthand notation in the above expressions for the three energy terms have the following definitions:

<math>~\mathfrak{f}_{WM}</math>	<math>~\equiv</math>	<math>~\frac{\nu^2}{q} \cdot f \, ,</math>
<math>~s_\mathrm{core}</math>	<math>~\equiv</math>	<math> 1 + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q^4} \cdot \lambda_{ic} \, , </math>
<math>~s_\mathrm{env}</math>	<math>~\equiv</math>	<math> 1 + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q(1-q^3)} \cdot \lambda_{ie} \mathfrak{F} \, ,</math>

Hence, if all the interface pressures are equal — that is, if <math>~P_i = P_{ic} = P_{ie}</math> and, hence also, <math>~\lambda_{i} = \lambda_{ic} = \lambda_{ie}</math> — then the total thermal energy is,

<math>~S_\mathrm{tot} = S_\mathrm{core} + S_\mathrm{env}</math>

<math> 2\pi P_{i} R^3 \biggl[ 1 + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q} \cdot \lambda_{i} (1+\mathfrak{F}) \biggr] \, ; </math>

and the virial is,

<math>~2S_\mathrm{tot} + W_\mathrm{grav}</math>

<math> 4\pi P_{i} R^3 \biggl[ 1 + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q} \cdot \lambda_{i} (1+\mathfrak{F} - f ) \biggr] \, . </math>

The virial should sum to zero in equilibrium, which means,

<math>~\frac{1}{\lambda_i} \biggr\|_\mathrm{eq}</math>	<math>~=</math>	<math> \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q} \cdot (f - 1- \mathfrak{F} ) </math>
<math>\Rightarrow ~~~~ \biggl[ \biggl( \frac{2^2\cdot 5\pi}{3} \biggr) \frac{q}{\nu^2} \biggr] \frac{R_\mathrm{eq}^4 P_i}{GM_\mathrm{tot}^2} </math>	<math>~=</math>	<math> f - 1- \mathfrak{F} </math>
<math>\Rightarrow ~~~~ \biggl( \frac{\rho_e}{\rho_c} \biggr)^{-1} \biggl[ \biggl( \frac{2^3\pi}{3} \biggr) \frac{q^6}{\nu^2} \biggr] \frac{R_\mathrm{eq}^4 P_i}{GM_\mathrm{tot}^2} </math>	<math>~=</math>	<math> \biggl[ (q^3 - q^5 )+ \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl( \frac{2}{5} - q^3 + \frac{3}{5}q^5 \biggr) \biggr] - \biggl[ (-2q^2 + 3q^3 - q^5) + \biggl( \frac{\rho_e}{\rho_c}\biggr) (-\frac{3}{5} +3q^2 - 3q^3 + \frac{3}{5} q^5) \biggr] </math>
	<math>~=</math>	<math> 2q^2(1-q) + \biggl( \frac{\rho_e}{\rho_c}\biggr) (1 -3q^2 + 2q^3 ) </math>
	<math>~=</math>	<math> q^2 \biggl( \frac{\rho_e}{\rho_c} \biggr)^{-1} (g^2-1) </math>
<math>\Rightarrow~~~~ \frac{1}{\lambda_i} \biggr\|_\mathrm{eq}</math>	<math>~=</math>	<math> \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^2 q^2 (g^2-1) \, . </math>

Shift to Central Pressure Normalization

Let's rework the definition of <math>~\lambda_i</math> in two ways: (1) Normalize <math>~R_\mathrm{eq}</math> to <math>~R_\mathrm{norm}</math> and normalize the pressure to <math>~P_\mathrm{norm}</math>; (2) shift the referenced pressure from the pressure at the interface <math>~(P_i)</math> to the central pressure <math>~(P_0)</math>, because it is <math>~P_0</math> that is directly related to <math>~K_c</math> and <math>~\rho_c</math>; specifically, <math>P_0 = K_c \rho_c^{\gamma_c}</math>. Appreciating that, in equilibrium,

<math>~P_0 - q^2 \Pi_\mathrm{eq} = K_c \rho_c^{\gamma_c} - \frac{3}{2^3 \pi} \biggl( \frac{GM_\mathrm{tot}^2}{R_\mathrm{eq}^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) q^2 \, ,</math>

the left-hand-side of the last expression, above, can be rewritten as,

<math>~\frac{1}{\lambda_i} \biggr\|_\mathrm{eq}</math>	<math>~\equiv</math>	<math>~\frac{R_\mathrm{eq}^4 P_i}{GM_\mathrm{tot}^2} </math>
	<math>~=</math>	<math>~\frac{R_\mathrm{eq}^4}{GM_\mathrm{tot}^2} \biggl[ P_0 - \frac{3}{2^3 \pi} \biggl( \frac{GM_\mathrm{tot}^2}{R_\mathrm{eq}^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) q^2\biggr] </math>
	<math>~=</math>	<math>~\frac{R_\mathrm{eq}^4 P_0}{GM_\mathrm{tot}^2} - \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^2 q^2 \, .</math>

Hence, the virial equilibrium condition gives,

<math> \frac{R_\mathrm{eq}^4 P_0}{GM_\mathrm{tot}^2} - \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^2 q^2 </math>	<math>~=</math>	<math> \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^2 q^2 (g^2-1) </math>
<math> \Rightarrow ~~~~ \frac{R_\mathrm{eq}^4 P_0}{GM_\mathrm{tot}^2} </math>	<math>~=</math>	<math> \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{2} q^2 g^2 \, . </math>

This result precisely matches the result obtained via the detailed force-balanced conditions imposed through hydrostatic equilibrium.

Adopting our new variable normalizations and realizing, in particular, that,

<math>~R_\mathrm{norm}^4 P_\mathrm{norm}</math>

<math>~GM_\mathrm{tot}^2 \, ,</math>

the expression alternatively can be rewritten as,

<math>~\frac{1}{\lambda_i} \biggr\|_\mathrm{eq}</math>	<math>~\equiv</math>	<math>~\frac{R_\mathrm{eq}^4 P_i}{GM_\mathrm{tot}^2} = \chi_\mathrm{eq}^4 \biggl( \frac{P_i}{P_\mathrm{norm}} \biggr) </math>
	<math>~=</math>	<math> \chi_\mathrm{eq}^4 \biggl\{ \frac{K_c \rho_c^{\gamma_c}}{P_\mathrm{norm}} - \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3}\biggr]^2 q^2 \biggl( \frac{GM_\mathrm{tot}^2}{R_\mathrm{eq}^4P_\mathrm{norm}} \biggr) \biggr\} </math>
	<math>~=</math>	<math> \chi_\mathrm{eq}^4 \biggl\{ \frac{K_c }{P_\mathrm{norm}} \biggl[ \frac{\rho_c}{\bar\rho} \biggl( \frac{3M_\mathrm{tot}}{4\pi R_\mathrm{norm}^3} \biggr) \chi_\mathrm{eq}^{-3} \biggr]^{\gamma_c} - \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3}\biggr]^2 q^2 \chi_\mathrm{eq}^{-4} \biggr\} </math>
	<math>~=</math>	<math> \chi_\mathrm{eq}^{4-3\gamma_c} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} \frac{K_c }{P_\mathrm{norm}} \biggl( \frac{M_\mathrm{tot}^{\gamma_c}}{R_\mathrm{norm}^{3\gamma_c}} \biggr) - \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3}\biggr]^2 q^2 </math>
	<math>~=</math>	<math> \chi_\mathrm{eq}^{4-3\gamma_c} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} - \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3}\biggr]^2 q^2 \, . </math>

Normalized in this manner, the virial equilibrium (as well as the hydrostatic balance) condition gives,

<math> \chi_\mathrm{eq}^{4-3\gamma_c} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} - \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3}\biggr]^2 q^2 </math>	<math>~=</math>	<math> \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^2 q^2 (g^2-1) </math>
<math> \Rightarrow ~~~~ \chi_\mathrm{eq}^{4-3\gamma_c} </math>	<math>~=</math>	<math> \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{2-\gamma_c} q^2 g^2 \, . </math>

Free-Energy Coefficients

Therefore, for an <math>~(n_c, n_e) = (0, 0)</math> bipolytrope, the coefficients in the normalized free-energy function are,

<math>~\mathcal{A}</math>	<math>~=</math>	<math> ~\frac{\nu^2}{5q} \cdot f = \frac{1}{5} \biggl( \frac{\nu}{q^3} \biggr)^2 \biggl[ q^5 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) q^3 (1 - q^2 ) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl( 1 - \frac{5}{2} q^3 + \frac{3}{2}q^5 \biggr) \biggr] \, ,</math>
<math>~\mathcal{B}</math>	<math>~\equiv</math>	<math> \biggl( \frac{4\pi }{3} \biggr) q^3 s_\mathrm{core} \biggl[ \frac{P_{ic} }{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_c} =\biggl( \frac{4\pi }{3} \biggr) q^3 s_\mathrm{core} \biggl[ \frac{1}{\lambda_{ic}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_c-4} </math>
	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ 1 + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q^4} \cdot \lambda_{ic} \biggr] \biggl[ \frac{1}{\lambda_{ic}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_c-4} =\biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ \frac{1}{\lambda_{ic}} + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q^4} \biggr] \chi_\mathrm{eq}^{3\gamma_c-4} </math>
	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl\{ \chi_\mathrm{eq}^{4-3\gamma_c} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} - \frac{2\pi}{3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3}\biggr]^2 q^2 + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q^4} \biggr\} \chi_\mathrm{eq}^{3\gamma_c-4} </math>
	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl\{ \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} + \chi_\mathrm{eq}^{3\gamma_c-4} \biggl[\frac{3}{2^2\cdot 5\pi} - \frac{3}{2^3\pi} \biggr]\frac{\nu^2}{q^4} \biggr\} </math>
	<math>~=</math>	<math> \biggl\{ \nu \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c-1} - \chi_\mathrm{eq}^{3\gamma_c-4} \biggl( \frac{3^2}{2^3\cdot 5\pi} \biggr) \frac{\nu^2}{q^4} \biggl( \frac{4\pi }{3} \biggr) q^3 \biggr\} = \nu \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c-1} - \chi_\mathrm{eq}^{3\gamma_c-4} \biggl( \frac{3}{10} \biggr) \frac{\nu^2}{q} </math>
<math>~\mathcal{C}</math>	<math>~\equiv</math>	<math>\biggl( \frac{4\pi }{3} \biggr) (1-q^3) s_\mathrm{env} \biggl[ \frac{P_{ie} \chi^{3\gamma_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \, .</math>

Note that, because <math>~P_{ie} = P_{ic}</math> in equilibrium, the ratio of coefficients,

<math>~\frac{\mathcal{C}}{\mathcal{B}}</math>	<math>~=</math>	<math>~\chi_\mathrm{eq}^{3(\gamma_e - \gamma_c)}\biggl\{\frac{(1-q^3) s_\mathrm{env}}{q^3 s_\mathrm{core}} \biggr\}</math>
<math>\Rightarrow ~~~~ \chi_\mathrm{eq}^{3(\gamma_c - \gamma_e)} \biggl( \frac{\mathcal{C}}{\mathcal{B}} \biggr)</math>	<math>~=</math>	<math> \frac{ 20\pi q(1-q^3) \lambda_i^{-1} + 3\nu^2 \mathfrak{F} }{ 20\pi q^4 \lambda_i^{-1} + 3\nu^2 } \, . </math>

The equilibrium condition is,

<math>\frac{\mathcal{A}}{\mathcal{B} + \mathcal{C}^'} = \chi_\mathrm{eq}^{4-3\gamma_c} \, ,</math>

where,

<math> \mathcal{C}^' \equiv \mathcal{C} \chi_\mathrm{eq}^{3(\gamma_c-\gamma_e)} \, . </math>

More General Derivation of Free-Energy Coefficients B and C

Keep in mind that, generally,

<math>GM_\mathrm{tot}^2</math>	<math>~=</math>	<math>~R_\mathrm{norm}^4 P_\mathrm{norm} = E_\mathrm{norm} R_\mathrm{norm} \, ;</math>
<math>~\frac{1}{\lambda_i}</math>	<math>~\equiv</math>	<math> \frac{R^4 P_i}{GM_\mathrm{tot}^2} = \biggl( \frac{P_i}{P_\mathrm{norm}} \biggr) \chi^4 </math> … and, note that … <math> \frac{1}{\Lambda} = \biggl( \frac{3\cdot 5}{2^2\pi} \biggr) \frac{1}{\lambda_i} \cdot \frac{1}{q^2 \sigma^2} \, ; </math>
<math>~\frac{\Pi}{P_\mathrm{norm}}</math>	<math>~=</math>	<math> \frac{3}{2^3\pi} \biggl( \frac{GM_\mathrm{tot}^2}{P_\mathrm{norm} R^4} \biggr) \frac{\nu^2}{q^6} =\biggl( \frac{2\pi}{3} \biggr) \sigma^2 \chi^{-4} \, ; </math>
<math>~\frac{K_c \rho_c^{\gamma_c} }{P_\mathrm{norm}}</math>	<math>~=</math>	<math> \frac{K_c}{P_\mathrm{norm}} \biggl( \frac{\rho_c}{\bar\rho} \biggr)^{\gamma_c} \biggl[ \frac{3M_\mathrm{tot}}{4\pi R^3} \biggr]^{\gamma_c} = \frac{K_c M_\mathrm{tot}^{\gamma_c} }{R_\mathrm{norm}^{3\gamma_c} P_\mathrm{norm}} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} \chi^{-3\gamma_c} = \sigma^{\gamma_c} \chi^{-3\gamma_c} \, , </math>

where we have introduced the notation,

<math> \sigma \equiv \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \, . </math>

So, the free-energy coefficient,

<math>~\mathcal{B}</math>

<math> \biggl( \frac{4\pi }{3} \biggr) q^3 s_\mathrm{core} \biggl[ \frac{P_{ic} }{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_c} =\biggl( \frac{4\pi }{3} \biggr) q^3 s_\mathrm{core} \biggl[ \frac{1}{\lambda_{ic}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_c-4} </math>

<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ 1 + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q^4} \cdot \lambda_{ic} \biggr]_\mathrm{eq} \biggl[ \frac{1}{\lambda_{ic}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_c-4} =\biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ \frac{1}{\lambda_{ic}} + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q^4} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_c-4} </math>

<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ \chi_\mathrm{eq}^4 \biggl( \frac{P_{ic}}{P_\mathrm{norm}} \biggr)_\mathrm{eq} +

\biggl( \frac{4\pi}{3\cdot 5} \biggr) q^2 \sigma^2 \biggr] \chi_\mathrm{eq}^{3\gamma_c-4} \, .

</math>

And the free-energy coefficient,

<math>~\mathcal{C}</math>	<math>~\equiv</math>	<math> \biggl( \frac{4\pi }{3} \biggr) (1-q^3) s_\mathrm{env} \biggl[ \frac{P_{ie} \chi^{3\gamma_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} = \biggl( \frac{4\pi }{3} \biggr) (1-q^3) s_\mathrm{env} \biggl[ \frac{1 }{\lambda_{ie}} \biggr]_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_e-4} </math>
	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) (1-q^3) \biggl\{ \frac{1 }{\lambda_{ie}} + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q(1-q^3)} \cdot \mathfrak{F} \biggr\}_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_e-4} </math>
	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) \biggl\{ (1-q^3) \chi_\mathrm{eq}^4 \biggl( \frac{P_{ie}}{P_\mathrm{norm}} \biggr)_\mathrm{eq}+ \biggl( \frac{2\pi}{3} \biggr) \sigma^2 \biggl[ \frac{2}{5} q^5 \mathfrak{F} \biggr] \biggr\}_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_e-4} \, . </math>

OLD DERIVATION

<math>P_{ic} = K_c \rho_c^{\gamma_c}</math>

NEW DERIVATION

<math>P_{ic} = P_0 - q^2\Pi = K_c \rho_c^{\gamma_c} - q^2\Pi</math>

… therefore …

<math>~\mathcal{B}_\mathrm{OLD}</math>

<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ \sigma^{\gamma_c} \chi_\mathrm{eq}^{4 - 3\gamma_c} +

\biggl( \frac{4\pi}{3\cdot 5} \biggr) q^2 \sigma^2 \biggr] \chi_\mathrm{eq}^{3\gamma_c-4}

</math>

<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ \sigma^{\gamma_c} +

\biggl( \frac{4\pi}{3\cdot 5} \biggr) q^2 \sigma^2 \chi_\mathrm{eq}^{3\gamma_c-4} \biggr]

</math>

<math>~\mathcal{B}_\mathrm{NEW}</math>	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ \sigma^{\gamma_c} \chi_\mathrm{eq}^{4 - 3\gamma_c} - \biggl( \frac{2\pi}{3} \biggr) q^2 \sigma^2 + \biggl( \frac{4\pi}{3\cdot 5} \biggr) q^2 \sigma^2 \biggr] \chi_\mathrm{eq}^{3\gamma_c-4} </math>
	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) q^3 \biggl[ \sigma^{\gamma_c} - \biggl( \frac{2\pi}{5} \biggr) q^2 \sigma^2 \chi_\mathrm{eq}^{3\gamma_c-4} \biggr] </math>

… and, enforcing in equilibrium <math>~P_{ie} = P_{ic}</math> …

<math>~\mathcal{C}_\mathrm{OLD}</math>	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) \biggl\{ (1-q^3) \biggl[ \sigma^{\gamma_c} \chi_\mathrm{eq}^{4-3\gamma_c} \biggr] + \biggl( \frac{2\pi}{3} \biggr) \sigma^2 \biggl[ \frac{2}{5} q^5 \mathfrak{F} \biggr] \biggr\}_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_e-4} </math>
	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) \biggl[ (1-q^3) \sigma^{\gamma_c} + \biggl( \frac{2\pi}{3} \biggr) \sigma^2 \biggl( \frac{2}{5} q^5 \mathfrak{F} \biggr) \chi_\mathrm{eq}^{3\gamma_e-4}\biggr] </math>

<math>~\mathcal{C}_\mathrm{NEW}</math>	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) \biggl\{ (1-q^3) \biggl[ \sigma^{\gamma_c} \chi_\mathrm{eq}^{4-3\gamma_c} - \biggl( \frac{2\pi}{3} \biggr) q^2 \sigma^2 \biggr] + \biggl( \frac{2\pi}{3} \biggr) \sigma^2 \biggl[ \frac{2}{5} q^5 \mathfrak{F} \biggr] \biggr\}_\mathrm{eq} \chi_\mathrm{eq}^{3\gamma_e-4} </math>
	<math>~=</math>	<math> \biggl( \frac{4\pi }{3} \biggr) \biggl\{ (1-q^3) \sigma^{\gamma_c} \chi_\mathrm{eq}^{3\gamma_e - 3\gamma_c} + \biggl( \frac{2\pi}{3} \biggr) \sigma^2 \biggl[ \biggl( \frac{2}{5} q^5 \mathfrak{F} \biggr) - q^2 (1-q^3) \biggr] \chi_\mathrm{eq}^{3\gamma_e-4} \biggr\} </math>

… and, also …

<math>~\frac{1}{\lambda_i} \biggr\|_\mathrm{eq}</math>	<math>~=</math>	<math> \sigma^{\gamma_c} \chi_\mathrm{eq}^{4 - 3\gamma_c} </math>
<math>\Rightarrow ~~~~ \frac{1}{\Lambda_\mathrm{eq}} </math>	<math>~=</math>	<math> \frac{1}{q^2 } \biggl( \frac{3\cdot 5}{2^2\pi} \biggr) \sigma^{\gamma_c - 2} \chi_\mathrm{eq}^{4 - 3\gamma_c} </math>
	<math>~=</math>	<math> \biggl( \frac{5q}{\nu} \biggr) \sigma^{\gamma_c - 1} \chi_\mathrm{eq}^{4 - 3\gamma_c} </math>

<math>~\frac{1}{\lambda_i} \biggr\|_\mathrm{eq}</math>	<math>~=</math>	<math> \sigma^{\gamma_c} \chi_\mathrm{eq}^{4 - 3\gamma_c} - \biggl( \frac{2\pi}{3} \biggr) q^2 \sigma^2 </math>
<math>\Rightarrow ~~~~ \frac{1}{\Lambda_\mathrm{eq}} </math>	<math>~=</math>	<math> \frac{1}{q^2 \sigma^2} \biggl( \frac{3\cdot 5}{2^2\pi} \biggr) \biggl[ \sigma^{\gamma_c} \chi_\mathrm{eq}^{4 - 3\gamma_c} - \biggl( \frac{2\pi}{3} \biggr) q^2 \sigma^2 \biggr] </math>
	<math>~=</math>	<math> \biggl( \frac{5q}{\nu} \biggr) \sigma^{\gamma_c - 1} \chi_\mathrm{eq}^{4 - 3\gamma_c} - \frac{5}{2} </math>

Extrema

Extrema in the free energy occur when,

<math>~\mathcal{A}</math>

<math>~\mathcal{B} \chi_\mathrm{eq}^{4-3\gamma_c} + \mathcal{C} \chi_\mathrm{eq}^{4-3\gamma_e} \, .</math>

Also, as stated above, because <math>~P_{ie} = P_{ic}</math> in equilibrium, the ratio of coefficients,

<math>~\frac{\mathcal{C}}{\mathcal{B}}</math>

<math>~\chi_\mathrm{eq}^{3(\gamma_e - \gamma_c)}\biggl[ \frac{(1-q^3) s_\mathrm{env}}{q^3 s_\mathrm{core}} \biggr] \, .</math>

When put together, these two relations imply,

<math>~\mathcal{A}</math>	<math>~=</math>	<math>~\mathcal{B} \chi_\mathrm{eq}^{4-3\gamma_c} + \chi_\mathrm{eq}^{4-3\gamma_c} \mathcal{B} \biggl[ \frac{(1-q^3) s_\mathrm{env}}{q^3 s_\mathrm{core}} \biggr] </math>
	<math>~=</math>	<math>~\mathcal{B} \chi_\mathrm{eq}^{4-3\gamma_c} \biggl[ 1+ \frac{(1-q^3) s_\mathrm{env}}{q^3 s_\mathrm{core}} \biggr] \, .</math>

But the definition of <math>~\mathcal{B}</math> gives,

<math>~\mathcal{B} \chi_\mathrm{eq}^{4-3\gamma_c}</math>

<math> \biggl( \frac{4\pi }{3} \biggr) q^3 s_\mathrm{core} \biggl[ \frac{1}{\lambda_{ic}} \biggr]_\mathrm{eq} \, . </math>

Hence, extrema occur when,

<math>~\mathcal{A}</math>	<math>~=</math>	<math>~\biggl( \frac{4\pi }{3} \biggr) q^3 s_\mathrm{core} \biggl[ \frac{1}{\lambda_{ic}} \biggr]_\mathrm{eq} \biggl[ 1+ \frac{(1-q^3) s_\mathrm{env}}{q^3 s_\mathrm{core}} \biggr] </math>
<math>\Rightarrow ~~~~ \biggl( \frac{3}{2^2 \cdot 5\pi } \biggr) \frac{\nu^2}{q} \cdot f</math>	<math>~=</math>	<math> \biggl[ \frac{1}{\lambda_{ic}} \biggr]_\mathrm{eq} \biggl[ q^3 s_\mathrm{core} + (1-q^3) s_\mathrm{env} \biggr] </math>
	<math>~=</math>	<math> \frac{q^3}{[\lambda_i]_\mathrm{eq}} + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q} + \frac{(1-q^3)}{[\lambda_i]_\mathrm{eq}} + \biggl( \frac{3}{2^2\cdot 5\pi} \biggr) \frac{\nu^2}{q} \cdot \mathfrak{F} </math>
<math>\Rightarrow ~~~~\biggl[ \frac{1}{\lambda_{ic}} \biggr]_\mathrm{eq} </math>	<math>~=</math>	<math> \biggl( \frac{3}{2^2 \cdot 5\pi } \biggr) \frac{\nu^2}{q} \cdot (f - 1 - \mathfrak{F}) </math>
	<math>~=</math>	<math> \biggl( \frac{2\pi}{3 } \biggr) \sigma^{2} q^2 (g^2 - 1 ) \, . </math>

In what follows, keep in mind that,

<math>~\chi_\mathrm{eq}^{4-3\gamma_c}</math>	<math>~=</math>	<math> \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{norm}} \biggr)^{4-3\gamma_c} = R_\mathrm{eq}^{4-3\gamma_c} \biggl( \frac{K_c}{G} \biggr) M_\mathrm{tot}^{\gamma_c-2} \, ; </math>
<math>~K_c \rho_c^{\gamma_c}</math>	<math>~=</math>	<math> K_c \biggl( \frac{\rho_c}{\bar\rho} \biggr)^{\gamma_c} \biggl[ \frac{3M_\mathrm{tot}}{4\pi R^3} \biggr]^{\gamma_c} = K_c \sigma^{\gamma_c} M_\mathrm{tot}^{\gamma_c} R^{-3\gamma_c} \, ; </math>
<math>~\Pi</math>	<math>~=</math>	<math> \frac{3}{2^3 \pi} \biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \frac{\nu^2}{q^6} = \frac{2\pi}{3} \biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \sigma^2 \, . </math>

OLD DERIVATION

<math>P_{i} = K_c \rho_c^{\gamma_c}</math>

<math>\Rightarrow ~~~~ K_c = P_{i} \sigma^{-\gamma_c} M_\mathrm{tot}^{-\gamma_c} R^{+3\gamma_c} </math>

NEW DERIVATION

<math>P_0 = K_c \rho_c^{\gamma_c} </math>

<math>\Rightarrow ~~~~ K_c = P_0 \sigma^{-\gamma_c} M_\mathrm{tot}^{-\gamma_c} R^{+3\gamma_c} </math>

… hence, as derived in the above table …

<math>~\frac{1}{\lambda_i} \biggr|_\mathrm{eq}</math>

<math> \sigma^{\gamma_c} \chi_\mathrm{eq}^{4 - 3\gamma_c} </math>

<math>~\frac{1}{\lambda_i} \biggr|_\mathrm{eq}</math>

<math> \sigma^{\gamma_c} \chi_\mathrm{eq}^{4 - 3\gamma_c} - \biggl( \frac{2\pi}{3} \biggr) q^2 \sigma^2 </math>

… which, when combined with the condition that identifies extrema, gives …

<math>~\chi_\mathrm{eq}^{4 - 3\gamma_c}</math>	<math>~=</math>	<math> \biggl( \frac{2\pi}{3 } \biggr) \sigma^{2-\gamma_c} q^2 (g^2 - 1 ) </math>
<math>~\Rightarrow ~~~~ R_\mathrm{eq}^{4-3\gamma_c} \biggl( \frac{K_c}{G} \biggr) M_\mathrm{tot}^{\gamma_c-2} </math>	<math>~=</math>	<math> \biggl( \frac{2\pi}{3 } \biggr) \sigma^{2-\gamma_c} q^2 (g^2 - 1 ) </math>
<math>\Rightarrow ~~~~ \frac{ R_\mathrm{eq}^{4} P_i }{GM_\mathrm{tot}^{2} } </math>	<math>~=</math>	<math> \biggl( \frac{2\pi}{3 } \biggr) \sigma^{2} q^2 (g^2-1) </math>

<math>~\sigma^{\gamma_c}\chi_\mathrm{eq}^{4 - 3\gamma_c} - \biggl( \frac{2\pi}{3} \biggr) q^2 \sigma^2</math>	<math>~=</math>	<math> \biggl( \frac{2\pi}{3 } \biggr) \sigma^2 q^2 (g^2 - 1 ) </math>
<math>\Rightarrow ~~~~ \chi_\mathrm{eq}^{4 - 3\gamma_c} </math>	<math>~=</math>	<math> \biggl( \frac{2\pi}{3 } \biggr) \sigma^{2-\gamma_c} q^2 g^2 </math>
<math>\Rightarrow ~~~~ R_\mathrm{eq}^{4-3\gamma_c} \biggl( \frac{K_c}{G} \biggr) M_\mathrm{tot}^{\gamma_c-2} </math>	<math>~=</math>	<math> \biggl( \frac{2\pi}{3 } \biggr) \sigma^{2-\gamma_c} q^2 g^2 </math>
<math>\Rightarrow ~~~~ \frac{ R_\mathrm{eq}^{4} P_0 }{GM_\mathrm{tot}^{2} } </math>	<math>~=</math>	<math> \biggl( \frac{2\pi}{3 } \biggr) \sigma^{2} q^2 g^2 </math>

These are consistent results because they result in the detailed force-balance relation, <math>P_0 - P_i = q^2 \Pi_\mathrm{eq} \, .</math>

Examples

Identification of Local Extrema in Free Energy
<math>~\nu</math>	<math>~q</math>	<math>~ \frac{\rho_e}{\rho_c} </math>	<math>~f\biggl(q, \frac{\rho_e}{\rho_c} \biggr)</math>	<math>~g^2\biggl(q, \frac{\rho_e}{\rho_c} \biggr)</math>	<math>~\Lambda_\mathrm{eq}</math>	<math>~\chi_\mathrm{eq}</math>	<math>~\mathcal{A}</math>	<math>~\mathcal{B}</math>	<math>~\mathcal{C}</math>	MIN/MAX
<math>~0.2</math>	<math>~9^{-1/3} = 0.48075</math>	<math>~0.5</math>	<math>~12.5644</math>	<math>~2.091312</math>	<math>~0.366531</math>	<math>~0.037453</math>	<math>~0.2090801</math>	<math>~0.2308269</math>	<math>~2.06252 \times 10^{-4}</math>	MIN
<math>~0.4</math>	<math>~4^{-1/3} = 0.62996</math>	<math>~0.5</math>	<math>~4.21974</math>	<math>~1.56498</math>	<math>~0.707989</math>	<math>~0.0220475</math>	<math>~0.2143496</math>	<math>~0.5635746</math>	<math>~4.4626 \times 10^{-5}</math>	MIN
<math>~0.473473</math>	<math>~0.681838</math>	<math>~0.516107</math>			<math>~0.462927</math>	<math>~0.08255</math>	"	"	"	MAX
<math>~0.5</math>	<math>~3^{-1/3} = 0.693361</math>	<math>~0.5</math>	<math>~2.985115</math>	<math>~1.42334</math>	<math>~0.9448663</math>	<math>~0.0152116</math>	<math>~0.2152641</math>	<math>~0.791882</math>	<math>~1.5464 \times 10^{-5}</math>	MIN
<math>~0.559839</math>	<math>~0.729581</math>	<math>~0.499188</math>			<math>~0.75089</math>	<math>~0.032196</math>	"	"	"	MAX

Free Energy Extrema when: <math>~~~~~\frac{\rho_e}{\rho_c} = \frac{1}{2} ~~~~\Rightarrow~~~~ q^3 = \frac{\nu}{2-\nu}</math>
<math>~\nu</math>	<math>~q</math>	<math>~ \frac{\rho_e}{\rho_c} </math>	<math>~f\biggl(q, \frac{\rho_e}{\rho_c} \biggr)</math>	<math>~g^2\biggl(q, \frac{\rho_e}{\rho_c} \biggr)</math>	<math>~\frac{1}{\lambda_{i}}\biggr\|_\mathrm{eq}</math>	<math>~\chi_\mathrm{eq}</math>	<math>~\mathcal{A}</math>	<math>~\mathcal{B}_\mathrm{NEW}</math>	<math>~\mathcal{C}_\mathrm{NEW}</math>	<math>~G^*</math>	Stability	MIN/MAX
<math>~0.5</math>	<math>~\biggl( \frac{1}{3} \biggr)^{1/3}</math>	<math>~0.5</math>	<math>~2.985115</math>	<math>~1.423340</math>	<math>~0.05466039</math>	<math>~0.3152983</math>	<math>~0.21526406</math>	<math>~0.23552725</math>	<math>~6.643899 \times 10^{-3}</math>	<math>~+ 0.5176146</math>	<math>~+0.429245</math>	MIN
…	…	…				<math>~0.6674</math>	"	"	"	<math>~+0.55572115</math>		MAX
<math>~0.6</math>	<math>~\biggl( \frac{3}{7} \biggr)^{1/3}</math>	<math>~0.5</math>	<math>~2.2507129</math>	<math>~1.31282895</math>	<math>~0.04160318</math>	<math>~0.3411545</math>	<math>~0.21493717</math>	<math>~0.26165939</math>	<math>~5.208750 \times 10^{-3}</math>	<math>~+ 0.73532249</math>	<math>~+0.0935217</math>	MIN
…	…	…				<math>~0.431745</math>	"	"	"	<math>~+0.7367797</math>		MAX
<math>~0.7</math>	<math>~\biggl( \frac{7}{13} \biggr)^{1/3}</math>	<math>~0.5</math>	<math>~1.7707809</math>	<math>~1.2209446</math>	<math>~0.029500</math>	<math>~0.3589388</math>	<math>~0.21330744</math>	<math>~0.28172532</math>	<math>~3.389793 \times 10^{-3}</math>	<math>~+ 0.8953395</math>	<math>~- 0.0767108</math>	MAX
…	…	…				<math>~0.270615</math>	"	"	"	<math>~+0.89227216</math>		MIN
System should be stable (with free energy minimum) if: <math>~~~~~\frac{( \gamma_e - \frac{4}{3})}{(\gamma_e - \gamma_c)} f - \biggl[1 + \frac{5}{2} (g^2-1) \biggr] ~>~ 0</math>

Solution Strategy

For a given set of free-energy coefficients, <math>~\mathcal{A}, \mathcal{B},</math> and <math>~\mathcal{C}</math>, along with a choice of the two adiabatic exponents <math>~(\gamma_c, \gamma_e)</math>, here's how to determine all of the physical parameters that are detailed in the above example table.

Step 1: Guess a value of <math>~0 < q < 1</math>.

Step 2: Given the pair of parameter values, <math>~(\mathcal{A}, q)</math>, determine the interface-density ratio, <math>~\rho_e/\rho_c</math>, by finding the appropriate root of the expression that defines the function, <math>~\mathcal{A}(q, \rho_e/\rho_c)</math>. This can be straightforwardly accomplished because, as demonstrated below, the relevant expression can be written as a quadratic function of <math>~(\rho_e/\rho_c)</math>.

Step 3: Given the pair of parameter values, <math>~(q, \rho_e/\rho_c)</math>, determine the value of the core-to-total mass ratio, <math>~\nu</math>, from the expression that was obtained from an integration over the mass, namely,

<math>~1 + \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl( \frac{1}{q^3} - 1 \biggr) \, .</math>

Step 4: Given the value of <math>~\mathcal{B}</math> along with the pair of parameter values, <math>~(q, \nu)</math>, the above expression that defines <math>~\mathcal{B}</math> can be solved to give the relevant value of the dimensionless parameter, <math>~\Lambda_\mathrm{eq}.</math>

Step 5: The value of <math>~\mathcal{C}^'</math> — the coefficient that appears on the right-hand-side of the above expression that defines <math>~\mathcal{C}</math> — can be determined, given the values of parameter triplet, <math>~(q, \nu, \Lambda_\mathrm{eq})</math>.

Step 6: Given the value of <math>~\mathcal{C}</math> and the just-determined value of the coefficient <math>~\mathcal{C}^'</math>, the normalized equilibrium radius, <math>~\chi_\mathrm{eq},</math> that corresponds to the value of <math>~q</math> that was guessed in Step #1 can be determined from the above definition of <math>~\mathcal{C}</math>, specifically,

<math>~\chi_\mathrm{eq}\biggr|_\mathrm{guess} </math>

<math>~\biggl( \frac{\mathcal{C}}{\mathcal{C}^'} \biggr)^{1/(3\gamma_e - 3\gamma_c)} \, .</math>

Step 7: But, independent of this guessed value of <math>~\chi_\mathrm{eq},</math> the condition for virial equilibrium — which identifies extrema in the free-energy function — gives the following expression for the normalized equilibrium radius:

<math>~\chi_\mathrm{eq}\biggr|_\mathrm{virial} </math>

<math>~\biggl[ \frac{\mathcal{A}}{\mathcal{B} + \mathcal{C}^'} \biggr]^{1/(4 - 3\gamma_c)} \, .</math>

Step 8: If <math>~\chi_\mathrm{eq}|_\mathrm{guess} \ne \chi_\mathrm{eq}|_\mathrm{virial}</math>, return to Step #1 and guess a different value of <math>~q</math>. Repeat Steps #1 through #7 until the two independently derived values of the normalized radius match, to a desired level of precision.

Keep in mind: (A) A graphical representation of the free-energy function, <math>~\mathfrak{G}(\chi)</math>, can also be used to identify the "correct" value of <math>~\chi_\mathrm{eq}</math> and, ultimately, the above-described iteration loop should converge on this value. (B) The free-energy function may exhibit more than one (or, actually, no) extrema, in which case more than one (or no) value of <math>~q</math> should lead to convergence of the above-described iteration loop.

Material that appears after this point in our presentation is under development and therefore
may contain incorrect mathematical equations and/or physical misinterpretations.
| Go Home |

Detailed Derivations

Dividing the free-energy expression through by <math>~E_\mathrm{norm}</math> generates,

<math>~\mathfrak{G}^* \equiv \frac{\mathfrak{G}}{E_\mathrm{norm}}</math>	<math>~=</math>	<math> - \frac{3}{5} \biggl( \frac{GM_\mathrm{tot}^2}{E_\mathrm{norm}} \biggr) \biggl( \frac{1}{R} \biggr) \cdot \mathfrak{f}_{WM} + \biggl[ \frac{\nu s_\mathrm{core} }{(\gamma_c - 1)} \biggr] \biggl[ \frac{M_\mathrm{tot} K_c \rho_{ic}^{\gamma_c-1} }{E_\mathrm{norm}} \biggr] </math>
		<math> ~+ \biggl[ \frac{(1-\nu) s_\mathrm{env} }{(\gamma_e - 1)} \biggr] \biggl[ \frac{M_\mathrm{tot} K_e \rho_{ie}^{\gamma_e-1} }{E_\mathrm{norm}} \biggr] </math>
	<math>~=</math>	<math> - \frac{3\cdot \mathfrak{f}_{WM}}{5} \biggl[ \frac{K_c G^{(3\gamma_c -4)}M_\mathrm{tot}^{2(3\gamma_c -4)}}{G^{3\gamma_c-3} M_\mathrm{tot}^{5\gamma_c-6}} \biggr]^{1/(3\gamma_c -4)} \biggl( \frac{1}{R} \biggr) </math>
		<math> + \biggl[ \frac{\nu s_\mathrm{core} }{(\gamma_c - 1)} \biggr] \biggl[ \frac{K_c M_\mathrm{tot}^{3\gamma_c -4} K_c^{3\gamma_c -4} }{G^{3\gamma_c-3} M_\mathrm{tot}^{5\gamma_c-6}} \biggr]^{1/(3\gamma_c -4)} \biggl( \frac{\rho_{ic}}{\bar\rho} \biggr)^{\gamma_c-1} \biggl[ \frac{3M_\mathrm{tot}}{4\pi R^3} \biggr]^{\gamma_c-1} </math>
		<math> ~+ \biggl[ \frac{(1-\nu) s_\mathrm{env} }{(\gamma_e - 1)} \biggr] \biggl[ \frac{K_c M_\mathrm{tot}^{3\gamma_c -4} K_e^{3\gamma_c -4} }{G^{3\gamma_c-3} M_\mathrm{tot}^{5\gamma_c-6}} \biggr]^{1/(3\gamma_c -4)} \biggl( \frac{\rho_{ie}}{\bar\rho} \biggr)^{\gamma_e-1} \biggl[ \frac{3M_\mathrm{tot}}{4\pi R^3} \biggr]^{\gamma_e-1} </math>
	<math>~=</math>	<math> - \frac{3\cdot \mathfrak{f}_{WM}}{5} \biggl( \frac{R_\mathrm{norm}}{R} \biggr) </math>
		<math> + \biggl[ \frac{\nu s_\mathrm{core} }{(\gamma_c - 1)} \biggr] \biggl( \frac{3M_\mathrm{tot}}{4\pi} \biggr)^{\gamma_c-1} \biggl[ \frac{ K_c^{3\gamma_c -3} }{G^{3\gamma_c-3} M_\mathrm{tot}^{2\gamma_c-2}} \biggr]^{1/(3\gamma_c -4)} \biggl[ \frac{R_\mathrm{norm}}{R} \biggr]^{3(\gamma_c-1)} \biggl[ \biggl( \frac{K_c}{G} \biggr) M_\mathrm{tot}^{\gamma_c-2} \biggr]^{-3(\gamma_c-1)/(3\gamma_c -4)} \biggl( \frac{\rho_{ic}}{\bar\rho} \biggr)^{\gamma_c-1} </math>
		<math> ~+ \biggl[ \frac{(1-\nu) s_\mathrm{env} }{(\gamma_e - 1)} \biggr] \biggl( \frac{3M_\mathrm{tot}}{4\pi} \biggr)^{\gamma_e-1} \biggl[ \frac{K_c^{3\gamma_c - 3} (K_e/K_c)^{3\gamma_c -4} }{G^{3\gamma_c-3} M_\mathrm{tot}^{2\gamma_c-2}} \biggr]^{1/(3\gamma_c -4)} \biggl[ \frac{R_\mathrm{norm}}{R} \biggr]^{3(\gamma_e-1)} \biggl[ \biggl( \frac{K_c}{G} \biggr) M_\mathrm{tot}^{\gamma_c-2} \biggr]^{-3(\gamma_e-1)/(3\gamma_c -4)} \biggl( \frac{\rho_{ie}}{\bar\rho} \biggr)^{\gamma_e-1} </math>
	<math>~=</math>	<math> - \frac{3\cdot \mathfrak{f}_{WM}}{5} \biggl( \frac{R_\mathrm{norm}}{R} \biggr) </math>
		<math> + \biggl[ \frac{\nu s_\mathrm{core} }{(\gamma_c - 1)} \biggr] \biggl( \frac{3}{4\pi} \biggr)^{\gamma_c-1} \biggl[ M_\mathrm{tot}^{3\gamma_c-4} \biggr]^{(\gamma_c-1)/(3\gamma_c-4)} \biggl[ M_\mathrm{tot}^{-2} \biggr]^{(\gamma_c-1)/(3\gamma_c -4)} \biggl[ M_\mathrm{tot}^{-3\gamma_c+6} \biggr]^{(\gamma_c-1)/(3\gamma_c -4)} \biggl[ \frac{R_\mathrm{norm}}{R} \biggr]^{3(\gamma_c-1)} \biggl( \frac{\rho_{ic}}{\bar\rho} \biggr)^{\gamma_c-1} </math>
		<math> ~+ \biggl[ \frac{(1-\nu) s_\mathrm{env} }{(\gamma_e - 1)} \biggr] \biggl( \frac{3}{4\pi} \biggr)^{\gamma_e-1} \biggl( \frac{K_e}{K_c} \biggr) \biggl[ M_\mathrm{tot}^{2(\gamma_e-1)-2(\gamma_c-1)}\biggr]^{1/(3\gamma_c-4)} \biggl[ \frac{R_\mathrm{norm}}{R} \biggr]^{3(\gamma_e-1)} \biggl[ \biggl( \frac{K_c}{G} \biggr)^{(\gamma_c-1)-(\gamma_e-1)} \biggr]^{3/(3\gamma_c -4)} \biggl( \frac{\rho_{ie}}{\bar\rho} \biggr)^{\gamma_e-1} </math>
	<math>~=</math>	<math> - \frac{3\cdot \mathfrak{f}_{WM}}{5} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-1} + \frac{\nu s_\mathrm{core} }{(\gamma_c - 1)} \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c-1} \biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{-3(\gamma_c-1)} </math>
		<math> ~+ \frac{(1-\nu) s_\mathrm{env} }{(\gamma_e - 1)} \biggl( \frac{K_e}{K_c} \biggr) \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2} \biggr]^{(\gamma_c-\gamma_e)/(3\gamma_c -4)} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ie}}{\bar\rho} \biggr]^{\gamma_e-1} \biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{-3(\gamma_e-1)} \, . </math>

We also want to ensure that envelope pressure matches the core pressure at the interface. This means,

<math>~K_e \rho_{ie}^{\gamma_e}</math>	<math>~=</math>	<math>~K_c \rho_{ic}^{\gamma_c}</math>
<math>\Rightarrow ~~~~\frac{K_e}{K_c} </math>	<math>~=</math>	<math>~\rho_{ic}^{\gamma_c} \rho_{ie}^{-\gamma_e} </math>
	<math>~=</math>	<math>~\biggl[ \frac{\rho_{ic}}{\rho_\mathrm{norm}} \biggr]^{\gamma_c} \biggl[ \frac{\rho_{ie}}{\rho_\mathrm{norm}} \biggr]^{-\gamma_e} \rho_\mathrm{norm}^{\gamma_c - \gamma_e}</math>
	<math>~=</math>	<math>~\biggl[ \frac{\rho_{ic}}{\rho_\mathrm{norm}} \biggr]^{\gamma_c} \biggl[ \frac{\rho_{ie}}{\rho_\mathrm{norm}} \biggr]^{-\gamma_e} \biggl\{ \frac{3}{4\pi} \biggl[ \frac{G^3 M_\mathrm{tot}^2}{K_c^3} \biggr]^{1/(3\gamma_c -4)} \biggr\}^{\gamma_c - \gamma_e}</math>
<math>\Rightarrow ~~~~\frac{K_e}{K_c} \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2} \biggr]^{(\gamma_c - \gamma_e)/(3\gamma_c -4)} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ie}}{\bar\rho} \biggr]^{\gamma_e-1}</math>	<math>~=</math>	<math>~\biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ic}}{\rho_\mathrm{norm}} \biggr]^{\gamma_c} \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ie}}{\rho_\mathrm{norm}} \biggr]^{-\gamma_e} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ie}}{\bar\rho} \biggr]^{\gamma_e-1} </math>
	<math>~=</math>	<math> \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c-1} \biggl( \frac{\rho_{ic}}{\rho_{ie}} \biggr) \biggl( \frac{\rho_\mathrm{norm}}{ \bar\rho } \biggr)^{\gamma_e - \gamma_c} </math>
	<math>~=</math>	<math> \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c-1} \biggl( \frac{\rho_{ic}}{\rho_{ie}} \biggr) \biggl( \frac{ R}{R_\mathrm{norm}} \biggr)^{3(\gamma_e - \gamma_c)} </math>

Hence, we can write the normalized (and dimensionless) free energy as,

<math>~\mathfrak{G}^*</math>

<math> - \frac{3\cdot \mathfrak{f}_{WM}}{5} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-1} + \biggl\{ \frac{\nu s_\mathrm{core} }{(\gamma_c - 1)} + \frac{(1-\nu) s_\mathrm{env} }{(\gamma_e - 1)} \biggl( \frac{ \rho_{ic} }{ \rho_{ie} } \biggr)\biggr\} \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c-1} \biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{-3(\gamma_c-1)} \, . </math>

Keep in mind that, if the envelope and core both have uniform (but different) densities, then <math>~\rho_{ic} = \rho_c</math>, <math>~\rho_{ie} = \rho_e</math>, and

Free Energy and Its Derivatives

Now, the free energy can be written as,

<math>~\mathfrak{G}</math>	<math>~=~</math>	<math>~U_\mathrm{tot} + W</math>
	<math>~=~</math>	<math>~\biggl[ \frac{2}{3(\gamma_c - 1)} \biggr] S_\mathrm{core} + \biggl[ \frac{2}{3(\gamma_e - 1)} \biggr] S_\mathrm{env} + W</math>
	<math>~=~</math>	<math>~\biggl[ \frac{2}{3(\gamma_c - 1)} \biggr] C_\mathrm{core} R^{3-3\gamma_c} + \biggl[ \frac{2}{3(\gamma_e - 1)} \biggr] C_\mathrm{env} R^{3-3\gamma_e} - A R^{-1} \, .</math>

The first derivative of the free energy with respect to radius is, then,

<math>~\frac{d\mathfrak{G}}{dR}</math>

<math>~ -2 C_\mathrm{core} R^{2-3\gamma_c} -2 C_\mathrm{env} R^{2-3\gamma_e} + A R^{-2} \, .</math>

And the second derivative is,

<math>~\frac{d^2\mathfrak{G}}{dR^2}</math>	<math>~=~</math>	<math>~ -2 (2-3\gamma_c) C_\mathrm{core} R^{1-3\gamma_c} -2 (2-3\gamma_e) C_\mathrm{env} R^{1-3\gamma_e} - 2A R^{-3} \, .</math>
	<math>~=~</math>	<math>~ \frac{2}{R^2} \biggl[(3\gamma_c-2) C_\mathrm{core} R^{3-3\gamma_c} + (3\gamma_e-2) C_\mathrm{env} R^{3-3\gamma_e} - A R^{-1} \biggr]</math>
	<math>~=~</math>	<math>~ \frac{2}{R^2} \biggl[(3\gamma_c-2) S_\mathrm{core} + (3\gamma_e-2) S_\mathrm{env} +W \biggr] \, .</math>

Equilibrium

The radius, <math>~R_\mathrm{eq}</math>, of the equilibrium configuration(s) is determined by setting the first derivative of the free energy to zero. Hence,

<math>~0 </math>	<math>~=~</math>	<math>~ 2 C_\mathrm{core} R_\mathrm{eq}^{2-3\gamma_c} + 2 C_\mathrm{env} R_\mathrm{eq}^{2-3\gamma_e} - A R_\mathrm{eq}^{-2} </math>
	<math>~=~</math>	<math>~ R_\mathrm{eq}^{-1} \biggl[ 2 C_\mathrm{core} R_\mathrm{eq}^{3-3\gamma_c} + 2 C_\mathrm{env} R_\mathrm{eq}^{3-3\gamma_e} - A R_\mathrm{eq}^{-1} \biggr]</math>
	<math>~=~</math>	<math>~ R_\mathrm{eq}^{-1} \biggl[ 2 S_\mathrm{core} + 2 S_\mathrm{env} +W \biggr]</math>
<math>\Rightarrow ~~~~ 2 S_\mathrm{tot} + W </math>	<math>~=~</math>	<math>~0 \, .</math>

This is the familiar statement of virial equilibrium. From it we should always be able to derive the radius of equilibrium configurations.

Stability

To assess the relative stability of an equilibrium configuration, we need to determine the sign of the second derivative of the free energy, evaluated at the equilibrium radius. If the sign of the second derivative is positive, the system is dynamically stable; if the sign is negative, he system is dynamically unstable. Using the above statement of virial equilibrium, that is, setting,

<math>~2 S_\mathrm{tot} + W</math>	<math>~=~</math>	<math>~0 \, ,</math>
<math>\Rightarrow ~~~~ S_\mathrm{env} </math>	<math>~=~</math>	<math>~- S_\mathrm{core} - \frac{W}{2} \, ,</math>

we obtain,

<math>~\frac{d^2\mathfrak{G}}{dR^2}\biggr\|_\mathrm{eq}</math>	<math>~=~</math>	<math>~ \frac{2}{R_\mathrm{eq}^2} \biggl[ (3\gamma_c-2) S_\mathrm{core} +W - (3\gamma_e-2)\biggl( S_\mathrm{core} + \frac{W}{2}\biggr) \biggr]_\mathrm{eq} </math>
	<math>~=~</math>	<math>~ \frac{2}{R_\mathrm{eq}^2} \biggl[ 3(\gamma_c-\gamma_e) S_\mathrm{core} + \biggl(2 - \frac{3}{2}\gamma_e\biggr)W \biggr]_\mathrm{eq} </math>
	<math>~=~</math>	<math>~ \frac{6}{R_\mathrm{eq}^2} \biggl[ (\gamma_c-\gamma_e) S_\mathrm{core} + \frac{1}{2}\biggl(\frac{4}{3} - \gamma_e\biggr)W \biggr]_\mathrm{eq} </math>
	<math>~=~</math>	<math>~ \frac{6}{R_\mathrm{eq}^2} \biggl[ -\frac{W}{2}\biggl( \gamma_e - \frac{4}{3}\biggr) - (\gamma_e-\gamma_c) S_\mathrm{core} \biggr]_\mathrm{eq} \, .</math>

So, if when evaluated at the equilibrium state, the expression inside of the square brackets of this last expression is negative, the equilibrium configuration will be dynamically unstable. We have chosen to write the expression in this particular final form because we generally will be interested in bipolytropes for which the adiabatic exponent of the envelope is greater than <math>~4/3</math> and the adiabatic exponent of the core is less than or equal to <math>~4/3</math> — that is, <math>~\gamma_e > 4/3 \ge \gamma_c</math>. Hence, because the gravitational potential energy, <math>~W</math>, is intrinsically negative, the system will be dynamically unstable only if the second term (involving <math>~S_\mathrm{core}</math>) is greater in magnitude than the first term (involving <math>~W</math>).

It is worth noting that, instead of drawing upon <math>~S_\mathrm{core}</math> and <math>~W</math> to define the stability condition, we could have used an appropriate combination of <math>~S_\mathrm{env}</math> and <math>~W</math>, or the <math>~S_\mathrm{core}</math> and <math>~S_\mathrm{env}</math> pair. Also, for example, because the virial equilibrium condition is <math>~S_\mathrm{tot} = -W/2</math>, it is easy to see that the following inequality also equivalently defines stability:

<math>~ S_\mathrm{tot}\biggl( \gamma_e - \frac{4}{3}\biggr) - (\gamma_e-\gamma_c) S_\mathrm{core} </math>

Examples

(0, 0) Bipolytropes

Review

In an accompanying discussion we have derived analytic expressions describing the equilibrium structure and the stability of bipolytropes in which both the core and the envelope have uniform densities, that is, bipolytropes with <math>~(n_c, n_e) = (0, 0)</math>. From this work, we find that integrals over the mass and pressure distributions give:

<math>~ \frac{W}{R_\mathrm{eq}^3 P_i} = - \frac{A}{R_\mathrm{eq}^4 P_i} </math>	<math>~=~</math>	<math>- ~ \frac{3}{5} \biggl[ \frac{GM_\mathrm{tot}^2}{R^4P_i} \biggr] \biggl( \frac{\nu^2}{q} \biggr) f </math>
	<math>~=~</math>	<math>- ~4\pi q^3 \Lambda f \, ,</math>
<math>~\frac{S_\mathrm{core}}{R_\mathrm{eq}^3 P_i} = B_\mathrm{core}</math>	<math>~=~</math>	<math> ~2\pi q^3 (1 + \Lambda) \, ,</math>
<math>~\frac{S_\mathrm{env}}{R_\mathrm{eq}^3 P_i} = B_\mathrm{env}</math>	<math>~=~</math>	<math> 2\pi \biggl[ (1-q^3) + \frac{5}{2} \Lambda \biggl( \frac{\rho_e}{\rho_0}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \Lambda \biggl( \frac{\rho_e}{\rho_0}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] \, ,</math>

where,

<math>~\Lambda</math>	<math>~\equiv~</math>	<math> \frac{3}{2^2 \cdot 5\pi} \biggl( \frac{GM_\mathrm{tot}^2}{R_\mathrm{eq}^4 P_i} \biggr) \frac{\nu^2}{q^4} \, ,</math>
<math>~f(q,\rho_e/\rho_c)</math>	<math>~\equiv~</math>	<math>1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^2} - 1 \biggr) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \biggl(\frac{1}{q^5} - 1 \biggr) - \frac{5}{2}\biggl(\frac{1}{q^2} - 1 \biggr) \biggr] </math>
	<math>~=~</math>	<math>1 + \frac{5}{2q^2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^2) + \frac{1}{2q^5} \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 (2 - 5q^3 + 3q^5) \, ,</math>
<math>~g^2(q,\rho_e/\rho_c)</math>	<math>~\equiv~</math>	<math>1 + \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( 1-q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] </math>
	<math>~\equiv~</math>	<math>1 + \biggl[ 2\biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q) + \frac{1}{q^2} \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 (1 - 3q^2 + 2q^3 ) \biggr] \, , </math>

Renormalize

Let's renormalize these energy terms in order to more readily relate them to the generalized expressions derived above.

<math>~R^3 P_i</math>	<math>~=</math>	<math> ~R^3 K_c \biggl(\frac{\rho_{ic}}{\bar\rho} \biggr)^{\gamma_c} \biggl[ \biggl( \frac{3}{4\pi}\biggr) \frac{M_\mathrm{tot}}{R^3} \biggr]^{\gamma_c} </math>
	<math>~=</math>	<math>~\biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{3-3\gamma_c} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c} K_c M_\mathrm{tot}^{\gamma_c} R_\mathrm{norm}^{3-3\gamma_c}</math>
	<math>~=</math>	<math>~\biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{3-3\gamma_c} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c} \biggl\{ K_c^{3\gamma_c -4} M_\mathrm{tot}^{\gamma_c(3\gamma_c -4)} \biggl[ \biggl( \frac{K_c}{G} \biggr) M_\mathrm{tot}^{\gamma_c-2} \biggr]^{3-3\gamma_c} \biggr\}^{1/(3\gamma_c -4)}</math>
	<math>~=</math>	<math>~\biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{3-3\gamma_c} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c} \biggl\{ \frac{G^{3\gamma_c-3} M_\mathrm{tot}^{5\gamma_c-6} }{K_c} \biggr\}^{1/(3\gamma_c -4)}</math>
	<math>~=</math>	<math>~\biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{3-3\gamma_c} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c} E_\mathrm{norm} \, .</math>

Also,

<math>~\biggl[ \frac{GM_\mathrm{tot}^2}{R} \biggr]^{3\gamma_c -4}</math>	<math>~=</math>	<math> ~\biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-(3\gamma_c-4)} G^{3\gamma_c -4} M_\mathrm{tot}^{6\gamma_c -8} \biggl( \frac{G}{K_c} \biggr) M_\mathrm{tot}^{2-\gamma_c} </math>
	<math>~=</math>	<math> ~\biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-(3\gamma_c-4)} \biggl[ \frac{ G^{3\gamma_c -3} M_\mathrm{tot}^{5\gamma_c -6} }{K_c} \biggr] </math>
	<math>~=</math>	<math> ~\biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-(3\gamma_c-4)} E_\mathrm{norm}^{3\gamma_c-4} \, . </math>
<math> \Rightarrow ~~~~\frac{GM_\mathrm{tot}^2}{R^4 P_i} </math>	<math>~=</math>	<math> ~\biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-1} \biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{3\gamma_c - 3} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{-\gamma_c} </math>
	<math>~=</math>	<math> \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{3\gamma_c - 4} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{-\gamma_c} \, . </math>

Hence,

<math>~\Lambda</math>

<math>~\equiv~</math>

<math>

\frac{3}{2^2 \cdot 5\pi} \frac{\nu^2}{q^4} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{-\gamma_c}

\biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{3\gamma_c - 4} \, .</math>

Given that <math>~\rho_{ic}/\bar\rho = \nu/q^3</math> for the <math>~(n_c, n_e) = (0, 0)</math> bipolytrope, we can finally write,

<math>~\frac{R^3 P_i}{E_\mathrm{norm}}</math>

<math>~\biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{3-3\gamma_c} \, ,</math>

and,

<math>~\Lambda</math>

<math>~\equiv~</math>

<math> \frac{3}{2^2 \cdot 5\pi} \frac{\nu^2}{q^4} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{-\gamma_c} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{3\gamma_c - 4} = \frac{1}{5} \frac{\nu}{q} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{1-\gamma_c} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{3\gamma_c - 4} \, .</math>

Hence the renormalized gravitational potential energy becomes,

<math> \frac{W_\mathrm{grav}}{E_\mathrm{norm}} </math>

<math> - \biggl( \frac{3}{5} \biggr) \frac{\nu^2}{q} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-1} \cdot f \, ;</math>

and the two, renormalized contributions to the thermal energy become,

<math>~\frac{U_\mathrm{core}}{ E_\mathrm{norm} } = \frac{2}{3(\gamma_c-1)} \biggl[ \frac{S_\mathrm{core}}{ E_\mathrm{norm} } \biggr]</math>	<math>~=~</math>	<math> \frac{4\pi q^3 (1 + \Lambda) }{3(\gamma_c-1)} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{3-3\gamma_c} \, ,</math>
<math>~\frac{U_\mathrm{env}}{ E_\mathrm{norm} } = \frac{2}{3(\gamma_e-1)} \biggl[ \frac{S_\mathrm{env}}{ E_\mathrm{norm} } \biggr]</math>	<math>~=~</math>	<math> \frac{4\pi}{3(\gamma_e-1)} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e} \biggl( \frac{K_e}{K_c} \biggr) \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2} \biggr]^{(3\gamma_c-3\gamma_e)/(3\gamma_c-4)} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{3-3\gamma_e} </math>
		<math> \times \biggl[ (1-q^3) + \frac{5}{2} \Lambda \biggl( \frac{\rho_e}{\rho_0}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \Lambda \biggl( \frac{\rho_e}{\rho_0}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] \, ,</math>

Finally, then, we can state that,

<math>~\mathfrak{f}_{WM}</math>	<math>~\equiv</math>	<math>~\frac{\nu^2}{q} \cdot f \, ,</math>
<math>~s_\mathrm{core}</math>	<math>~\equiv</math>	<math> 1 + \Lambda \, , </math>
<math>~(1-q^3) s_\mathrm{env}</math>	<math>~\equiv</math>	<math> (1-q^3) + \Lambda\biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_0}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_0}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] \, .</math>

Virial Equilibrium and Stability Evaluation

With these expressions in hand, we can deduce the equilibrium radius and relativity stability of <math>~(n_c, n_e) = (0, 0)</math> bipolytropes using the generalized expressions provided above. For example, from the statement of virial equilibrium <math>~(2S_\mathrm{tot} = - W )</math> we obtain,

<math>~q^3 (1 + \Lambda) + (1-q^3) + \frac{5}{2} \Lambda \biggl( \frac{\rho_e}{\rho_0}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \Lambda \biggl( \frac{\rho_e}{\rho_0}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) </math>

<math>~q^3 \Lambda \biggl[ 1 + \frac{5}{2q^2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^2) + \frac{1}{2q^5} \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 (2 - 5q^3 + 3q^5) \biggr] </math>

<math>\Rightarrow ~~~~ \frac{1}{\Lambda}</math>	<math>~=~</math>	<math>\frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (q-q^3) + \frac{1}{2q^2} \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 (2 - 5q^3 + 3q^5) - \biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_0}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_0}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] </math>
	<math>~=~</math>	<math>\frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (q-q^3 + 2 -3q +q^3) + \frac{1}{2q^2} \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 (2 - 5q^3 + 3q^5 +3 - 15q^2+15q^3 -3q^5) </math>
	<math>~=~</math>	<math>\frac{5}{2}\biggl[ 2\biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q) + \frac{1}{q^2} \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 (1 - 3q^2 + 2q^3 ) \biggr] </math>
	<math>~=~</math>	<math>\frac{5}{2}(g^2-1) </math>
<math>\Rightarrow ~~~~ \biggl[ \frac{P_i}{GM_\mathrm{tot}^2} \biggr] R_\mathrm{eq}^4</math>	<math>~=~</math>	<math>\biggl( \frac{3}{2^3 \pi } \biggr) \frac{\nu^2}{q^4} (g^2-1) \, . </math>

Or, given the above renormalization, this expression can be written as,

<math> \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{4-3\gamma_c } \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c} </math>	<math>~=~</math>	<math>\biggl( \frac{3}{2^3 \pi } \biggr) \frac{\nu^2}{q^4} (g^2-1) </math>
<math> \Rightarrow ~~~~ \frac{R}{R_\mathrm{norm}} </math>	<math>~=~</math>	<math> \biggl\{ \biggl( \frac{3}{2^3 \pi } \biggr) \frac{\nu^2}{q^4} (g^2-1) \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c} \biggr\}^{1/(4-3\gamma_c)} \, . </math>

And the condition for dynamical stability is,

<math>-\frac{W}{2}\biggl( \gamma_e - \frac{4}{3}\biggr) - (\gamma_e-\gamma_c) S_\mathrm{core} </math>	<math>~>~</math>	<math>~0 \, .</math>
<math>\Rightarrow ~~~~ 2\pi q^3 \Lambda \biggl[ \biggl( \gamma_e - \frac{4}{3}\biggr) f - (\gamma_e-\gamma_c) \biggl( 1 + \frac{1}{\Lambda}\biggr) \biggr] </math>	<math>~>~</math>	<math>~0 \, .</math>
<math>~\biggl( \gamma_e - \frac{4}{3} \biggr)f - (\gamma_e - \gamma_c) \biggl[1 + \frac{5}{2}(g^2-1) \biggr]</math>	<math>~>~</math>	<math>~0 \, .</math>

(5, 1) Bipolytropes

In another accompanying discussion we have derived analytic expressions describing the equilibrium structure of bipolytropes with <math>~(n_c, n_e) = (5, 1)</math>. Can we perform a similar stability analysis of these configurations? Work in progress!

Best of the Best

One Derivation of Free Energy

<math>~\mathfrak{G}^*</math>	<math>~=</math>	<math> - \frac{3\cdot \mathfrak{f}_{WM}}{5} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-1} + \frac{\nu s_\mathrm{core} }{(\gamma_c - 1)} \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\nu}{q^3} \biggr]^{\gamma_c-1} \biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{-3(\gamma_c-1)} </math>
		<math> ~+ \frac{(1-\nu) s_\mathrm{env} }{(\gamma_e - 1)} \biggl( \frac{K_e}{K_c} \biggr) \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2} \biggr]^{(\gamma_c-\gamma_e)/(3\gamma_c -4)} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e-1} \biggl[ \frac{R}{R_\mathrm{norm}} \biggr]^{-3(\gamma_e-1)} \, . </math>

Another Derivation of Free Energy

Hence the renormalized gravitational potential energy becomes,

<math> \frac{W_\mathrm{grav}}{E_\mathrm{norm}} </math>

<math> - \biggl( \frac{3}{5} \biggr) \frac{\nu^2}{q} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{-1} \cdot f \, ;</math>

and the two, renormalized contributions to the thermal energy become,

<math>~\frac{U_\mathrm{core}}{ E_\mathrm{norm} } = \frac{2}{3(\gamma_c-1)} \biggl[ \frac{S_\mathrm{core}}{ E_\mathrm{norm} } \biggr]</math>	<math>~=~</math>	<math> \frac{4\pi q^3 (1 + \Lambda) }{3(\gamma_c-1)} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c} \biggl( \frac{R}{R_\mathrm{norm}} \biggr)^{3-3\gamma_c} </math>
	<math>~=~</math>	<math> \frac{\nu (1 + \Lambda) }{(\gamma_c-1)} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{\gamma_c-1} \chi^{3-3\gamma_c} \, ,</math>
<math>~\frac{U_\mathrm{env}}{ E_\mathrm{norm} } = \frac{2}{3(\gamma_e-1)} \biggl[ \frac{S_\mathrm{env}}{ E_\mathrm{norm} } \biggr]</math>	<math>~=~</math>	<math> \frac{2 (2\pi) }{3(\gamma_e-1)} \biggl[ \frac{R^3 P_{ie}}{ E_\mathrm{norm} } \biggr] \biggl[ (1-q^3) + \frac{5}{2} \Lambda \biggl( \frac{\rho_e}{\rho_0}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \Lambda \biggl( \frac{\rho_e}{\rho_0}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] </math>
	<math>~=~</math>	<math> \frac{2 (2\pi) }{3(\gamma_e-1)} \biggl[ \frac{\mathrm{BigTerm}}{E_\mathrm{norm}} \biggr] R^3 K_e \rho_{ie}^{\gamma_e} </math>
	<math>~=~</math>	<math> \frac{2 (2\pi) }{3(\gamma_e-1)} \biggl[ \frac{\mathrm{BigTerm}}{E_\mathrm{norm}} \biggr] R^3 K_e \rho_\mathrm{norm}^{\gamma_e} \biggl( \frac{\rho_{ie}}{\bar\rho} \biggr)^{\gamma_e} \biggl( \frac{\bar\rho}{\rho_\mathrm{norm}} \biggr)^{\gamma_e} </math>
	<math>~=~</math>	<math> \frac{2 (2\pi) }{3(\gamma_e-1)} \biggl[ \frac{\mathrm{BigTerm}}{E_\mathrm{norm}} \biggr] ( \rho_\mathrm{norm} R_\mathrm{norm}^3) K_e \rho_\mathrm{norm}^{\gamma_e-1} \biggl[ \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e} \chi^{3-3\gamma_e} </math>
	<math>~=~</math>	<math> \frac{2 (2\pi) }{3(\gamma_e-1)} \biggl[ \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e} \chi^{3-3\gamma_e} \biggl[ \mathrm{BigTerm}\biggr] \biggl( \frac{3M_\mathrm{tot}}{4\pi} \biggr) \frac{K_e}{E_\mathrm{norm}} \biggl( \frac{3}{4\pi} \biggr)^{\gamma_e-1} \biggl[ \frac{G^3 M_\mathrm{tot}^2 }{K_c^3}\biggr]^{(\gamma_e-1)/(3\gamma_c-4)} </math>
	<math>~=~</math>	<math> \frac{(1-\nu)}{ (1-q^3) (\gamma_e-1)} \biggl[ \frac{3}{4\pi} \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e-1} \biggl( \frac{K_e}{K_c} \biggr) \chi^{3-3\gamma_e} \biggl[ \mathrm{BigTerm}\biggr] \frac{K_c M_\mathrm{tot}}{E_\mathrm{norm}} \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2}\biggr]^{(1-\gamma_e)/(3\gamma_c-4)} </math>
	<math>~=~</math>	<math> \frac{(1-\nu)}{ (1-q^3) (\gamma_e-1)} \biggl[ \frac{3}{4\pi} \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e-1} \biggl( \frac{K_e}{K_c} \biggr) \chi^{3-3\gamma_e} \biggl[ \mathrm{BigTerm}\biggr] \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2}\biggr]^{(1-\gamma_e)/(3\gamma_c-4)} \biggl[ \frac{K_c^{3\gamma_c-4} M_\mathrm{tot}^{3\gamma_c-4}}{ G^{3\gamma_c-3} M_\mathrm{tot}^{5\gamma_c-6} K_c^{-1}} \biggr]^{1/(3\gamma_c-4)} </math>
	<math>~=~</math>	<math> \frac{(1-\nu)}{ (1-q^3) (\gamma_e-1)} \biggl[ \frac{3}{4\pi} \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e-1} \biggl( \frac{K_e}{K_c} \biggr) \chi^{3-3\gamma_e} \biggl[ \mathrm{BigTerm}\biggr] \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2}\biggr]^{(1-\gamma_e)/(3\gamma_c-4)} \biggl[ \frac{K_c^{3}}{G^{3} M_\mathrm{tot}^{2}} \biggr]^{(\gamma_c-1)/(3\gamma_c-4)} </math>
	<math>~=~</math>	<math> \frac{(1-\nu)}{ (1-q^3) (\gamma_e-1)} \biggl[ \frac{3}{4\pi} \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e-1} \biggl( \frac{K_e}{K_c} \biggr) \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2}\biggr]^{(\gamma_c-\gamma_e)/(3\gamma_c-4)} \biggl[ \mathrm{BigTerm}\biggr] \chi^{3-3\gamma_e} </math>

Finally, then, we can state that,

<math>~\mathfrak{f}_{WM}</math>	<math>~\equiv</math>	<math>~\frac{\nu^2}{q} \cdot f \, ,</math>
<math>~s_\mathrm{core}</math>	<math>~\equiv</math>	<math> 1 + \Lambda \, , </math>
<math>~(1-q^3) s_\mathrm{env}</math>	<math>~\equiv</math>	<math> (1-q^3) + \Lambda\biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_0}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_0}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] \, .</math>

Note,

<math>~\Lambda</math>

<math>~\equiv~</math>

<math>

\frac{3}{2^2 \cdot 5\pi} \frac{\nu^2}{q^4} \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{-\gamma_c}

\biggl( \frac{R_\mathrm{eq}}{R_\mathrm{norm}} \biggr)^{3\gamma_c - 4} = \frac{1}{5} \biggl( \frac{\nu}{q} \biggr) \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{1-\gamma_c} \chi_\mathrm{eq}^{3\gamma_c - 4} \, .</math>

We also want to ensure that envelope pressure matches the core pressure at the interface. This means,

<math>~K_e \rho_{ie}^{\gamma_e}</math>	<math>~=</math>	<math>~K_c \rho_{ic}^{\gamma_c}</math>
<math>\Rightarrow ~~~~\frac{K_e}{K_c} </math>	<math>~=</math>	<math>~\rho_{ic}^{\gamma_c} \rho_{ie}^{-\gamma_e} </math>
	<math>~=</math>	<math>~\biggl[ \frac{\rho_{ic}}{\rho_\mathrm{norm}} \biggr]^{\gamma_c} \biggl[ \frac{\rho_{ie}}{\rho_\mathrm{norm}} \biggr]^{-\gamma_e} \rho_\mathrm{norm}^{\gamma_c - \gamma_e}</math>
	<math>~=</math>	<math>~\biggl[ \frac{\rho_{ic}}{\rho_\mathrm{norm}} \biggr]^{\gamma_c} \biggl[ \frac{\rho_{ie}}{\rho_\mathrm{norm}} \biggr]^{-\gamma_e} \biggl\{ \frac{3}{4\pi} \biggl[ \frac{G^3 M_\mathrm{tot}^2}{K_c^3} \biggr]^{1/(3\gamma_c -4)} \biggr\}^{\gamma_c - \gamma_e}</math>
<math>\Rightarrow ~~~~\frac{K_e}{K_c} \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2} \biggr]^{(\gamma_c - \gamma_e)/(3\gamma_c -4)} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ie}}{\bar\rho} \biggr]^{\gamma_e-1}</math>	<math>~=</math>	<math>~\biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ic}}{\rho_\mathrm{norm}} \biggr]^{\gamma_c} \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ie}}{\rho_\mathrm{norm}} \biggr]^{-\gamma_e} \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{\rho_{ie}}{\bar\rho} \biggr]^{\gamma_e-1} </math>
	<math>~=</math>	<math> \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c-1} \biggl( \frac{\rho_{ic}}{\rho_{ie}} \biggr) \biggl( \frac{\rho_\mathrm{norm}}{ \bar\rho } \biggr)^{\gamma_e - \gamma_c} </math>
	<math>~=</math>	<math> \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\rho_{ic}}{\bar\rho} \biggr]^{\gamma_c-1} \biggl( \frac{\rho_{ic}}{\rho_{ie}} \biggr) \biggl( \frac{ R}{R_\mathrm{norm}} \biggr)^{3(\gamma_e - \gamma_c)} </math>

Keep in mind that, if the envelope and core both have uniform (but different) densities, then <math>~\rho_{ic} = \rho_c</math>, <math>~\rho_{ie} = \rho_e</math>, and

Summary

Understanding Free-Energy Behavior

Step 1: Pick values for the separate coefficients, <math>\mathcal{A}, \mathcal{B},</math> and <math>\mathcal{C},</math> of the three terms in the normalized free-energy expression,

<math>~\mathfrak{G}^*</math>

<math>~ -~ 3\mathcal{A} \chi^{-1} - \frac{\mathcal{B}}{(1-\gamma_c)} ~\chi^{3-3\gamma_c} - \frac{\mathcal{C}}{(1-\gamma_e)} ~\chi^{3-3\gamma_e} </math>

then plot the function, <math>\mathfrak{G}^*(\chi)</math>, and identify the value(s) of <math>~\chi_\mathrm{eq}</math> at which the function has an extremum (or multiple extrema).

Step 2: Note that,

<math>~\mathcal{A}</math>	<math>~\equiv</math>	<math>~\frac{\nu^2}{5q} \biggl\{ 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^2} - 1 \biggr) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \biggl(\frac{1}{q^5} - 1 \biggr) - \frac{5}{2}\biggl(\frac{1}{q^2} - 1 \biggr) \biggr] \biggr\}</math>
	<math>~=</math>	<math>~\frac{1}{5} \biggl( \frac{\nu}{q^3} \biggr)^2 \biggl[ q^5 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (1 - q^2 )q^3 + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl( 1 - \frac{5}{2}q^3 + \frac{3}{2}q^5\biggr) \biggr]</math>
<math>~\mathcal{B}</math>	<math>~\equiv</math>	<math>~\nu \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\nu}{q^3} \biggr]^{\gamma_c-1} \biggl[ 1+\Lambda_\mathrm{eq} \biggr]</math>
<math>~\mathcal{C}</math>	<math>~\equiv</math>	<math>~(1-\nu)\biggl( \frac{K_e}{K_c} \biggr)^* \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{(1-\nu)}{(1-q^3)} \biggr]^{\gamma_e-1} \biggl\{ 1 + \frac{\Lambda_\mathrm{eq}}{(1-q^3)}\biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] \biggr\}</math>
	<math>~\equiv</math>	<math>~ \nu \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\nu}{q^3} \biggr]^{\gamma_c-1} \biggl\{ \frac{(1-q^3)}{q^3} + \frac{\Lambda_\mathrm{eq}}{q^3} \biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] \biggr\} \chi_\mathrm{eq}^{3(\gamma_e - \gamma_c)}</math>

where (see, for example, in the context of its original definition),

<math>~\Lambda_\mathrm{eq} \equiv

\frac{3}{2^2\pi \cdot 5} \biggl( \frac{GM_\mathrm{tot}^2}{R_\mathrm{eq}^4 P_i} \biggr) \frac{\nu^2}{q^4}

</math>

<math>

\frac{1}{5} \biggl( \frac{\nu}{q} \biggr) \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{1-\gamma_c}

\chi_\mathrm{eq}^{3\gamma_c - 4} </math>

and, where,

<math>\biggl( \frac{K_e}{K_c} \biggr)^* \equiv \frac{K_e}{K_c} \biggl[ \frac{K_c^3}{G^3 M_\mathrm{tot}^2} \biggr]^{(\gamma_c - \gamma_e)/(3\gamma_c -4)} </math>

<math> \biggl[\biggl( \frac{3}{4\pi} \biggr) \frac{1-\nu}{1-q^3} \biggr]^{-\gamma_e} \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\nu}{q^3} \biggr]^{\gamma_c} \chi_\mathrm{eq}^{3(\gamma_e - \gamma_c)} \, . </math>

Also, keep in mind that, if the envelope and core both have uniform (but different) densities, then <math>~\rho_{ic} = \rho_c</math>, <math>~\rho_{ie} = \rho_e</math>, and

<math> \frac{\rho_c}{\bar\rho} = \frac{\nu}{q^3} \, ; ~~~~~ \frac{\rho_e}{\bar\rho} = \frac{1-\nu}{1-q^3} \, ; ~~~~~ \frac{\rho_e}{\rho_c} = \frac{q^3(1-\nu)}{\nu (1-q^3)} ~~\Rightarrow ~~~ \frac{q^3}{\nu} = \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3 \, . </math>

Step 3: An analytic evaluation tells us that the following should happen. Using the numerically derived value for <math>~\chi_\mathrm{eq}</math>, define,

<math>~\mathcal{C}^' \equiv \mathcal{C} \chi_\mathrm{eq}^{3(\gamma_c - \gamma_e)} \, .</math>

We should then discover that,

<math>\frac{\mathcal{A}}{\mathcal{B} + \mathcal{C}^'} = \chi_\mathrm{eq}^{4-3\gamma_c} = \frac{1}{\Lambda_\mathrm{eq}} \cdot \frac{1}{5} \biggl( \frac{\nu}{q} \biggr) \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{1-\gamma_c} \, .</math>

Check It

<math>~\mathcal{B} + \mathcal{C}^'</math>

<math>~ \nu \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\nu}{q^3} \biggr]^{\gamma_c-1} \biggl\{ \biggl[ 1+\Lambda_\mathrm{eq} \biggr] + \frac{(1-q^3)}{q^3} + \frac{\Lambda_\mathrm{eq}}{q^3} \biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] \biggr\} </math>

<math>\Rightarrow~~~~\mathcal{A} \biggl[ \Lambda_\mathrm{eq} \cdot 5\biggl( \frac{q}{\nu^2} \biggr) \biggr]</math>

<math>~ 1+\Lambda_\mathrm{eq} + \frac{(1-q^3)}{q^3} + \frac{\Lambda_\mathrm{eq}}{q^3} \biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] </math>

<math> \Rightarrow~~~~\Lambda_\mathrm{eq} \biggl\{ 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^2} - 1 \biggr) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \biggl(\frac{1}{q^5} - 1 \biggr) - \frac{5}{2}\biggl(\frac{1}{q^2} - 1 \biggr) \biggr] \biggr\} </math>

<math>~ 1+\Lambda_\mathrm{eq} + \frac{(1-q^3)}{q^3} + \frac{\Lambda_\mathrm{eq}}{q^3} \biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] </math>

<math> \Rightarrow~~~~\Lambda_\mathrm{eq} \biggl\{ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^2} - 1 \biggr) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \biggl(\frac{1}{q^5} - 1 \biggr) - \frac{5}{2}\biggl(\frac{1}{q^2} - 1 \biggr) \biggr] \biggr\} </math>

<math>~ \frac{1}{q^3} + \frac{\Lambda_\mathrm{eq}}{q^3} \biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] </math>

<math>\Rightarrow~~~~ \frac{1}{\Lambda_\mathrm{eq}} </math>	<math>~=</math>	<math> q^3 \biggl\{ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^2} - 1 \biggr) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \biggl(\frac{1}{q^5} - 1 \biggr) - \frac{5}{2}\biggl(\frac{1}{q^2} - 1 \biggr) \biggr] \biggr\} - \biggl[ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) + \frac{3}{2q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) \biggr] </math>
<math>\Rightarrow~~~~ \frac{2}{\Lambda_\mathrm{eq}} </math>	<math>~=</math>	<math> 5\biggl( \frac{\rho_e}{\rho_c} \biggr) ( q - q^3 ) + \frac{2}{q^2}\biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (1 - q^5 ) - \frac{5}{2}( q^3 - q^5 ) \biggr] - 5 \biggl( \frac{\rho_e}{\rho_c}\biggr) (-2 + 3q - q^3) - \frac{3}{q^2} \biggl( \frac{\rho_e}{\rho_c}\biggr)^2 (-1 +5q^2 - 5q^3 + q^5) </math>
	<math>~=</math>	<math> 5\biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl[ ( q - q^3 ) - (-2 + 3q - q^3) \biggr] + \frac{1}{q^2}\biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ 2(1 - q^5 ) - 5( q^3 - q^5 ) - 3(-1 +5q^2 - 5q^3 + q^5) \biggr] </math>
	<math>~=</math>	<math> 10\biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl[ 1-q \biggr] + \frac{5}{q^2}\biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ 1 - 3q^2 + 2q^3 \biggr] </math>
<math>\Rightarrow~~~~ \frac{1}{\Lambda_\mathrm{eq}} \biggl[ \frac{2q^2}{5} \biggl( \frac{\rho_e}{\rho_c} \biggr)^{-1} \biggr] </math>	<math>~=</math>	<math> 2q^2 (1-q )+ \biggl( \frac{\rho_e}{\rho_c} \biggr) ( 1 - 3q^2 + 2q^3 ) </math>

Fortunately, this precisely matches our earlier derivation, which states that,

<math>~\frac{1}{\Lambda}</math>

<math>\frac{5}{2}(g^2-1) = \frac{5}{2}\biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( 1-q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] \, .</math>

Playing With One Example

By setting,

<math>~\mathcal{A}</math>	<math>~\mathcal{B}</math>	<math>~\mathcal{C}</math>
<math>~\gamma_c = 6/5; ~~~ \gamma_e = 2</math>
2.5	1.0	2.0

a plot of <math>~\mathfrak{G}^*</math> versus <math>~\chi</math> exhibits the following, two extrema:

extremum	<math>~\chi_\mathrm{eq}</math>	<math>~\mathfrak{G}^*</math>		<math>~\chi_\mathrm{eq}^{3(\gamma_c - \gamma_e)}</math>	<math>~\mathcal{C}^'</math>	<math>~\chi_\mathrm{eq}^{4-3\gamma_c}</math>	<math>~\frac{\mathcal{A}}{\mathcal{B} +\mathcal{C}^'}</math>
MIN	<math>1.1824</math>	<math>-0.611367</math>	<math>~~~~\Rightarrow~~~~</math>	<math>0.66891</math>	<math>1.3378</math>	<math>1.0693</math>	<math>1.0694</math>
MAX	<math>9.6722</math>	<math>+0.508104</math>	<math>~\Rightarrow</math>	<math>0.004313</math>	<math>0.008625</math>	<math>2.4786</math>	<math>2.4786</math>

The last two columns of this table confirm the internal consistency of the relationships presented in Step 3, above. But what does this mean in terms of the values of <math>~\nu</math>, <math>~q</math>, and the related ratio of densities at the interface, <math>~\rho_e/\rho_c</math>?

Let's assume that what we're trying to display and examine is the behavior of the free-energy surface for a fixed value of the ratio of densities at the interface. Once the value of <math>~\rho_e/\rho_c</math> has been specified, it is clear that the value of <math>~q</math> (and, hence, also <math>~\nu</math>) is set because <math>~\mathcal{A}</math> has also been specified. But our specification of <math>~\mathcal{B}</math> along with <math>~\rho_e/\rho_c</math> also forces a particular value of <math>~q</math>. It is unlikely that these two values of <math>~q</math> will be the same. In reality, once <math>~\mathcal{A}</math> and <math>~\mathcal{B}</math> have both been specified, they force a particular <math>~(\nu, q)</math> pair. How do we (easily) figure out what this pair is?

Let's begin by rewriting the expressions for <math>~\mathcal{A}</math> and <math>~\mathcal{B}</math> in terms of just <math>~q</math> and the ratio, <math>~\rho_e/\rho_c</math>, keeping in mind that, for the case of a uniform-density core (of density, <math>~\rho_c</math>) and a uniform-density envelope (of density, <math>~\rho_e</math>),

hence,

<math>~\biggl[ 1 + \biggl( \frac{\rho_e}{\rho_c} \biggr) \frac{(1-q^3)}{q^3} \biggr]^{-1}</math>

and

<math>~\biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3 \biggr] \, .</math>

Putting the expression for <math>~\mathcal{A}</math> in the desired form is simple because <math>~\nu</math> only appears as a leading factor. Specifically, we have,

<math>~\mathcal{A}</math>	<math>~=</math>	<math>~\frac{\pi q^5}{5} \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3 \biggr]^{-2} \biggl\{ 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^2} - 1 \biggr) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \biggl(\frac{1}{q^5} - 1 \biggr) - \frac{5}{2}\biggl(\frac{1}{q^2} - 1 \biggr) \biggr] \biggr\}</math>
	<math>~=</math>	<math>~\frac{\pi }{5} \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3 \biggr]^{-2} \biggl\{ q^5 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (q^3 - q^5 ) + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[1-\frac{5}{2} q^3 + \frac{3}{2}q^5 \biggr] \biggr\} \, .</math>

The expression for <math>~\mathcal{B}</math> can be written in the form,

<math>~\mathcal{B}</math>	<math>~=</math>	<math>~\nu \biggl[ \biggl( \frac{3}{4\pi} \biggr)\frac{\nu}{q^3} \biggr]^{\gamma_c-1} \biggl\{ 1+\frac{\pi}{5} \biggl( \frac{\nu}{q} \biggr) \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{1-\gamma_c} \chi_\mathrm{eq}^{3\gamma_c - 4} \biggr\}</math>
	<math>~=</math>	<math>~\nu \biggl[ \biggl( \frac{4\pi}{3} \biggr)\frac{q^3}{\nu} \biggr]^{1-\gamma_c} +\frac{\pi q^5}{5} \biggl( \frac{\nu^2}{q^6} \biggr) \chi_\mathrm{eq}^{3\gamma_c - 4} </math>
	<math>~=</math>	<math>~q^3 \biggl( \frac{4\pi}{3} \biggr)^{1-\gamma_c} \biggl[ \frac{q^3}{\nu} \biggr]^{-\gamma_c} +\frac{\pi q^5}{5} \biggl( \frac{q^3}{\nu} \biggr)^{-2} \chi_\mathrm{eq}^{3\gamma_c - 4} </math>
	<math>~=</math>	<math>~q^3 \biggl( \frac{4\pi}{3} \biggr)^{1-\gamma_c} \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3 \biggr]^{-\gamma_c} +\frac{\pi q^5}{5} \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3 \biggr]^{-2} \chi_\mathrm{eq}^{3\gamma_c - 4} \, . </math>

Generally speaking, the equilibrium radius, <math>~\chi_\mathrm{eq}</math>, which appears in the expression for <math>~\mathcal{B}</math>, is not known ahead of time. Indeed, as is illustrated in our simple example immediately above, the normal path is to pick values for the coefficients, <math>~\mathcal{A}</math>, <math>~\mathcal{B}</math>, and <math>~\mathcal{C}</math>, and determine the equilibrium radius by looking for extrema in the free-energy function. And because <math>~\chi_\mathrm{eq}</math> is not known ahead of time, it isn't clear how to (easily) figure out what pair of physical parameter values, <math>~(\nu, q)</math>, give self-consistent values for the coefficient pair, <math>~(\mathcal{A}, \mathcal{B})</math>.

Because we are using a uniform density core and uniform density envelope as our base model, however, we do know the analytic solution for <math>~\chi_\mathrm{eq}</math>. As stated above, it is,

<math>~\chi_\mathrm{eq}^{4-3\gamma_c}</math>	<math>~=</math>	<math>~\frac{1}{\Lambda_\mathrm{eq}} \cdot \frac{\pi}{5} \biggl( \frac{\nu}{q} \biggr) \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{\nu}{q^3} \biggr]^{1-\gamma_c} </math>
	<math>~=</math>	<math>~\frac{\pi q^2}{2} \biggl( \frac{3}{4\pi} \biggr)^{1-\gamma_c}\biggl( \frac{\nu}{q^3} \biggr) \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_c} \biggr) \biggl( 1-q \biggr) + \frac{\rho_e}{\rho_c} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] \cdot \biggl[ \frac{q^3}{\nu} \biggr]^{\gamma_c-1} </math>
	<math>~=</math>	<math>~\frac{\pi}{2} \biggl( \frac{3}{4\pi} \biggr)^{1-\gamma_c} \biggl(\frac{\rho_e}{\rho_c}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_c} \biggr) ( q^2-q^3 ) + \frac{\rho_e}{\rho_c} (1 - q^2) \biggr] \cdot \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3\biggr]^{\gamma_c-2} </math>

Combining this expression with the one for <math>~\mathcal{B}</math> gives us the desired result — although, strictly speaking, it is cheating! We can now methodically choose <math>~(\nu, q)</math> pairs and map them into the corresponding values of <math>~\mathcal{A}</math> and <math>~\mathcal{B}</math>. And, via an analogous "cheat," the choice of <math>~(\nu, q)</math> also gives us the self-consistent value of <math>~\mathcal{C}</math>. In this manner, we should be able to map out the free-energy surface for any desired set of physical parameters.

Second Example

Explain Logic

The figure presented here, on the right, shows a plot of the free energy, as a function of the dimensionless radius, <math>~\mathfrak{G}^*(\chi)</math>, where,

<math>~\mathfrak{G}^*</math>

<math>~ -~ 3\mathcal{A} \chi^{-1} - \frac{\mathcal{B}}{(1-\gamma_c)} ~\chi^{3-3\gamma_c} - \frac{\mathcal{C}}{(1-\gamma_e)} ~\chi^{3-3\gamma_e} \, ,</math>

and, where we have used the parameter values,

<math>~\mathcal{A}</math>	<math>~\mathcal{B}</math>	<math>~\mathcal{C}</math>
<math>~\gamma_c = 6/5; ~~~ \gamma_e = 2</math>
0.201707	0.0896	0.002484

Directly from this plot we deduce that this free-energy function exhibits a minimum at <math>~\chi_\mathrm{eq} = 0.1235</math> and that, at this equilibrium radius, the configuration has a free-energy value, <math>~\mathfrak{G}^*(\chi_\mathrm{eq} ) = -2.0097</math>. Via the steps described below, we demonstrate that this identified equilibrium radius is appropriate for an <math>~(n_c, n_e) = (0, 0)</math> bipolytrope (with the just-specified core and envelope adiabatic indexes) that has the following physical properties:

Fractional core mass, <math>~\nu = 0.1</math>;
Core-envelope interface located at <math>~r_i/R = q = 0.435</math>;
Density jump at the core-envelope interface, <math>~\rho_e/\rho_c = 0.8</math>.

Step 1: Because the ratio, <math>~q^3/\nu</math>, is a linear function of the density ratio, <math>~\rho_e/\rho_c</math>, the full definition of the free-energy coefficient, <math>~\mathcal{A}</math>, can be restructured into a quadratic equation that gives the density ratio for any choice of the parameter pair, <math>~(q, \mathcal{A})</math>. Specifically,

<math>~5 \biggl( \frac{q^3}{\nu} \biggr)^2 \mathcal{A}</math>	<math>~=</math>	<math> q^5 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (1 - q^2 )q^3 + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl( 1 - \frac{5}{2}q^3 + \frac{3}{2}q^5\biggr) </math>
<math>~\Rightarrow ~~~~ 5\mathcal{A} \biggl[ \biggl( \frac{\rho_e}{\rho_c} \biggr) (1-q^3) + q^3 \biggr]^2 </math>	<math>~=</math>	<math> q^5 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (1 - q^2 )q^3 + \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl( 1 - \frac{5}{2}q^3 + \frac{3}{2}q^5\biggr) \, , </math>

and this can be written in the form,

<math> \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 a + \biggl( \frac{\rho_e}{\rho_c} \biggr) b + c</math>

where,

<math>~a</math>	<math>~\equiv</math>	<math>~5\mathcal{A} (1-q^3)^2 - 1 + \frac{5}{2}q^3 - \frac{3}{2}q^5 \, ,</math>
<math>~b</math>	<math>~\equiv</math>	<math>~10\mathcal{A} q^3(1-q^3) - \frac{5}{2}q^3 (1-q^2) \, ,</math>
<math>~c</math>	<math>~\equiv</math>	<math>~5\mathcal{A} q^6 - q^5 \, .</math>

Hence,

<math>~\frac{1}{2a} \biggl[\pm ( b^2 - 4ac)^{1/2} - b \biggr] \, .</math>

(For our physical problem it appears as though only the positive root is relevant.) For the purposes of this example, we set <math>~\mathcal{A} = 0.2017</math> and examined a range of values of <math>~q</math> to find a physically interesting value for the density ratio. We picked:

<math>~\mathcal{A}</math>	<math>~q</math>		<math>~a</math>	<math>~b</math>	<math>~c</math>		<math>~\frac{\rho_e}{\rho_c}</math>		<math>~\nu</math>
<math>~0.2017</math>	<math>~0.435</math>	<math>\Rightarrow</math>	<math>~0.03173</math>	<math>~-0.01448</math>	<math>~-0.008743</math>	<math>\Rightarrow</math>	<math>~0.80068</math>	<math>\Rightarrow</math>	<math>~0.10074</math>

Step 2: Next, we chose the parameter pair,

<math> ~\biggl(q, \frac{\rho_e}{\rho_c} \biggr) = (0.43500, 0.80000) </math>

and determined the following parameter values from the known analytic solution:

<math>~\nu</math>	<math>~f\biggl(q, \frac{\rho_e}{\rho_c} \biggr)</math>	<math>~g^2\biggl(q, \frac{\rho_e}{\rho_c} \biggr)</math>	<math>~\Lambda_\mathrm{eq}</math>	<math>~\chi_\mathrm{eq}</math>	<math>~\mathcal{A}</math>	<math>~\mathcal{B}</math>	<math>~\mathcal{C}</math>
<math>~0.100816</math>	<math>~43.16365</math>	<math>~3.923017</math>	<math>~0.13684</math>	<math>~0.12349</math>	<math>~0.201707</math>	<math>~0.089625</math>	<math>~0.002484</math>

Construction Multiple Curves to Define a Free-Energy Surface

Okay. Now that we have the hang of this, let's construct a sequence of curves that represent physical evolution at a fixed interface-density ratio, <math>~\rho_e/\rho_c</math>, but for steadily increasing core-to-total mass ratio, <math>~\nu</math>. Specifically, we choose,

From the known analytic solution, here are parameters defining several different equilibrium models:

Identification of Local Minimum in Free Energy
<math>~\nu</math>	<math>~q</math>	<math>~f\biggl(q, \frac{\rho_e}{\rho_c} \biggr)</math>	<math>~g^2\biggl(q, \frac{\rho_e}{\rho_c} \biggr)</math>	<math>~\Lambda_\mathrm{eq}</math>	<math>~\chi_\mathrm{eq}</math>	<math>~\mathcal{A}</math>	<math>~\mathcal{B}</math>	<math>~\mathcal{C}</math>
<math>~0.2</math>	<math>~9^{-1/3} = 0.48075</math>	<math>~12.5644</math>	<math>~2.091312</math>	<math>~0.366531</math>	<math>~0.037453</math>	<math>~0.2090801</math>	<math>~0.2308269</math>	<math>~2.06252 \times 10^{-4}</math>
<math>~0.4</math>	<math>~4^{-1/3} = 0.62996</math>	<math>~4.21974</math>	<math>~1.56498</math>	<math>~0.707989</math>	<math>~0.0220475</math>	<math>~0.2143496</math>	<math>~0.5635746</math>	<math>~4.4626 \times 10^{-5}</math>
<math>~0.5</math>	<math>~3^{-1/3} = 0.693361</math>	<math>~2.985115</math>	<math>~1.42334</math>	<math>~0.9448663</math>	<math>~0.0152116</math>	<math>~0.2152641</math>	<math>~0.791882</math>	<math>~1.5464 \times 10^{-5}</math>
Here we are examining the behavior of the free-energy function for bipolytropic models having <math>~(n_c, n_e) = (0, 0)</math>, <math>~(\gamma_c, \gamma_e) = (6/5, 2)</math>, and a density ratio at the core-envelope interface of <math>~\rho_e/\rho_c = 1/2</math>. The figure shown here, on the right, displays the three separate free-energy curves, <math>~\mathfrak{G}^*(\chi)</math> — where, <math>~\chi \equiv R/R_\mathrm{norm}</math> is the normalized configuration radius — that correspond to the three values of <math>~\nu \equiv M_\mathrm{core}/M_\mathrm{tot}</math> given in the first column of the above table. Along each curve, the local free-energy minimum corresponds to the the equilibrium radius, <math>~\chi_\mathrm{eq}</math>, recorded in column 6 of the above table.

Each of the free-energy curves shown above has been entirely defined by our specification of the three coefficients in the free-energy function, <math>~\mathcal{A}, \mathcal{B}</math>, and <math>~\mathcal{C}</math>. In each case, the values of these three coefficients was judiciously chosen to produce a curve with a local minimum at the correct value of <math>~\chi_\mathrm{eq}</math> corresponding to an equilibrium configuration having the desired <math>~(\nu, \rho_e/\rho_c)</math> model parameters. Upon plotting these three curves, we noticed that two of the curves — curves for <math>~\nu = 0.4</math> and <math>~\nu = 0.5</math> — also display a local maximum. Presumably, these maxima also identify equilibrium configurations, albeit unstable ones. From a careful inspection of the plotted curves, we have identified the value of <math>~\chi_\mathrm{eq}</math> that corresponds to the two newly discovered (unstable) equilibrium models; these values are recorded in the table that immediately follows this paragraph. By construction, we also know what values of <math>~\mathcal{A}, \mathcal{B}</math>, and <math>~\mathcal{C}</math> are associated with these two identified equilibria; these values also have been recorded in the table. But it is not immediately obvious what the values are of the <math>~(\nu, \rho_e/\rho_c)</math> model parameters that correspond to these two equilibrium models.

Subsequently Identified Local Energy Maxima
<math>~\chi_\mathrm{eq}</math>	<math>~\mathfrak{G}^*</math>	<math>~\chi_\mathrm{eq}^{4-3\gamma_c}</math>	<math>~\therefore</math>	<math>~\mathcal{A}</math>	<math>~\mathcal{B}</math>	<math>~\mathcal{C}</math>	<math>~\mathcal{C}^' = \mathcal{A} \chi_\mathrm{eq}^{3\gamma_c-4} - \mathcal{B}</math>	<math>~\biggl( \frac{\mathcal{C}}{\mathcal{C}^'} \biggr)^{1/(3\gamma_e - 3\gamma_c)}</math>
<math>~0.08255</math>	<math>~+ 4.87562</math>	<math>~0.368715</math>	<math>~0</math>	<math>~0.2143496</math>	<math>~0.5635746</math>	<math>~4.4626 \times 10^{-5}</math>	<math>~1.7768 \times 10^{-2}</math>	<math>~0.08254</math>
<math>~0.032196</math>	<math>~+11.5187</math>	<math>~0.25300</math>	<math>~0</math>	<math>~0.2152641</math>	<math>~0.791882</math>	<math>~1.5464 \times 10^{-5}</math>	<math>~5.8964 \times 10^{-2}</math>	<math>~0.032196</math>

Related Discussions

Analytic solution with <math>~n_c = 0</math> and <math>~n_e=0</math>.
Analytic solution with <math>~n_c = 5</math> and <math>~n_e=1</math>.

User:Tohline/SSC/BipolytropeGeneralization

Bipolytrope Generalization

Setup

Review of Isolated Polytrope

Bipolytrope

Partitioning the Mass

Separate Contributions to Gravitational Potential Energy

Separate Thermodynamic Energy Reservoirs

Sanity Check

Old Stuff

Simplest Bipolytrope

Familiar Setup

Cleaner Virial Presentation

Shift to Central Pressure Normalization

Free-Energy Coefficients

More General Derivation of Free-Energy Coefficients B and C

Extrema

Examples

Solution Strategy

Detailed Derivations

Free Energy and Its Derivatives

Equilibrium

Stability

Examples

(0, 0) Bipolytropes

Review

Renormalize

Virial Equilibrium and Stability Evaluation

(5, 1) Bipolytropes

Best of the Best

One Derivation of Free Energy

Another Derivation of Free Energy

Summary

Understanding Free-Energy Behavior

Check It

Playing With One Example

Second Example

Explain Logic

Construction Multiple Curves to Define a Free-Energy Surface

Related Discussions

Navigation menu

Search