Editing SSC/BipolytropeGeneralizationVersion2 (section)

===Bipolytrope===
When considering an isolated, spherically symmetric ''bi''polytropic configuration, each energy term will be made up of separate contributions coming from the core and envelope, that is,
<div align="center">
<math>
\mathfrak{G} = (W_\mathrm{grav} + \mathfrak{S}_A)_\mathrm{core} + (W_\mathrm{grav} + \mathfrak{S}_A)_\mathrm{env}  \, .
</math>
</div> 

====Partitioning the Mass====
The core will be principally defined in terms of two dimensionless parameters &#8212; <math>~q</math> and <math>~\nu</math> &#8212; which are, respectively, the core's radius relative to the bipolytrope's total radius, and the core's mass relative to the total mass of the bipolytropic configuration, specifically,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>q</math>
  </td>
  <td align="center">
<math>\equiv</math>
  </td>
  <td align="left">
<math>x_i = \frac{r_i}{R_\mathrm{edge}} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>\nu</math>
  </td>
  <td align="center">
<math>\equiv</math>
  </td>
  <td align="left">
<math>\frac{M_\mathrm{core}}{M_\mathrm{tot}} \, .</math>
  </td>
</tr>
</table>
</div>
Given the separate (equilibrium) density profiles of the core and the envelope while sticking to the notation used in our [[SSCpt1/Virial#Normalizations|introductory discussion]], we can write,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>(\mathrm{For}~0 \leq r^* \leq r_i^*)</math> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>M_r </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{4\pi}{3} \biggr) R_\mathrm{norm}^3 \rho_\mathrm{norm} \int_0^{r^*}  3 (r^*)^2 \rho_\mathrm{core}^* dr^* 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
M_\mathrm{tot} \cdot \chi^3 \int_0^{x}  3 \biggl[ \frac{\rho_\mathrm{core}(x)}{{\bar\rho}_\mathrm{core}} \biggr] 
\biggl[ \frac{M_\mathrm{core}/(x_i R_\mathrm{edge})^3}{M_\mathrm{tot}/R^3_\mathrm{norm}} \biggr] x^2 dx 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
M_\mathrm{tot} \biggl( \frac{\nu}{q^3} \biggr) \int_0^{x}  3
\biggl[ \frac{\rho(x)}{\bar\rho} \biggr]_\mathrm{core}  
x^2 dx \, ; 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>(\mathrm{For}~r_i^* \leq r^* \leq R_\mathrm{edge}^*)</math> &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>~M_r </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math> M_\mathrm{core} + 
\biggl( \frac{4\pi}{3} \biggr) R_\mathrm{norm}^3 \rho_\mathrm{norm} \int_{r_i^*}^{r^*}  3 (r^*)^2 \rho_\mathrm{env}^* dr^* 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>M_\mathrm{core} +
M_\mathrm{tot} \cdot \chi^3 \int_{x_i}^{x}  3 \biggl[ \frac{\rho_\mathrm{env}(x)}{{\bar\rho}_\mathrm{env}} \biggr] 
\biggl\{ \frac{M_\mathrm{env}/[ (1-x_i^3) R_\mathrm{edge}^3]}{M_\mathrm{tot}/R^3_\mathrm{norm}} \biggr\} x^2 dx 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
M_\mathrm{tot} \biggl\{\nu + \biggl( \frac{1-\nu}{1-q^3} \biggr) \int_{x_i}^{x}  3
\biggl[ \frac{\rho(x)}{\bar\rho} \biggr]_\mathrm{env}  
x^2 dx \biggr\} \, .
</math>
  </td>
</tr>
</table>
</div>

'''EXAMPLES'''
* [[SSC/Structure/BiPolytropes/FreeEnergy00#Mass_Profile|Mass profile of bipolytrope]] with <math>(n_c, n_e) = (0, 0)</math>.
* [[SSC/Structure/BiPolytropes/FreeEnergy51#Mass_Profile|Mass profile of bipolytrope]] with <math>(n_c, n_e) = (5, 1)</math>.

====Separate Contributions to Gravitational Potential Energy====

Given the separate (equilibrium) density and <math>~M_r</math> profiles of the core and the envelope while sticking to the notation used in our [[SSCpt1/Virial#Normalizations|introductory discussion]], we can write,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~W_\mathrm{grav}\biggr|_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
- E_\mathrm{norm} \int_0^{r_i^*} 3\biggl[\frac{M_r(r^*)}{M_\mathrm{tot}} \biggr]  r^* \rho^* dr^* 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
 <td align="left">
<math>
- E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{\nu}{q^3} \biggr)
\int_0^{x_i} 3x \biggl[\frac{M_r(x)}{M_\mathrm{tot}} \biggr]_\mathrm{core}  \biggl[ \frac{\rho(x)}{\bar\rho} \biggr]_\mathrm{core}  dx 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~W_\mathrm{grav}\biggr|_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
- E_\mathrm{norm} \int_{r_i^*}^{\chi = R_\mathrm{edge}^*} 3\biggl[\frac{M_r(r^*)}{M_\mathrm{tot}} \biggr]  r^* \rho^* dr^* 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
 <td align="left">
<math>
- E_\mathrm{norm} \cdot \chi^{-1} \biggl( \frac{1-\nu}{1-q^3} \biggr)
\int_{x_i}^{1} 3x \biggl[\frac{M_r(x)}{M_\mathrm{tot}} \biggr]_\mathrm{env}  \biggl[ \frac{\rho(x)}{\bar\rho} \biggr]_\mathrm{env}  dx 
</math>
  </td>
</tr>

</table>
</div>


'''EXAMPLES'''
* [[SSC/Structure/BiPolytropes/FreeEnergy00#Gravitational_Potential_Energy|Gravitational Potential Energy of bipolytrope]] with <math>(n_c, n_e) = (0, 0)</math>.
* [[SSC/Structure/BiPolytropes/FreeEnergy51#Gravitational_Potential_Energy|Gravitational Potential Energy of bipolytrope]] with <math>(n_c, n_e) = (5, 1)</math>.

====Separate Thermodynamic Energy Reservoirs====
In our [[SSCpt1/Virial#Looking_Ahead_to_Bipolytropes|introductory discussion of the free energy function for spherically symmetric configurations]], we developed expressions that define the separate contributions to the thermodynamic energy reservoir that pertain to the core and envelope of bipolytropic configurations.  In that discussion we pointed out that, in general for the core, the pressure drops monotonically from a value of <math>~P_0</math> at the center of the configuration according to an expression of the form,
<div align="center">
<math>~P_\mathrm{core}(x) = P_0 [1 - p_c(x)]</math> &nbsp; &nbsp;&nbsp; for &nbsp; &nbsp;&nbsp; <math>~0 \leq x \leq q \, ,</math>
</div>
and that, for the envelope, the pressure drops monotonically from a value of <math>~P_{ie}</math> at the interface according to an expression of the form,
<div align="center">
<math>~P_\mathrm{env}(x) = P_{ie} [1 - p_e(x)]</math> &nbsp; &nbsp;&nbsp; for &nbsp; &nbsp;&nbsp; <math>~q \leq x \leq 1 \, ,</math>
</div>
where <math>~p_c(x)</math> and <math>~p_e(x)</math> are both dimensionless functions that will depend on the equations of state that are chosen for the core and envelope, respectively.  By prescription, the pressure in the envelope must drop to zero at the surface of the bipolytropic configuration, hence, we should expect that <math>~p_e(1) = 1</math>.  Furthermore, by prescription, the pressure in the core will drop to a value, <math>~P_{ic}</math>, at the interface, so we can write,
<div align="center">
<math>~P_{ic} = P_0 [1 - p_c(q)] \, .</math>
</div>

In equilibrium &#8212; that is, when <math>~R_\mathrm{edge} = R_\mathrm{eq}</math> &#8212; we will demand that the pressure at the interface be the same, whether it is referenced in the core or in the envelope, that is, we will demand that <math>~P_{ic} = P_{ie} \, .</math>  It is therefore strategically advantageous to rewrite the expression for the run of pressure through the core in terms of the pressure at the interface rather than in terms of the central pressure; specifically,
<div align="center">
<math>~P_\mathrm{core}(x) = P_{ic} \biggl[\frac{1 - p_c(x)}{1-p_c(q)} \biggr] \, .</math> 
</div>

<div align="center">
<table border="1" cellpadding="5" width="75%">
<tr><td align="left">
<div align="center">
<font color="red" size="+1"> KEY QUESTION:  How should piecewise pressure be normalized? </font>
</div>
Here we will operate under two constraints when relating the central pressure, <math>~P_0</math>, to the ''core'' pressure at the interface, <math>~P_{ic}</math>.  First, the pressure in the core will drop monotonically from its central value according to the following very general prescription,
<div align="center">
<math>~P_\mathrm{core}(x) = P_0 [1 - p_c(x)] \, .</math>
</div>
Second, it is the pressure ''at the interface'', <math>~P_{ic} = P_\mathrm{core}(q)</math>, that will be determined by the virial equilibrium condition.  There are two possible ways to determine the central pressure from knowledge of <math>~P_{ic}</math>, but which is the physically correct method to embrace?

<font color="darkgreen">'''Case A:'''</font>  Addition &hellip;
<div align="center">
<math>~P_0 = P_{ic} + \Delta P \, ,</math>&nbsp;&nbsp;&nbsp;&nbsp; where &nbsp;&nbsp;&nbsp;&nbsp;<math>\Delta P \equiv P_0 p_c(q) \, ,</math>
</div>
in which case, we should write that,
<div align="center">
<math>~P_\mathrm{core}(x) = ( P_{ic} +  \Delta P ) - P_0 p_c(x)  \, .</math>
</div>


<font color="darkgreen">'''Case B:'''</font>  Multiplication &hellip;
<div align="center">
<math>~P_{ic} = P_0 [1-p_c(q)] ~~~~\Rightarrow~~~~~ P_0 = \frac{P_{ic}}{[1-p_c(q)]} \, ,</math>
</div>
in which case, we should write that,
<div align="center">
<math>~P_\mathrm{core}(x) = P_{ic} \biggl[ \frac{ 1-p_c(x) }{ 1-p_c(q) } \biggr]  \, .</math>
</div>

Prior to August 2014, we have been naively implementing "Case A,"  effectively assuming that the quantity, <math>~\Delta P</math> (as well as <math>~P_{ic}</math>), is held fixed as we search for the equilibrium value of <math>~P_0</math>.  See, for example, the comment dated <font color="red">12 February 2014</font> in connection with my [[SSC/VirialStability#Discussion_with_Kundan_Kadam|discussions with Kundan Kadam]], or even [[SSC/BipolytropeGeneralization#DerivationTable|the "new derivation" summarized in the table below]], where we have set,
<div align="center">
<math>~\Delta P = \Pi q^2 \, .</math>
</div>

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

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


With these generic expressions for the pressure profile in hand, the separate components of the thermodynamic energy reservoir derived in our [[SSCpt1/Virial#Looking_Ahead_to_Bipolytropes|introductory discussion]] are,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\biggl( \frac{\mathfrak{S}_A}{E_\mathrm{norm}} \biggr)_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\frac{4\pi }{({\gamma_c}-1)}  \biggl[ \frac{P_{ic} \chi^{3\gamma_c}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_c} 
\int_0^q  \biggl[\frac{1 - p_c(x)}{1-p_c(q)} \biggr]  x^2 dx
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\biggl( \frac{\mathfrak{S}_A}{E_\mathrm{norm}} \biggr)_\mathrm{env}</math>
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>
\frac{4\pi }{({\gamma_e}-1)}  \biggl[ \frac{P_{ie} \chi^{3\gamma_e}}{P_\mathrm{norm}} \biggr]_\mathrm{eq} \chi^{3-3\gamma_e} 
\int_q^1  \biggl[1 - p_e(x) \biggr]  x^2 dx \, .
</math>
  </td>
</tr>
</table>
</div>


'''EXAMPLES'''
* [[SSC/Structure/BiPolytropes/FreeEnergy00#Thermodynamic_Energy_Reservoir|Thermodynamic Energy Reservoir of bipolytrope]] with <math>(n_c, n_e) = (0, 0)</math>.
* [[SSC/Structure/BiPolytropes/FreeEnergy51#Thermodynamic_Energy_Reservoir|Thermodynamic Energy Reservoir of bipolytrope]] with <math>(n_c, n_e) = (5, 1)</math>.