View source for SSC/Structure/BiPolytropes/Analytic51

__FORCETOC__
=BiPolytrope with n<sub>c</sub> = 5 and n<sub>e</sub> = 1=
{| class="PGEclass" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
|- 
! style="height: 125px; width: 125px; background-color:white;" |<font size="-1">[[H_BookTiledMenu#MoreModels|<b>Eggleton, Faulkner<br />&amp; Cannon (1998)<br /><br />Analytic</b>]]<br />(n<sub>c</sub>, n<sub>e</sub>) = (5, 1)</font>
|}
[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 30 March 2013:  As far as I have been able to determine, this analytic structural model has not previously been published in a refereed, archival journal.       Subsequent comment by J. E. Tohline on 23 June 2013:  Last night I stumbled upon an article by Eagleton, Faulkner, and Cannon (1998) in which this identical analytically definable bipolytrope has been presented.  Insight drawn from this article is presented in an additional subsection, below.]]
Here we construct a [[SSC/Structure/BiPolytropes#BiPolytropes|bipolytrope]] in which the core has an <math>~n_c=5</math> polytropic index and the envelope has an <math>~n_e=1</math> polytropic index.  This system is particularly interesting because the entire structure can be described by closed-form, analytic expressions.  In deriving the properties of this model, we will follow the [[SSC/Structure/BiPolytropes#Solution_Steps|general solution steps for constructing a bipolytrope]] that we have outlined elsewhere.   
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />

==Steps 2 &amp; 3==
Based on the discussion [[SSC/Structure/Polytropes#n_.3D_5_Polytrope|presented elsewhere of the structure of an isolated <math>n=5</math> polytrope]], the core of this bipolytrope will have the following properties:
<div align="center">
<math>
\theta(\xi) = \biggl[ 1 + \frac{1}{3}\xi^2 \biggr]^{-1/2} ~~~~\Rightarrow ~~~~  \theta_i = \biggl[ 1 + \frac{1}{3}\xi_i^2 \biggr]^{-1/2} ;
</math>

<math>
\frac{d\theta}{d\xi} = - \frac{\xi}{3}\biggl[ 1 + \frac{1}{3}\xi^2 \biggr]^{-3/2} ~~~~\Rightarrow ~~~~  \biggl(\frac{d\theta}{d\xi}\biggr)_i = - \frac{\xi_i}{3}\biggl[ 1 + \frac{1}{3}\xi_i^2 \biggr]^{-3/2} \, .
</math>
</div>
The first zero of the function <math>~\theta(\xi)</math> and, hence, the surface of the corresponding isolated <math>~n=5</math> polytrope is located at <math>~\xi_s = \infty</math>.  Hence, the interface between the core and the envelope can be positioned anywhere within the range, <math>~0 < \xi_i < \infty</math>.

==Step 4:  Throughout the core (0 &le; &xi; &le; &xi;<sub>i</sub>)==
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="center" colspan="3">
Specify:  <math>K_c</math> and <math>\rho_0 ~\Rightarrow</math>
  </td>
  <td colspan="2">
&nbsp;
  </td>
</tr>
<tr>
  <td align="right">
<math>\rho</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0 \theta^{n_c}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-5/2}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>P</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_c \rho_0^{1+1/n_c} \theta^{n_c + 1}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_c \rho_0^{6/5} \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>r</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{(n_c + 1)K_c}{4\pi G} \biggr]^{1/2} \rho_0^{(1-n_c)/(2n_c)} \xi</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c}{G\rho_0^{4/5}} \biggr]^{1/2} \biggl(\frac{3}{2\pi}\biggr)^{1/2} \xi</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>M_r</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>4\pi \biggl[ \frac{(n_c + 1)K_c}{4\pi G} \biggr]^{3/2} \rho_0^{(3-n_c)/(2n_c)} \biggl(-\xi^2 \frac{d\theta}{d\xi} \biggr)</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5} } \biggr]^{1/2} \biggl( \frac{2\cdot 3}{\pi } \biggr)^{1/2} \biggl[ \xi^3 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3/2} \biggr]</math>
  </td>
</tr>
</table>

</div>

<table border="1" align="center" cellpadding="5" width="80%"><tr><td align="left">

Equations (A2) from {{ EFC98 }} present the same relations but adopt the following notations:<font color="red" size="+2"><b>*</b></font>

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

<tr>
  <td align="right"><math>P</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>p_0 \theta^6</math></td>
<td align="center">&nbsp; &nbsp; &nbsp; &nbsp; <math>\Rightarrow~~~</math></td>
  <td align="right"><math>K_c</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>p_0 \rho_0^{-6/5} \, ;</math></td>
</tr>

<tr>
  <td align="right"><math>\rho</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\rho_0 \theta^5</math></td>
<td align="center">&nbsp; &nbsp; &nbsp; where: </td>
  <td align="right"><math>\theta</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>a_\mathrm{efc} \biggl[a_\mathrm{efc}^2 + r^2 \biggr]^{-1 / 2} 
= \biggl[1 + \biggl(\frac{r}{a_\mathrm{efc}}\biggr)^2 \biggr]^{-1 / 2}</math></td>
</tr>

<tr>
  <td align="right"><math>a_\mathrm{efc}^2 </math></td>
  <td align="center"><math>\equiv</math></td>
  <td align="left"><math>\frac{18p_0}{4\pi G \rho_0^2}</math></td>
<td align="center">&nbsp; &nbsp; &nbsp; &nbsp; <math>\Rightarrow~~~</math></td>
  <td align="right"><math>a_\mathrm{efc}^2</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\frac{18K_c \rho_0^{6/5}}{4\pi G \rho_0^2}
=
3 \biggl[ \frac{K_c }{G \rho_0^{4/5}} \cdot \biggl( \frac{3}{2\pi}\biggr) \biggr]</math></td>
</tr>

<tr>
  <td align="center" colspan="3">&nbsp;</td>
<td align="center">&nbsp; &nbsp; &nbsp; &nbsp; <math>\Rightarrow~~~</math></td>
  <td align="right"><math>a_\mathrm{efc}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>
\sqrt{3} \biggl[ \frac{K_c }{G \rho_0^{4/5}} \biggr]^{1 / 2} \biggl( \frac{3}{2\pi}\biggr)^{1 / 2}
=
\sqrt{3}\biggl(\frac{r}{\xi}\biggr)
</math></td>
</tr>

<tr>
  <td align="center" colspan="3">&nbsp;</td>
<td align="center">&nbsp; &nbsp; &nbsp; &nbsp; <math>\Rightarrow~~~</math></td>
  <td align="right"><math>\frac{r}{a_\mathrm{efc}}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\frac{\xi}{\sqrt{3}} \, ;</math></td>
</tr>
</table>
Hence,

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

<tr>
  <td align="right"><math>\theta</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\biggl[ 1 + \frac{\xi^2}{3}\biggr]^{-1 / 2} \, ,</math></td>
</tr>
</table>
which matches our expression for the core's polytrope function, <math>\theta</math>.

----

Now, look at the {{ EFC98hereafter }} expression for the core's integrated mass.
<table border="0" align="center" cellpadding="5">

<tr>
  <td align="right"><math>M_r</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\frac{4\pi \rho_0 a_\mathrm{efc}^3}{3} \cdot \frac{r^3}{(a_\mathrm{efc}^2 + r^2)^{3/ 2}}</math></td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\frac{4\pi \rho_0 }{3} \biggl[ a_\mathrm{efc}^3 \biggr] 
(r/a_\mathrm{efc})^3 \biggl[1 + (r/a_\mathrm{efc})^2 \biggr]^{-3/ 2}</math></td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\frac{4\pi \rho_0 }{3} \biggl[ \frac{K_c }{G \rho_0^{4/5}} \cdot \biggl( \frac{3^2}{2\pi}\biggr)  \biggr]^{3 / 2} 
\frac{\xi^3}{3^{3 / 2}} \biggl[1 + \frac{\xi^2}{3} \biggr]^{-3/ 2}</math></td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\frac{4\pi }{3} \biggl( \frac{3}{2\pi}\biggr)^{3 / 2} \biggl[ \frac{K_c^3}{G^3}\cdot \frac{\rho_0^2}{ \rho_0^{12/5}}  \biggr]^{1 / 2} 
\xi^3 \biggl[1 + \frac{\xi^2}{3} \biggr]^{-3/ 2}</math></td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\biggl[ \frac{K_c^3}{G^3\rho_0^{2/5}} \biggr]\biggl( \frac{2\cdot 3}{\pi} \biggr)^{1 / 2} 
\xi^3 \biggl(1 + \frac{\xi^2}{3} \biggr)^{-3/ 2}</math></td>
</tr>
</table>
This expression matches ours.
</td></tr></table>

==Step 5: Interface Conditions==

<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td colspan="3">
&nbsp;
  </td>
  <td align="left" colspan="2">
Setting <math>n_c=5</math>, <math>n_e=1</math>, and <math>\phi_i = 1 ~~~~\Rightarrow</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\frac{\rho_e}{\rho_0}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{n_c}_i \phi_i^{-n_e}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{5}_i </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\biggl( \frac{K_e}{K_c} \biggr) </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0^{1/n_c - 1/n_e}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-(1+1/n_e)} \theta^{1 - n_c/n_e}_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0^{-4/5}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-4}_i</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\frac{\eta_i}{\xi_i}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{n_c + 1}{n_e+1} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c}\biggr) \theta_i^{(n_c-1)/2} \phi_i^{(1-n_e)/2}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>3^{1/2} \biggl( \frac{\mu_e}{\mu_c}\biggr) \theta_i^{2}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\biggl( \frac{d\phi}{d\eta} \biggr)_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{n_c + 1}{n_e + 1} \biggr]^{1/2} \theta_i^{- (n_c + 1)/2} \phi_i^{(n_e+1)/2} \biggl( \frac{d\theta}{d\xi} \biggr)_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>3^{1/2} \theta_i^{- 3} \biggl( \frac{d\theta}{d\xi} \biggr)_i</math>
  </td>
</tr>
</table>
</div>

==Step 6:  Envelope Solution==

Adopting equation (8) of {{ Beech88full }}, the most general solution to the <math>n=1</math> Lane-Emden equation can be written in the form,
<div align="center">
<math>
\phi = A \biggl[ \frac{\sin(\eta - B)}{\eta} \biggr] \, ,
</math>
</div>
where <math>A</math> and <math>B</math> are constants.  The first derivative of this function is,
<div align="center">
<math>
\frac{d\phi}{d\eta} = \frac{A}{\eta^2} \biggl[ \eta\cos(\eta-B) - \sin(\eta-B) \biggr] \, .
</math>
</div>


<table border="1" align="center" cellpadding="5" width="80%"><tr><td align="left">

Equations (A2) from {{ EFC98 }} present the same relations but adopt the following notations:<font color="red" size="+2"><b>*</b></font>

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

<tr>
  <td align="right"><math>P</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>p_0 \theta^6</math></td>
<td align="center">&nbsp; &nbsp; &nbsp; &nbsp; <math>\Rightarrow~~~</math></td>
  <td align="right"><math>K_c</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>p_0 \rho_0^{-6/5} \, ;</math></td>
</tr>
</table>

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


From Step 5, above, we know the value of the function, <math>\phi</math> and its first derivative at the interface; specifically,
<div align="center">
<math>
\phi_i = 1~~~~\mathrm{and}
~~~~\biggl( \frac{d\phi}{d\eta}\biggr)_i =3^{1/2} \theta_i^{- 3} \biggl( \frac{d\theta}{d\xi} \biggr)_i~~~~
\mathrm{at}~~~~\eta_i =3^{1/2} \xi_i \biggl( \frac{\mu_e}{\mu_c}\biggr) \theta_i^{2} \, .</math>
</div>
From this information we can determine the constants <math>A</math> and <math>B</math>; specifically,
<div align="center">
<math>
\eta_i - B = \tan^{-1}(\Lambda_i^{-1}) = \frac{\pi}{2}- \tan^{-1}(\Lambda_i) \, ,
</math>

<math>
A = \frac{\phi_i \eta_i}{\sin(\eta_i - B)} = \phi_i \eta_i (1 + \Lambda_i^2)^{1/2} \, ,
</math>
</div>
where,
<div align="center">
<math>
\Lambda_i = \frac{1}{\eta_i} + \frac{1}{\phi_i} \biggl(\frac{d\phi}{d\eta}\biggr)_i \, .
</math>
</div>

==Step 7==
The surface will be defined by the location, <math>\eta_s</math>, at which the function <math>\phi(\eta)</math> first goes to zero, that is, 
<div align="center">
<math>
\eta_s = \pi + B =  \frac{\pi}{2} + \eta_i + \tan^{-1}(\Lambda_i) \, .
</math>
</div>

<table border="1" align="center" cellpadding="5" width="80%"><tr><td align="left">

Equations (A2) from {{ EFC98 }} present the same relations but adopt the following notations:<font color="red" size="+2"><b>*</b></font>

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

<tr>
  <td align="right"><math>P</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>p_0 \theta^6</math></td>
<td align="center">&nbsp; &nbsp; &nbsp; &nbsp; <math>\Rightarrow~~~</math></td>
  <td align="right"><math>K_c</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>p_0 \rho_0^{-6/5} \, ;</math></td>
</tr>
</table>

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

==Step 8:  Throughout the envelope (&eta;<sub>i</sub> &le; &eta; &le; &eta;<sub>s</sub>)==
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="center" colspan="5">
&nbsp;
  </td>
  <td align="left" colspan="2">
Knowing:  <math>K_e/K_c</math> and <math>\rho_e/\rho_0</math> from Step 5 &nbsp; <math>\Rightarrow</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\rho</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_e \phi^{n_e}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0 \biggl(\frac{\rho_e}{\rho_0}\biggr) \phi</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0 \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{5}_i \phi</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>P</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_e \rho_e^{1+1/n_e} \phi^{n_e + 1}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_c \rho_0^{6/5} \biggl(\frac{K_e \rho_0^{4/5}}{K_c}\biggr) \biggl(\frac{\rho_e}{\rho_0}\biggr)^{2} \phi^{2}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_c \rho_0^{6/5}  \theta^{6}_i \phi^{2}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>r</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{(n_e + 1)K_e}{4\pi G} \biggr]^{1/2} \rho_e^{(1-n_e)/(2n_e)} \eta</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c}{G \rho_0^{4/5}} \biggr]^{1/2} \biggl( \frac{K_e \rho_0^{4/5}}{K_c} \biggr)^{1/2} (2\pi)^{-1/2}\eta</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c}{G \rho_0^{4/5}} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \theta^{-2}_i (2\pi)^{-1/2}\eta</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>M_r</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>4\pi \biggl[ \frac{(n_e + 1)K_e}{4\pi G} \biggr]^{3/2} \rho_e^{(3-n_e)/(2n_e)} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5}} \biggr]^{1/2} \biggl( \frac{K_e \rho_0^{4/5}}{K_c} \biggr)^{3/2} \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl( \frac{2}{\pi} \biggr)^{1/2} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5}} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-1}_i \biggl( \frac{2}{\pi} \biggr)^{1/2} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>
</tr>
</table>

</div>

<table border="1" align="center" cellpadding="5" width="80%"><tr><td align="left">

An examination of their equations (A3) reveals that {{ EFC98hereafter }} continue to use <math>\theta</math> to represent the polytropic function throughout the envelope &#8212; for clarity, we will attach the subscript <math>\theta_\mathrm{efc} </math> &#8212; whereas we use <math>\phi</math>. Henceforth we will assume that these functions are interchangeable, that is, <math>\theta_\mathrm{efc} \leftrightarrow \phi</math>, and examine whether or not their various physical parameter expressions match ours. 

[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline:  As detailed in the text, there appears to be a type-setting error in both of these expressions; as published by EFC98, the exponent on the coefficient of theta_i should be 6 and 5, respectively, whereas it appears as 4.]]In their Eqs. (A3), {{ EFC98hereafter }} state that, throughout the envelope,

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

<tr>
  <td align="right"><math>P</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>p_0 \theta_c^4 \theta_\mathrm{efc}^2 = K_c \rho_0^{6 / 5} \theta_i^4 \phi^2</math></td>
<td align="center">&nbsp; &nbsp; &nbsp; &nbsp; and,</td>
  <td align="right"><math>\rho</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\frac{\rho_0}{\alpha} \cdot \theta_c^4 \theta_\mathrm{efc} 
= \rho_0\biggl(\frac{\mu_e}{\mu_c}\biggr)  \theta_i^4 \phi \, ,</math></td>
</tr>
</table>

where, in both expressions, we have replaced their label for the value of the polytropic function at the core-envelope interface, <math>\theta_c</math>, with our label, <math>\theta_i</math>. Both of their expressions match ours <font color="red"><b>EXCEPT &hellip; NOTE:</b></font> &nbsp; in both of their expressions, <math>\theta_i</math> is raised to the 4<sup>th</sup> power whereas, according to our derivation, this interface value should be raised to the 6<sup>th</sup> power in the expression for pressure and it should be raised to the 5<sup>th</sup> power in the expression for the density.  We are confident that our expressions are the correct ones and therefore presume that type-setting errors are present in both of the {{ EFC98hereafter }} expressions.

----

We state that the envelope's polytropic function has the form,

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

<tr>
  <td align="right"><math>\phi</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>-\frac{A}{\eta} \cdot \sin(B-\eta) \, ,</math></td>
</tr>
</table>
where,

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

<tr>
  <td align="right"><math>\eta</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>\biggl[ \frac{G \rho_0^{4/5}}{K_c} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{2}_i (2\pi)^{1 / 2}~r \, .</math>
  </td>
</tr>
</table>

{{ EFC98hereafter }} state that,

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

<tr>
  <td align="right"><math>\theta_\mathrm{efc}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left"><math>\frac{B_\mathrm{efc}}{r} \cdot \sin[\beta(r_s - r)] \, ,</math></td>
</tr>
</table>
where,

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

<tr>
  <td align="right"><math>\beta ~r</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>\frac{3\theta_c^2}{\alpha } \biggl[a_\mathrm{efc}\biggr]^{-1} r</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>3 \theta_c^2 \biggl(\frac{\mu_e}{\mu_c}\biggr) 
\biggl\{\sqrt{3} \biggl[ \frac{K_c }{G \rho_0^{4/5}} \biggr]^{1 / 2} \biggl( \frac{3}{2\pi}\biggr)^{1 / 2}\biggr\}^{-1} r
</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>\biggl[ \frac{G \rho_0^{4/5}}{K_c} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{2}_i (2\pi)^{1 / 2}~r \, .</math>
  </td>
</tr>
</table>
We therefore conclude that <math>\beta r = \eta</math>, and  <math>\beta r_s = B</math>.  If, as we assume to be the case, <math>\theta_\mathrm{efc} = \phi</math>, it must also be the case that,

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

<tr>
  <td align="right"><math>\frac{B_\mathrm{efc}}{r}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>-\frac{A}{\eta } </math>
  </td>
</tr>

<tr>
  <td align="right"><math>\Rightarrow ~~~ B_\mathrm{efc}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>-\frac{A}{\beta } \, .</math>
  </td>
</tr>
</table>

----

Our expression for the integrated mass throughout the envelope is,

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

<tr>
  <td align="right"><math>M_r</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>\biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5}} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-1}_i \biggl( \frac{2}{\pi} \biggr)^{1/2} \biggl\{- A \biggl[ \eta\cos(\eta-B) - \sin(\eta-B) \biggr]\biggr\}</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>-A \biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5}} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-1}_i \biggl( \frac{2}{\pi} \biggr)^{1/2} \biggl[\sin(B - \eta) + \eta\cos(B - \eta) \biggr] \, .</math>
  </td>
</tr>
</table>

According to {{ EFC98hereafter }},
<table border="0" align="center" cellpadding="5">

<tr>
  <td align="right"><math>M_r</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{4\pi \rho_0 \alpha}{9} \biggl[ B_\mathrm{efc} a_\mathrm{efc}^2 \biggr] 
\biggl\{
\sin[\beta(r_s - r)] + \beta r \cos[\beta(r_s - r)]
\biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{4\pi \rho_0 }{9}\biggl(\frac{\mu_e}{\mu_c}\biggr)^{-1}\biggl[ -\frac{A a_\mathrm{efc}^2}{\beta} \biggr] 
\biggl[
\sin(B-\eta) + \eta \cos(B-\eta)
\biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
-A \cdot
\frac{4\pi \rho_0 }{9}\biggl(\frac{\mu_e}{\mu_c}\biggr)^{-1} 
\biggl[
\sin(B-\eta) + \eta \cos(B-\eta)
\biggr]
\biggl\{ \biggl[ \frac{G \rho_0^{4/5}}{K_c} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{2}_i (2\pi)^{1 / 2}  \biggr\}^{-1}
\biggl\{ 3 \biggl[ \frac{K_c }{G \rho_0^{4/5}} \cdot \biggl( \frac{3}{2\pi}\biggr) \biggr] \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">&nbsp;</td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
-A \cdot
\biggl(\frac{2 }{\pi}\biggr)^{1 / 2}\biggl(\frac{\mu_e}{\mu_c}\biggr)^{-2} \theta^{-2}_i 
\biggl[
\sin(B-\eta) + \eta \cos(B-\eta)
\biggr]
\biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5}} \biggr]^{1/2}
</math>
  </td>
</tr>
</table>
</td></tr></table>

=Examples=

==Normalization==
The dimensionless variables used in Tables 1 &amp; 2 are defined as follows:
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>\rho^*</math>
  </td>
  <td align="center">
<math>\equiv</math>
  </td>
  <td align="left">
<math>\frac{\rho}{\rho_0}</math>
  </td>

  <td align="center">; &nbsp;&nbsp;&nbsp;</td>

  <td align="right">
<math>r^*</math>
  </td>
  <td align="center">
<math>\equiv</math>
  </td>
  <td align="left">
<math>\frac{r}{[K_c^{1/2}/(G^{1/2}\rho_0^{2/5})]}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>P^*</math>
  </td>
  <td align="center">
<math>\equiv</math>
  </td>
  <td align="left">
<math>\frac{P}{K_c\rho_0^{6/5}}</math>
  </td>

  <td align="center">; &nbsp;&nbsp;&nbsp;</td>

  <td align="right">
<math>M_r^*</math>
  </td>
  <td align="center">
<math>\equiv</math>
  </td>
  <td align="left">
<math>\frac{M_r}{[K_c^{3/2}/(G^{3/2}\rho_0^{1/5})]}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>H^*</math>
  </td>
  <td align="center">
<math>\equiv</math>
  </td>
  <td align="left">
<math>\frac{H}{K_c\rho_0^{1/5}}</math>
  </td>

  <td align="center">. &nbsp;&nbsp;&nbsp;</td>

  <td align="right" colspan="3">
&nbsp;
  </td>
</tr>

</table>
</div>

==Parameter Values==

The <math>2^\mathrm{nd}</math> column of Table 1 catalogues the analytic expressions that define various parameters and physical properties (as identified, respectively, in column 1) of the <math>(n_c, n_e) = (5, 1)</math> bipolytrope.  We have evaluated these expressions for various choices of the dimensionless interface radius, <math>\xi_i</math>, and have tabulated the results in the last few columns of the table.  The tabulated values have been derived assuming <math>\mu_e/\mu_c = 1</math>, that is, assuming that the core and the envelope have the same mean molecular weights.  

<!--  BEGIN TABLE OF PARAMETERS --->
<div align="center">
<b>Table 1:  Properties of <math>(n_c, n_3) = (5, 1)</math> BiPolytrope Having Various Interface Locations, <math>\xi_i</math></b><br>
[[File:BiPolytropeParametersV01.xml|Accompanying spreadsheet with parameter values]]

<table border="1" cellpadding="5" width="80%">
<tr>
  <td align="center">
Parameter
  </td>
  <td align="center">
<math>\xi_i</math>
  </td>
  <td align="center">
0.5
  </td>
  <td align="center">
1.0
  </td>
  <td align="center">
3.0
  </td>
  <td rowspan="20">
[[File:Bipolytrope51Boundaries02.png|500px|Examples]]<br />
<table border="0" align="left" cellpadding="10"><tr><td align="left">
For bipolytropic models having <math>\mu_e/\mu_c = 1.0</math>, this figure shows how the interface location, <math>\eta_i</math> (solid purple curve), the surface radius, <math>\eta_s</math> (green circular markers), and the parameter, <math>\tan^{-1}\Lambda_i</math> (orange circular markers), vary with <math>\xi_i</math> (ordinate) over the range, <math>0 \le \xi_i \le 12</math>.  The three horizontal, red-dashed line segments identify the values of <math>\xi_i</math> for which numerical values of these (and other) parameters have been listed in the table shown here on the left.</td></tr></table>
  </td>
</tr>

<tr>
  <td align="center">
<math>\theta_i</math>
  </td>
  <td align="center">
<math>\biggl( 1+\frac{1}{3}\xi_i^2 \biggr)^{-1/2}</math>
  </td>
  <td align="center">
0.96077
  </td>

  <td align="center">
0.86603
  </td>
  <td align="center">
0.50000
  </td>
</tr>

<tr>
  <td align="center">
<math>-\biggl(\frac{d\theta_i}{d\xi}\biggr)_i</math>
  </td>
  <td align="center">
<math>\frac{1}{3} \xi_i \biggl( 1+\frac{1}{3}\xi_i^2 \biggr)^{-3/2}</math>
  </td>
  <td align="center">
0.14781
  </td>

  <td align="center">
0.21651
  </td>
  <td align="center">
0.12500
  </td>
</tr>

<tr>
  <td align="center">
<math>r^*_\mathrm{core} \equiv r^*_i</math>
  </td>
  <td align="center">
<math>\biggl( \frac{3}{2\pi} \biggr)^{1/2} \xi_i</math>
  </td>
  <td align="center">
0.34549
  </td>

  <td align="center">
0.69099
  </td>
  <td align="center">
2.07297
  </td>
</tr>

<tr>
  <td align="center">
<math>\rho^*_i \biggr|_c = \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \rho^*_i \biggr|_e</math>
  </td>
  <td align="center">
<math>\biggl( 1+\frac{1}{3}\xi_i^2 \biggr)^{-5/2}</math>
  </td>
  <td align="center">
0.81864
  </td>

  <td align="center">
0.48714
  </td>
  <td align="center">
0.03125
  </td>
</tr>

<tr>
  <td align="center">
<math>P^*_i</math>
  </td>
  <td align="center">
<math>\biggl( 1+\frac{1}{3}\xi_i^2 \biggr)^{-3}</math>
  </td>
  <td align="center">
0.78653
  </td>

  <td align="center">
0.42188
  </td>
  <td align="center">
0.01563
  </td>
</tr>

<tr>
  <td align="center">
<math>H^*_i \biggr|_c = \frac{n_c+1}{n_e+1} \biggl( \frac{\mu_e}{\mu_c} \biggr) H^*_i \biggr|_e</math>
  </td>
  <td align="center">
<math>6 \biggl( 1+\frac{1}{3}\xi_i^2 \biggr)^{-1/2}</math>
  </td>
  <td align="center">
5.76461
  </td>

  <td align="center">
5.19615
  </td>
  <td align="center">
3.00000
  </td>
</tr>

<tr>
  <td align="center">
<math>M^*_\mathrm{core}</math>
  </td>
  <td align="center">
<math>\biggl( \frac{6}{\pi}\biggr)^{1/2} (\xi_i \theta_i)^3</math>
  </td>
  <td align="center">
0.15320
  </td>

  <td align="center">
0.89762
  </td>
  <td align="center">
4.66417
  </td>
</tr>

<tr>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1}\eta_i</math>
  </td>
  <td align="center">
<math>\sqrt{3} ~\theta_i^2 \xi_i</math>
  </td>
  <td align="center">
0.79941
  </td>

  <td align="center">
1.29904
  </td>
  <td align="center">
1.29904
  </td>
</tr>

<tr>
  <td align="center">
<math>-\biggl( \frac{d\phi}{d\eta} \biggr)_i</math>
  </td>
  <td align="center">
<math>\sqrt{3} ~\theta_i^{-3} \biggl( - \frac{d\theta}{d\xi} \biggr)_i = \frac{\xi_i}{\sqrt{3}}</math>
  </td>
  <td align="center">
0.28868
  </td>

  <td align="center">
0.57735
  </td>
  <td align="center">
1.73205
  </td>
</tr>

<tr>
  <td align="center">
<math>\Lambda_i</math>
  </td>
  <td align="center">
<math>\frac{1}{\eta_i} + \biggl( \frac{d\phi}{d\eta} \biggr)_i</math>
  </td>
  <td align="center">
0.96225
  </td>

  <td align="center">
0.19245
  </td>
  <td align="center">
-0.96225
  </td>
</tr>

<tr>
  <td align="center">
<math>A</math>
  </td>
  <td align="center">
<math>\eta_i (1 + \Lambda_i^2)^{1/2}</math>
  </td>
  <td align="center">
1.10940
  </td>

  <td align="center">
1.32288
  </td>
  <td align="center">
1.80278
  </td>
</tr>

<tr>
  <td align="center">
<math>B</math>
  </td>
  <td align="center">
<math>\eta_i - \frac{\pi}{2} + \tan^{-1}( \Lambda_i)</math>
  </td>
  <td align="center">
- 0.00523
  </td>

  <td align="center">
-0.08163
  </td>
  <td align="center">
-1.03792
  </td>
</tr>

<tr>
  <td align="center">
<math>\eta_s</math>
  </td>
  <td align="center">
<math>\pi + B</math>
  </td>
  <td align="center">
3.13637
  </td>

  <td align="center">
3.05996
  </td>
  <td align="center">
2.10367
  </td>
</tr>

<tr>
  <td align="center">
<math>- \biggl( \frac{d\phi}{d\eta} \biggr)_s</math>
  </td>
  <td align="center">
<math>\frac{A}{\eta_s}</math>
  </td>
  <td align="center">
0.35372
  </td>

  <td align="center">
0.43232
  </td>
  <td align="center">
0.85697
  </td>
</tr>

<tr>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \cdot \biggl[ R^* \equiv r^*_s \biggr]</math>
  </td>
  <td align="center">
<math>\frac{\eta_s}{\sqrt{2\pi} ~\theta_i^2}</math>
  </td>
  <td align="center">
1.35550
  </td>

  <td align="center">
1.62766
  </td>
  <td align="center">
3.35697
  </td>
</tr>

<tr>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^2 M^*_\mathrm{tot}</math>
  </td>
  <td align="center">
<math>\biggl(\frac{2}{\pi}\biggr)^{1/2} \theta_i^{-1} \biggl( -\eta^2 \frac{d\phi}{d\eta} \biggr)_s = \biggl(\frac{2}{\pi}\biggr)^{1/2} \frac{A\eta_s}{\theta_i}</math>
  </td>
  <td align="center">
2.88959
  </td>

  <td align="center">
3.72945
  </td>
  <td align="center">
6.05187
  </td>
</tr>

<tr>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \frac{\rho_c}{\bar\rho}</math>
  </td>
  <td align="center">
<math>~\frac{\eta_s^2}{3A\theta_i^5}</math>
  </td>
  <td align="center">
3.61035
  </td>

  <td align="center">
4.84326
  </td>
  <td align="center">
26.1844
  </td>
</tr>

<tr>
  <td align="center">
<math>\biggl(\frac{\mu_e}{\mu_c}\biggr)^{-2} \cdot \biggl[ \nu \equiv \frac{M_\mathrm{core}}{M_\mathrm{tot}} \biggr]</math>
  </td>
  <td align="center">
<math>\sqrt{3} ~\biggl( \frac{\xi_i^3 \theta_i^4}{A\eta_s} \biggr)</math>
  </td>
  <td align="center">
0.05302
  </td>

  <td align="center">
0.24068
  </td>
  <td align="center">
0.77070
  </td>
</tr>

<tr>
  <td align="center">
<math>\biggl(\frac{\mu_e}{\mu_c}\biggr)^{-1} \cdot \biggl[ q \equiv \frac{r_\mathrm{core}}{R} \biggr]</math>
  </td>
  <td align="center">
<math>\sqrt{3}~\biggl[\frac{\xi_i \theta_i^2}{\eta_s}\biggr]</math>
  </td>
  <td align="center">
0.25488
  </td>

  <td align="center">
0.42453
  </td>
  <td align="center">
0.61751
  </td>
</tr>

</table>
</div>
<!--  END TABLE OF PARAMETERS --->


Alternatively, if given <math>\mu_e/\mu_c</math> and the value of the parameter, <math>\eta_i</math>, then we have,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\biggl(\frac{\mu_e}{\mu_c} \biggr)^{-1}\eta_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{3^{3/2}\xi_i}{3 + \xi_i^2}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~ 0</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\xi_i^2 - \biggl[ \biggl(\frac{\mu_e}{\mu_c} \biggr) \frac{3^{3/2}}{\eta_i }\biggr] \xi_i + 3</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Rightarrow ~~~\xi_i </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\sqrt{3}
\biggl(\frac{\mu_e}{\mu_c} \biggr) \frac{3}{2\eta_i } \biggl\{ 1 \pm  
\sqrt{1 - \biggl[ \biggl(\frac{\mu_e}{\mu_c} \biggr)^{-1} \frac{2 \eta_i }{3}\biggr]^2  }
\biggr\} \, .
</math>
  </td>
</tr>
</table>

It must be understood, therefore, that the interface location is restricted to the range,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~\le \eta_i \le</math>
  </td>
  <td align="left">
<math>~\frac{3}{2}\biggl(\frac{\mu_e}{\mu_c} \biggr)\, ,</math>
  </td>
</tr>
</table>

and that this upper limit on <math>\eta_i</math> is associated with a model whose core radius is, <math>\xi_i = \sqrt{3}</math>.  Also,

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

<tr>
  <td align="right">
<math>\Lambda_i </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\frac{1}{\eta_i} -
\biggl(\frac{\mu_e}{\mu_c} \biggr) \frac{3}{2\eta_i } \biggl\{ 1 \pm  
\sqrt{1 - \biggl[ \biggl(\frac{\mu_e}{\mu_c} \biggr)^{-1} \frac{2 \eta_i }{3}\biggr]^2  }
\biggr\} \, .
</math>
  </td>
</tr>
</table>

==Profile==

Once the values of the key set of parameters have been determined as illustrated in Table 1, the radial profile of various physical variables can be determined throughout the bipolytrope as detailed in [[#Steps_2_.26_3|step &#35;4]] and [[#Step_7|step &#35;8]], above.  Table 2 summarizes the mathematical expressions that define the profile throughout the core (column 2) and throughout the envelope (column 3) of the normalized mass density, <math>\rho^*(r^*)</math>, the normalized gas pressure, <math>P^*(r^*)</math>, and the normalized mass interior to <math>r^*</math>, <math>M_r^*(r^*)</math>.  For all profiles, the relevant normalized radial coordinate is <math>r^*</math>, as defined in the <math>2^\mathrm{nd}</math> row of Table 2.  Graphical illustrations of these resulting profiles can be viewed by clicking on the thumbnail images posted in the last few columns of Table 2.  

<div align="center">
<b>Table 2:  Radial Profile of Various Physical Variables</b>
<table border="1" cellpadding="6">
<tr>
  <td align="center" rowspan="2">
Variable
  </td>
  <td align="center" rowspan="2">
Throughout the Core<br>
<math>0 \le \xi \le \xi_i</math>
  </td>
  <td align="center" rowspan="2">
Throughout the Envelope<sup>&dagger;</sup><br>
<math>\eta_i \le \eta \le \eta_s</math>
  </td>

  <td align="center" colspan="3">
Plotted Profiles
  </td>
</tr>

<tr>
  <td align="center">
<math>\xi_i = 0.5</math>
  </td>
  <td align="center">
<math>\xi_i = 1.0</math>
  </td>
  <td align="center">
<math>\xi_i = 3.0</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>r^*</math>
  </td>
  <td align="center">
<math>\biggl( \frac{3}{2\pi} \biggr)^{1/2} \xi</math>
  </td>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \theta^{-2}_i (2\pi)^{-1/2}\eta</math>
  </td>

  <td align="center" colspan="3">
&nbsp;
  </td>
</tr>

<tr>
  <td align="center">
<math>\rho^*</math>
  </td>
  <td align="center">
<math>\biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-5/2}</math>
  </td>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{5}_i \phi(\eta)</math>
  </td>

  <td align="center">
<!--  [[File:PlotDensity_xi_0.5.jpg|thumb|75px]] -->
[[Image:DenXi05.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:DenXi10.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:DenXi30.jpg|thumb|75px]]
  </td>
</tr>

<tr>
  <td align="center">
<math>P^*</math>
  </td>
  <td align="center">
<math>\biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3}</math>
  </td>
  <td align="center">
<math>\theta^{6}_i [\phi(\eta)]^{2}</math>
  </td>

  <td align="center">
<!-- [[File:PlotPressure_xi_0.5.jpg|thumb|75px]] -->
[[Image:PresXi05.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:PresXi10.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:PresXi30.jpg|thumb|75px]]
  </td>
</tr>

<tr>
  <td align="center">
<math>M_r^*</math>
  </td>
  <td align="center">
<math>\biggl( \frac{2\cdot 3}{\pi } \biggr)^{1/2} \biggl[ \xi^3 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3/2} \biggr]</math>
  </td>
  <td align="center">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-1}_i \biggl( \frac{2}{\pi} \biggr)^{1/2} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>

  <td align="center">
<!-- [[File:PlotPressure_xi_0.5.jpg|thumb|75px]] -->
[[Image:MassXi05.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:MassXi10.jpg|thumb|75px]]
  </td>
  <td align="center">
[[Image:MassXi30.jpg|thumb|75px]]
  </td>
</tr>

<tr>
  <td align="left" colspan="6">
<sup>&dagger;</sup>In order to obtain the various envelope profiles, it is necessary to evaluate <math>\phi(\eta)</math> and its first derivative using the information presented in Step 6, above.
  </td>
</tr>
</table>
</div>

[As of 28 April 2013] For the interface locations <math>\xi_i = 0.5, 1.0,\mathrm{and}~3.0</math>, Table 2 provides profiles for three values of the molecular weight ratio:  <math>\mu_e/\mu_c = 1.0, 1/2,\mathrm{and}~1/4</math>.  In all nine graphs, blue diamonds trace the structure of the <math>n_c=5</math> core; the core extends to a radius, <math>r^*_\mathrm{core}</math>, that is independent of molecular weight ratio but varies in direct proportion to the choice of <math>\xi_i</math>.  Specifically, as tabulated in the fourth row of Table 1, <math>r^*_\mathrm{core} = 0.34549, ~0.69099, \mathrm{and} ~2.07297</math> for, respectively, <math>\xi_i = 0.5,~1,~\mathrm{and}~3</math>.  Notice that, while the pressure profile and mass profile are continuous at the interface for all choices of the molecular weight ratio, the density profile exhibits a discontinuous jump that is in direct proportion to the chosen value of <math>\mu_e/\mu_c</math>.  

Throughout the <math>n_e = 1</math> envelope, the profile of all physical variables varies with the choice of the molecular weight ratio.  In the Table 2 graphs, red squares trace the envelope profile for <math>\mu_e/\mu_c = 1.0</math>; green triangles trace the envelope profile for <math>\mu_e/\mu_c = 1/2</math>; and purple crosses trace the envelope profile for <math>\mu_e/\mu_c = 1/4</math>.  The surface of the bipolytropic configuration is defined by the (normalized) radius, <math>R^*</math>, at which the envelope density and pressure drop to zero; the values tabulated in row 16 of Table 1 &#8212; <math>1.35550, ~1.62766, ~\mathrm{and} ~3.35697</math> for, respectively, <math>\xi_i = 0.5,~1,~\mathrm{and}~3</math> &#8212; correspond to a molecular weight ratio of unity and, hence also, to the envelope profiles traced by red squares in the Table 2 graphs.  As the molecular weight ratio is decreased from unity to <math>1/2</math> and, then, <math>1/4</math> for a given choice of <math>\xi_i</math>, the (normalized) radius of the bipolytrope increases roughly in inverse proportion to <math>\mu_e/\mu_c</math> as suggested by the formula for <math>R^*</math> shown in Table 1.  This proportional relation is not exact, however, because the parameter <math>\eta_s</math>, which also appears in the formula for <math>R^*</math>, contains an implicit dependence on the chosen value of the molecular weight ratio through the parameter <math>\eta_i</math>.

For a given choice of the interface parameter, <math>\xi_i</math>, the (normalized) mass that is contained in the core is independent of the choice of the molecular weight ratio.  However, the (normalized) total mass, <math>M_\mathrm{tot}^*</math>, varies significantly with the choice of <math>\mu_e/\mu_c</math>; as suggested by the expression provided in row 17 of Table 1, the variation is in rough proportion to <math>(\mu_e/\mu_c)^{-2}</math> but, as with <math>R^*</math>, this proportional relation is not exact because the parameters <math>\eta_s</math> and <math>A</math> which also appear in the formula for <math>M_\mathrm{tot}^*</math> harbor an implicit dependence on the molecular weight ratio.

==Model Sequences==


For a given choice of <math>\mu_e/\mu_c</math> a physically relevant sequence of models can be constructed by steadily increasing the value of <math>\xi_i</math> from zero to infinity &#8212; or at least to some value, <math>\xi_i \gg 1</math>.  Figure 1 shows how the fractional core mass, <math>\nu \equiv M_\mathrm{core}/M_\mathrm{tot}</math>, varies with the fractional core radius, <math>q \equiv r_\mathrm{core}/R</math>, along sequences having six different values of <math>\mu_e/\mu_c</math>, as detailed in the figure caption.  The natural expectation is that an increase in <math>\xi_i</math> along a given sequence will correspond to an increase in the relative size &#8212; both the radius and the mass &#8212; of the core. This expectation is realized along the sequences marked by blue diamonds (<math>\mu_e/\mu_c = 1</math>) and by red squares (<math>\mu_e/\mu_c = </math>&frac12;).  But the behavior is different along the other four illustrated sequences.  For sufficiently large <math>\xi_i</math>, the relative radius of the core begins to decrease; then, as <math>\xi_i</math> is pushed to even larger values, eventually the relative core mass begins to decrease.  Additional properties of these equilibrium sequences are discussed in [[SSC/FreeEnergy/PolytropesEmbedded#Behavior_of_Equilibrium_Sequence|an accompanying chapter]].

<div align="center">
<table border="0" cellpadding="5" width="85%">

<tr>
  <td align="center" colspan="2" bgcolor="white">
[[Image:PlotSequencesBest02.png|500px|center]]
  </td>
</tr>

<tr>
  <td align="left" colspan="2">
'''Figure 1:'''  Analytically determined plot of fractional core mass (<math>\nu</math>) versus fractional core radius (<math>q</math>) for <math>(n_c, n_e) = (5, 1)</math> bipolytrope model sequences having six different values of <math>\mu_e/\mu_c</math>: 1 (blue diamonds), &frac12; (red squares), 0.345 (dark purple crosses), &#8531; (pink triangles), 0.309 (light green dashes),  and &frac14; (purple asterisks).  Along each of the model sequences, points marked by solid-colored circles correspond to models whose interface parameter, <math>\xi_i</math>, has one of three values:  0.5 (green circles), 1 (dark blue circles), or 3 (orange circles); the images linked to Table 2 provide plots of the density, pressure and mass profiles for nine of these identified models.
  </td>
</tr>
</table>
</div>

The variation of <math>\nu</math> with <math>q</math> for a seventh analytically determined model sequence &#8212; one for which <math>\mu_e/\mu_c = 1/5</math> &#8212; is mapped out by a string of blue diamond symbols in the left-hand side of Figure 2. It behaves in an analogous fashion to the <math>\mu_e/\mu_c = </math>&frac14; (purple asterisks) sequence displayed in Figure 1.  It also quantitatively, as well as qualitatively, resembles the sequence that was numerically constructed by {{ SC42full }} for models with an isothermal core (<math>n_c = \infty</math>) and an <math>n_e=3/2</math> envelope; Fig. 1 from their paper has been reproduced here on the right-hand side of Figure 2.

<div align="center">
<table border="1" cellpadding="5" width="80%">
<tr>
  <td align="center" colspan="3">
<b>
Figure 2: Relationship to Sch&ouml;nberg-Chandrasekhar Mass Limit
</b>
  </td>
</tr>

<tr>
  <td align="center">
Analytic BiPolytrope with <math>n_c=5</math>, <math>n_e = 1</math>, and <math>\mu_e/\mu_c = 1/5</math>
  </td>
  <td align="center">
Edited excerpt from [http://adsabs.harvard.edu/abs/1942ApJ....96..161S Sch&ouml;nberg &amp; Chandrasekhar (1942)]
  </td>
  <td align="center">
Figure from [http://adsabs.harvard.edu/abs/1941ApJ....94..525H Henrich &amp; Chandraskhar (1941)]
  </td>
</tr>
<tr>
  <td align="center">
[[Image:SC_42Comparison.jpg|400px|center]]
  </td>
  <td align="center" rowspan="2">
[[Image:SC42_Fig1.jpg|400px|center]]
<!-- [[Image:AAAwaiting01.png|400px|center]] -->
  </td>
  <td align="center" rowspan="2">
[[Image:HenrichChandra41b.jpg|200px|center]]
<!-- [[Image:AAAwaiting01.png|200px|center]] -->
  </td>
</tr>

<tr>
  <td align="left">
(''Above'') Plot of fractional core mass (<math>\nu</math>) versus fractional core radius (<math>q</math>) for the analytic bipolytrope having <math>\mu_e/\mu_c = 1/5</math>.  The behavior of this analytically defined model sequence resembles the behavior of the numerically constructed isothermal core models presented by (''center'') {{ SC42 }} and by (''far right'') {{ HC41full }}.
  </td>
</tr>
</table>
</div>

=Limiting Mass=
==Background==
As early as 1941, Chandraskhar and his collaborators realized that the shape of the model sequence in a <math>\nu</math> versus <math>q</math> diagram, as displayed in Figure 2 above, implies that equilibrium structures can exist only if the fractional core mass lies below some limiting value.  This realization is documented, for example, by the following excerpt from &sect;5 of [http://adsabs.harvard.edu/abs/1941ApJ....94..525H Henrich &amp; Chandraskhar (1941)].

<div align="center">
<table border="1" cellpadding="5" width="60%">
<tr>
  <td align="center" colspan="1">
Text excerpt from &sect;5 (pp. 532 - 533) of<br />{{ HC41figure }}
  </td>
</tr>

<tr>
  <td align="left" colspan="1">
<!-- [[Image:HenrichChandra41a.jpg|600px|center|HenrichChandra1941]] -->
<!-- [[Image:AAAwaiting01.png|600px|center]] -->
<font color="darkgreen">"&hellip; at a fixed central temperature, the fraction of the total mass, <math>\nu</math>, contained in the core increases slowly at first and soon very rapidly as <math>q</math> approaches <math>q_\mathrm{max}</math>.  However, this increase of <math>\nu</math> does not continue indefinitely; <math>\nu</math> soon attains a maximum value <math>\nu_\mathrm{max}</math>.  There exists, therefore, an upper limit to the mass which can be contained in the isothermal core."</font>
  </td>
</tr>

</table>
</div>

Given that our bipolytropic sequence has been defined analytically, it may be possible to analytically determine the limiting core mass of our model.  In order to accomplish this, we need to identify the point along the sequence &#8212; in particular, the value of the dimensionless interface location &#8212; at which <math>d\nu/dq = 0</math> or, equivalently, <math>d\nu/d\xi_i = 0</math>.

Before carrying out the desired differentiations, we will find it useful to rewrite the relevant expressions in terms of the parameters,
<div align="center">
<math>
\ell_i \equiv \frac{\xi_i}{\sqrt{3}}  \, ;
</math>
&nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; 
<math>
m_3 \equiv 3 \biggl( \frac{\mu_e}{\mu_c} \biggr) \, .
</math>
</div>

We obtain,

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

<tr>
  <td align="right">
<math>\eta_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>m_3 \biggl( \frac{\ell_i}{1 + \ell_i^2} \biggr) \, ;</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\Lambda_i </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\frac{1}{m_3\ell_i} [ 1 + (1-m_3)\ell_i^2]  ~~~ \Rightarrow ~~~ 
</math>
&nbsp;&nbsp;'''<font color="red">Believe it or not &hellip; &nbsp;</font>'''
<math>
(1 + \Lambda^2) = \frac{(1+\ell_i^2)}{m_3^2 \ell_i^2} \biggl[ 1 + (1-m_3)^2 \ell_i^2 \biggr] \, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>A</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{1 + (1-m_3)^2 \ell_i^2}{1 + \ell_i^2} \biggr]^{1/2} \, ;</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\biggl( \frac{\pi}{2} \biggr)^{1/2} \frac{M_\mathrm{core}}{9}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{\ell_i^3}{(1 + \ell_i^2)^{3/2}} \, ;</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\biggl( \frac{\pi}{2} \biggr)^{1/2} \frac{M_\mathrm{tot}}{9}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\frac{1}{m_3^2} \biggl[ 1 + (1-m_3)^2 \ell_i^2 \biggr]^{1/2}
\biggl\{  \biggl(\frac{\pi}{2} + \tan^{-1} \Lambda_i \biggr) + m_3 \ell_i (1 + \ell_i^2)^{-1} \biggr\} \, .
</math>
  </td>
</tr>

</table>
</div>

Hence,
<div align="center">
<math>
\nu \equiv \frac{M_\mathrm{core}}{M_\mathrm{tot}} 
= (m_3^2 \ell_i^3) (1 + \ell_i^2)^{-1/2} [1 + (1-m_3)^2 \ell_i^2]^{-1/2}  \biggl[ m_3\ell_i + (1+\ell_i^2) \biggl(\frac{\pi}{2} + \tan^{-1} \Lambda_i \biggr) \biggr]^{-1}
</math>
</div>

<table border="1" align="center" width="80%" cellpadding="5"><tr><td align="left">
An interesting limiting case is <math>m_3 = 1</math>, in which case,
<table border="0" align="center" cellpadding="5">
<tr>
  <td align="right">
<math>
\nu
</math>
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
(\ell_i^3) (1 + \ell_i^2)^{-1/2}  \biggl[ \ell_i + (1+\ell_i^2) \biggl(\frac{\pi}{2} + \tan^{-1} \biggl(\frac{1}{\ell_i}\biggr) \biggr) \biggr]^{-1} \, ,</math>
  </td>
</tr>
</table>
and the maximum value of <math>\nu</math> along this sequence arises when <math>\ell_i \rightarrow \infty</math>, in which case,
<table border="0" align="center" cellpadding="5">
<tr>
  <td align="right">
<math>
\nu
</math>
  </td>
  <td align="center">
<math>
\rightarrow
</math>
  </td>
  <td align="left">
<math>
\ell_i^2   \biggl[ \ell_i + (1+\ell_i^2) \biggl(\frac{\pi}{2}  \biggr) \biggr]^{-1} 
\rightarrow
\frac{2}{\pi}
\, .</math>
  </td>
</tr>
</table>
</td></tr></table>


The condition, <math>d\nu/d\xi_i = 0</math>, also will be satisfied if the condition,
<div align="center">
<math>
\frac{d\ln\nu}{d\ln\ell_i} = 0 \, ,
</math>
</div>
is met.

==Derivation==
My manual derivation gives,

<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
(1+\ell_i^2)  \biggl[ \frac{\pi}{2} + \tan^{-1} \Lambda_i \biggr]
\biggl\{ 3  - \frac{(1-m_3)^2 \ell_i^2 (1+\ell_i^2)}{[ 1 + (1-m_3)^2 \ell_i^2 ]} \biggr\} 
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math> (1+\ell_i^2)
\frac{\partial \tan^{-1}\Lambda_i}{\partial \ln \ell_i} 
-m_3\ell_i \biggl\{ \ell_i^2 + 2  - \frac{(1-m_3)^2 \ell_i^2 (1+\ell_i^2)}{[ 1 + (1-m_3)^2 \ell_i^2 ]} \biggr\}
</math> 
  </td>
</tr>
</table>
</div>

where,

<div align="center">
<math>
\frac{\partial \tan^{-1}\Lambda_i}{\partial \ln \ell_i} =  
\frac{[(1-m_3)\ell_i^2 - 1 ] }{m_3\ell_i (1 + \Lambda_i^2)}
= \frac{m_3 \ell_i [(1-m_3)\ell_i^2 - 1 ]}{(1 + \ell_i^2) [ 1 + (1-m_3)^2 \ell_i^2]}  \, .
</math>
</div>

Upon rearrangement, this gives,

<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
(1+\ell_i^2)  \biggl[ \frac{\pi}{2} + \tan^{-1} \Lambda_i \biggr]
\biggl\{ 3[ 1 + (1-m_3)^2 \ell_i^2 ]  - (1-m_3)^2 \ell_i^2 (1+\ell_i^2) \biggr\} 
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
 m_3\ell_i \biggl\{[(1-m_3)\ell_i^2 - 1 ] -(\ell_i^2 + 2)[ 1 + (1-m_3)^2 \ell_i^2 ]  + (1-m_3)^2 \ell_i^2 (1+\ell_i^2) \biggr\} \, ,
</math> 
  </td>
</tr>
</table>
</div>


and further simplification <font color="red">[completed on 19 May 2013]</font> gives,

<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>
\underbrace{\biggl(\frac{\pi}{2} + \tan^{-1} \Lambda_i\biggr) (1+\ell_i^2) [ 3 + (1-m_3)^2(2-\ell_i^2)\ell_i^2]}_{\mathrm{LHS}} 
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\underbrace{m_3 \ell_i [(1-m_3)\ell_i^4 - (m_3^2 - m_3 +2)\ell_i^2 - 3]}_{\mathrm{RHS}} \, .
</math> 
  </td>
</tr>
</table>
</div>

<table border="1" align="center" cellpadding="8">
<tr>
  <td align="center" colspan="12">
<b>Maximum Fractional Core Mass, <math>\nu = M_\mathrm{core}/M_\mathrm{tot}</math> (solid green circular markers)<br />for Equilibrium Sequences having Various Values of <math>\mu_e/\mu_c</math>
  </td>
</tr>
<tr>
  <td align="center">
<math>\frac{\mu_e}{\mu_c}</math>
  </td>
  <td align="center">
<math>\xi_i</math>
  </td>
  <td align="center">
<math>\theta_i</math>
  </td>
  <td align="center">
<math>\eta_i</math>
  </td>
  <td align="center">
<math>\Lambda_i</math>
  </td>
  <td align="center">
<math>A</math>
  </td>
  <td align="center">
<math>\eta_s</math>
  </td>
  <td align="center">
LHS
  </td>
  <td align="center">
RHS
  </td>
  <td align="center">
<math>q \equiv \frac{r_\mathrm{core}}{R}</math>
  </td>
  <td align="center">
<math>\nu \equiv \frac{M_\mathrm{core}}{M_\mathrm{tot}}</math>
  </td>
  <td align="center" rowspan="7">[[File:TurningPoints51Bipolytropes.png|450px|Extrema along Various Equilibrium Sequences]]</td>
</tr>

<tr>
  <td align="center">
<math>\frac{1}{3}</math>
  </td>
  <td align="center">
<math>\infty</math> </td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">---</td>
  <td align="center">0.0 </td>
  <td align="center">
<math>\frac{2}{\pi}</math> </td>
</tr>

<tr>
  <td align="center">
0.33
  </td>
  <td align="right">
24.00496  </td>
  <td align="right">
0.0719668  </td>
  <td align="right">
0.0710624  </td>
  <td align="right">
0.2128753  </td>
  <td align="right">
0.0726547  </td>
  <td align="right">
1.8516032  </td>
  <td align="right">
-223.8157  </td>
  <td align="right">
-223.8159  </td>
  <td align="right">
0.038378833  </td>
  <td align="right">
0.52024552  </td>
</tr>

<tr>
  <td align="center">
0.316943
  </td>
  <td align="right">
10.744571  </td>
  <td align="right">
0.1591479  </td>
  <td align="right">
0.1493938  </td>
  <td align="right">
0.4903393  </td>
  <td align="right">
0.1663869  </td>
  <td align="right">
2.1760793  </td>
  <td align="right">
-31.55254  </td>
  <td align="right">
-31.55254  </td>
  <td align="right">
0.068652714  </td>
  <td align="right">
0.382383875  </td>
</tr>

<tr>
  <td align="center">
0.3090
  </td>
  <td align="right">
8.8301772  </td>
  <td align="right">
0.1924833  </td>
  <td align="right">
0.1750954  </td>
  <td align="right">
0.6130669  </td>
  <td align="right">
0.2053811  </td>
  <td align="right">
2.2958639  </td>
  <td align="right">
-18.47809  </td>
  <td align="right">
-18.47808  </td>
  <td align="right">
0.076265588  </td>
  <td align="right">
0.331475715  </td>
</tr>

<tr>
  <td align="center">
<math>\frac{1}{4}</math>
  </td>
  <td align="right">
4.9379256  </td>
  <td align="right">
0.3309933  </td>
  <td align="right">
0.2342522  </td>
  <td align="right">
1.4179907  </td>
  <td align="right">
0.4064595  </td>
  <td align="right">
2.761622  </td>
  <td align="right">
-2.601255  </td>
  <td align="right">
-2.601257  </td>
  <td align="right">
0.084824137  </td>
  <td align="right">
0.139370157  </td>
</tr>

<tr>
  <td align="left" colspan="11">
Recall that,
<div align="center">
<math>
\ell_i \equiv \frac{\xi_i}{\sqrt{3}}  \, ;
</math>
&nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; 
<math>
m_3 \equiv 3 \biggl( \frac{\mu_e}{\mu_c} \biggr) \, .
</math>
</div>
  </td>
</tr>
</table>

==Limit when m<sub>3</sub> = 0==
It is instructive to examine the root of this equation in the limit where <math>m_3 = 0</math> &#8212; that is, when <math>\mu_e/\mu_c = 0</math>.  First, we note that,
<div align="center">
<math>\Lambda_i\biggr|_{m_3 \rightarrow 0} = \biggl\{ \frac{1}{m_3\ell_i} [ 1 + (1-m_3)\ell_i^2] \biggr\}_{m_3 \rightarrow 0} = \infty \, .</math> 
</div>
Hence,
<div align="center">
<math>\biggl[\tan^{-1}\Lambda_i\biggr]_{m_3 \rightarrow 0} = \frac{\pi}{2} \, ,</math> 
</div>
and the limiting relation becomes,
<div align="center">
<math>
\pi (1+\ell_i^2) [ 3 + (2-\ell_i^2)\ell_i^2] = 0 \, ,
</math>
</div>
or, more simply,
<div align="center">
<math>
\ell_i^4 - 2\ell_i^2 - 3 = 0 \, .
</math>
</div>
The real root is,
<div align="center">
<math>\ell_i^2 = \frac{1}{2} \biggl[ 2 + \sqrt{4 + 12} \biggr] = 3 ~~~~ \Rightarrow ~~~~ \xi_i = 3 \, .</math>
</div>
For <math>\xi_i = 3</math>, the radius of the core, the mass of the core, and the pressure at the edge of the core are, respectively,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>r^*_\mathrm{core}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{3^3}{2\pi}\biggr)^{1/2} </math>
  </td>

<td align="center">&nbsp; &nbsp; <math>\Rightarrow</math>  &nbsp; &nbsp; </td>

  <td align="right">
<math>r_\mathrm{core}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{3^3}{2\pi}\biggr)^{1/2} \biggl[ \frac{K_c^{1/2}}{G^{1/2} \rho_0^{2/5}} \biggr] \, ;</math>
  </td>

</tr>

<tr>
  <td align="right">
<math>M^*_\mathrm{core}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{3^7}{2^5\pi}\biggr)^{1/2} </math>
  </td>

<td align="center">&nbsp; &nbsp; <math>\Rightarrow</math>  &nbsp; &nbsp; </td>

  <td align="right">
<math>M_\mathrm{core}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{3^7}{2^5\pi}\biggr)^{1/2}  \biggl[ \frac{K_c^{3/2}}{G^{3/2} \rho_0^{1/5}} \biggr] \, ;</math>
  </td>

</tr>

<tr>
  <td align="right">
<math>P^*_i</math>
  </td>
  <td align="right">
<math>=</math>
  </td>
  <td align="left">
<math>2^{-6} </math>
  </td>

<td align="center">&nbsp; &nbsp; <math>\Rightarrow</math>  &nbsp; &nbsp; </td>

  <td align="right">
<math>P_i</math>
  </td>
  <td align="right">
<math>=</math>
  </td>
  <td align="left">
<math>2^{-6} [ K_c \rho_0^{6/5}] \, .</math>
  </td>

</tr>
</table>
</div>

If we invert the middle expression to obtain <math>\rho_0</math> in terms of <math>M_\mathrm{core}</math>, specifically,
<div align="center">
<math>\rho_0^{1/5} = \biggl(\frac{3^7}{2^5\pi}\biggr)^{1/2}  \biggl[ \frac{K_c^{3/2}}{G^{3/2} M_\mathrm{core}} \biggr] \, ,</math> 
</div>
then we can rewrite <math>r_\mathrm{core}</math> and <math>P_i</math> in terms of, respectively, the ''reference'' radius,
<math>R_\mathrm{rf}</math>, and reference pressure, <math>P_\mathrm{rf}</math>, as defined in 
[[SSC/Structure/PolytropesEmbedded#Extension_to_Bounded_Sphere_2|our discussion of isolated <math>n=5</math>
polytropes embedded in an external medium]].  Specifically, we obtain,
<div align="center">
<table border="1" cellpadding="5">
<tr>

  <td align="right">
<math>r_\mathrm{core}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{3^3}{2\pi}\biggr)^{1/2} \biggl[ \frac{K_c^{1/2}}{G^{1/2} } \biggr] \biggl(\frac{3^7}{2^5\pi}\biggr)^{-1}  \biggl[ \frac{K_c^{3/2}}{G^{3/2} M_\mathrm{core}} \biggr]^{-2}</math>
  </td>

  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{2^9 \pi}{3^{11}}\biggr)^{1/2} \biggl[ \frac{G^{5/2} M^2_\mathrm{core}}{K_c^{5/2}} \biggr]</math>
  </td>

  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{2^9 \pi}{3^{11}}\biggr)^{1/2} \frac{3^3}{2^6} \biggl( \frac{5^5}{\pi} \biggr)^{1/2} R_\mathrm{rf} \biggr|_{n=5}</math>
  </td>

  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{5^5}{2^3 \cdot 3^5} \biggr)^{1/2} R_\mathrm{rf} \biggr|_{n=5}</math>
  </td>

</tr>
<tr>

  <td align="right">
<math>P_i</math>
  </td>
  <td align="right">
<math>=</math>
  </td>
  <td align="left">
<math>2^{-6} [ K_c ] \biggl(\frac{3^7}{2^5\pi}\biggr)^{3}  \biggl[ \frac{K_c^{3/2}}{G^{3/2} M_\mathrm{core}} \biggr]^6 </math>
  </td>

  <td align="right">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{3^7}{2^{7}\pi}\biggr)^{3}  \biggl[ \frac{K_c^{10}}{G^{9} M^6_\mathrm{core}} \biggr] </math>
  </td>

  <td align="right">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl(\frac{3^7}{2^7\pi}\biggr)^{3}  \biggl( \frac{2^{26} \pi^3}{3^{12} 5^9} \biggr) P_\mathrm{rf} \biggr|_{n=5}</math>
  </td>

  <td align="right">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{2^{5}\cdot 3^9 }{5^9} \biggr) P_\mathrm{rf} \biggr|_{n=5}</math>
  </td>

</tr>
</table>
</div>

['''<font color="red">26 May 2013</font>''' with further elaboration on '''<font color="red">28 May 2013</font>'''] This is the same result that was obtained when we
[[SSC/Structure/PolytropesEmbedded#Extension_to_Bounded_Sphere_2|embedded an isolated <math>n=5</math> polytrope in an external medium]].  Apparently, therefore, the physics that leads to the mass limit for a [[SSC/Structure/BonnorEbert#Pressure-Bounded_Isothermal_Sphere|Bonnor-Ebert sphere]] is the same physics that sets the {{ SC42 }} mass limit.

=Derivation by Eggleton, Faulkner, and Cannon (1998)=

The analytically prescribable sequence of bipolytropic models having <math>(n_c, n_e) = (5, 1)</math> displays an interesting behavior that extends beyond identification of a Sch&ouml;nberg-Chandrasekhar-like mass limit. After reaching a maximum value of <math>q</math> but before reaching the maximum value of <math>\nu</math>, the sequence bends back on itself.  This means that, even though the fraction of mass enclosed in the core is steadily increasing, the total radius of the configuration is increasing faster than the radius of the core.  Qualitatively, at least, this mimics the behavior exhibited by normal stars as they evolve off the main sequence and up the red giant branch.

As I pondered whether or not to probe this analogy in more depth, I recalled &#8212; even dating back to my years as a graduate student at Lick Observatory &#8212; hearing John Faulkner profess that he finally understood why stars become red giants.  I also recalled the following passage from  [<b>[[Appendix/References#HK94|<font color="red">HK94</font>]]</b>]'s textbook on ''Stellar Interiors'':
<table align="center" width="75%" border="1" cellpadding="10">
<tr>
  <td align="center">
Excerpt from &sect;2.3, p. 55 of [<b>[[Appendix/References#HK94|<font color="red">HK94</font>]]</b>]
  </td>
</tr>
<tr>
  <td align="left">
<font color="darkgreen">"The enormous increase in radius that accompanies hydrogen exhaustion moves the star into the red giant region of the H-R diagram.  While the transition to red giant dimensions is a fundamental result of all evolutionary calculations, a convincing yet intuitively satisfactory explanation of this dramatic transformation has not been formulated.  Our discussion of this phenomenon follows that of [http://adsabs.harvard.edu/abs/1984PhR...105..329I Iben and Renzini (1984)] although we must state that it is not the whole story."</font><sup>&dagger;</sup>
<br />_____________<br />
<sup>&dagger;</sup><font color="darkgreen">"Other attempts include:  {{ EF81 }}; [http://adsabs.harvard.edu/abs/1983A%26A...127..411W Weiss (1983)]; [http://adsabs.harvard.edu/abs/1985ApJ...296..554Y Yahil &amp; Van den Horn (1985)]; [http://adsabs.harvard.edu/abs/1988ApJ...329..803A Applegate (1988)]; [http://adsabs.harvard.edu/abs/1989MNRAS.236..505W Whitworth (1989)]; [http://adsabs.harvard.edu/abs/1992ApJ...400..280R Renzini et al. (1992)]. [http://adsabs.harvard.edu/abs/1991ApJ...372..592B Bhaskar &amp; Nigam (1991)] use an interesting set of dimensional arguments plus notions from polytrope theory.  We suspect the answers may lie in their paper but someone has yet to come along and translate the mathematics into an easily comprehensible physical picture."</font>
  </td>
</tr>
</table>

While examining the set of authors who more recently have cited the work by {{ EF81 }}, I discovered a paper by {{ EFC98full }} with the following abstract:

<!--
<table align="center" border="1" cellpadding="5" width="65%">
<tr>
  <td align="center" colspan="1">
[[Image:EagletonFaulknerCannon98.jpg|600px|center|Eggleton, Faulkner, &amp; Cannon (1998, MNRAS, 298, 831)]]
  </td>
</tr>
</table>
-->
<table border="1" cellpadding="10" align="center" width="75%">
<tr><td align="left">
<!-- [[Image:Tohline1985_Eq9.png|500px|center]] -->
<div align="center">{{ EFC98figure }}</div>
Abstract: &nbsp;<font color="darkgreen">"We present a simple analytic model of a composite polytropic star, which exhibits a limiting Sch&ouml;nberg-Chandrasekhar core mass fraction strongly analogous to the classic numerical result for an isothermal core, a radiative envelope and a &mu;-jump (i.e. a molecular weight jump) at the interface.  Our model consists of an n<sub>c</sub> = 5 core, an n<sub>e</sub> = 1 envelope and a &mu;-jump by a factor &ge; 3; the core mass fraction cannot exceed 2/&pi;.  We use the classic ''U, V'' plane to show that composite models will exhibit a Sch&ouml;nberg-Chandrasekhar limit only if the core is 'soft', i.e. has n<sub>c</sub> &ge; 5, and the envelope is 'hard', i.e. has n<sub>c</sub> < 5; in the critical case (n<sub>c</sub> = 5), the limit only exists if the &mu;-jump is sufficiently large, &ge; 6/(n<sub>e</sub> + 1)."</font>
</td></tr>
</table>
This paper uses analytic techniques to derive precisely the same sequence of <math>(n_c, n_e) = (5, 1)</math> bipolytropic models that we have presented above.

=Free Energy=
Here we use this bipolytrope's free energy function to probe the relative dynamical stability of various equilibrium models.  This derivation for <math>(n_c, n_e) = (5, 1)</math> bipolytropes is similar to the one that has been [[SSC/Structure/BiPolytropes/Analytic00#Free_Energy|presented elsewhere in the context of <math>(n_c, n_e) = (0, 0)</math> bipolytropes]] and follows the analysis outline provided in our [[SSC/BipolytropeGeneralization#Bipolytrope_Generalization|discussion of the stability of generalized bipolytropes]].

==Expression for Free Energy==
In order to construct the free energy function, we need mathematical expressions for the gravitational potential energy, <math>W</math>, and for the thermal energy content, <math>S</math>, of the models; and it will be natural to break both energy expressions into separate components derived for the <math>n_c=5</math> core and for the <math>n_e = 1</math> envelope.  Consistent with the above equilibrium model derivations, we will work with dimensionless variables.  Specifically, we define,

<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>~W^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>\frac{W}{[K_c^5/G^3]^{1/2}}</math>
  </td>

  <td align="center">; &nbsp;&nbsp;&nbsp;</td>

  <td align="right">
<math>~S^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>\frac{S}{[K_c^5/G^3]^{1/2}} \, .</math>
  </td>
</tr>
</table>
</div>

Drawing on the various functional expressions that are provided in the above derivations, including the [[SSC/Structure/BiPolytropes/Analytic51#Parameter_Values|Table of Parameters]], integrals over the material in the core give us,

<div align="center">
<table border="0" cellpadding="4">
<tr>
  <td align="right">
<math>~S^*_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ \frac{3}{2} \int_0^{r_i} \biggl(\frac{P^*}{\rho^*}\biggr)_\mathrm{core} (4\pi \rho^*)_\mathrm{core} (r^*)^2 dr^*</math>
  </td>

  <td align="center" rowspan="5" width="8%">
&nbsp;
  </td>
  <td align="center" rowspan="5">
[[File:Mathematica01.png|275px|center|Mathematica Integral]]
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ 6\pi \biggl( \frac{3}{2\pi} \biggr)^{3/2} \int_0^{\xi_i} \biggl(1+\frac{1}{3}\xi^2\biggr)^{-3} \xi^2 d\xi</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ 6\pi \biggl( \frac{3^2}{2\pi} \biggr)^{3/2} \int_0^{x_i} \biggl(1+x^2\biggr)^{-3} x^2 dx</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ \biggl( \frac{3^8}{2^7\pi} \biggr)^{1/2} \biggl[ \frac{x_i}{(1+x_i^2)} - \frac{2x_i}{(1+x_i^2)^2} + \tan^{-1}(x_i) \biggr] </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ \frac{1}{2} \biggl( \frac{3^8}{2^5\pi} \biggr)^{1/2} \biggl[ x_i (x_i^4 - 1 )(1+x_i^2)^{-3}   + \tan^{-1}(x_i) \biggr] \, ,</math>
  </td>
</tr>
</table>
</div>
where, in order to streamline the integral for Mathematica, we have used the substitution, <math>~x \equiv \xi/\sqrt{3}</math>; and,

<div align="center">
<table border="0" cellpadding="4">
<tr>
  <td align="right">
<math>~W^*_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ - \int_0^{r_i} (4\pi M_r^* \rho^*)_\mathrm{core} (r^*) dr^*</math>
  </td>

  <td align="center" rowspan="5" width="4%">
&nbsp;
  </td>
  <td align="center" rowspan="5">
[[File:Mathematica02.png|275px|center|Mathematica Integral]]
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ - 4\pi \int_0^{\xi_i} \biggl( \frac{2\cdot 3}{\pi } \biggr)^{1/2} \biggl[ \xi^3 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3/2} \biggr] 
\biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-5/2} \biggl( \frac{3}{2\pi} \biggr) \xi d\xi</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ - \biggl( \frac{2^3\cdot 3^8}{\pi } \biggr)^{1/2}\int_0^{x_i}  \biggl( 1 + x^2 \biggr)^{-4} x^4 dx</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ - \biggl( \frac{2^3\cdot 3^8}{\pi } \biggr)^{1/2}
\biggl[ 3\tan^{-1}(x_i) + \frac{x_i(3x_i^4 -8x_i^2 -3)}{(1+x_i^2)^3} \biggr] \biggl( \frac{1}{2^4 \cdot 3} \biggr) </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ - \biggl( \frac{3^8}{2^5\pi } \biggr)^{1/2}
\biggl[ x_i \biggl(x_i^4 - \frac{8}{3} x_i^2 -1 \biggr) (1 + x_i^2)^{-3}  + \tan^{-1}(x_i) \biggr] \, .</math>
  </td>
</tr>
</table>
</div>
(Apology:  The parameter <math>~x_i</math> introduced here is identical to the parameter <math>~\ell_i</math> that was introduced earlier in the context of [[SSC/Structure/BiPolytropes/Analytic51#Limiting_Mass|our discussion of the "Limiting Mass"]] of these models.  Sorry for the unnecessary duplication of parameters and possible confusion!)  

<div id="Buchdahl1978">
<table border="1" align="center" width="85%" cellpadding="8">
<tr><td align="left">
While our aim, here, has been to determine an expression for the gravitational potential energy of a ''truncated'' <math>n = 5</math> polytropic sphere, our derived expression can also give the gravitational potential of an ''isolated'' <math>n = 5</math> polytrope by evaluating the expression in the limit <math>x_i \rightarrow \infty</math>.  In this limit, the first term inside the square brackets goes to zero, while the second term,
<div align="center">
<math>\lim_{x_i \to \infty}\tan^{-1}(x_i) = \frac{\pi}{2} \, .</math>
</div>

We see, therefore, that,
<div align="center">
<math>
W^* \biggr|_\mathrm{tot} =
\lim_{x_i \to \infty}W^* = - \biggl( \frac{3^8}{2^5\pi } \biggr)^{1/2}\frac{\pi}{2} 
= - \biggl( \frac{3^8 \pi}{2^7} \biggr)^{1/2} \, .
</math>
</div>

Taking into account our adopted energy normalization, this can be rewritten with the dimensions of energy as,

<div align="center" id="twoSplusWcore">
<table border="0" cellpadding="4">
<tr>
  <td align="right">
<math>W_\mathrm{grav} \biggr|_\mathrm{tot}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math> 
- \biggl( \frac{3^8 \pi}{2^7} \biggr)^{1/2} \biggl( \frac{K_c^5}{G^3} \biggr)^{1/2} 
= - \biggl( \frac{3^8 \pi}{2^7} \biggr)^{1/2} \biggl[ \frac{\pi}{2^3\cdot 3^7} \biggr]^{1/2} \frac{GM^2 }{a_5} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math> 
- \biggl( \frac{3 \pi^2}{2^{10}} \biggr)^{1/2} \frac{GM^2}{a_5} \, ,
</math>
  </td>
</tr>
</table>
</div>

where we have elected to write the total gravitational potential energy in terms of the natural scale length for <math>~n = 5</math> polytropes, which, [[SSC/Structure/Polytropes#Primary_E-Type_Solution_2|as documented elsewhere]], is, 

<div align="center" id="twoSplusWcore">
<table border="0" cellpadding="4">
<tr>
  <td align="right">
<math>a_{5}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math> 
\biggl[ \frac{3K}{2\pi G} \biggr]^{1/2} \rho_c^{-2/5} 
= \biggl[ \frac{3K}{2\pi G} \biggr]^{1/2} \biggl[ \frac{\pi M^2 G^3}{2\cdot 3^4} \biggr] K^{-3}
= GM^2 \biggl[ \frac{\pi G^3}{2^3\cdot 3^7 K^5} \biggr]^{1/2} \, .
</math>
  </td>
</tr>
</table>
</div>

As can be seen from the following, boxed-in equation excerpt, our derived expression for the total gravitational potential energy of an ''isolated'' <math>n=5</math> polytrope exactly matches the result derived by {{ Buchdahl78full }}.  The primary purpose of Buchdahl's ''short communication'' was to point out that, despite the fact that its radius extends to infinity, "the gravitational potential energy of [an ''isolated''] polytrope of index 5 is finite."

<div align="center">
<table border="1" align="center" cellpadding="8" width="75%">
<tr><td align="center">
Equation excerpt from p. 116 of<br />{{ Buchdahl78figure }}
<!--[http://adsabs.harvard.edu/abs/1978AuJPh..31..115B H. A. Buchdahl (1978, Astralian J. Phys., 31, 115)]-->
</td></tr>
<tr><td align="center">
<!-- [[File:Buchdahl1978.png|350px|center|Buchdahl (1978, Australian J. Phys., 31, 115)]] -->
<!-- [[Image:AAAwaiting01.png|350px|center|Buchdahl (1978, Australian J. Phys., 31, 115)]] -->
<math>\Omega = -(\pi\sqrt{3}/32)GM^2/\alpha \, .</math>
</td></tr>
<tr><td align="left">Note that a comparison between Buchdahl's derived expression and our expression in the limit <math>x_i \rightarrow \infty</math> requires the parameter substitutions,
<div align="center"><math>\Omega \rightarrow W_\mathrm{grav}|_\mathrm{tot}</math> &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; <math>\alpha \rightarrow a_5</math></div>
</td></tr>
</table>
</div>

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

Notice that these two terms combine to give, for the core,

<div align="center" id="twoSplusWcore">
<table border="0" cellpadding="4">
<tr>
  <td align="right">
<math>\biggl( 2S + W \biggr)_\mathrm{core}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{2 \cdot 3^6}{\pi } \biggr)^{1/2} \frac{x_i^3}{(1 + x_i^2)^3} 
= \biggl( \frac{2}{\pi } \biggr)^{1/2} 3^{3/2} \xi_i^3 \biggl( 1 + \frac{1}{3}\xi_i^2 \biggr)^{-3}\, .</math>
  </td>
</tr>
</table>
</div>

Similarly, integrals over the material in the envelope give us,

<div align="center">
<table border="0" cellpadding="4">
<tr>
  <td align="right">
<math>S^*_\mathrm{env}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{3}{2} \int_{r_i}^R \biggl(\frac{P^*}{\rho^*}\biggr)_\mathrm{env} (4\pi \rho^*)_\mathrm{env} (r^*)^2 dr^*</math>
  </td>

  <td align="center" rowspan="5" width="8%">
&nbsp;
  </td>
  <td align="center" rowspan="5">
[[File:Mathematica03.png|300px|center|Mathematica Integral]]
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>6\pi \int_{\eta_i}^{\eta_s} [\theta^{6}_i \phi^{2}] \biggl[ \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \theta^{-2}_i (2\pi)^{-1/2} \biggr]^3 \eta^2 d\eta</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{3^2}{2\pi} \biggr)^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2 \int_{\eta_i}^{\eta_s} [\sin(\eta - B)]^2 d\eta</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{3^2}{2^5\pi} \biggr)^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2 
\biggl\{  2(\eta - B) - \sin[2(\eta-B)] \biggr\}_{\eta_i}^{\eta_s} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2^5\pi} \biggr)^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2 
\biggl\{  6(\eta - B) - 3\sin[2(\eta-B)] \biggr\}_{\eta_i}^{\eta_s} \, ; </math>
  </td>
</tr>
</table>
</div>

and,

<div align="center">
<table border="0" cellpadding="4">
<tr>
  <td align="right">
<math>W^*_\mathrm{env}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>- \int_{r_i}^{R} (4\pi M_r^* \rho^*)_\mathrm{env} (r^*) dr^*</math>
  </td>

  <td align="center" rowspan="6" width="4%">
&nbsp;
  </td>
  <td align="center" rowspan="6">
[[File:Mathematica04.png|300px|center|Mathematica Integral]]
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>- 4\pi \int_{\eta_i}^{\eta_s} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-1}_i \biggl( \frac{2}{\pi} \biggr)^{1/2} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)
\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{5}_i \phi
\biggl[ \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \theta^{-2}_i (2\pi)^{-1/2} \biggr]^2 \eta d\eta</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>-\biggl( \frac{2^3}{\pi} \biggr)^{1/2}  \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} 
\int_{\eta_i}^{\eta_s}  \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr) \phi  \eta d\eta
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>-\biggl( \frac{2^3}{\pi} \biggr)^{1/2}  \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2
\int_{\eta_i}^{\eta_s}  [ \sin(\eta-B) - \eta\cos(\eta-B) ] \sin(\eta - B)  d\eta
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>-\biggl( \frac{1}{2^3\pi} \biggr)^{1/2}  \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2
\biggl\{ - 3\sin[2(\eta - B)] +2\eta \cos[2(\eta - B)] + 4(\eta - B) + 2B \biggr\}_{\eta_i}^{\eta_s}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>-\biggl( \frac{1}{2^3\pi} \biggr)^{1/2}  \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2
\biggl\{6(\eta-B) - 3\sin[2(\eta - B)] -4\eta\sin^2(\eta-B) + 4B \biggr\}_{\eta_i}^{\eta_s} \, .
</math>
  </td>
</tr>
</table>
</div>

In this case, the two terms combine to give, for the envelope,

<div align="center">
<table border="0" cellpadding="4">

<tr>
  <td align="right">
<math>\biggl( 2S + W \biggr)_\mathrm{env}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{1}{2^3\pi} \biggr)^{1/2}  \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2
\biggl[4\eta\sin^2(\eta-B) + 4B \biggr]_{\eta_i}^{\eta_s}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{2}{\pi} \biggr)^{1/2}  \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2
\biggl[\eta_s\sin^2(\eta_s-B) - \eta_i\sin^2(\eta_i-B) \biggr] \, .</math>
  </td>
</tr>
</table>
</div>

==Equilibrium Condition==
===Global===

Recognizing from the above [[SSC/Structure/BiPolytropes/Analytic51#Parameter_Values|Table of Parameters]] that,
<div align="center">
<table border="0" cellpadding="3">

<tr>
  <td align="right">
<math>~A</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~\frac{\eta_i}{\sin(\eta_i - B)} \, ,</math>
  </td>
  <td align="left">
[because <math>~\phi_i = 1</math>]
  </td>
</tr>

<tr>
  <td align="right">
<math>~(\eta_s - B)</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~\pi \, ,</math>
  </td>
  <td align="left">
[hence, <math>~\sin^2(\eta_s - B) = 0</math>]
  </td>
</tr>

<tr>
  <td align="right">
<math>~\eta_i</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~3^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr) \xi_i \biggl( 1 + \frac{1}{3} \xi_i^2 \biggr)^{-1}\, ,</math>
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>
</table>
</div>
we can rewrite this last "envelope virial" expression as,

<div align="center">
<table border="0" cellpadding="4">

<tr>
  <td align="right">
<math>~\biggl( 2S + W \biggr)_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>- ~ \biggl( \frac{2}{\pi} \biggr)^{1/2}  \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} \eta_i^3</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>- ~ \biggl( \frac{2}{\pi} \biggr)^{1/2} 3^{3/2}  \xi_i^3 \biggl( 1 + \frac{1}{3} \xi_i^2 \biggr)^{-3} \, .</math>
  </td>
</tr>
</table>
</div>
This expression is equal in magnitude, but opposite in sign to the "core virial" expression derived earlier.   Hence, putting the core and envelope contributions together, we find,
<div align="center">
<table border="0" cellpadding="4">

<tr>
  <td align="right">
<math>~\biggl( 2S + W \biggr)_\mathrm{tot} ~=~ 2(S_\mathrm{core} + S_\mathrm{env}) + (W_\mathrm{core} + W_\mathrm{env})</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>~ 0 \, .</math>
  </td>
</tr>
</table>
</div>
This demonstrates that the detailed force-balanced models of <math>~(n_c, n_e) = (5,1)</math> bipolytropes derived above are also all in virial equilibrium, as should be the case.  More importantly, showing that these four separate energy integrals sum to zero helps provide confirmation that the four energy integrals have been derived correctly.  This allows us to confidently proceed to an evaluation of the relative dynamical stability of the models.

===In Parts===
In section <b><font color="maroon" size="+1">&#x2469;</font></b> of our ''[[SSC/SynopsisStyleSheet#Bipolytropes|Tabular Overview]]'', we speculated that, in bipolytropic equilibrium structures, the statements

<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~2S_\mathrm{core} + W_\mathrm{core} = 3P_i V_\mathrm{core}</math>
  </td>
  <td align="center">
&nbsp; &nbsp; and &nbsp; &nbsp;
  </td>
  <td align="left">
<math>~2S_\mathrm{env} + W_\mathrm{env} = - 3P_i V_\mathrm{core}  \, ,</math>
  </td>
</tr>
</table>
hold separately.  Let's evaluate the "PV" term.  We find that,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~
3P_i V_\mathrm{core} = 4\pi P_i r_i^3
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
4\pi \biggl( 1 + \frac{\xi_i^2}{3} \biggr)^{- 3}\biggl(\frac{3}{2\pi}\biggr)^{3 / 2} \xi_i^3 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( 1 + \frac{\xi_i^2}{3} \biggr)^{- 3}\biggl(\frac{2 \cdot 3^3 }{\pi} \biggr)^{1 / 2} \xi_i^3 \, .
</math>
  </td>
</tr>
</table>
This is '''precisely''' the "extra term" that shows up (with opposite signs) in the above-derived expressions for the separate quantities, <math>~(2S + W)_\mathrm{core}</math> and <math>~(2S + W)_\mathrm{env}</math>.  Hence our speculation has been shown to be correct, at least for the case of bipolytropes with, <math>~(\gamma_c, \gamma_e) = (\tfrac{6}{5}, 2)</math>.

==Stability Condition==

According to the accompanying [[SSC/BipolytropeGeneralization#Stability|free-energy based, generalized formulation of stability in bipolytropes]], our above derived <math>~(n_c, n_e) = (5,1)</math> bipolytropes &#8212; or, equivalently, <math>~(\gamma_c, \gamma_e) = (6/5, 2)</math> bipolytropes &#8212; will be dynamically stable only if,

<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>~ - (W_\mathrm{core} + W_\mathrm{env}) \biggl( \gamma_e - \frac{4}{3}\biggr) 
</math>
  </td>
  <td align="center">
<math>~>~</math>
  </td>
  <td align="left">
<math> ~ 2(\gamma_e-\gamma_c) S_\mathrm{core}  \, .</math>
  </td>
</tr>
</table>
</div>
Otherwise, they will be dynamically unstable toward radial perturbations.  For various values of the <math>~\mu_e/\mu_c</math> ratio, Table 3 identifies the value of <math>~\xi_i</math> &#8212; and the corresponding values of <math>~q</math> and <math>~\nu</math> &#8212; at which the left-hand side of this stability relation equals the right-hand side.  The locus of points provided by Table 3 defines the curve that separates stable from unstable regions of the <math>~q-\nu</math> parameter space. The red-dashed curve drawn in Figure 3 graphically depicts this demarcation:  the region below the curve identifies bipolytrope models that are dynamically stable while the region above the curve identifies unstable models.

<div align="center">
<table border="0" cellpadding="5" width="85%">

<tr>
  <td align="center" rowspan="2">

<table border="1" cellpadding="5">
<tr>
  <th align="center" colspan="9">
Table 3:  Points Defining Stability Curve
  </th>
</tr>

<tr>
  <td align="center" width="25%">
<math>~\mu_e/\mu_c</math>
  </td>
  <td align="center" width="25%">
<math>~\xi_i</math>
  </td>
  <td align="center" width="25%">
<math>~q </math>
  </td>
  <td align="center">
<math>~\nu </math>
  </td>

  <td align="center" rowspan="12" bgcolor="#EEEEEE">
&nbsp;
  </td>

  <td align="center" width="25%">
<math>~\mu_e/\mu_c</math>
  </td>
  <td align="center" width="25%">
<math>~\xi_i</math>
  </td>
  <td align="center" width="25%">
<math>~q </math>
  </td>
  <td align="center">
<math>~\nu </math>
  </td>
</tr>

<tr>
  <td align="center">
1
  </td>
  <td align="center">
2.416
  </td>
  <td align="center">
0.5952
  </td>
  <td align="center">
0.6830
  </td>

  <td align="center">
0.375
  </td>
  <td align="center">
6.259
  </td>
  <td align="center">
0.1695
  </td>
  <td align="center">
0.6054
  </td>
</tr>

<tr>
  <td align="center">
0.95
  </td>
  <td align="center">
2.500
  </td>
  <td align="center">
0.5805
  </td>
  <td align="center">
0.6884
  </td>

  <td align="center">
0.350
  </td>
  <td align="center">
7.341
  </td>
  <td align="center">
0.1284
  </td>
  <td align="center">
0.5439
  </td>
</tr>

<tr>
  <td align="center">
0.90
  </td>
  <td align="center">
2.594
  </td>
  <td align="center">
0.5642
  </td>
  <td align="center">
0.6937
  </td>

  <td align="center">
0.340
  </td>
  <td align="center">
7.991
  </td>
  <td align="center">
0.1109
  </td>
  <td align="center">
0.5081
  </td>
</tr>

<tr>
  <td align="center">
0.80
  </td>
  <td align="center">
2.816
  </td>
  <td align="center">
0.5255
  </td>
  <td align="center">
0.7031
  </td>

  <td align="center">
&#8531;
  </td>
  <td align="center">
8.548
  </td>
  <td align="center">
0.0990
  </td>
  <td align="center">
0.4790
  </td>
</tr>

<tr>
  <td align="center">
0.70
  </td>
  <td align="center">
3.109
  </td>
  <td align="center">
0.4775
  </td>
  <td align="center">
0.7104
  </td>

  <td align="center">
0.32
  </td>
  <td align="center">
10.2
  </td>
  <td align="center">
0.0744
  </td>
  <td align="center">
0.4038
  </td>
</tr>

<tr>
  <td align="center">
0.65
  </td>
  <td align="center">
3.296
  </td>
  <td align="center">
0.4481
  </td>
  <td align="center">
0.7124
  </td>

  <td align="center">
0.31
  </td>
  <td align="center">
12.4
  </td>
  <td align="center">
0.05536
  </td>
  <td align="center">
0.3264
  </td>
</tr>

<tr>
  <td align="center">
0.60
  </td>
  <td align="center">
3.523
  </td>
  <td align="center">
0.4142
  </td>
  <td align="center">
0.7125
  </td>

  <td align="center">
0.305
  </td>
  <td align="center">
14.4
  </td>
  <td align="center">
0.04494
  </td>
  <td align="center">
0.2772
  </td>
</tr>

<tr>
  <td align="center">
0.55
  </td>
  <td align="center">
3.809
  </td>
  <td align="center">
0.3748
  </td>
  <td align="center">
0.7096
  </td>

  <td align="center">
0.3
  </td>
  <td align="center">
17.733
  </td>
  <td align="center">
0.03412
  </td>
  <td align="center">
0.2186
  </td>
</tr>

<tr>
  <td align="center">
&frac12;
  </td>
  <td align="center">
4.186
  </td>
  <td align="center">
0.3284
  </td>
  <td align="center">
0.7014
  </td>

  <td align="center">
0.295
  </td>
  <td align="center">
25.737
  </td>
  <td align="center">
0.02165
  </td>
  <td align="center">
0.14347
  </td>
</tr>

<tr>
  <td align="center">
0.45
  </td>
  <td align="center">
4.719
  </td>
  <td align="center">
0.2733
  </td>
  <td align="center">
0.6830
  </td>

  <td align="center">
0.291
  </td>
  <td align="center">
75.510
  </td>
  <td align="center">
0.00666
  </td>
  <td align="center">
0.0450
  </td>
</tr>

<tr>
  <td align="center">
0.40
  </td>
  <td align="center">
5.574
  </td>
  <td align="center">
0.2073
  </td>
  <td align="center">
0.6429
  </td>

  <td align="center">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
</tr>
</table>

  </td>
   <td align="center" colspan="1" bgcolor="white">
[[Image:PlotStabilityBest02.png|500px|center]]
  </td>
</tr>

<tr>
 <td align="left" colspan="1">
'''Figure 3:'''  Largely the same as Figure 1, above, but a red-dashed curve has been added that separates the <math>~q - \nu</math> domain into regions that contain stable models (lying below the curve) from dynamically unstable models (lying above the curve), as determined by the virial stability analysis presented here.
  </td>
</tr>
</table>
</div>

=Related Discussions=
* [[SSC/Structure/PolytropesEmbedded#Embedded_Polytropic_Spheres|Polytropes emdeded in an external medium]]
* [[SSC/Structure/BiPolytropes#BiPolytropes|Constructing BiPolytropes]]
* [[SSC/Structure/BonnorEbert#Pressure-Bounded_Isothermal_Sphere|Bonnor-Ebert spheres]]
** [http://en.wikipedia.org/wiki/Bonnor-Ebert_mass Bonnor-Ebert Mass] according to Wikipedia
** <font color="red">Link has disappeared:</font>  [http://www.astro.umd.edu/~cychen/MATLAB/ASTR320/matlabFrom320spring2011/Bonnor-EbertSphere/html/BonnorEbert.html A MATLAB script to determine the Bonnor-Ebert Mass coefficient] developed by [http://www.astro.umd.edu/people/cychen.html Che-Yu Chen] as a graduate student in the University of Maryland Department of Astronomy
* [[SSC/Structure/LimitingMasses#Sch.C3.B6nberg-Chandrasekhar_Mass|Sch&ouml;nberg-Chandrasekhar limiting mass]]
* [[SSC/Structure/LimitingMasses#Relationship_Between_the_Bonnor-Ebert_and_Sch.C3.B6nberg-Chandrasekhar_Critical_Masses|Relationship between Bonnor-Ebert and Sch&ouml;nberg-Chandrasekhar limiting masses]]

{{ SGFfooter }}