Editing SSC/VirialStability

__FORCETOC__ <!-- will force the creation of a Table of Contents -->
<!-- __NOTOC__ will force TOC off -->

=Virial Stability of BiPolytropes=

==Discussion with Kundan Kadam==
See particularly &hellip; [https://ui.adsabs.harvard.edu/abs/2016MNRAS.462.2237K/abstract K. Kadam ''et al.'' (2016, MNRAS, 462, pp. 2237-2245)]

The comments provided inside the following box are for Kundan Kadam, providing feedback on a derivation he sent to me ( J. E. Tohline) on Monday, 27 January 2014.

<table align="center" border="1" width="95%" cellpadding="5">
<tr><td align="left">
[<font color="red">31 January 2014</font>] I agree with the following derived expressions; they all already appear in my presentation elsewhere on this page:
<div align="center">
Throughout the envelope &hellip; &nbsp; &nbsp; &nbsp;
<math>
M_r = \frac{4\pi}{3} \biggl( r^3\rho_e + r_i^3\rho_c - r_i^3\rho_e\biggr) \, ;
</math>

<math>
M_\mathrm{tot} = \frac{4\pi}{3}\rho_c \biggl( \frac{q^3 R^3}{\nu} \biggr) \, ;
</math>

<math>
\nu = \biggl[1 + \biggl(\frac{\rho_e}{\rho_c}\biggr)(q^{-3} - 1) \biggr]^{-1} \, .
</math>
</div>
During the recent December holiday break, I also used the hydrostatic balance relation to derive expressions for the pressure throughout the two-component (uniform density core + uniform density envelope) model, as you have done.  I agree completely with your derivation of the expressions for &#8212; and relationships between &#8212; <math>P_{cc}</math>, <math>P_{ci}</math>, and <math>P_{ei}</math>.  In addition, I chose to normalize all the pressures to,
<div align="center">
<math>
P_\mathrm{norm} \equiv \biggl(G^3 \rho_c^4 M_\mathrm{tot}^2 \biggr)^{1/3} \, .
</math>
</div>
Letting an asterisk superscript denote normalized pressures, that is, <math>P^* \equiv P/P_\mathrm{norm}</math>, my derived expressions are:
<div align="center">
<math>
P^*_{cc} = P^*_{ci} + \biggr[\frac{\pi \nu^2}{6} \biggr]^{1/3} \, ,
</math>
</div>
and,
<div align="center">
<math>
P^*_{ei} = P^*_{ci} = \biggr[\frac{\pi \nu^2}{6} \biggr]^{1/3} \biggl( \frac{\rho_e}{\rho_c}\biggr) \biggl[ \frac{\rho_e}{\rho_c} 
\biggl(\frac{1}{q^2}-1\biggr) + 2 \biggl(1- \frac{\rho_e}{\rho_c}\biggr)(1-q) \biggr] \, .
</math>
</div>
<br />

<div align="center">
<hr width="85%" />
</div>

<br />
[<font color="red">4 February 2014</font>] In deriving the above analytic expression for the pressure, we have effectively determined the detailed [[SSC/Structure/BiPolytropes/Analytic00#BiPolytrope_with_nc_.3D_0_and_ne_.3D_0|equilibrium structure of a bipolytrope that has <math>n_c = 0</math> and <math>n_e = 0</math>]] without ever using the virial theorem.  This is an excellent accomplishment for a variety of reasons.  

<ins>First</ins>, we can now compare and contrast these newly derived expressions with the analogous expressions that have been derived earlier for a [[SSC/Structure/BiPolytropes/Analytic51#BiPolytrope_with_nc_.3D_5_and_ne_.3D_1|bipolytrope with  <math>n_c = 5</math> and <math>n_e = 1</math>]].  For example, at the end of [[SSC/Structure/BiPolytropes/Analytic00#Step_7:__Surface_Boundary_Condition|Step #7 in the new derivation]], we find the following expressions for the equilibrium radius and total mass of a bipolytrope with <math>(n_c, n_e) = (0, 0)</math> expressed in terms of the chosen central density, <math>\rho_0</math>, and central pressure, <math>P_0</math>, of the configuration:
<div align="center">
<table border="0" cellpadding="5">

<tr>
  <td align="right">
<math>
\biggl[ \frac{G\rho_0^2}{P_0} \biggr]^{1/2} R</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp; 
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2\pi} \biggr)^{1/2} \frac{\sqrt{3}}{q g} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\biggl[ \frac{G^3\rho_0^4}{P_0^3} \biggr]^{1/2} M_\mathrm{tot}</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp; 
  </td>
  <td align="left">
<math>\biggl( \frac{2}{\pi} \biggr)^{1/2} \frac{\sqrt{3}}{\nu g^3} \, ,</math>
  </td>
</tr>
</table>
</div>
where the function (<font color="red">note that, for a few days, there was a sign error in this definition of <math>g</math></font>),
<div align="center">
<math>
g(q, \rho_e/\rho_0) \equiv \biggl\{ 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] \biggr\}^{1/2}\, .
</math>
</div>
We will also find it useful to combine these two expressions to eliminate direct reference to the central density, <math>\rho_0</math>, obtaining,
<div align="center">
<table border="0">
<tr>
  <td align="right">
<math>
\biggl[ \frac{R^4}{GM_\mathrm{tot}^2} \biggr] P_0</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp; 
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2^3\pi} \biggr) \frac{3\nu^2 g^2}{q^4} \, .</math>
  </td>
</tr>
</table>
</div>

For comparison, referring back to the "Normalization" discussion and table of "Parameter Values" provided in our [[SSC/Structure/BiPolytropes/Analytic51#Examples|example bipolytrope with <math>(n_c, n_e) = (5, 1)</math>]], we have,
<div align="center">
<table border="0" cellpadding="5">

<tr>
  <td align="right">
<math>
\biggl( \frac{\mu_e}{\mu_c} \biggr) \biggl[ \frac{G\rho_0^2}{(K_c \rho_0^{6/5})} \biggr]^{1/2} R</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp; 
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2\pi} \biggr)^{1/2} \frac{\eta_s}{\theta_i^2} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\biggl( \frac{\mu_e}{\mu_c} \biggr) ^2 \biggl[ \frac{G^3\rho_0^4}{(K_c \rho_0^{6/5})^3} \biggr]^{1/2} M_\mathrm{tot}</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp; 
  </td>
  <td align="left">
<math>\biggl( \frac{2}{\pi} \biggr)^{1/2} \frac{A\eta_s}{\theta_i} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\biggl[ \frac{R^4}{GM^2_\mathrm{tot}} \biggr]^{1/2} K_c \rho_0^{6/5}</math>
  </td>
  <td align="center">
&nbsp; <math>~=</math>&nbsp; 
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2^3\pi} \biggr) \frac{\eta_s^2}{A^2\theta_i^6} \, ,</math>
  </td>
</tr>
</table>
</div>
where expressions for the various functions, <math>A, \eta_s, \theta_i</math>, are also provided in the table of "Parameter Values."  The dimensional normalizations are clearly the same in both types of bipolytropes because, in the latter case, <math>P_0 = K_c \rho_0^{6/5}</math>.

<ins>Second</ins>, a comparison between these two types of bipolytropes may help clarify how a discontinuous jump in the mean molecular weight should be handled in the more general virial analysis.  Usually <math>\rho_e/\rho_c</math> is set equal to <math>\mu_e/\mu_c</math> at the interface, but if the core and envelope are both treated as having uniform densities, setting the ratio of densities to the ratio of mean molecular weights seems to overconstrain the problem.  I haven't figured out yet how to handle this, but maybe this model comparison will help.

<ins>Third</ins>, when conducting the virial analysis, it should now be clear how to evaluate the system's free energy.  When the system is treated as being composed of a central uniform-density spherical core surrounded by a uniform-density envelope, we have already derived an analytic expression for the total gravitational potential energy, <math>W</math>.  The above derivation of the pressure distribution throughout such a configuration &#8212; that is, throughout a bipolytrope with <math>(n_c, n_e) = (0, 0)</math> &#8212; now lets us derive an analytic expression for the thermal energy of the core, <math>S_\mathrm{core}</math>, and an analytic expression for the thermal energy of the envelope, <math>S_\mathrm{env}</math>.  Whether focusing on the core or the envelope, start with an appreciation that,
<div align="center">
<math>
dS = \frac{3}{2} \biggl( \frac{P}{\rho} \biggr) dm = \frac{3}{2} \biggl( \frac{P}{\rho} \biggr) \biggl( \frac{4\pi}{3} \biggr) \rho r^2 dr
= 2\pi Pr^2 dr \, .
</math>
</div>
Hence,
<div align="center">
<math>
S_\mathrm{core} = 2\pi \int_0^{r_i} [P(r)]_\mathrm{core} r^2 dr; ~~~~~\mathrm{and}~~~~~
S_\mathrm{env} = 2\pi \int_{r_i}^R [P(r)]_\mathrm{env} r^2 dr \, .
</math>
</div>
Then, for either volume segment, the relationship between the internal energy, <math>U</math>, and the thermal energy, <math>S</math>, is,
<div align="center">
<math>
U = \biggl[ \frac{2}{3(\gamma - 1)} \biggr] S \, .
</math>
</div>

My initial derivation gives,
<table border="1" align="center" width="85%">
<tr><td align="left">
<font color="darkblue">
Do not use these expressions as they are flawed.  The correct expressions follow in the comments dated 12 February 2014.
</font>
</td></tr>
<tr><td>
<math>
S_\mathrm{core} = 2\pi P_0 R^3 \biggl[q^3 - q^5 \biggl( \frac{2\pi G\rho_0^2 R^2}{5P_0} \biggr) \biggr] \, ;
</math>
</td></tr>
<tr><td>
<math>
S_\mathrm{env} = \frac{4\pi}{3} R^3 \biggl\{ \frac{3P_i}{2} (1-q^3) + \pi G \rho_0^2 R^2 \biggl(\frac{\rho_e}{\rho_c} \biggr) 
\biggl[  \frac{1}{6}(6q^3 - q^2 - 5q^5) - \frac{1}{10}\biggl(\frac{\rho_e}{\rho_c} \biggr) (2 + 10q^3 - 5q^2-7q^5)      \biggr] \biggr\} \, .
</math>
</td></tr>
</table>

<br />

<div align="center">
<hr width="85%" />
</div>

<br />
[<font color="red">12 February 2014</font>] In an [[SSC/Structure/BiPolytropes/Analytic00#Thermal_Energy_Content|accompanying discussion of the thermal energy content of an <math>(n_c, n_e) = (0, 0)</math> bipolytrope]], I've re-derived and cross-checked the expressions for <math>S_\mathrm{core}</math> and <math>S_\mathrm{env}</math>.  I am quite confident that the correct expressions are:

<div align="center">
<table border="0">
<tr>
  <td align="right">
<math>~S_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>\biggl( \frac{4\pi}{5} \biggr) R^3 q^5 \biggl (\frac{5P_i}{2q^2} + \Pi \biggr)  \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~S_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math> 
\frac{\pi R^3}{5} \biggl\{ 10 P_i  (1-q^3)  
+ 10 \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \biggl[-2q^2 + 3q^3 - q^5  \biggr]
  + 6 \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr)^2 \biggl[ -1 + 5q^2 -5q^3 +  q^5   \biggr] \biggr\} \, ,
</math>
  </td>
</tr>
</table>
</div>

where, <math>P_i</math> is the pressure at the interface,
<div align="center">
<math>
\Pi \equiv \frac{3}{2^3 \pi} 
\biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) \, ,
</math>
</div>
and the relationship between the central pressure and the pressure at the interface is,
<div align="center">
<math>~P_0 = P_i + \Pi q^2 \, .</math>
</div>

<br />

<div align="center">
<hr width="85%" />
</div>

<br />
[<font color="red">7 February 2014</font>] More succinctly, this <math>(n_c, n_e) = (0,0)</math> bipolytrope is very interesting because it provides three independent ways of deriving the proper equilibrium structure.  If we work through the detailed derivation in all three cases, we should have a very firm foundation from which to perform the stability analysis of a broad range of bipolytropic configurations. The three derivations are &hellip;
<ol>
<li>Detailed force balance (Kundan has already done this):  Use the hydrostatic balance equation to determine how the pressure varies with radius throughout the two-component (core-envelope) structure.  For a given choice of <math>P_0, \rho_0, \nu,</math> and <math>q</math>, the radius and, hence, total mass of the configuration is determined by the surface boundary condition, that is, <math>P = 0</math> at <math>r = R</math>.  These results are summarized above.
<li>Virial equilibrium:  The structures derived from demanding detailed force balance should also be in virial equilibrium.  (The virial equilibrium condition is less demanding than detailed force balance because it is a global statement of energy balance and is not concerned with internal structural details. Virial equilibrium is a necessary but not sufficient condition; the force-balance solution gives the more detailed information but must, itself, also define a structure that has virial energy balance.)  It is well known that, for the simplest spherically symmetric stars, the statement of virial balance is,
<div align="center">
<math>
2S + W = 0 \, .
</math>
</div>
Since we have analytic expressions for both <math>S</math> and <math>W</math>, we should check to see if they sum appropriately to zero as expected from this simple statement of virial equilibrium.  [Note inserted <font color="red">12 February 2014</font>:  An [[SSC/Structure/BiPolytropes/Analytic00#Virial_Equilibrium|accompanying discussion of the virial equilibrium of this particular bipolytrope]] shows that the relation does hold precisely.]
</ol>
<ol start="3">
<li>Extrema in the Free Energy:  We should also expect to be able to define the free energy of this bipolytrope by the expression,
<div align="center">
<math>
\mathfrak{G} = U + W \, .
</math>
</div>
Then, for a given choice of <math>P_0, \rho_0, \nu,</math> and <math>q</math> (perhaps we should choose <math>M_\mathrm{tot}</math> instead of <math>P_0</math>), the equilibrium radius should be defined by the radius at which,
<div align="center">
<math>
\frac{\partial \mathfrak{G}}{\partial R} = 0 \, .
</math>
</div>
We should check to make sure that this is also true.
</ol>
Items #2 and #3 are the next things I will be working on.

<br />

<div align="center">
<hr width="85%" />
</div>

<br />
[<font color="red">17 February 2014; by J. E. Tohline</font>] Using derivatives of the free energy function, <math>~\mathfrak{G}</math>, I have finished deriving the equilibrium radius for bipolytropic configurations where the density of the core <math>~(\rho_0)</math> and the density of the envelope <math>~(\rho_e)</math> are both assumed to be uniform.  I have also finished deriving the condition for the dynamical stability of such configurations assuming that, as the equilibrium structure undergoes a radial perturbation, the core compresses along a <math>~\gamma_c</math> adiabat while the envelope compresses along a <math>~\gamma_e</math> adiabat.  The [[SSC/Structure/BiPolytropes/Analytic00#Equilibrium_Condition|equilibrium radius is derived here]] and the [[SSC/Structure/BiPolytropes/Analytic00#Stability_Condition|stability condition is described here]].

<br />

<div align="center">
<hr width="85%" />
</div>

<br />
[<font color="red">25 February 2014; by J. E. Tohline</font>] I have extended the free-energy-based stability analysis to models of <math>~(n_c, n_e) = (5, 1)</math> bipolytropes.  As is [[SSC/Structure/BiPolytropes/Analytic51#Expression_for_Free_Energy|detailed here]], I have been able to derive analytic, closed form expressions for the gravitational potential energy and internal thermal energies of these configurations.  As is [[SSC/Structure/BiPolytropes/Analytic51#Equilibrium_Condition|shown here]], an examination of the first derivative of the free energy identifies equilibrium configurations that exactly match those derived from the earlier detailed force-balance derivation; and, finally, as [[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|shown here (see especially Figure 3)]], the second derivative of the free energy expression allows us to determine which equilibrium models are stable and which are dynamically unstable.  Hooray!
</td></tr>
</table>


{{ SGFworkInProgress }}

=BiPolytrope Structural Relations=
[Following a discussion that Tohline had with Kundan Kadam on <font color="red">3 July 2013</font>, we have decided to carry out a virial equilibrium and stability analysis of nonrotating bipolytropes.]  

We will adopt the following approach:
* Properties of the core <math>\cdots</math> 
** Uniform density, <math>\rho_c</math>;
** Polytropic constant, <math>K_c</math>, and polytropic index, <math>n_c</math>;
** Surface of the core at <math>r_i</math>;
* Properties of the envelope <math>\cdots</math> 
** Uniform density, <math>\rho_e</math>;
** Polytropic constant, <math>K_e</math>, and polytropic index, <math>n_e</math>;
** Base of the core at <math>r_i</math> and surface at <math>R</math>.

Use the dimensionless radius,
<div align="center">
<math>\xi \equiv \frac{r}{r_i}</math>.
</div>
Then, <math>\xi_i = 1</math> and <math>\xi_s \equiv R/r_i</math>.

==Expressions for Mass==
Inside the core, the expression for the mass interior to any radius, <math>0 \le \xi \le 1</math>, is,
<div align="center">
<math>M_\xi = \frac{4\pi}{3} \rho_c r_i^3 \xi^3</math> .
</div>
The expression for the mass interior to any position within the envelope, <math>1 \le \xi \le \xi_s</math>, is,
<div align="center">
<math>M_\xi = \frac{4\pi}{3} r_i^3 \biggl[\rho_c  + \rho_e(\xi^3 - 1) \biggr]</math> .
</div>
Hence, the mass of the core, the mass of the envelope, and the total mass are, respectively,
<div align="center">
<math>
M_\mathrm{core} = \frac{4\pi}{3} \rho_c r_i^3
</math> ;

<math>
M_\mathrm{env} = \frac{4\pi}{3} r_i^3 \biggl[\rho_e (\xi_s^3 - 1) \biggr] 
</math> ;

<math>
M_\mathrm{tot} = \frac{4\pi}{3} r_i^3 \biggl[\rho_c  + \rho_e(\xi_s^3 - 1) \biggr] 
</math> .
</div>

Following the work of [http://adsabs.harvard.edu/abs/1942ApJ....96..161S Sch&ouml;nberg &amp; Chandrasekhar (1942)] &#8212; see [[SSC/Structure/LimitingMasses#Sch.C3.B6nberg-Chandrasekhar_Mass|our accompanying discussion]]  &#8212; we are seeking equilibrium configurations in the <math>\nu - q</math> plane where,
<table align="center" border="0" cellpadding="10">
<tr>
  <td align="center">
<math>\nu \equiv \frac{M_\mathrm{core}}{M_\mathrm{tot}} </math> &nbsp; &nbsp; and &nbsp; &nbsp;
<math>q \equiv \frac{r_i}{R} = \frac{1}{\xi_s}</math>.
  </td>
</tr>
</table>
So we can combine the above expressions to obtain,
<div align="center">
<math>\frac{\rho_e}{\rho_c} =  \frac{M_\mathrm{env}}{M_\mathrm{core}} (\xi_s^3 - 1)^{-1} 
= \biggl[ \frac{1-\nu}{\nu}\biggr] (\xi_s^3 - 1)^{-1}
= \frac{q^3}{\nu}\biggl( \frac{1 - \nu}{1- q^3} \biggr) \, ,
</math>
</div>
or,
<div align="center">
<math>\nu = \biggl[ 1 + \biggl( \frac{\rho_e}{\rho_c} \biggr) \biggl(\frac{1}{q^3} - 1\biggr) \biggr]^{-1}</math> .
</div>

<table align="center" border="3" cellpadding="10" width="65%">
<tr>
  <td align="left">
It is worth noting that exactly the same result arises from an examination of the analytically definable, structural properties of [[SSC/Structure/BiPolytropes/Analytic51#Parameter_Values|bipolytropes having <math>n_c = 5</math> and <math>n_e = 1</math>]].  That is, the ratio of the average density in the envelope to the average density in the core is,
<div align="center">
<math>
\frac{\bar{\rho}_e}{\bar{\rho}_c} = 
\frac{(M^*_\mathrm{tot} - M^*_\mathrm{core}) / [(R^*)^3 - (r^*_\mathrm{core})^3]}{M^*_\mathrm{core} / (r^*_\mathrm{core})^3 } =
\frac{q^3 (1- \nu)}{\nu (1-q^3)} \, .
</math>
</div>
This is, of course, at it should be.
  </td>
</tr>
</table>


It is worth noting that, because <math>\bar\rho \equiv 3M_\mathrm{tot}/(4\pi R^3) \,</math>, we can write,
<div align="center">
<math>\frac{\rho_c}{\bar\rho} = \frac{\nu}{q^3} \, ,</math> &nbsp;&nbsp;&nbsp;and &nbsp;&nbsp;&nbsp; 
<math>\frac{\rho_e}{\bar\rho} = \frac{1-\nu}{1-q^3} \, ,</math>
</div>
which is consistent with the above expression for the ratio, <math>\rho_e/\rho_c \, .</math>  The following figure shows how <math>\nu \,</math> varies with <math>q \,</math> for various choices of the mass density ratio, <math>\rho_e/\rho_c \,</math>.  It illustrates that, for a given core-to-total mass ratio, <math>\nu \,</math>, the relative location of the interface radius, <math>q \,</math>, can vary between zero and one, but each value of <math>q \,</math> reflects a different ratio of envelope-to-core mass density.

<div align="center">
<table border="1">
<tr>
  <td align="center">
[[File:NuVersusQ.jpg|600px|center|Nu versus Q]]
  </td>
</tr>
</table>
</div>

==Energy Expressions==

The gravitational potential energy of the bipolytropic configuration is obtained by integrating over the following differential energy contribution,
<div align="center">
<math>dW = - \biggl( \frac{GM_r}{r} \biggr) dm</math> .
</div>

Hence,
<table border="0" align="center">

<tr>
  <td align="right">
<math>W = W_\mathrm{core} + W_\mathrm{env}</math>
  </td>
  <td align="left">
<math>
= - G \biggl\{ \int_0^{r_i} \biggl( \frac{M_r}{r} \biggr) 4\pi r^2 \rho_c dr + \int^R_{r_i} \biggl( \frac{M_r}{r} \biggr) 4\pi r^2 \rho_e dr \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="left">
<math>
= - G \biggl\{ \int_0^1 \biggl( \frac{4\pi }{3} \rho_c r_i^3 \xi^3 \biggr) 4\pi r_i^2 \rho_c \xi d\xi 
+ \int_1^{\xi_s} \frac{4\pi}{3} \rho_c r_i^3 \biggl[ 1 + \frac{\rho_e}{\rho_c}(\xi^3 - 1) \biggr] 4\pi r_i^2 \rho_e \xi d\xi \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="left">
<math>
= - \frac{3GM^2_\mathrm{core}}{r_i} \biggl\{ \int_0^1  \xi^4 d\xi 
+ \int_1^{\xi_s} \biggl[ 1 + \frac{\rho_e}{\rho_c}(\xi^3 - 1) \biggr] \biggl( \frac{\rho_e}{\rho_c} \biggr) \xi d\xi \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="left">
<math>
= - \frac{3GM^2_\mathrm{core}}{r_i} \biggl\{ \frac{1}{5} 
+ \biggl( \frac{\rho_e}{\rho_c} \biggr) \int_1^{\xi_s} \xi d\xi 
+ \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \int_1^{\xi_s}  (\xi^3 - 1) \xi d\xi \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="left">
<math>
= - \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr) 3\nu^2 \xi_s \biggl\{ \frac{1}{5} 
+ \frac{1}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1)
+ \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ \frac{1}{5}(\xi_s^5 - 1) - \frac{1}{2}(\xi_s^2-1) \biggr] \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="left">
<math>
= - \biggl( \frac{3GM^2_\mathrm{tot}}{5R_0} \biggr) \biggl( \frac{R}{R_0}\biggr)^{-1} \nu^2 \xi_s f(\nu,q) \, ,
</math>
  </td>
</tr>

</table>
where <math>R_0</math> is an, as yet unspecified, normalization radius, and
<div align="center">
<math>
f(\nu,q) \equiv 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1)
+ \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (\xi_s^5 - 1) - \frac{5}{2}(\xi_s^2-1) \biggr] \, .
</math>
</div>

I like the form of this expression.  The leading term, which scales as <math>R^{-1}</math>, encapsulates the behavior of the gravitational potential energy for a given choice of the internal structure, namely, a given choice of <math>\xi_s</math>, <math>\nu</math>, and density ratio <math>(\rho_e/\rho_c)</math>.  Actually, only two internal structural parameters need to be specified &#8212; <math>\nu</math> and <math>\xi_s</math> (or, <math>q</math>).  From these two, the expression shown above allows the determination of <math>(\rho_e/\rho_c)</math>.

<table border="3" cellpadding="10" align="center">
<tr>
  <td align="left">

'''Sanity Check:'''  Uniform Density Configuration
<div align="center">
<math>
M_\mathrm{core} = \frac{4\pi}{3} r_i^3 \rho_c \, ;
</math>

<math>
M_\mathrm{env} = \frac{4\pi}{3} r_i^3 \rho_e (\xi_s^3 -1) \, ;
</math>

<math>
M_\mathrm{tot} = \frac{4\pi}{3} r_i^3 \rho_c \biggl[ 1 + \frac{\rho_e}{\rho_c} (\xi_s^3 -1) \biggr] \, .
</math>
</div>
If <math>\rho_e/\rho_c = 1</math>, then,
<div align="center">
<math>
\frac{M_\mathrm{env}}{M_\mathrm{core}} = \biggl(\frac{1}{\nu} - 1\biggr) = (\xi_s^3 -1) = \biggl(\frac{1}{q^3}-1\biggr) 
~~~~\Rightarrow ~~~~ \nu = q^3 \, .
</math>
</div>

The gravitational potential energy is,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>W\biggr|_{\rho_e/\rho_c = 1}</math>
  </td>
  <td align="left">
<math>
= - \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr) 3\nu^2 \xi_s \biggl\{ \frac{1}{5} 
+ \frac{1}{2}  (\xi_s^2 - 1)
+ \biggl[ \frac{1}{5}(\xi_s^5 - 1) - \frac{1}{2}(\xi_s^2-1) \biggr] \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="left">
<math>
= - \frac{3}{5} \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr)\nu^2 \xi_s \biggl\{ \frac{1}{5} 
+ \biggl[ \frac{1}{5}(\xi_s^5 - 1)  \biggr] \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="left">
<math>
= - \frac{3}{5} \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr)\nu^2 \xi_s^6 
= - \frac{3}{5} \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr) \biggl( \frac{\nu}{q^3}\biggr)^2 =  - \frac{3}{5} \biggl( \frac{GM^2_\mathrm{tot}}{R} \biggr)  \, .
</math>
  </td>
</tr>
</table>
</div>

  </td>
</tr>
</table>

Drawing on expressions developed in our [[VE#Virial_Equation|introductory discussion of the virial equation]], the internal energy of the bipolytropic configuration is,
<table border="0" align="center">
<tr>
  <td align="right">
<math>
U = U_\mathrm{core} + U_\mathrm{env}
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\biggl\{ M_\mathrm{core} \biggl[ (1 - \delta_{\infty n_c}) n_c K_c \rho_c^{1/n_c}  + \delta_{\infty n_c} c_s^2 \ln(\rho_c/\rho_0) \biggr] 
+ n_e M_\mathrm{env} K_e \rho_e^{1/n_e}  \biggr\} \, ,
</math>
  </td>
</tr>
</table>
where <math>\rho_0</math> is an, as yet unspecified, normalization density, and we have allowed for either an isothermal (<math>\delta_{\infty n_c} = 1</math>) or an adiabatic (<math>\delta_{\infty n_c} = 0</math>) core.  


=Strategy=
==Scaling Parameters==
We want to vary the total radius, <math>R \,</math>, of the configuration and look for extrema in the free energy, while holding the following parameters fixed:  <math>M_\mathrm{tot} \,</math>, <math>c_s^2 \,</math> (or <math>K_c \,</math>), <math>\nu \,</math>, and <math>q \,</math>.  So we need to rewrite the expressions for <math>W \,</math> and <math>U \,</math> such that everything is constant except for <math>R \,</math>, or <math>\chi \equiv R/R_0 \,</math>.  And, in order to put everything explicitly in terms of the fixed parameters just specified, let's go ahead and define the length, time, and density against which all dimensional quantities will ultimately be scaled.

<table border="1" cellpadding="10" align="center">
<tr>
  <td colspan="3" align="center">
<font size="+1"><b>Chosen Scaling Parameters</b></font>
  </td>
</tr>

<tr>
  <th>
&nbsp;
  </th>
  <th align="center">
Polytropic Core
  </th>
  <th align="center">
Isothermal Core
  </th>
</tr>

<tr>
   <td align="center">
<math>R_0 \,</math>
  </td>
   <td align="center">
<math>\biggl[ \biggl( \frac{K_c}{G} \biggr)^{n_c} M_\mathrm{tot}^{1-n_c} \biggr]^{1/(3-n_c)} \,</math>
  </td>
   <td align="center">
<math>\frac{GM_\mathrm{tot}}{c_s^2} \,</math>
  </td>
</tr>

<tr>
   <td align="center">
<math>t_0 \,</math>
  </td>
   <td align="center">
<math>\biggl[ \biggl( \frac{K_c^3}{M_\mathrm{tot}^2} \biggr)^{n_c} G^{-(2n_c+3)} \biggr]^{1/(6-2n_c)} \,</math>
 </td>
   <td align="center">
<math>\frac{GM_\mathrm{tot}}{c_s^3} \,</math>
  </td>
</tr>

<tr>
   <td align="center">
<math>\rho_0 \equiv \frac{3M_\mathrm{tot}}{4\pi R_0^3} \,</math>
  </td>
   <td align="center">
<math>\frac{3}{4\pi} \biggl[ \frac{G^3 M_\mathrm{tot}^2}{K_c^3} \biggr]^{n_c/(3-n_c)} \,</math>
  </td>
   <td align="center">
<math>\frac{3}{4\pi} \biggl[ \frac{c_s^6}{G^3 M_\mathrm{tot}^2} \biggr] \,</math>
  </td>
</tr>

</table>


==Rescaled Energy Expressions==
===Isothermal Core===
In the case of an isothermal core, the expression for <math>W</math> shown above needs only a slight modification to put it in the appropriate form, namely,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>W \,</math> 
  </td>
  <td align="left">
<math>
= - \frac{3}{5}  M_\mathrm{tot} c_s^2 \chi^{-1} \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, ,
</math>
  </td>
</tr>

</table>

In the case of an isothermal core, the expression for the total internal energy may be rewritten as,
<table border="0" align="center" cellpadding="5">
<tr>
  <td align="right">
<math>
\frac{U}{M_\mathrm{tot}} \,
</math>
  </td>
  <td align="center">
<math>= \,</math>
  </td>
  <td align="left">
<math>
\nu c_s^2 \ln(\rho_c/\rho_0)  + (1-\nu) n_e K_e \rho_e^{1/n_e}  
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\nu c_s^2 [\ln(\rho_c/\bar{\rho}) + \ln(\bar\rho/\rho_0)]   + (1-\nu) n_e K_e \rho_0^{1/n_e} \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr) \biggr]^{1/n_e} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\nu c_s^2 [\ln(\nu/q^3) -3 \ln\chi]   + (1-\nu) n_e K_e \rho_0^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \chi^{-3/n_e} \, .
</math>
  </td>
</tr>
</table>
Hence (in the case of an isothermal core),

<table border="0" align="center" cellpadding="5">
<tr>
  <td align="right">
<math>\mathfrak{G}^* \equiv \frac{\mathfrak{G}}{M_\mathrm{tot} c_s^2}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
 - \frac{3}{5}  \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \chi^{-1} 
+ \nu [\ln(\nu/q^3) -3 \ln\chi]   + (1-\nu) n_e \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \chi^{-3/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\mathfrak{G}_0^* - A \chi^{-1} - C \ln\chi  + D \chi^{-3/n_e} \, ,
</math>
  </td>
</tr>

</table>
where,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\mathfrak{G}_0^* \,</math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
\nu \ln(\nu/q^3)   \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>A \,</math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
\frac{3}{5} \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>C \,</math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
3 \nu  \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>D \, </math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
(1-\nu) n_e \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} 
= (1-\nu) n_e \biggl( \frac{K_e}{c_s^2} \biggr) \biggl[  \frac{3}{4\pi} \biggl(\frac{c_s^6}{G^3M^2_\mathrm{tot}}\biggr)   \biggr]^{1/n_e} \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} \, .
</math>
  </td>
</tr>
</table>

===Polytropic Core===
In the case of a polytropic core,

<table border="0" align="center" cellpadding="5">
<tr>
  <td align="right">
<math>\mathfrak{G}^* \equiv \frac{\mathfrak{G}}{M_\mathrm{tot} K_c \rho_0^{1/n_c}}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
- A \chi^{-1}  + B_c \chi^{-3/n_c} + B_e \chi^{-3/n_e} \, ,
</math>
  </td>
</tr>

</table>
where,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>A \,</math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
\frac{3}{5} \biggl[ \frac{GM_\mathrm{tot}}{R_0 K_c \rho_0^{1/n_c}} \biggr] \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) =
\frac{3}{5} \biggl( \frac{4\pi}{3} \biggr)^{1/n_c} \biggl( \frac{\nu^2}{q} \biggr) f(\nu,q) \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>B_c \,</math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
n_c \nu \biggl( \frac{\nu}{q^3}\biggr)^{1/n_c} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>B_e \, </math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
n_e (1-\nu) \biggl[ \biggl( \frac{K_e }{K_c}\biggr)\rho_0^{1/n_e - 1/n_c} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} 
\, .
</math>
  </td>
</tr>
</table>

==Temperature Across the Interface==
In order to ensure that the temperature of the envelope is the same as the temperature of the core when there is a drop in the mean molecular weight at the interface, we need to have,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>\biggl[ \frac{P_e}{(\rho_e/\mu_e)}\biggr]_E</math>
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\biggl[ \frac{P_c}{(\rho_c/\mu_c)} \biggr]_E \, .</math>
  </td>
</tr>
</table>
Note that this reflects the same physical condition as the constraint that is placed on the enthalpy at the interface when we analyze the detailed structure of <math>n_c = 5</math>, <math>n_e=1</math> bipolytropes (see the [[SSC/Structure/BiPolytropes/Analytic51#Parameter_Values|middle of Table 1 in the accompanying discussion]]), namely,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>\biggl[ \frac{\mu_e H_e}{(n_e+1)}\biggr]_i</math>
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\biggl[ \frac{\mu_c H_c}{(n_c+1)} \biggr]_i \, .</math>
  </td>
</tr>
</table>


In the case of an isothermal core, this implies,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e}{ c_s^2}</math>
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\biggl[\rho_e^{-1/n_e}\biggr]_E</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~~ \biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e \rho_0^{1/n_e}}{ c_s^2}</math>
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math> \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E \biggr]^{-1/n_e}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>
\biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E^{-1/n_e}
= \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \chi_E^{3/n_e}
</math>
  </td>
</tr>
</table>

<table border="1" cellpadding="5" align="center">
<tr>
  <td colspan="3" align="center">
Not Necessarily Useful
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~~ \chi_E^{3/n_e}</math>
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e}{ c_s^2} \biggl[  \frac{3}{4\pi} \biggl(\frac{c_s^6}{G^3M^2_\mathrm{tot}}\biggr)   \biggr]^{1/n_e} 
\biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \biggl[  \frac{D}{(1-\nu)n_e M_\mathrm{tot} c_s^2}  \biggr]</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \biggl[  \frac{3\nu}{(1-\nu)n_e}  \biggl(\frac{D}{C}\biggr)\biggr]</math>
  </td>
</tr>
</table>
Therefore, the coefficient <math>D\,</math> becomes,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>D \, </math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
n_e (1-\nu) \biggl[ \frac{K_e \rho_0^{1/n_e}}{c_s^2} \biggr] \biggl( \frac{1-\nu}{1-q^3} \biggr)^{1/n_e} 
= n_e (1-\nu) \biggl( \frac{\mu_c }{\mu_e}\biggr) \chi_E^{3/n_e} \, .
</math>
  </td>
</tr>

</table>


In the case of a polytropic core, this implies,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>\biggl( \frac{\mu_e}{\mu_c}\biggr) \frac{K_e}{ K_c}</math>
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>\rho_0^{1/n_c} \biggl[\rho_e^{-1/n_e}\biggr]_E</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~~ \biggl( \frac{\mu_e}{\mu_c}\biggr) \biggl( \frac{K_e }{ K_c} \biggr) \rho_0^{(1/n_e-1/n_c)}</math>
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math> \biggl[\biggl( \frac{\rho_e}{\bar\rho} \biggr) \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E \biggr]^{-1/n_e}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
   <td align="center">
<math>= \, </math>
  </td>
  <td align="left">
<math>
\biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \biggl( \frac{\bar\rho}{\rho_0} \biggr)_E^{-1/n_e}
= \biggl( \frac{1-\nu}{1-q^3} \biggr)^{-1/n_e} \chi_E^{3/n_e}
</math>
  </td>
</tr>

</table>
Therefore, the coefficient <math>B_e\,</math> becomes,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>B_e \, </math>
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
n_e (1-\nu) \biggl( \frac{\mu_c }{\mu_e}\biggr) \chi_E^{3/n_e} \, .
</math>
  </td>
</tr>

</table>


==Pressure Across the Interface==

We will relate <math>K_e</math> to <math>K_c</math> by demanding that initially the pressure is identical in both layers.  
The relevant algebraic relation will depend on whether the core is isothermal (<math>n_c = \infty</math>), or whether it
has a finite polytropic index and therefore adjusts adiabatically to compressions or expansions.

===Isothermal Core===
Guided by the interface conditions presented in [[SSC/Structure/BiPolytropes#Interface_Conditions|Table 2 of our accompanying discussion of the structure of bipolytropes]], the condition for pressure balance in the case of an isothermal core should be,
<div align="center">
<math>\frac{c_s^2}{K_e \rho_0^{1/n_e}} = \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e} \biggl( \frac{\rho_e}{\rho_c} \biggr)^{1+1/n_e} 
= \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \, .</math>
</div>

Actually, in the full structural solution,
<div align="center">
<table border="0" cellpadding="3">

<tr>
  <td align="right">
<math>\frac{c_s^2}{K_e \rho_0^{1/n_e}}</math>
  </td>
  <td align="center">
<math>= </math>
  </td>
  <td align="left">
<math>
\biggl( \frac{\rho_c}{\rho_0} \biggr)^{1/n_e} e^{\psi_i}
\biggl( \frac{\rho_e}{\rho_c} \biggr)^{1+1/n_e} \phi_i^{1+n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>= </math>
  </td>
  <td align="left">
<math>
\biggl( \frac{\rho_c}{\rho_0} \biggr)^{1/n_e} e^{\psi_i}
\biggl[ \frac{\mu_e}{\mu_c} e^{-\psi_i} \phi_i^{-{n_e}}\biggr]^{1+1/n_e} \phi_i^{1+n_e} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>= </math>
  </td>
  <td align="left">
<math>
\biggl( \frac{\rho_c}{\rho_0} e^{-\psi_i} \biggr)^{1/n_e} 
\biggl[ \frac{\mu_e}{\mu_c} \biggr]^{1+1/n_e} \, .
</math>
  </td>
</tr>
</table>
</div>

===Adiabatic Core===
Guided by the interface conditions presented in [[SSC/Structure/BiPolytropes#Interface_Conditions|Table 2 of our accompanying discussion of the structure of bipolytropes]], the condition for pressure balance in the case of polytropic core should be,
<div align="center">
<math>\biggl(\frac{K_c}{K_e}\biggr) \rho_0^{1/n_c -1/n_e} = \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e - 1/n_c} \biggl( \frac{\rho_e}{\rho_c} \biggr)^{1+1/n_e} 
= \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e - 1/n_c} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \, .</math>
</div>

Actually, in the full structural solution,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>\biggl(\frac{K_c}{K_e}\biggr) \rho_0^{1/n_c -1/n_e} </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math> 
\rho_0^{1/n_c -1/n_e} \rho_c^{-(1+1/n_c)} \theta_i^{-(1+n_c)} \rho_e^{1+1/n_e} \phi_i^{1+n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math> 
\biggl( \frac{\rho_c}{\rho_0} \biggr)^{1/n_e - 1/n_c} \theta_i^{-(1+n_c)} \biggl( \frac{\rho_e}{\rho_c}\biggr)^{1+1/n_e} \phi_i^{1+n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math> 
\biggl( \frac{\rho_c}{\rho_0} \biggr)^{1/n_e - 1/n_c} \biggl[ \frac{\mu_e}{\mu_c} \theta_i^{n_c} 
 \biggr]^{1+1/n_e} \theta_i^{-(1+n_c)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math> 
\biggl( \frac{\rho_c}{\rho_0} \theta_i^{n_c}\biggr)^{1/n_e - 1/n_c} \biggl[ \frac{\mu_e}{\mu_c}  
 \biggr]^{1+1/n_e} \, .
</math>
  </td>
</tr>

</table>
</div>

=Virial Analysis=

==Isothermal Core Equilibrium Condition==
For a given set of fixed coefficients in the free energy expression, equilibria are identified by setting <math>\partial\mathfrak{G}/\partial\chi = 0</math>.  Generally,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>\frac{\partial\mathfrak{G}}{\partial\chi}</math> 
  </td>
 <td align="center">
<math>= \,</math>
  </td>
 <td align="left">
<math>
A\chi^{-2} - C\chi^{-1} - \biggl(\frac{3}{n_e}\biggr)D\chi^{-(1+3/n_e)} \, .
</math>
  </td>
</tr>
</table>
So, for the case of <math>n_e = 3/2</math>, the equilibrium radius is given by the condition,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>A\chi_E^{-2}</math> 
  </td>
 <td align="center">
<math>= \,</math>
  </td>
 <td align="left">
<math>
 C\chi_E^{-1} + 2D\chi_E^{-3}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~ C\chi_E^{2} - A\chi_E + 2D
</math> 
  </td>
 <td align="center">
<math>= \,</math>
  </td>
 <td align="left">
<math>0 \,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~ \chi_E
</math> 
  </td>
 <td align="center">
<math>= \,</math>
  </td>
 <td align="left">
<math>
\frac{1}{2C} \biggl[ A \pm \biggl( A^2 - 8DC \biggr)^{1/2}\biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~ \chi_E
</math> 
  </td>
 <td align="center">
<math>= \,</math>
  </td>
 <td align="left">
<math>
\frac{A}{2C} \biggl[ 1 \pm \biggl( 1 - \frac{8DC}{A^2} \biggr)^{1/2}\biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp; 
  </td>
 <td align="center">
<math>= \,</math>
  </td>
 <td align="left">
<math>
\chi_0 [ 1 \pm \sqrt{1-\Pi} ] \, ,
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\chi_0 \,</math> 
  </td>
 <td align="center">
<math>\equiv \,</math>
  </td>
 <td align="left">
<math>
\frac{A}{2C}  \, ,
</math>
  </td>
  <td align="center">
&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;
  </td>

  <td align="right">
<math>\Pi \,</math> 
  </td>
 <td align="center">
<math>\equiv \,</math>
  </td>
 <td align="left">
<math>
\frac{8DC}{A^2}  \, .
</math>
  </td>
</tr>
</table>

==Combined Equilibrium Constraints==
But, from the condition on the temperature at the interface we also need,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\chi_E \,</math> 
  </td>
 <td align="center">
<math>=\,</math>
  </td>
 <td align="left">
<math>
\biggl(\frac{\mu_e}{\mu_c} \biggr)^{1/2} \biggl[ \frac{2\nu}{(1-\nu)} \biggl( \frac{D}{C} \biggr)\biggr]^{1/2}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp; 
  </td>
 <td align="center">
<math>=\,</math>
  </td>
 <td align="left">
<math>
\chi_0 \biggl(\frac{\mu_e}{\mu_c} \biggr)^{1/2} \biggl[ \biggl(\frac{\nu}{1-\nu}\biggr) \Pi \biggr]^{1/2} \, .
</math>
  </td>
</tr>
</table>

Hence, the two conditions combined imply,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>1 \pm \sqrt{1-\Pi}</math> 
  </td>
 <td align="center">
<math>=\,</math>
  </td>
 <td align="left">
<math>
\lambda^{1/2} \Pi^{1/2} \, ,
</math>
  </td>
</tr>
</table>
where,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>\lambda \,</math> 
  </td>
 <td align="center">
<math>\equiv\,</math>
  </td>
 <td align="left">
<math>
\biggl(\frac{\mu_e}{\mu_c} \biggr) \biggl(\frac{\nu}{1-\nu}\biggr) \, .
</math>
  </td>
</tr>
</table>
This, in turn, implies,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>1-\Pi \, </math> 
  </td>
 <td align="center">
<math>=\,</math>
  </td>
 <td align="left">
<math>
(\lambda^{1/2} \Pi^{1/2} -1 )^2\, 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp; 
  </td>
 <td align="center">
<math>=\,</math>
  </td>
 <td align="left">
<math>
\lambda\Pi -2(\lambda\Pi)^{1/2} + 1 \, 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~2(\lambda\Pi)^{1/2} \,</math> 
  </td>
 <td align="center">
<math>=\,</math>
  </td>
 <td align="left">
<math>
\Pi(1+\lambda) \, 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~\Pi  \,</math> 
  </td>
 <td align="center">
<math>=\,</math>
  </td>
 <td align="left">
<math>
\frac{4\lambda}{(1+\lambda)^2} \, .
</math>
  </td>
</tr>
</table>

But notice, as well, that in order to have real roots of the equilibrium condition, we need <math>\Pi \le 1</math>.  This means,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>4\lambda \,</math> 
  </td>
 <td align="center">
<math>\le \,</math>
  </td>
 <td align="left">
<math>
(1+\lambda)^2 \, 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~ \lambda^2 - 2\lambda + 1  \,</math> 
  </td>
 <td align="center">
<math>\ge \,</math>
  </td>
 <td align="left">
<math>
0 \, .
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~ (\lambda - 1)^2  \,</math> 
  </td>
 <td align="center">
<math>\ge \,</math>
  </td>
 <td align="left">
<math>
0 \, 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~ \lambda  \,</math> 
  </td>
 <td align="center">
<math>\ge \,</math>
  </td>
 <td align="left">
<math>
1 \,  
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~ \frac{\mu_e}{\mu_c}  \,</math> 
  </td>
 <td align="center">
<math>\ge \,</math>
  </td>
 <td align="left">
<math>
\frac{1}{\nu} - 1 \,  
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~~ \nu  \,</math> 
  </td>
 <td align="center">
<math>\ge \,</math>
  </td>
 <td align="left">
<math>
\biggl(\frac{\mu_e}{\mu_c} + 1\biggr)^{-1} \,  .
</math>
  </td>
</tr>
</table>

==Free Energy Expression==
To within an additive constant, the free energy may now be written as,
<div align="center">
<math>
\mathfrak{G} = W + U  = - A \chi^{-1} + B_e \chi^{-3/n_e} + (1-\delta_{\infty n_c}) B_c \chi^{-3/n_c} - \delta_{\infty n_c} \biggl[ B_I \ln\chi  - \frac{1}{3}B_I \ln \biggl( \frac{\rho_c|_0}{\rho_0} \biggr) \biggr] \, ,
</math>
</div>
where, <math>\chi \equiv R/R_0</math> and,
<table border="0" align="center" cellpadding="5">

<tr>
  <td align="right">
<math>
A \,
</math>  
  </td>
  <td align="center">
<math>=\,</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{3GM^2_\mathrm{tot}}{5R_0} \biggr) \nu^2 \xi_s \biggl\{ 1 
+ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1)
+ \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (\xi_s^5 - 1) - \frac{5}{2}(\xi_s^2-1) \biggr] \biggr\} 
= \biggl( \frac{3GM^2_\mathrm{tot}}{5R_0} \biggr) \frac{\nu^2}{q}  f(\nu, q) \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>B_e \,</math>  
</td>
  <td align="center">
<math>=\,</math>
  </td>
  <td align="left">
<math>
M_\mathrm{tot}(1-\nu) n_eK_e (\rho_e|_0)^{1/n_e} = M_\mathrm{tot} (1-\nu) K_e\rho_0^{1/n_e} n_e \biggl( \frac{\rho_e|_0}{\rho_0} \biggr)^{1/n_e} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>B_c \,</math>  
</td>
  <td align="center">
<math>= \,</math>
  </td>
  <td align="left">
<math>
M_\mathrm{tot} \nu n_c K_c (\rho_c|_0)^{1/n_c} = M_\mathrm{tot} K_e\rho_0^{1/n_e} \nu n_c \biggl[\biggl(\frac{K_c}{K_e}\biggr) \rho_0^{1/n_c - 1/n_e} \biggr] \biggl( \frac{\rho_c|_0}{\rho_0}\biggr)^{1/n_c} 
= M_\mathrm{tot} K_e\rho_0^{1/n_e} \nu n_c \biggl[ \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e - 1/n_c} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \biggr] \biggl( \frac{\rho_c|_0}{\rho_0}\biggr)^{1/n_c}  \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>B_I\,</math>  
</td>
  <td align="center">
<math>=\,</math>
  </td>
  <td align="left">
<math>
3 M_\mathrm{tot} c_s^2 \nu = M_\mathrm{tot} K_e\rho_0^{1/n_e} \biggl[ \frac{c_s^2}{K_e\rho_0^{1/n_e}  } \biggr] 3\nu 
= M_\mathrm{tot} K_e\rho_0^{1/n_e} \biggl[ \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e}\biggr] 3\nu \, .
</math>
  </td>
</tr>
</table>

==Derivatives of Free Energy==
<div align="center">
<math>
\frac{\partial\mathfrak{G}}{\partial \chi} = A \chi^{-2} -\frac{3}{n_e} B_e \chi^{-(1+3/n_e)} - (1-\delta_{\infty n_c}) \frac{3}{n_c} B_c \chi^{-(1+3/n_c)} 
- \delta_{\infty n_c} B_I \chi^{-1}  \, ;
</math>

<math>
\frac{\partial^2\mathfrak{G}}{\partial \chi^2} = -2 A \chi^{-3} + \frac{3}{n_e} \biggl(1+\frac{3}{n_e}\biggr) B_e \chi^{-(2+3/n_e)} + (1-\delta_{\infty n_c}) \frac{3}{n_c} \biggl(1+\frac{3}{n_c}\biggr) B_c \chi^{-(2+3/n_c)} 
+ \delta_{\infty n_c} B_I \chi^{-2} \, .
</math>
</div>

==Equilibrium Condition==
We obtain the equilibrium radius, <math>\chi_E</math>, when <math>\partial\mathfrak{G}/\partial\chi = 0</math>.  Hence,
the relation governing the equilibrium radius is,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
A \chi_E^{-2} 
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
(1-\delta_{\infty n_c}) \frac{3}{n_c} B_c \chi_E^{-1- 3/n_c} 
+\delta_{\infty n_c} B_I \chi_E^{-1} +\frac{3}{n_e} B_e \chi_E^{-1-3/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~~ \frac{n_e A}{3B_e}
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
(1-\delta_{\infty n_c}) \frac{n_e B_c}{n_c B_e} \chi_E^{1- 3/n_c} 
+\delta_{\infty n_c} \frac{n_e B_I}{3B_e} \chi_E + \chi_E^{1-3/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~~ \chi_E^{1-3/n_e}
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\alpha - (1-\delta_{\infty n_c}) \beta \chi_E^{1- 3/n_c} - \delta_{\infty n_c} \beta_I \chi_E \, ,
</math>
  </td>
</tr>

</table>
</div>
where,
<div align="center">
<table border="0" cellpadding="5">

<tr>
  <td align="right">
<math>\alpha \,</math> 
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
\frac{n_e A}{3B_e} = \biggl[ \frac{GM_\mathrm{tot}}{5R_0 K_e (\rho_e|_0)^{1/n_e} } \biggr]  \frac{\nu^2 f(\nu, q)}{q} 
= \biggl[ \frac{GM_\mathrm{tot}}{5R_0 K_e \rho_0^{1/n_e} } \biggr]\biggl( \frac{\rho_e|_0}{\rho_0} \biggr)^{-1/n_e}  \frac{\nu^2 f(\nu, q)}{q} 
 \, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\beta \,</math> 
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>\frac{n_e B_c}{n_c B_e} = \nu \biggl[ \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e - 1/n_c}  \biggl( \frac{\rho_c|_0}{\rho_0}\biggr)^{1/n_c} 
\biggl( \frac{\rho_e|_0}{\rho_0} \biggr)^{-1/n_e} \biggr] \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e} \, ; </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\beta_I \,</math> 
  </td>
  <td align="center">
<math>\equiv \,</math>
  </td>
  <td align="left">
<math>
\frac{n_e B_I}{3B_e} = \frac{c_s^2 \nu}{(1-\nu)K_e (\rho_e|_0)^{1/n_e}} = 
\nu \biggl[ \biggl(\frac{\rho_c}{\rho_0}\biggr)^{1/n_e} \biggl( \frac{\rho_e|_0}{\rho_0} \biggr)^{-1/n_e}
\biggr] \biggl( \frac{\mu_e}{\mu_c} \biggr)^{1+1/n_e}
\, .
</math>
  </td>
</tr>
</table>
</div>

===Isothermal Core===
In the case of an isothermal core (<math>\delta_{\infty n_c} = 1</math>),

<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
A \chi_E^{-2} 
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
B_I \chi_E^{-1} +\frac{3}{n_e} B_e \chi_E^{-1-3/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~ B_I \chi_E^{3/n_e} -  A \chi_E^{3/n_e-1} +\frac{3}{n_e} B_e 
</math>
  </td>
  <td align="center"><math>=\,</math></td>
  <td align="left">
<math>
0 \, .
</math>
  </td>
</tr>

</table>
</div>

Or, alternatively, via a dimensionless treatment,

<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\chi_E^{1-3/n_e}
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\alpha - \beta_I \chi_E \, .
</math>
  </td>
</tr>

</table>
</div>
===Adiabatic Core===
In the case of an adiabatice core (<math>\delta_{\infty n_c} = 0</math>),

<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\chi_E^{1-3/n_e}
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\alpha - \beta \chi_E^{1- 3/n_c}  \, .
</math>
  </td>
</tr>

</table>
</div>

==Stability==
At this equilibrium radius, the second derivative of the free energy has the value,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\chi_E^3 \biggl( \frac{n_e}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E
</math>
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
-2 A\biggl( \frac{n_e}{3B_e} \biggr)  + (1-\delta_{\infty n_c}) \biggl( \frac{n_e}{3B_e} \biggr) \frac{3}{n_c} \biggl(1+\frac{3}{n_c}\biggr) B_c \chi_E^{1-3/n_c} 
+ \delta_{\infty n_c} \biggl( \frac{n_e}{3B_e} \biggr) B_I \chi_E + \frac{3}{n_e}\biggl( \frac{n_e}{3B_e} \biggr)  \biggl(1+\frac{3}{n_e}\biggr) B_e \chi_E^{1-3/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
-2 \alpha + (1-\delta_{\infty n_c}) \beta \biggl(1+\frac{3}{n_c}\biggr) \chi_E^{1-3/n_c} 
+ \delta_{\infty n_c} \beta_I \chi_E + \biggl(1+\frac{3}{n_e}\biggr) \chi_E^{1-3/n_e} \, ,
</math>
  </td>
</tr>

</table>
</div>
which, when combined with the condition for equilibrium gives,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\chi_E^3 \biggl( \frac{n_e}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E
</math>
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
-2 \alpha + (1-\delta_{\infty n_c}) \beta \biggl(1+\frac{3}{n_c}\biggr) \chi_E^{1-3/n_c} 
+ \delta_{\infty n_c} \beta_I \chi_E + \biggl(1+\frac{3}{n_e}\biggr) \biggl[  \alpha - (1-\delta_{\infty n_c}) \beta \chi_E^{1- 3/n_c} - \delta_{\infty n_c} \beta_I \chi_E \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\alpha \biggl(\frac{3}{n_e}-1\biggr) + (1-\delta_{\infty n_c}) \beta \biggl(\frac{3}{n_c}-\frac{3}{n_e}\biggr) \chi_E^{1-3/n_c} 
-  \delta_{\infty n_c} \beta_I \biggl(\frac{3}{n_e}\biggr)\chi_E 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~ \chi_E^3 \biggl( \frac{n_e^2}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E
</math>
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\alpha (3-n_e) - (1-\delta_{\infty n_c}) 3\beta \biggl(1 - \frac{n_e}{n_c} \biggr) \chi_E^{1-3/n_c} -  \delta_{\infty n_c} 3\beta_I \chi_E \, .
</math>
  </td>
</tr>

</table>
</div>

The equilibrium configuration is stable as long as this second derivative is positive.  

===Isothermal Core===
Hence, for a bipolytrope with an isothermal core (<math>\delta_{\infty n_c} = 1</math>), the configuration is stable as long as,
<div align="center">
<math>
\chi_E <  \frac{\alpha (3-n_e)}{3\beta_I} \, .
</math>
</div>

===Adiabatic Core===
In the adiabatic case (<math>\delta_{\infty n_c} = 0</math>), the configuration is stable as long as,
<div align="center">
<math>
\chi_E^{1-3/n_c} <  \frac{\alpha n_c (3-n_e)}{3\beta (n_c-n_e)} \, .
</math>
</div>


==Examples==
===Isothermal Core with <math>n_e=3/2</math>===
Consider the case examined by [http://adsabs.harvard.edu/abs/1942ApJ....96..161S Sch&ouml;nberg &amp; Chandrasekhar (1942)], that is, the case of an isothermal core and an envelope with <math>n_e = 3/2</math>.  The equilibrium radius is given by the expression,

<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\chi_E^{-1}
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\alpha - \beta_I \chi_E \,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~~ \beta_I \chi_E^{2} -\alpha \chi_E + 1
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
0 \,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~~ \chi_E
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{1}{2\beta_I} \biggl[\alpha \pm \biggr(\alpha^2-4\beta_I \biggr)^{1/2}\biggr] \,
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{\alpha}{2\beta_I} \biggl[1 \pm \biggr(1-\frac{4\beta_I}{\alpha^2} \biggr)^{1/2}\biggr] \, .
</math>
  </td>
</tr>

</table>
</div>
And the system is stable when,
<div align="center">
<math>
\chi_E < \chi_0 \equiv  \frac{\alpha }{2\beta_I} \, .
</math>
</div>

A couple of physical attributes are now clear:
*Physical configurations only exist for <math>(4\beta_I/\alpha^2) \le 1</math>.
*For each value of <math>(4\beta_I/\alpha^2) < 1 \,</math>, there are two equilibrium configurations, given by the <math>\pm</math> roots of the quadratic equation for <math>\chi_E</math>; the "negative" branch is stable but the "positive" branch is unstable.

Note that,
<div align="center">
<math>
\chi_0 \equiv \frac{\alpha }{2\beta_I}  = \biggl( \frac{GM_\mathrm{tot}}{10R_0 c_s^2} \biggr) \frac{\nu f(\nu,q)}{q} ~~~\Rightarrow ~~~ 
\biggl( \frac{10R_0 c_s^2}{GM_\mathrm{tot}} \biggr) = \frac{\nu f(\nu,q)}{q \chi_0}  \, ,
</math>
</div>
and,
<div align="center">
<math>
\frac{4\beta_I}{\alpha^2 }  = \biggl( \frac{10 R_0 c_s^2 }{GM_\mathrm{tot}} \biggr)^2 \biggl[ \frac{K_e (\rho_e|_0)^{1/n_e}}{c_s^2} \biggr] \frac{q^2(1-\nu)}{\nu^3 f^2(\nu,q)}  
= \biggl[ \frac{K_e (\rho_e|_0)^{1/n_e}}{c_s^2} \biggr] \frac{q^2(1-\nu)}{\nu^3 f^2(\nu,q)} \biggl[\frac{\nu f(\nu,q)}{q \chi_0} \biggr]^2 
= \biggl[ \frac{\mu_c}{\mu_e} \biggr] \biggl( \frac{1}{\nu} - 1 \biggr) \, .
</math>
</div>
But, this last expression must be less than or equal to unity, which implies,
<div align="center">
<math>
\frac{1}{\nu} \le 1 + \frac{\mu_e}{\mu_c} ~~~\Rightarrow ~~~ \nu \ge \biggl(1 + \frac{\mu_e}{\mu_c}  \biggr)^{-1}
</math>
This doesn't seem to have the correct behavior, for example, the smaller values of <math>\nu</math> should be the stable ones, so there must be a mistake in the derivation.
</div>

===Adiabatic Core with <math>n_c = 5</math> and <math>n_e=1</math>===
Consider the case with an analytical structure derived by  [http://adsabs.harvard.edu/abs/1998MNRAS.298..831E Eagleton, Faulkner, and Cannon] (1998, MNRAS, 298, 831), that is, the case of an adiabatic core having <math>n_c=5</math> and an envelope with <math>n_e = 1</math>.  The equilibrium radius is,

=Old and Probably Irrelevant Discussion=

==Summary Expressions (New)==
In the above derivations, we have adopted the notation, 
<div align="center">
<math>
\rho_\mathrm{norm} \equiv \frac{3M_\mathrm{tot}}{4\pi R_0^3} \, .
</math>
</div>
Now, guided by the [[SSCpt1/Virial#Nonrotating_Configuration_Embedded_in_an_External_Medium|earlier discussion of pressure-bounded isothermal spheres]], we choose the following normalization energy and radius:
<div align="center">
<math>
E_0 = 3M_\mathrm{tot} c_s^2 
</math>
&nbsp; &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; &nbsp;
<math>
R_0 = \frac{GM_\mathrm{tot}}{5c_s^2} \, .
</math>
</div>
Also, by analogy, it is useful to define the dimensionless parameter,
<div align="center">
<math>
\Pi_I \equiv \frac{K_e \rho_\mathrm{norm}^{1/n_e}}{c_s^2}
= \frac{K_e}{c_s^2} \biggl[ \frac{3M_\mathrm{tot}}{4\pi R_0^3} \biggr]^{1/n_e}
= \biggl( \frac{3\cdot 5^3}{4\pi} \biggr)^{1/n_e} \frac{K_e}{c_s^2} \biggl[ \frac{c_s^6}{G^3 M_\mathrm{tot}^2} \biggr]^{1/n_e} \, .
</math>
</div>
(It is worth noting that if we set <math>n_e = -1</math>, the dimensionless parameter <math>\Pi_I</math> becomes identical to the parameter <math>\Pi</math> as defined [[SSCpt1/Virial#P-V_Diagram|in the context of our discussion of the Bonnor-Ebert sphere]].  But in order to complete the analogy with the Bonnor-Ebert sphere discussion, we would also need to change the sign  on the last term in the above expression for the free energy because in the earlier discussion the external pressure was an external, confining condition whereas here it is included as an internal energy of the system.)


<table border="1" cellpadding="5" align="center">
<tr>
  <th align="center" colspan="3">
Relevant Expressions for Isothermal Core
  </th>
</tr>

<tr>
  <td align="center">
<math>
\frac{\rho_e}{\rho_c}= \frac{\mu_e}{\mu_c}
</math>
  </td>
  <td align="center">
<math>
\frac{q^3}{\nu}\biggl( \frac{1 - \nu}{1- q^3} \biggr) 
</math>
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>

<tr>
  <td align="center">
<math>
\chi \equiv \frac{R}{R_0} 
</math>
  </td>
  <td align="center">
<math>
\Pi_I^{n_e/3}  [\nu^{-n_e}\ (1-\nu)^{n_e+1} ]^{1/3} q^{n_e} (1-q^3)^{-(n_e+1)/3} 
</math>
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>

<tr>
  <td align="center">
<math>
\chi^3 = \biggl(\frac{R}{R_0} \biggr)^{3}
</math>
  </td>
  <td align="center">
<math>
\Pi_I^{n_e}  \biggl( \frac{\rho_e}{\rho_c}\biggr)^{n_e + 1} \biggl[ q^3 + \biggl( \frac{\rho_e}{\rho_c} \biggr)(1-q^3)  \biggr]^{-1} 
</math>
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>

<tr>
  <td align="center">
<math>
\frac{A}{E_0} 
</math>
  </td>
  <td align="center">
<math>
\nu^2 \xi_s \biggl\{ 1 + \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1)
+ \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (\xi_s^5 - 1) - \frac{5}{2}(\xi_s^2-1) \biggr] \biggr\}
</math>
  </td>
  <td align="left">
<math>= \biggl( \frac{\nu^2}{q} \biggr) \biggl\{ f_A(\nu,q) \biggr\}</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>
\frac{B_e}{E_0} 
</math>
  </td>
  <td align="center">
<math>
\Pi_I \biggl( \frac{n_e}{3} \biggr) (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e}  (\xi_s^3 - 1)^{-1/n_e}
</math>
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>

<tr>
  <td align="center">
<math>
\frac{B_I}{E_0} 
</math>
  </td>
  <td align="center">
<math>\nu </math>
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>

<tr>
  <td align="center">
<math>
\frac{\mathfrak{G}}{E_0} 
</math>
  </td>
  <td align="center">
<math>
- \frac{A}{E_0} \chi^{-1}  - \frac{B_I}{E_0} \ln\chi + \frac{B_e}{E_0} \chi^{-3/n_e} 
</math>
  </td>
  <td align="left">
<math>
= \nu \biggl[ - \biggl( \frac{\nu}{q} \biggr) \frac{f_A(\nu,q)}{\chi} - \ln\chi + \frac{n_e}{3} \biggl(\frac{1}{q^3}-1\biggr)\biggr]
</math>
  </td>
</tr>
</table>

[On <font color="red">8 November 2013</font>, J. E. Tohline wrote: I just confirmed that the simpler expression for the normalized total free energy, <math>\mathfrak{G}/E_0</math>, matches the more complicated version.  I don't like the result because the third term in the free energy -- the one contributed by the internal energy of the envelope -- is independent of the radius of the configuration, <math>\chi</math>; it works out this way because my expression for <math>\Pi_I</math> has a radial dependence that exactly cancels out the explicit radial dependence that appears in the more complicated expression.  But maybe it's okay after all because this expression is intended to show how the free energy varies across the <math>(q,\nu)</math> plane, and the effect of <math>\Pi_I</math> appears implicitly through the specification of <math>\chi</math>, or visa versa.]


==Summary Expressions (Old)==
In the above derivations, we have adopted the notation, 
<div align="center">
<math>
\rho_\mathrm{norm} \equiv \frac{3M_\mathrm{tot}}{4\pi R_0^3} \, .
</math>
</div>
Now, guided by the dimensional aspects of the various coefficients in the free energy expression, we choose the following normalization energy and radius:
<div align="center">
<math>
E_0 = M_\mathrm{tot} K_e \rho_\mathrm{norm}^{1/n_e} 
</math>
&nbsp; &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; &nbsp;
<math>
R_0 = \frac{GM_\mathrm{tot}}{K_e \rho_\mathrm{norm}^{1/n_e}} \, .
</math>
</div>
When combined with the expression for <math>\rho_\mathrm{norm}</math>, these become,
<div align="center">
<math>
E_0 = \biggl[ \biggl( \frac{4\pi}{3} \biggr) G^3 M_\mathrm{tot}^{5-n_e} K_e^{-n_e}\biggr]^{1/(3-n_e)} 
</math>
&nbsp; &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; &nbsp;
<math>
R_0 = \biggl[ \frac{3M_\mathrm{tot}}{4\pi} \biggl( \frac{K_e}{GM_\mathrm{tot}} \biggr)^{n_e} \biggr]^{1/(3-n_e)}  \, .
</math>
</div>
So, the primary scales are determined after specifying two parameters:  <math>M_\mathrm{tot}</math> and <math>K_e</math>.
We also obtain,
<div align="center">
<math>
\kappa_I \equiv \frac{c_s^2}{K_e \rho_\mathrm{norm}^{1/n_e}}  = \frac{M_\mathrm{tot} c_s^2}{E_0} = 
c_s^2 \biggl[ \biggl( \frac{3}{4\pi} \biggr) \frac{K_e^{n_e}}{G^3 M_\mathrm{tot}^2} \biggr]^{1/(3-n_e)}  \, .
</math>
</div>


<table border="1" cellpadding="5" align="center">
<tr>
  <th align="center" colspan="2">
Relevant Expressions for Isothermal Core
  </th>
</tr>

<tr>
  <td align="center">
<math>
\frac{\rho_e}{\rho_c}
</math>
  </td>
  <td align="center">
<math>
\frac{q^3}{\nu}\biggl( \frac{1 - \nu}{1- q^3} \biggr) 
</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>
\chi \equiv \frac{R}{R_0} 
</math>
  </td>
  <td align="center">
<math>
q^{n_e} (1-q^3)^{-(n_e+1)/3} \biggl[\nu^{-n_e}\ (1-\nu)^{n_e+1} \kappa_I^{-n_e} \biggr]^{1/3} 
</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>
\frac{A}{E_0} 
</math>
  </td>
  <td align="center">
<math>
\biggl( \frac{3}{5} \biggr)\nu^2 \xi_s \biggl\{ 1 
+ \frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) (\xi_s^2 - 1)
+ \biggl( \frac{\rho_e}{\rho_c} \biggr)^2 \biggl[ (\xi_s^5 - 1) - \frac{5}{2}(\xi_s^2-1) \biggr] \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>
\frac{B_e}{E_0} 
</math>
  </td>
  <td align="center">
<math>
n_e(1 - \nu)^{1+1/n_e} \xi_s^{3/n_e}  (\xi_s^3 - 1)^{-1/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>
\frac{B_I}{E_0} 
</math>
  </td>
  <td align="center">
<math>
3 \kappa_I \nu 
</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>
\frac{\mathfrak{G}}{E_0} 
</math>
  </td>
  <td align="center">
<math>
- \frac{A}{E_0} \chi^{-1}  - \frac{B_I}{E_0} \ln\chi + \frac{B_e}{E_0} \chi^{-3/n_e} 
</math>
  </td>
</tr>
</table>

Subsequently, we will also find it useful to have expressions for the following coefficient ratios:

<table border="0" align="center">

<tr>
  <td align="right">
<math>\frac{n_e A}{3 B_e}</math>  
</td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{GM^2_\mathrm{tot}}{5R_0} \biggr) \nu^2 \xi_s \biggl\{ f(\rho_e/\rho_c , \xi_s) \biggr\}
\biggl[K_e M_\mathrm{tot}  \rho_\mathrm{norm}^{1/n_e} (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e}  (\xi_s^3 - 1)^{-1/n_e}\biggr]^{-1}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{GM_\mathrm{tot}}{5K_e \rho_\mathrm{norm}^{1/n_e} R_0} \biggr) \nu^2 \xi_s \biggl\{ f(\rho_e/\rho_c , \xi_s) \biggr\}
(1 - \nu)^{-(1+1/n_e)} q^{3/n_e}  \biggl(\frac{1}{q^3} - 1 \biggr)^{1/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{GM_\mathrm{tot}}{5K_e \rho_\mathrm{norm}^{1/n_e} R_0} \biggr) 
\nu^2 q^{-1} [ (1 - \nu)^{-(n_e+1)} (1 - q^3 ) ]^{1/n_e} \biggl\{ f(\rho_e/\rho_c , \xi_s) \biggr\} \, ;
</math>
  </td>
</tr>
</table>


<table border="0" align="center">
<tr>
  <td align="right">
<math>\frac{n_e B_I}{3 B_e}</math>  
</td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
M_\mathrm{tot} c_s^2 \nu \biggl[K_e M_\mathrm{tot}  \rho_\mathrm{norm}^{1/n_e} (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e}  (\xi_s^3 - 1)^{-1/n_e}\biggr]^{-1}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
 </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\kappa_I  \nu (1 - \nu)^{-(1+1/n_e)} q^{3/n_e}  \biggl(\frac{1}{q^3} - 1 \biggr)^{1/n_e} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
 </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\kappa_I  \nu \biggl[ (1 - \nu)^{-(n_e+1)} (1-q^{3})  \biggr]^{1/n_e} \, ;
</math>
  </td>
</tr>
</table>


<table border="0" align="center">

<tr>
  <td align="right">
<math>\frac{n_e B_c}{n_c B_e}</math>  
</td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\biggl[K_c M_\mathrm{tot}  \rho_\mathrm{norm}^{1/n_c} \nu^{1+1/n_c}\xi_s^{3/n_c}  \biggr]
\biggl[K_e M_\mathrm{tot}  \rho_\mathrm{norm}^{1/n_e} (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e}  (\xi_s^3 - 1)^{-1/n_e}\biggr]^{-1}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\kappa_A \nu^{1+1/n_c}\xi_s^{3/n_c}  \biggl[ (1 - \nu)^{1+1/n_e} \xi_s^{3/n_e}  (\xi_s^3 - 1)^{-1/n_e}\biggr]^{-1}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\kappa_A \biggl[\nu^{1+1/n_c} (1 - \nu)^{-(1+1/n_e)}  \biggr]  \biggl[ \xi_s^{3(1/n_c-1/n_e)}   (\xi_s^3 - 1)^{1/n_e}\biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\kappa_A \biggl[\nu^{1+1/n_c} (1 - \nu)^{-(1+1/n_e)}  \biggr]  \biggl[ q^{3(1/n_e-1/n_c)}   \biggl( \frac{1}{q^3} - 1\biggr)^{1/n_e}\biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\kappa_A \biggl[\nu^{1+1/n_c} (1 - \nu)^{-(1+1/n_e)}  \biggr]  \biggl[ q^{-3n_e/n_c}   (1-q^3) \biggr]^{1/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\kappa_A \biggl[\nu^{1+1/n_c} q^{-3/n_c}   \biggr]  \biggl[  (1 - \nu)^{-(n_e+1)} (1-q^3) \biggr]^{1/n_e} \, .
</math>
  </td>
</tr>
</table>

==Derivatives of Free Energy==
<div align="center">
<math>
\frac{\partial\mathfrak{G}}{\partial \chi} = A \chi^{-2} -(1-\delta_{\infty n_c}) \frac{3}{n_c} B_c \chi^{-(1+3/n_c)} 
- \delta_{\infty n_c} B_I \chi^{-1} -\frac{3}{n_e} B_e \chi^{-(1+3/n_e)} \, ;
</math>

<math>
\frac{\partial^2\mathfrak{G}}{\partial \chi^2} = -2 A \chi^{-3} + (1-\delta_{\infty n_c}) \frac{3}{n_c} \biggl(1+\frac{3}{n_c}\biggr) B_c \chi^{-(2+3/n_c)} 
+ \delta_{\infty n_c} B_I \chi^{-2} + \frac{3}{n_e} \biggl(1+\frac{3}{n_e}\biggr) B_e \chi^{-(2+3/n_e)} \, .
</math>
</div>

==Equilibrium Condition==
We obtain the equilibrium radius, <math>\chi_E</math>, when <math>\partial\mathfrak{G}/\partial\chi = 0</math>.  Hence,
the relation governing the equilibrium radius is,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
A \chi_E^{-2} 
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
(1-\delta_{\infty n_c}) \frac{3}{n_c} B_c \chi_E^{-1- 3/n_c)} 
+\delta_{\infty n_c} B_I \chi_E^{-1} +\frac{3}{n_e} B_e \chi_E^{-1-3/n_e)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~~ \frac{n_e A}{3B_e}
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
(1-\delta_{\infty n_c}) \frac{n_e B_c}{n_c B_e} \chi_E^{1- 3/n_c} 
+\delta_{\infty n_c} \frac{n_e B_I}{3B_e} \chi_E + \chi_E^{1-3/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~~ \chi_E^{1-3/n_e}
</math>
  </td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\alpha - (1-\delta_{\infty n_c}) \beta \chi_E^{1- 3/n_c} - \delta_{\infty n_c} \beta_I \chi_E \, ,
</math>
  </td>
</tr>

</table>
</div>
where,
<div align="center">
<math>\alpha \equiv \frac{n_e A}{3B_e} \, ; ~~~ \beta \equiv \frac{n_e B_c}{n_c B_e} \, ; ~~~ \beta_I \equiv \frac{n_e B_I}{3B_e} \, .</math> 
</div>

==Stability==
At this equilibrium radius, the second derivative of the free energy has the value,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\chi_E^3 \biggl( \frac{n_e}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E
</math>
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
-2 A\biggl( \frac{n_e}{3B_e} \biggr)  + (1-\delta_{\infty n_c}) \biggl( \frac{n_e}{3B_e} \biggr) \frac{3}{n_c} \biggl(1+\frac{3}{n_c}\biggr) B_c \chi_E^{1-3/n_c} 
+ \delta_{\infty n_c} \biggl( \frac{n_e}{3B_e} \biggr) B_I \chi_E + \frac{3}{n_e}\biggl( \frac{n_e}{3B_e} \biggr)  \biggl(1+\frac{3}{n_e}\biggr) B_e \chi_E^{1-3/n_e}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
-2 \alpha + (1-\delta_{\infty n_c}) \beta \biggl(1+\frac{3}{n_c}\biggr) \chi_E^{1-3/n_c} 
+ \delta_{\infty n_c} \beta_I \chi_E + \biggl(1+\frac{3}{n_e}\biggr) \chi_E^{1-3/n_e} \, ,
</math>
  </td>
</tr>

</table>
</div>
which, when combined with the condition for equilibrium gives,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\chi_E^3 \biggl( \frac{n_e}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E
</math>
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
-2 \alpha + (1-\delta_{\infty n_c}) \beta \biggl(1+\frac{3}{n_c}\biggr) \chi_E^{1-3/n_c} 
+ \delta_{\infty n_c} \beta_I \chi_E + \biggl(1+\frac{3}{n_e}\biggr) \biggl[  \alpha - (1-\delta_{\infty n_c}) \beta \chi_E^{1- 3/n_c} - \delta_{\infty n_c} \beta_I \chi_E \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\alpha \biggl(\frac{3}{n_e}-1\biggr) + (1-\delta_{\infty n_c}) \beta \biggl(\frac{3}{n_c}-\frac{3}{n_e}\biggr) \chi_E^{1-3/n_c} 
-  \delta_{\infty n_c} \beta_I \biggl(\frac{3}{n_e}\biggr)\chi_E 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~~ \chi_E^3 \biggl( \frac{n_e^2}{3B_e} \biggr) \frac{\partial^2\mathfrak{G}}{\partial \chi^2} \biggr|_E
</math>
  </td>
  <td align="right">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\alpha (3-n_e) - (1-\delta_{\infty n_c}) 3\beta \biggl(1 - \frac{n_e}{n_c} \biggr) \chi_E^{1-3/n_c} -  \delta_{\infty n_c} 3\beta_I \chi_E \, .
</math>
  </td>
</tr>

</table>
</div>

Finally, the equilibrium configuration is stable as long as this second derivative is positive.  Hence, for a bipolytrope with an 
isothermal core (<math>\delta_{\infty n_c} = 1</math>), the configuration is stable as long as,
<div align="center">
<math>
\chi_E <  \frac{\alpha (3-n_e)}{3\beta_I} \, .
</math>
</div>
In the adiabatic case (<math>\delta_{\infty n_c} = 0</math>), the configuration is stable as long as,
<div align="center">
<math>
\chi_E^{1-3/n_c} <  \frac{\alpha n_c (3-n_e)}{3\beta (n_c-n_e)} \, .
</math>
</div>


==Examples==
===Isothermal Core with <math>n=3/2</math> Envelope===
When the core is isothermal and <math>n_e = 3/2</math>, the equilibrium condition is:
<div align="center">
<math>
\chi_E^{-1} =  \alpha - \beta_I \chi_E \, ,
</math>

<math>
\Rightarrow ~~~~ \beta_I \chi_E^2 -  \alpha \chi_E + 1 = 0 \, ,
</math>

<math>
\Rightarrow ~~~~ \chi_E = \frac{1}{2\beta_I} \biggl[ \alpha \pm \sqrt{\alpha^2 - 4\beta_I} \biggr] = \frac{\alpha}{2\beta_I} \biggl[ 1 \pm \sqrt{1 - \frac{4\beta_I}{\alpha^2}} \biggr] \, .
</math>
</div>

At the same time, the condition for stability is,
<div align="center">
<math>
\chi_E <  \frac{\alpha}{2\beta_I} \, .
</math>
</div>

===Isothermal Core with <math>n=1</math> Envelope===
When the core is isothermal and <math>n_e = 1</math>, the equilibrium condition is:
<div align="center">
<math>
\chi_E^{-2} =  \alpha - \beta_I \chi_E \, ,
</math>

<math>
\Rightarrow ~~~~ \chi_E^3 -  \frac{\alpha}{\beta_I} \chi_E^2 + \frac{1}{\beta_I} = 0 \, .
</math>
</div>
(We need to solve this cubic equation.)

At the same time, the condition for stability is,
<div align="center">
<math>
\chi_E <  \frac{2\alpha}{3\beta_I} \, .
</math>
</div>

===Old (and probably incorrect) cases===
====Envelope with <math>n=3/2</math>====
If we choose an <math>n_e = 3/2</math> envelope, we obtain stability for,
<div align="center">
<math>
\chi_E <  \frac{2\beta}{\alpha}\, .
</math>
</div>
In this case, the equilibrium radius condition is,
<div align="center">
<math>
\chi_E^2 - \alpha \chi_E + \beta =0 
</math>

<math>
\Rightarrow ~~~~ \chi_E = \frac{1}{2}\biggl[\alpha \pm \biggl(  \alpha^2 -4\beta \biggr)^{1/2}   \biggr] = 
\frac{\alpha}{2}\biggl[1 \pm \biggl(  1 -\frac{4\beta}{\alpha^2} \biggr)^{1/2}   \biggr]
</math>

</div>


====Envelope with <math>n=1</math>====
If, instead, we choose an <math>n_e = 1</math> envelope, we obtain stability for,
<div align="center">
<math>
\chi_E <  \sqrt{\frac{3\beta}{\alpha} }\, .
</math>
</div>
In this case, the equilibrium radius condition is,
<div align="center">
<math>
\alpha = \chi_E + \beta \chi_E^{-2} \, ,
</math>

<math>
\Rightarrow ~~~~ \chi_E^3 - \alpha \chi_E^2 + \beta = 0 \, .
</math>

</div>

=See Also=


{{ SGFfooter }}