Editing SSC/FreeEnergy/PolytropesEmbedded/Pt3C

=Free Energy of Embedded Polytropes=
<table border="1" align="center" width="100%" colspan="8">
<tr>
  <td align="center" rowspan="2" bgcolor="lightblue" width="25%"><br />[[SSC/FreeEnergy/PolytropesEmbedded/Pt1|Part I: &nbsp; Synopsis]]<br />&nbsp;</td>
  <td align="center" rowspan="2" bgcolor="lightblue" width="25%"><br />[[SSC/FreeEnergy/PolytropesEmbedded/Pt2|Part II:&nbsp; Truncated Polytropes]]<br />&nbsp;</td>
  <td align="center" rowspan="1" colspan="3"bgcolor="lightblue"><br />Part III: &nbsp;Free-Energy of Bipolytropes<br />&nbsp;</td>
</tr>
<tr>
  <td align="center" rowspan="1" bgcolor="lightblue" width="17%"><br />[[SSC/FreeEnergy/PolytropesEmbedded/Pt3A|IIIA:&nbsp; Focus on (5, 1) Bipolytropes]]<br />&nbsp;</td>
  <td align="center" rowspan="1" bgcolor="lightblue" width="17%"><br />[[SSC/FreeEnergy/PolytropesEmbedded/Pt3B|IIIB:&nbsp; Focus on (0, 0) Bipolytropes]]<br />&nbsp;</td>
  <td align="center" rowspan="1" bgcolor="lightblue"><br />[[SSC/FreeEnergy/PolytropesEmbedded/Pt3C|IIIC:&nbsp; Overview]]<br />&nbsp;</td>
</tr>
</table>

=Free-Energy of Bipolytropes=

In this case, the Gibbs-like free energy is given by the sum of four separate energies,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{G}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[W_\mathrm{grav} + \mathfrak{S}_\mathrm{therm}\biggr]_\mathrm{core} + \biggl[W_\mathrm{grav} + \mathfrak{S}_\mathrm{therm}\biggr]_\mathrm{env} \, .
</math>
  </td>
</tr>

</table>
</div>

In addition to specifying (generally) separate polytropic indexes for the core, <math>~n_c</math>, and envelope, <math>~n_e</math>, and an envelope-to-core mean molecular weight ratio, <math>~\mu_e/\mu_c</math>, we will assume that the system is fully defined via specification of the following five physical parameters:
* Total mass, <math>~M_\mathrm{tot}</math>;
* Total radius, <math>~R</math>; 
* Interface radius, <math>~R_i</math>, and associated dimensionless interface marker, <math>~q \equiv R_i/R</math>;
* Core mass, <math>~M_c</math>, and associated dimensionless mass fraction, <math>~\nu \equiv M_c/M_\mathrm{tot}</math>;
* Polytropic constant in the core, <math>~K_c</math>.

In general, the warped free-energy surface drapes across a five-dimensional parameter "plane" such that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{G}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\mathfrak{G}(R, K_c, M_\mathrm{tot}, q, \nu) \, .</math>
  </td>
</tr>
</table>
</div>
==Overview==

===BiPolytrope51===

====Key Analytic Expressions====

<div align="center" id="FreeEnergy51">
<table border="1" cellpadding="5" align="center">
<tr>
<th align="center">
<font size="+1">Out-of-Equilibrium, Free-Energy Expression for BiPolytropes with</font> <math>~(n_c, n_e) = (5, 1)</math>
</th>
</tr>
<tr><td align="center">

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

<tr>
  <td align="right">
<math>~\mathfrak{G}^*_{51} \equiv 2^4\biggl( \frac{q}{\nu^2}\biggr) \chi_\mathrm{eq} \biggl[\frac{\mathfrak{G}_{51}}{E_\mathrm{norm}} \biggr]</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{1}{\ell_i^2} \biggl[
\Chi^{-3/5} (5 \mathfrak{L}_i)
+\Chi^{-3} (4\mathfrak{K}_i)
-\Chi^{-1} (3\mathfrak{L}_i  +12\mathfrak{K}_i ) \biggr]
</math>
  </td>
</tr>
</table>
</td></tr>
</table>
</div>
where,

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

<tr>
  <td align="right">
<math>~\mathfrak{L}_i</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{(\ell_i^4-1)}{\ell_i^2} + \frac{(1+\ell_i^2)^3}{\ell_i^3} \cdot \tan^{-1}\ell_i \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathfrak{K}_i</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{\Lambda_i}{\eta_i}  + \frac{(1+\Lambda_i^2)}{\eta_i} \biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i\biggr] \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Lambda_i</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{1}{\eta_i} - \ell_i \, ,</math>
  </td>
</tr>

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

From the [[SSC/Structure/BiPolytropes/Analytic51#Parameter_Values|accompanying Table 1 parameter values]], we also can write,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{1}{q}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\eta_s}{\eta_i}
= 1 + \frac{1}{\eta_i}\biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr] \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\nu</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\ell_i q}{(1+\Lambda_i^2)^{1/2}} \, .
</math>
  </td>
</tr>
</table>
</div>

Consistent with our [[SSC/BipolytropeGeneralization#Free_Energy_and_Its_Derivatives|generic discussion of the stability of bipolytropes]] and the ''specific'' discussion of [[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|the stability of bipolytropes having]] <math>~(n_c, n_e) = (5, 1)</math>, it can straightforwardly be shown that <math>~\partial \mathfrak{G}^*_{51}/\partial \chi = 0</math> is satisfied by setting <math>~\Chi = 1</math>; that is, the equilibrium condition is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\chi = \chi_\mathrm{eq}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(\frac{\pi}{2^3}\biggr)^{1/2} \frac{\nu^2}{q} \cdot \frac{(1+\ell_i^2)^3}{3^3\ell_i^5} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl\{ \biggl(\frac{\pi}{3}\biggr) 2^{2-n_c} \nu^{n_c-1} q^{3-n_c} \biggl[\frac{(1+\ell_i^2)^{6/5}}{(3\ell_i^2)} \biggr]^{n_c} \biggr\}^{1/(n_c-3)} \, ,</math>
  </td>
</tr>
</table>
</div>

where the last expression has been cast into a form that more clearly highlights overlap with the expression, below, for the equilibrium radius for zero-zero bipolytropes.  Furthermore, the equilibrium configuration is unstable whenever,
<div align="center">
<math>~\biggl[ \frac{\partial^2 \mathfrak{G}^*_{51}}{\partial \chi^2}\biggr]_{\Chi=1} < 0 \, ,</math>
</div>
that is, it is unstable whenever,

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

<tr>
  <td align="right">
<math>~\frac{ \mathfrak{L}_i}{\mathfrak{K}_i}</math>
  </td>
  <td align="center">
<math>~></math>
  </td>
  <td align="left">
<math>~20 \, .</math>
  </td>
</tr>
</table>
</div>
[[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|Table 1 of an accompanying chapter]] &#8212; and the red-dashed curve in the figure adjacent to that table &#8212; identifies some key properties of the model that marks the transition from stable to unstable configurations along equilibrium sequences that have various values of the mean-molecular weight ratio, <math>~\mu_e/\mu_c</math>.

====Behavior of Equilibrium Sequence====

Here we reprint Figure 1 from an [[SSC/Structure/BiPolytropes/Analytic51#Model_Sequences|accompanying chapter wherein the structure of five-one bipolytropes has been derived]]. It displays detailed force-balance sequences in the <math>~q - \nu</math> plane for a variety of choices of the ratio of mean-molecular-weights, <math>~\mu_e/\mu_c</math>, as labeled.

[[File:PlotSequencesBest02.png|450px|center|Five-One Bipolytropic Equilibrium Sequences for Various ratios of the mean molecular weight]]

=====Limiting Values=====
Each sequence begins <math>~(\ell_i = 0)</math> at the origin, that is, at <math>~(q,\nu) = (0,0)</math>.  As <math>~\ell_i \rightarrow \infty</math>, however, the sequences terminate at different coordinate locations, depending on the value of <math>~m_3 \equiv 3(\mu_e/\mu_c)</math>.  In deriving the various limits, it will be useful to note that,
<div align="center">
<table border="0" cellpadding="5" align="center">

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

<tr>
  <td align="right">
<math>~\Lambda_i</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{(1+\ell_i^2)}{m_3\ell_i}-\ell_i</math>
  </td>
</tr>

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

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

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

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

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

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{m_3^2\ell_i^2}\biggl\{ 1 + (2 -2m_3 + m_3^2) \ell_i^2 + (1 -m_3)^2 \ell_i^4 \biggr\}</math>
  </td>
</tr>
</table>
</div>

Examining the three relevant parameter regimes, we see that:

* For <math>~\mu_e/\mu_c < \tfrac{1}{3}</math>, that is, <math>~m_3 < 1</math> &hellip;
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\tan^{-1} \Lambda_i \biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx </math>
  </td>
  <td align="left">
<math>~\tan^{-1} \biggl[\frac{(1 -m_3)}{m_3} \biggr]\ell_i</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx </math>
  </td>
  <td align="left">
<math>~\frac{\pi}{2} -  \biggl[\frac{m_3}{(1 -m_3)\ell_i} \biggr]</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{1}{q}\biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ 
1 + \frac{(1 + \ell_i^2)}{m_3 \ell_i}\biggl[\pi - \frac{m_3}{(1 -m_3)\ell_i} \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ 
\frac{m_3 + \pi \ell_i}{m_3}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ q\biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ 
\frac{1}{1 + (\pi \ell_i/m_3)} \rightarrow 0 \, .
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
and
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>

<tr>
  <td align="right">
<math>~\biggl(\frac{\nu}{q}\biggr)^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_i^2}{1 + \Lambda_i^2}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~m_3^2\ell_i^4 \biggl\{ 1 + (2 -2m_3 + m_3^2) \ell_i^2 + (1 -m_3)^2 \ell_i^4 \biggr\}^{-1}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~m_3^2\biggl\{ \ell_i^{-4} + (2 -2m_3 + m_3^2) \ell_i^{-2} + (1 -m_3)^2 \biggr\}^{-1}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{\nu}{q}\biggr|_{\ell_i\rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~\frac{m_3}{1-m_3} </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \nu \biggr|_{\ell_i\rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~\biggl[\frac{m_3}{1-m_3}\biggr]\frac{1}{1 + (\pi \ell_i/m_3)} \rightarrow 0 \, . </math>
  </td>
</tr>
</table>
</div>

*  For <math>~\mu_e/\mu_c = \tfrac{1}{3}</math>, that is, <math>~m_3 = 1</math> &hellip;
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\tan^{-1} \Lambda_i </math>
  </td>
  <td align="center">
<math>~= </math>
  </td>
  <td align="left">
<math>~\tan^{-1} \biggl(\frac{1}{\ell_i}\biggr)</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \tan^{-1} \Lambda_i \biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx </math>
  </td>
  <td align="left">
<math>~\frac{1}{\ell_i}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{1}{q}\biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~
1 + \frac{(1 + \ell_i^2)}{\ell_i}\biggl[\frac{\pi}{2} + \frac{1}{\ell_i }\biggr] </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{\pi}{2}\biggr)\ell_i</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ q\biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~
\frac{2}{\pi \ell_i} \rightarrow 0</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
and
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>

<tr>
  <td align="right">
<math>~\biggl(\frac{\nu}{q}\biggr)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_i}{(1 + 1/\ell_i^2)^{1/2}}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \nu \biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~\ell_i \biggl(\frac{2}{\pi \ell_i} \biggr) = \frac{2}{\pi} \approx 0.63662</math>
  </td>
</tr>

</table>
</div>


*  For <math>~\mu_e/\mu_c > \tfrac{1}{3}</math>, that is, <math>~m_3 > 1</math> &hellip;
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\tan^{-1} \Lambda_i \biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx </math>
  </td>
  <td align="left">
<math>~\tan^{-1} \biggl[-\biggl(\frac{m_3-1}{m_3} \biggr)\ell_i\biggr]</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx </math>
  </td>
  <td align="left">
<math>~-\frac{\pi}{2} +  \biggl[\frac{m_3}{(m_3-1)\ell_i} \biggr]</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{1}{q}\biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ 
1 + \frac{(1 + \ell_i^2)}{m_3 \ell_i}\biggl[ \frac{m_3}{(m_3-1)\ell_i} \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ 
1 + \frac{(1 + 1/\ell_i^2)}{(m_3-1) } 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ 
1 + \frac{1}{(m_3-1) } = \frac{m_3}{(m_3-1)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ q\biggr|_{\ell_i \rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ 
\frac{(m_3-1)}{m_3}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
and
  </td>
  <td align="left">
&nbsp;
  </td>
</tr>

<tr>
  <td align="right">
<math>~\biggl(\frac{\nu}{q}\biggr)^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_i^2}{1 + \Lambda_i^2}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~m_3^2\ell_i^4 \biggl\{ 1 + (2 -2m_3 + m_3^2) \ell_i^2 + (m_3-1)^2 \ell_i^4 \biggr\}^{-1}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~m_3^2\biggl\{ \ell_i^{-4} + (2 -2m_3 + m_3^2) \ell_i^{-2} + (m_3-1)^2 \biggr\}^{-1}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{\nu}{q}\biggr|_{\ell_i\rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~\frac{m_3}{m_3-1} </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \nu \biggr|_{\ell_i\rightarrow \infty}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{m_3}{m_3-1} \biggl[\frac{m_3 - 1}{m_3}\biggr] \rightarrow 1 \, . </math>
  </td>
</tr>

</table>
</div>

Summarizing:

* For <math>~\mu_e/\mu_c < \tfrac{1}{3}</math>, that is, <math>~m_3 < 1</math> &nbsp;  &nbsp;  &nbsp; &hellip; &nbsp;  &nbsp;  &nbsp;  <math>~(q,\nu)_{\ell_i \rightarrow \infty} = (0, 0) \, .</math>


*  For <math>~\mu_e/\mu_c = \tfrac{1}{3}</math>, that is, <math>~m_3 = 1</math> &nbsp;  &nbsp;  &nbsp; &hellip; &nbsp;  &nbsp;  &nbsp;  <math>~(q,\nu)_{\ell_i \rightarrow \infty} = (0, \tfrac{2}{\pi}) \, .</math>


*  For <math>~\mu_e/\mu_c > \tfrac{1}{3}</math>, that is, <math>~m_3 > 1</math> &nbsp;  &nbsp;  &nbsp; &hellip; &nbsp;  &nbsp;  &nbsp;  <math>~(q,\nu)_{\ell_i \rightarrow \infty} = [(m_3-1)/m_3, 1] \, .</math>

=====Turning Points=====

Let's identify the location of two turning points along the <math>~\nu(q)</math> sequence &#8212; one defines <math>~q_\mathrm{max}</math> and the other identifies <math>~\nu_\mathrm{max}</math>.  They occur, respectively, where,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{d\ln q}{d\ln \ell_i} = 0</math>
  </td>
  <td align="center">
&nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; 
  </td>
  <td align="left">
<math>~\frac{d\ln \nu}{d\ln \ell_i} = 0 \, .</math>
  </td>
</tr>
</table>
</div>

In deriving these expressions, we will use the relations,
<div align="center">
<table border="0" cellpadding="5" align="center">

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

<tr>
  <td align="right">
<math>~\frac{d\Lambda_i}{d\ell_i}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~- \frac{1}{m_3\ell_i^2} \biggl[ 1 - \ell_i^2(1-m_3) \biggr] \, ,</math>
  </td>
</tr>
</table>
</div>
where, 
<div align="center">
<math>~m_3 \equiv 3\biggl(\frac{\mu_e}{\mu_c}\biggr) \, .</math>
</div>


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

<tr>
  <td align="right">
<math>~q </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl\{ 1 + \frac{1}{\eta_i}\biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr] \biggr\}^{-1} \, ,</math>
  </td>
</tr>
</table>
</div>

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

<tr>
  <td align="right">
<math>~\frac{d\ln q}{d\ln \ell_i}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_i}{q} \cdot ( -q^2) \frac{d}{d\ell_i} \biggl\{
\frac{1}{\eta_i}\biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr]
\biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~-q\ell_i \biggl\{-\frac{1}{\eta_i^2}\biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr]\frac{d\eta_i}{d\ell_i}
+ \frac{1}{\eta_i(1+\Lambda_i^2)} \frac{d\Lambda_i}{d\ell_i}
\biggr\}
</math>
  </td>
</tr>

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

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{q}{m_3^2}{\ell_i^2}\biggl\{m_3 \ell_i (1-\ell_i^2) \biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr]
+ \frac{(1+\ell_i^2)}{(1+\Lambda_i^2)} \biggl[ 1 - \ell_i^2(1-m_3) \biggr] 
\biggr\} \, .
</math>
  </td>
</tr>
</table>
</div>

And, given that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\nu </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\ell_i q}{(1+\Lambda_i^2)^{1/2}} \, .
</math>
  </td>
</tr>
</table>
</div>

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

<tr>
  <td align="right">
<math>~\frac{d\ln \nu}{d\ln \ell_i}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\ell_i}{\nu} \biggl\{
\frac{q}{(1+\Lambda_i^2)^{1/2}} + \frac{q}{(1+\Lambda_i^2)^{1/2}} \frac{d\ln q}{d\ln \ell_i}
- \frac{\ell_i q \Lambda_i }{(1+\Lambda_i^2)^{3/2}} \frac{d\Lambda_i}{d\ell_i}
\biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{q \ell_i}{ \nu(1+\Lambda_i^2)^{1/2}}\biggl\{
1 + \frac{d\ln q}{d\ln \ell_i}
+ \frac{\Lambda_i }{m_3 \ell_i (1+\Lambda_i^2)} \biggl[ 1 - \ell_i^2(1-m_3) \biggr] 
\biggr\}
</math>
  </td>
</tr>
</table>
</div>

In summary, then, the <math>~q_\mathrm{max}</math> turning point occurs where,

<div align="center">
<table border="0" cellpadding="5" align="center">
<!--
<tr>
  <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ m_3 \ell_i (1-\ell_i^2) \biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr]
+ \frac{(1+\ell_i^2)}{(1+\Lambda_i^2)} \biggl[ 1 - \ell_i^2(1-m_3) \biggr] 
\, ;
</math>
  </td>
</tr>
-->
<tr>
  <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ (1+\Lambda_i^2)\biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr]
+ \frac{(1+\ell_i^2)}{m_3 \ell_i (1-\ell_i^2)} \biggl[ 1 - \ell_i^2(1-m_3) \biggr] 
\, ;
</math>
  </td>
</tr>
</table>
</div>
 
and the <math>~\nu_\mathrm{max}</math> turning point occurs where,

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

<tr>
  <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
1 + \frac{\Lambda_i }{m_3 \ell_i (1+\Lambda_i^2)} \biggl[ 1 - \ell_i^2(1-m_3) \biggr]  
+ \frac{q \ell_i^3  (1-\ell_i^2)}{m_3} \biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr]
+ \frac{q\ell_i^2}{m_3^2}\cdot \frac{(1+\ell_i^2)}{(1+\Lambda_i^2)} \biggl[ 1 - \ell_i^2(1-m_3) \biggr] 
</math>
  </td>
</tr>

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

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
1 
+ \frac{q \ell_i^3  (1-\ell_i^2)}{m_3} \biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr]
+ \frac{1}{m_3 \ell_i} \biggl[ \frac{\Lambda_i }{(1+\Lambda_i^2)}  
+ \frac{q\ell_i^3}{m_3}\cdot \frac{(1+\ell_i^2)}{(1+\Lambda_i^2)} \biggr] \cdot \biggl[ 1 - \ell_i^2(1-m_3) \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>

<table border="0" width="80%" cellpadding="5" align="center"><tr><td align="left">
<font color="red"><b>NOTE:</b></font>&nbsp; As we show [[#Limiting_Values|above]], for the special case of <math>~m_3 = 1</math> &#8212; that is, when <math>~\mu_e/\mu_c = \tfrac{1}{3}</math>, precisely &#8212; the equilibrium sequence (as <math>~\ell_i \rightarrow \infty</math>) intersects the <math>~q = 0</math> axis at precisely the value, <math>~\nu = 2/\pi</math>.  As is illustrated graphically in [[SSC/Structure/BiPolytropes/Analytic51#Model_Sequences|Figure 1 of an accompanying chapter]], no <math>~\nu_\mathrm{max}</math> turning point exists for values of  <math>~m_3 > 1</math>.  
</td></tr>
</table>

For the record, we repeat, as well, that the transition from stable to dynamically unstable configurations occurs along the sequence when,

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

<tr>
  <td align="right">
<math>~\frac{(\ell_i^4-1)}{\ell_i^2} + \frac{(1+\ell_i^2)^3}{\ell_i^3} \cdot \tan^{-1}\ell_i </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~20 \biggl\{ \frac{\Lambda_i}{\eta_i}  + \frac{(1+\Lambda_i^2)}{\eta_i} \biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i\biggr] \biggr\} </math>
  </td>
</tr>

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

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ m_3 \ell_i (\ell_i^4-1) + m_3(1+\ell_i^2)^3\cdot \tan^{-1}\ell_i </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~20\ell_i^2 (1+\Lambda_i^2)(1+\ell_i^2) \biggl\{ \frac{\Lambda_i}{(1+\Lambda_i^2)}  + \biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i\biggr] \biggr\} </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{m_3 \ell_i (\ell_i^4-1) + m_3(1+\ell_i^2)^3\cdot \tan^{-1}\ell_i }{ 20\ell_i^2 (1+\ell_i^2)}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\Lambda_i  + (1+\Lambda_i^2)\biggl[\frac{\pi}{2} + \tan^{-1}\Lambda_i\biggr]  \, .</math>
  </td>
</tr>

</table>
</div>


In order to clarify what equilibrium sequences do not have any turning points, let's examine how the <math>~q_\mathrm{max}</math> turning-point expression behaves as <math>~\ell_i \rightarrow \infty</math>.

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

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

<tr>
  <td align="right">
<math>~ \Rightarrow ~~~ \frac{(1+\ell_i^2)}{ m_3 \ell_i (\ell_i^2-1)   } \biggl[ 1 + \ell_i^2(m_3-1) \biggr]</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \biggl\{ 1 + \frac{1}{m_3^2\ell_i^2}\biggl[ (m_3-1) \ell_i^2-1 \biggr]^2 \biggr\} \biggl\{ \frac{\pi}{2} 
+ \biggl[ -\frac{\pi}{2} - \frac{1}{\Lambda_i} + \frac{1 }{3\Lambda_i^3} + \mathcal{O}(\Lambda_i^{-5} )\biggr] \biggr\}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~ \Rightarrow ~~~ \frac{(1+\ell_i^2) \ell_i^2(m_3-1)}{ m_3 \ell_i (\ell_i^2-1)   } \biggl[ 1 + \frac{1}{\ell_i^2(m_3-1)} \biggr]</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \biggl\{ 1 + \frac{(m_3-1)^2 \ell_i^2}{m_3^2}\biggl[ 1 - \frac{1}{ (m_3-1) \ell_i^2 } \biggr]^2 \biggr\} \cdot  
\frac{1}{(-\Lambda_i)} \biggl[  1 - \frac{1 }{3\Lambda_i^2} + \cancelto{0}{\mathcal{O}(\Lambda_i^{-4} )}\biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~ \Rightarrow ~~~ \frac{(1+\ell_i^2)}{ \ell_i (\ell_i^2-1)   } \cdot \frac{m_3}{(m_3-1)} \biggl[ 1 + \frac{1}{\ell_i^2(m_3-1)} \biggr]</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \biggl[ 1 - \frac{2}{ (m_3-1) \ell_i^2 } + \frac{m_3^2}{(m_3-1)^2 \ell_i^2}  +  \frac{1}{ (m_3-1)^2 \ell_i^4 } \biggr] \cdot  \frac{m_3}{(m_3-1)\ell_i} \biggl[1 - \frac{1}{(m_3-1)\ell_i^2}\biggr]^{-1}
\biggl\{  1 - \frac{m_3^2}{3(m_3-1)^2\ell_i^2} \biggl[1 - \frac{1}{(m_3-1)\ell_i^2}\biggr]^{-2} \biggr\} 
</math>
  </td>
</tr>

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

</table>
</div>

The leading-order term is unity on both sides of this expression, so they cancel; let's see what results from keeping terms <math>~\propto \ell_i^{-2}</math>.

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

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

<tr>
  <td align="right">
<math>~\Rightarrow~~~ 2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ - \frac{2}{ (m_3-1) } + \frac{2m_3^2}{3(m_3-1)^2 }   
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~ 6(m_3-1)^2 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ - 6(m_3-1) + 2m_3^2   
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~ 6m_3^2-12m_3 + 6 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ - 6m_3+6 + 2m_3^2   
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~ m_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{3}{2} \, .  
</math>
  </td>
</tr>
</table>
</div>

We therefore conclude that the <math>~q_\mathrm{max}</math> turning point does not appear along any sequence for which,

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

<tr>
  <td align="right">
<math>~m_3</math>
  </td>
  <td align="center">
<math>~></math>
  </td>
  <td align="left">
<math>~\frac{3}{2}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{\mu_e}{\mu_c}</math>
  </td>
  <td align="center">
<math>~></math>
  </td>
  <td align="left">
<math>~\frac{1}{2}\, .</math>
  </td>
</tr>
</table>
</div>


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

<tr>
  <td align="right">
<math>~ \frac{1}{\ell_i^2} \biggl[ 2 + \frac{1}{\alpha m_3} \biggr]</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ \frac{1}{\ell_i^2}\biggl[
\frac{1}{\alpha^2 }  - \frac{2}{ m_3 \alpha  } -\frac{1}{3\alpha^2 } \biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \biggl[ 2\alpha^2 m_3 + \alpha \biggr]</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ \biggl[
\frac{2m_3}{3}  - 2\alpha  \biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ 2\alpha^2 m_3 </math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~ 
\frac{2m_3}{3}  - 3\alpha  
</math>
  </td>
</tr>
</table>
</div>
-->

<div align="center">
<table border="1" align="center" cellpadding="5">
<tr>
  <td align="center" colspan="2">Five-One Bipolytrope Equilibrium Sequences in <math>~q - \nu</math> Plane</td>
</tr>
<tr>
  <td align="center" width="50%">
Full Sequences for Various <math>~\frac{\mu_e}{\mu_c}</math>
  </td>
  <td align="center" width="50%">
Magnified View with Turning Points and Stability Transition-Points Identified
  </td>
</tr>
<tr>
  <td align="center" colspan="2">
[[File:Qvsnu51combined.png|750 px|Five-One Sequences]]
  </td>
</tr>
</table>
</div>

====Graphical Depiction of Free-Energy Surface====

<div align="center">
<table border="0" cellpadding="8" width="880px">
<tr><td align="center"><font size="+1"><b>
Figure 1: &nbsp; Free-Energy Surface for <math>~(n_c,n_e) = (5,1)</math> and <math>~\frac{\mu_e}{\mu_c} = 1</math></b></font></td></tr>
<tr><td align="center">

<table border="1" cellpadding="5">
<tr>
<td align="center" bgcolor="#CCFFFF">
[[File:FreeEnergy51Surface2.png|center|300px|Free-Energy surface for 5_1 bipolytrope]]
</td>
<td align="center" bgcolor="#CCFFFF">
[[File:Bipolytrope51Muratio1.gif|center|439px|Free-Energy surface for 5_1 bipolytrope]]
</td>
</tr>
<tr>
  <td align="left" colspan="2">
''Left Panel:''  The free energy (vertical, blue axis) is plotted as a function of the radial interface location, <math>~\xi_i</math> (red axis), and the normalized configuration radius, <math>~\Chi \equiv \chi/\chi_\mathrm{eq}</math> (green axis).  ''Right Panel:'' Same as the left panel, but animated in order to highlight undulations of the surface.  The value of the free energy is indicated by color as well as by the height of the warped surface &#8212; red identifies the lowest depicted energies while blue identifies the highest depicted energies; these same colors have been projected down onto the <math>~z = 0</math> plane to present a two-dimensional, color-contour plot.  A multi-colored line segment drawn parallel to the <math>~\xi_i</math> axis at the value, <math>~\Chi = 1</math>, identifies the configuration's ''equilibrium'' radius for each value of the interface location.  Equilibrium configurations marked in white lie at the bottom of the principal free-energy "valley" and are stable, while configurations marked in blue lie at the top of a free-energy "hill," indicating that they are unstable; the red dot identifies the location along this equilibrium sequence where the transition from stable to unstable configurations occurs.
  </td>
</tr>
</table>

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

For purposes of reproducibility, it is incumbent upon us to clarify how the values of the free energy were normalized in order to produce the free-energy surface displayed in Figure 1.  The normalization steps are explicitly detailed within the [[#Fortran_Code|fortran program, below]], that generated the data for plotting purposes; here we provide a brief summary.  We evaluated the normalized free energy, <math>~\mathfrak{G}^*_{51}</math>, across a <math>~200 \times 200</math> zone grid of uniform spacing, covering the following <math>~(x,y) = (\ell_i,\Chi)</math> domain:
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{1}{\sqrt{3}}</math>
  </td>
  <td align="center">
<math>~\le \ell_i \le</math>
  </td>
  <td align="left">
<math>~\frac{3}{\sqrt{3}}</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~0.469230769</math>
  </td>
  <td align="center">
<math>~\le \Chi\le</math>
  </td>
  <td align="left">
<math>~2.0</math>
  </td>
</tr>
</table>
</div>

(With this specific definition of the ''y''-coordinate grid, <math>~\Chi = 1</math> is associated with zone 70.)  After this evaluation, a constant, <math>~E_\mathrm{fudge} = -10</math>, was added to <math>~\mathfrak{G}^*</math> in order to ensure that the free energy was negative across the entire domain.  Then (inorm = 1), for each specified interface location, <math>~x = \ell_i</math>, employing the ''equilibrium'' value of the free energy,
<div align="center">
<math>~E_0 = \mathfrak{G}^*_{51}(\ell_i, \Chi = 1) + E_\mathrm{fudge} \, ,</math>
</div>
the free energy was normalized across all values of <math>~y = \Chi</math> via the expression,
<div align="center">
<math>~\mathrm{fe} = \frac{(\mathfrak{G}^*_{51} + E_\mathrm{fudge}) - (E_0)_i}{|E_0|_i} \, .</math>
</div>

Finally, for plotting purposes, at each grid cell vertex ("vertex") &#8212; as well as at each grid cell center ("cell") &#8212; the value of the free energy, <math>~\mathrm{fe}</math>, was renormalized as follows,
<div align="center">
<math>~\mathrm{vertex} =  \frac{\mathrm{fe} - \mathrm{min}(\mathrm{fe})}{ \mathrm{max}(\mathrm{fe}) - \mathrm{min}(\mathrm{fe})} \, .</math>
</div>
Via this last step, the minimum "vertex" energy across the entire domain was 0.0 while the maximum "vertex" energy was 1.0.


<div align="center">
<table border="1" cellpadding="8" align="center">
<tr>
  <th align="center" colspan="3">FORTRAN Program Documentation</th>
  <th align="center" colspan="3">Example Evaluations<p></p>(See also [[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|associated Table 1]])</th>
</tr>
<tr>
  <th align="center" colspan="1">Coord. Axis</th>
  <th align="center" colspan="1">Coord. Name</th>
  <th align="left">Associated Physical Quantity</th>
  <td align="center"><math>~\frac{\mu_e}{\mu_c} = 1</math></td>
  <td align="center" colspan="2"><math>~\frac{\mu_e}{\mu_c} = 0.305</math></td>
</tr>
<tr>
  <td align="center">x-axis</td>
  <td align="center">bsize</td>
  <td align="left"><math>~\ell_i \equiv \frac{\xi_i}{\sqrt{3}}</math></td>
  <td align="center"><math>~\frac{2.416}{\sqrt{3}} = 1.395</math></td>
  <td align="center"><math>~\frac{8.1938}{\sqrt{3}} = 4.7307</math></td>
  <td align="center"><math>~\frac{14.389}{\sqrt{3}} = 8.3076</math></td>
</tr>
<tr>
  <td align="center">y-axis</td>
  <td align="center">csize</td>
  <td align="left"><math>~\Chi \equiv \frac{\chi}{\chi_\mathrm{eq}}</math></td>
  <td align="center"><math>~1</math></td>
  <td align="center"><math>~1</math></td>
  <td align="center"><math>~1</math></td>
</tr>
<tr>
  <th align="center" colspan="3">Relevant Lines of Code</th>
  <td colspan="3" rowspan="3">&nbsp;</td>
</tr>
<tr>
  <td align="left" colspan="3">
<pre>
      eta = 3.0d0*muratio*bsize/(1.0d0+bsize**2)
      Gami = 1.0d0/eta-bsize
      frakL = (bsize**4-1.0d0)/bsize**2 + &
     &        DATAN(bsize)*((1.0d0+bsize**2)/bsize)**3
      frakK = Gami/eta + ((1.0d0+Gami**2)/eta)*(pii/2.0d0+DATAN(Gami))
      E0    = ((5.0d0*frakL) + (4.0d0*frakK)&
     &        - (3.0d0*frakL+12.0d0*frakK))/bsize**2+Efudge
          fescalar(j,k) = (csize**(-0.6d0)*(5.0d0*frakL)&
     &        + csize**(-3.0d0)*(4.0d0*frakK)&
     &   - (3.0d0*frakL+12.0d0*frakK)/csize)/bsize**2 + Efudge
          if(inorm.eq.1)fescalar(j,k)=fescalar(j,k)/DABS(E0) &
     &                  - E0/DABS(E0)
</pre>
  </td>
</tr>
<tr>
  <th align="center" colspan="1">Variable</th>
  <th align="center" colspan="1">Represents</th>
  <th align="left">Value calculated via the expression &hellip;</th>
</tr>
<tr>
  <td align="center">eta</td>
  <td align="center"><math>~\eta_i</math></td>
  <td align="left">
<math>~3 \biggl(\frac{\mu_e}{\mu_c}\biggr)\biggl[ \frac{\ell_i}{(1+\ell_i^2)} \biggr]</math>
  </td>
  <td align="center"><math>~1.421</math></td>
  <td align="center"><math>~0.1851</math></td>
  <td align="center"><math>~0.1086</math></td>
<tr>
  <td align="center">Gami</td>
  <td align="center"><math>~\Lambda_i</math></td>
  <td align="left"><math>~\frac{1}{\eta_i} - \ell_i</math></td>
  <td align="center"><math>~-0.691</math></td>
  <td align="center"><math>~0.6705</math></td>
  <td align="center"><math>~0.9033</math></td>
</tr>
<tr>
  <td align="center">frakL</td>
  <td align="center"><math>~\mathfrak{L}_i</math></td>
  <td align="left"><math>~
\frac{(\ell_i^4 - 1)}{\ell_i^2} + \biggl[ \frac{1+\ell_i^2}{\ell_i} \biggr]^3 \tan^{-1}\ell_i
</math></td>
  <td align="center"><math>~10.37</math></td>
  <td align="center"><math>~186.80</math></td>
  <td align="center"><math>~937.64</math></td>
</tr>
<tr>
  <td align="center">frakK</td>
  <td align="center"><math>~\mathfrak{K}_i</math></td>
  <td align="left"><math>~
\frac{\Lambda_i}{\eta_i} + \frac{(1 + \Lambda_i^2)}{\eta_i} \biggl[ \frac{\pi}{2} + \tan^{-1}\Lambda_i \biggr]
</math></td>
  <td align="center"><math>~0.518</math></td>
  <td align="center"><math>~20.544</math></td>
  <td align="center"><math>~46.882</math></td>
</tr>
<tr>
  <td align="center">&nbsp;</td>
  <td align="center"><math>~\frac{\mathfrak{L}_i}{\mathfrak{K}_i}</math></td>
  <td align="left">&nbsp;</td>
  <td align="center"><math>~20</math></td>
  <td align="center"><math>~9.093</math></td>
  <td align="center"><math>~20</math></td>
</tr>
<!-- OMIT
<tr>
  <td align="center">fescalar(x,y) - Efudge</td>
  <td align="center"><math>~\mathfrak{G}^*_{51}(\ell_i,\Chi)</math></td>
  <td align="left">
<math>~ \frac{1}{\ell_i^2} \biggl[
\Chi^{-3/5} (5 \mathfrak{L}_i)
+\Chi^{-3} (4\mathfrak{K}_i)
-\Chi^{-1} (3\mathfrak{L}_i  +12\mathfrak{K}_i ) \biggr]
</math>
</td>
  <td align="center"><math>~8.525</math></td>
  <td align="center">&hellip;</td>
  <td align="center">&hellip;</td>
</tr>
END OMIT-->
<tr>
  <td align="center">E0 - Efudge</td>
  <td align="center"><math>~\mathfrak{G}^*_{51}(\ell_i,\Chi=1)</math></td>
  <td align="left">
<math>~ 
\frac{1}{\ell_i^2} \biggl[ 5 \mathfrak{L}_i + 4\mathfrak{K}_i -  (3\mathfrak{L}_i  +12\mathfrak{K}_i ) \biggr]
= \frac{2(\mathfrak{L}_i - 4\mathfrak{K}_i)}{\ell_i^2} 
</math>
</td>
  <td align="center"><math>~8.525</math></td>
  <td align="center"><math>~9.3496</math></td>
  <td align="center"><math>~21.737</math></td>
</tr>
</table>
</div>

<div align="center">
<table border="0" cellpadding="8" width="880px">
<tr><td align="center"><font size="+1"><b>
Figure 2: &nbsp; Free-Energy Surface for <math>~(n_c,n_e) = (5,1)</math> and <math>~\frac{\mu_e}{\mu_c} = 0.305</math></b></font></td></tr>
<tr><td align="center">

<table border="1" cellpadding="5">
<tr>
<td align="center" bgcolor="#CCFFFF">
[[File:FEsurface01.png|center|400px|Free-Energy surface for 5_1 bipolytrope]]
</td>
</tr>
<tr>
  <td align="left" colspan="1">
''Left Panel:''  The free energy (vertical, blue axis) is plotted as a function of the radial interface location, <math>~\xi_i</math> (red axis), and the normalized configuration radius, <math>~\Chi \equiv \chi/\chi_\mathrm{eq}</math> (green axis).  ''Right Panel:'' Same as the left panel, but animated in order to highlight undulations of the surface.  The value of the free energy is indicated by color as well as by the height of the warped surface &#8212; red identifies the lowest depicted energies while blue identifies the highest depicted energies; these same colors have been projected down onto the <math>~z = 0</math> plane to present a two-dimensional, color-contour plot.  A multi-colored line segment drawn parallel to the <math>~\xi_i</math> axis at the value, <math>~\Chi = 1</math>, identifies the configuration's ''equilibrium'' radius for each value of the interface location.  Equilibrium configurations marked in white lie at the bottom of the principal free-energy "valley" and are stable, while configurations marked in blue lie at the top of a free-energy "hill," indicating that they are unstable; the red dot identifies the location along this equilibrium sequence where the transition from stable to unstable configurations occurs.
  </td>
</tr>
</table>

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

===BiPolytrope00===

<div align="center" id="FreeEnergy00">
<table border="1" cellpadding="5" align="center">
<tr>
<th align="center">
<font size="+1">Out-of-Equilibrium, Free-Energy Expression for BiPolytropes with ''Structural'' </font> <math>~(n_c, n_e) = (0, 0)</math>
</th>
</tr>
<tr><td align="center">

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

<tr>
  <td align="right">
<math>~\mathfrak{G}^*_{00} \equiv 5 \biggl(\frac{q}{\nu^2}\biggr) \chi_\mathrm{eq} 
\biggl[\frac{\mathfrak{G}_{00}}{E_\mathrm{norm}} \biggr]</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \frac{5}{2q^3} \biggl[
n_c A_2\Chi^{-3/n_c} + n_e B_2\Chi^{-3/n_e} - 3C_2\Chi^{-1} \biggr]
</math>
  </td>
</tr>
</table>
</td></tr>
</table>
</div>
where, 

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

<tr>
  <td align="right">
<math>~A_2 </math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>
\frac{2}{5}q^3  \biggl[1 + q^3 (f - 1-\mathfrak{F} ) \biggr]   \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~B_2 </math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>
\frac{2}{5} q^3  f - A_2 \, , 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~C_2 </math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
 <td align="left">
<math>
\frac{2}{5} q^3 f  \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~f</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
 <td align="left">
<math>
1+
\frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr)  \biggl(\frac{1}{q^2}-1 \biggr) +\biggl( \frac{\rho_e}{\rho_c}  
\biggr)^2 \biggl[ \frac{1}{q^5}-1 + \frac{5}{2} \biggl( 1-\frac{1}{q^2}\biggr)\biggr]  \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathfrak{F} </math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\frac{5}{2} \biggl( \frac{\rho_e}{\rho_c} \biggr) \frac{1}{q^5} \biggl[  (-2q^2 + 3q^3 - q^5) +
\frac{3}{5} \biggl( \frac{\rho_e}{\rho_c}\biggr) (-1 +5q^2 - 5q^3 + q^5) \biggr] \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\frac{\rho_e}{\rho_c} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{q^3(1-\nu)}{\nu(1-q^3)} \, .
</math>
  </td>
</tr>
</table>
</div>

The associated equilibrium radius is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\chi_\mathrm{eq}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl\{ \biggl(\frac{\pi}{3}\biggr) 2^{2-n_c} \nu^{n_c-1} q^{3-n_c} 
\biggl[1 + \frac{2}{5} q^3(f-1-\mathfrak{F})\biggr]^{n_c} \biggr\}^{1/(n_c-3)}
\, .
</math>
  </td>
</tr>
</table>
</div>

We have deduced that the system is unstable if,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{n_e}{3}\biggl[ \frac{3-n_e}{n_c-n_e} \biggr] </math>
  </td>
  <td align="center">
<math>~< </math>
  </td>
  <td align="left">
<math>~
\frac{A_2}{C_2}
= \frac{1}{f} \biggl[1 + q^3 (f - 1-\mathfrak{F} ) \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>

==Fortran Code==
This is the program that generated the free-energy data for the "five-one" bipolytrope that is displayed in the above, Figure 1 image/animation.

<pre>
      PROGRAM BiPolytrope
      real*8 pii
      real*8 bmin,bmax,cmin,cmax,db,dc
      real*8 c(200),b(200),chalf(199),bhalf(199)
      real*8 bsize,csize,emin,emax
      real*8 fepoint(200,200),fescalar(199,199)
      real*8 ell(200),ellhalf(199)
      real*8 muratio,eta,Gami,frakK,frakL
      real*8 eshift,ediff
      real   xx(200),yy(200),cell(199,199),vertex(200,200)
      real   valuemin,minlog,valufudge
      real*8 q,nu,chiEq,Enorm,E0,Efudge
      integer j,k,n,nmax,inorm
  101 format(4x,'bsize',7x,'csize',8x,'xi',10x,'A',12x,'B',12x,&
     &'fM',13x,'fW',11x,'fA',/)
! 102 format(1p8d12.4)
  103 format(2i5,1p3d14.6)
  201 format(5x,'valuemin = ',1pe15.5,//)
  205 format(5x,'For Cell-center ... emin, emax = ',1p2d14.6,/)
  206 format(5x,'For Cell-vertex ... emin, emax = ',1p2d14.6,/)
  701 format(5x,1p10d10.2)  
  700 format(8x,'xi',9x,'ell',8x,'eta',8x,'Lambda',5x,'frakK',&
     &       5x,'frakL',8x,'q',5x,'nu',5x,'chiEq',8x,'E0',/)
!!!!!!!!
!!!!!!!!
      inorm=1
      pii = 4.0d0*datan(1.0d0) 
      muratio = 1.0d0
      bsize = 0.0d0
      csize = 0.0d0
      Efudge = -10.0d0
      write(*,101)
!     write(*,102)bsize,csize,xival,coefA,coefB,fM,fW,fA
!!!!!!!!!!!
!
!  In this free-energy routine, c = X = chi/chi_eq and b = xi_i
!
!!!!!!!!!!!
      nmax = 200
      bmin = 1.0d0
      bmax = 3.0d0
      db   = (bmax-bmin)/dfloat(nmax-1)
      b(1) = bmin
      ell(1) = b(1)/dsqrt(3.0d0)
!  These values of cmin and cmax ensure that X=1 occurs at zone 70
      cmin = 0.469230769d0
      cmax = 2.00d0
      dc   = (cmax-cmin)/dfloat(nmax-1)
      c(1) = cmin
      do n=2,nmax
        b(n) = b(n-1)+db
        c(n) = c(n-1)+dc
        ell(n) = b(n)/dsqrt(3.0d0)
      enddo
      do n=1,nmax-1
        bhalf(n) = 0.5d0*(b(n)+b(n+1))
        chalf(n) = 0.5d0*(c(n)+c(n+1))
        ellhalf(n) = bhalf(n)/dsqrt(3.0d0)
      enddo
!
!  BEGIN LOOP to evaluate free energy (cell center)
!
      emin = 0.0d0
      emax = 0.0d0
      write(*,700)
      do j=1,nmax-1
        bsize = ellhalf(j)
      eta = 3.0d0*muratio*bsize/(1.0d0+bsize**2)
      Gami = 1.0d0/eta-bsize
      frakL = (bsize**4-1.0d0)/bsize**2 + &
     &        DATAN(bsize)*((1.0d0+bsize**2)/bsize)**3
      frakK = Gami/eta + ((1.0d0+Gami**2)/eta)*(pii/2.0d0+DATAN(Gami))
      q = 1.0d0/(1.0d0 + (0.5d0*pii+DATAN(Gami))/eta)
      nu = bsize*q/dsqrt(1.0d0+Gami**2)
      chiEq = dsqrt(pii/8.0d0)*(nu**2/(q*bsize**2))&
     &        *((1.0d0+bsize**2)/(3.0d0*bsize))**3
      Enorm = 16.0d0*(q/nu**2)*chiEq
      E0    = ((5.0d0*frakL) + (4.0d0*frakK)&
     &        - (3.0d0*frakL+12.0d0*frakK))/bsize**2+Efudge
      write(*,701)b(j),bsize,eta,Gami,frakK,frakL,q,nu,chiEq,E0
        do k=1,nmax-1
          csize=chalf(k)
          fescalar(j,k) = (csize**(-0.6d0)*(5.0d0*frakL)&
     &        + csize**(-3.0d0)*(4.0d0*frakK)&
     &   - (3.0d0*frakL+12.0d0*frakK)/csize)/bsize**2 + Efudge
          if(inorm.eq.1)fescalar(j,k)=fescalar(j,k)/DABS(E0) &
     &                  - E0/DABS(E0)
        if(fescalar(j,k).gt.0.5d0)fescalar(j,k)=0.5d0
          if(fescalar(j,k).lt.emin)emin=fescalar(j,k)
          if(fescalar(j,k).gt.emax)emax=fescalar(j,k)
!         write(*,103)j,k,bsize,csize,fescalar(j,k)
        enddo
      enddo
      write(*,205)emin,emax
!
!  BEGIN LOOP to evaluate free energy (cell vertex)
!
      emin = 0.0d0
      emax = 0.0d0
      do j=1,nmax
        bsize = ell(j)
      eta = 3.0d0*muratio*bsize/(1.0d0+bsize**2)
      Gami = 1.0d0/eta-bsize
      frakL = (bsize**4-1.0d0)/bsize**2 + &
     &        DATAN(bsize)*((1.0d0+bsize**2)/bsize)**3
      frakK = Gami/eta + ((1.0d0+Gami**2)/eta)*(pii/2.0d0+DATAN(Gami))
      q = 1.0d0/(1.0d0 + (0.5d0*pii+DATAN(Gami))/eta)
      nu = bsize*q/dsqrt(1.0d0+Gami**2)
      chiEq = dsqrt(pii/8.0d0)*(nu**2/(q*bsize**2))&
     &        *((1.0d0+bsize**2)/(3.0d0*bsize))**3
      Enorm = 16.0d0*(q/nu**2)*chiEq
      E0    = ((5.0d0*frakL) + (4.0d0*frakK)&
     &        - (3.0d0*frakL+12.0d0*frakK))/bsize**2 + Efudge
        do k=1,nmax
          csize=c(k)
          fepoint(j,k) = (csize**(-0.6d0)*(5.0d0*frakL)&
     &        + csize**(-3.0d0)*(4.0d0*frakK)&
     &   - (3.0d0*frakL+12.0d0*frakK)/csize)/bsize**2 + Efudge
          if(inorm.eq.1)fepoint(j,k)=fepoint(j,k)/DABS(E0) &
     &                  - E0/DABS(E0)
        if(fepoint(j,k).gt.0.5d0)fepoint(j,k)=0.5d0
          if(fepoint(j,k).lt.emin)emin=fepoint(j,k)
          if(fepoint(j,k).gt.emax)emax=fepoint(j,k)
!         write(*,103)j,k,bsize,csize,fepoint(j,k)
        enddo
      enddo
      write(*,206)emin,emax
!
!  Now fill single-precision arrays for plotting.
!
      do n=1,nmax
!       xx(n)=b(n)/b(nmax)
!       yy(n)=c(n)/c(nmax)
        xx(n)=b(n)-bmin
        yy(n)=c(n)-cmin
!       xx(n)=b(n)
!       yy(n)=c(n)
      enddo
      valuemin= -emin
      valufudge = 1.0d0/(emax-emin)
      do k=1,nmax
        do j=1,nmax
          vertex(j,k)=fepoint(j,k)+valuemin
          vertex(j,k)=vertex(j,k)*valufudge
        enddo
      enddo
      do k=1,nmax-1
        do j=1,nmax-1
          cell(j,k)=fescalar(j,k)+valuemin
          cell(j,k)=cell(j,k)*valufudge
        enddo
      enddo
      call XMLwriter01(nmax,xx,yy,cell,vertex)
      stop
      END PROGRAM BiPolytrope
      Subroutine XMLwriter01(imax,x,y,cell_scalar,point_scalar)
      real x(200),y(200),z(1)
      real cell_scalar(199,199),point_scalar(200,200)
      integer imax
      integer extentX,extentY,extentZ
      integer ix0,iy0,iz0
      integer norm(200,3)
!     imax=200
      ix0=0
      iy0=0
      iz0=0
      extentX=imax-1
      extentY=imax-1
      extentZ=0
      z(1) = 0.0
! Set normal vector 1D array
      do i=1,imax
        norm(i,1)=0
        norm(i,2)=0
        norm(i,3)=1
      enddo
  201 format('<?xml version="1.0"?>')
  202 format('<VTKFile type="RectilinearGrid" version="0.1" byte_order="LittleEndian">')
  302 format('</VTKFile>')
  203 format(2x,'<RectilinearGrid WholeExtent="',6I4,'">')
  303 format(2x,'</RectilinearGrid>')
  204 format(4x,'<Piece Extent="',6I4,'">')
  304 format(4x,'</Piece>')
  205 format(6x,'<CellData Scalars="cell_scalars" Normals="magnify">')
  305 format(6x,'</CellData>')
  206 format(8x,'<DataArray type="Float32" Name="magnify" NumberOfComponents="3" format="ascii">')
  207 format(8x,'<DataArray type="Float32" Name="cell_scalars" format="ascii">')
  399 format(8x,'</DataArray>')
  208 format(6x,'<PointData Scalars="colorful" Normals="direction">')
  308 format(6x,'</PointData>')
  209 format(8x,'<DataArray type="Float32" Name="colorful" format="ascii">')
  210 format(6x,'<Coordinates>')
  310 format(6x,'</Coordinates>')
  211 format(8x,'<DataArray type="Float32" format="ascii" RangeMin="0" RangeMax="5">')
  212 format(8x,'<DataArray type="Float32" format="ascii">')
  213 format(8x,'<DataArray type="Float32" Name="direction" NumberOfComponents="3" format="ascii">')
  501 format(10f9.5)
  502 format(10f9.5)
  503 format(5x,9(1x,3I2))
  504 format(10f9.5)
  505 format(5x,10(1x,3I2))
!!!!!
!
!  Begin writing out XML tags.
!
!!!!!
      write(*,201)                              !<?xml
      write(*,202)                              !VTKFile
        write(*,203)ix0,extentX,iy0,extentY,iz0,extentZ        !  RectilinearGrid
          write(*,204)ix0,extentX,iy0,extentY,iz0,extentZ      !    Piece
            write(*,205)                        !      CellData
              write(*,207)                      !        DataArray(cell_scalars)
      do j=1,imax-1
      write(*,501)(cell_scalar(i,j),i=1,imax-1)
      enddo
              write(*,399)                      !        /DataArray
              write(*,206)                      !        DataArray(cell_scalars)
      do j=1,imax-1
      write(*,503)(norm(i,1),norm(i,2),norm(i,3),i=1,imax-1)
      enddo
              write(*,399)                      !        /DataArray
            write(*,305)                        !      /CellData
            write(*,208)                        !      PointData
              write(*,209)                      !        DataArray(points)
      write(*,502)((point_scalar(i,j),i=1,imax),j=1,imax)
              write(*,399)                      !        /DataArray
              write(*,213)                      !        DataArray(cell_scalars)
      do j=1,imax
      write(*,505)(norm(i,1),norm(i,2),norm(i,3),i=1,imax)
      enddo
              write(*,399)                      !        /DataArray
            write(*,308)                        !      /PointData
            write(*,210)                        !      Coordinates
              write(*,212)                      !        DataArray(x-direction)
        write(*,504)(x(i),i=1,imax)
              write(*,399)                      !        /DataArray
              write(*,212)                      !        DataArray(y-direction)
        write(*,504)(y(i),i=1,imax)
              write(*,399)                      !        /DataArray
              write(*,212)                      !        DataArray(z-direction)
        write(*,504)z(1)
              write(*,399)                      !        /DataArray
            write(*,310)                        !      /Coordinates
          write(*,304)                          !    /Piece
        write(*,303)                            !  /RectilinearGrid
      write(*,302)                              !/VTKFile
      return
      end
</pre>



=Nonstandard Examination=
In our introductory remarks, [[#Free-Energy_of_Bipolytropes|above]], we said the warped free-energy surface drapes across a five-dimensional parameter "plane" such that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{G}_{51}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\mathfrak{G}(R, K_c, M_\mathrm{tot}, q, \nu) \, .</math>
  </td>
</tr>
</table>
</div>

From a more pragmatic point of view, we should have said that the "five-one" free-energy surface drapes across a five-dimensional parameter "plane" such that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{G}_{51}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\mathfrak{G}(R, K_c, M_\mathrm{tot}, \ell_i, \tfrac{\mu_e}{\mu_c}) \, .</math>
  </td>
</tr>
</table>
</div>
In our initial, standard examination of the structure of this warped free-energy surface, we held  three parameters fixed &#8212; namely, <math>~K_c, M_\mathrm{tot}, \tfrac{\mu_e}{\mu_c}</math> &#8212; and plotted <math>~\mathfrak{G}_{51}(\ell_i, \Chi\equiv R/R_\mathrm{eq})</math>.  Now, let's fix <math>~\Chi = 1</math> and plot <math>~\mathfrak{G}_{51}(\ell_i, \tfrac{\mu_e}{\mu_c})</math>.  The following plot shows how a portion of the <math>~(\ell_i, \mu_e/\mu_c)</math> grid maps onto the traditional <math>~(q, \nu)</math> plane.  The numerical labels identify lines of constant <math>~\xi_i = \sqrt{3}\ell_i</math> &#8212; 7 (light green), 9 (pink), and 12 (orange) &#8212; and lines of constant <math>~\mu_e/\mu_c</math> &#8212; 0.330 (purple), 0.315 (dark green), and 0.305 (white/blue).

[[File:GridOnNuQplot.png|center|500px|xi-ell grid drawn on q-nu grid]]

=See Also=
In <font color="red">October 2023</font>, this very long chapter was subdivided in order to more effectively accommodate edits.  Here is a list of the resulting set of shorter chapters:
<ol>
  <li>[[SSC/FreeEnergy/PolytropesEmbedded/Pt1|Free-Energy Synopsis]]</li>
  <li>[[SSC/FreeEnergy/PolytropesEmbedded/Pt2|Free-Energy of Truncated Polytropes]]</li>
  <li>Free-Energy of BiPolytropes
    <ul>
<li>[[SSC/FreeEnergy/PolytropesEmbedded/Pt3A|Focus on Five-One Free-Energy Expression]]</li>
<li>[[SSC/FreeEnergy/PolytropesEmbedded/Pt3B|Focus on Zero-Zero Free-Energy Expression]]</li>
<li>[[SSC/FreeEnergy/PolytropesEmbedded/Pt3C|Overview]]</li>
    </ul>
  </li>
</ol>


{{ SGFfooter }}