Editing SSC/Stability/BiPolytropes

__FORCETOC__ 
<!-- __NOTOC__ will force TOC off -->

=Organizational Index=
<table border="1" align="center" width="100%" colspan="8">
<tr>
  <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Stability/BiPolytropes|Part I: &nbsp; The Search]]<br />&nbsp;</td>
  <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Stability/BiPolytropes/Pt2|Part II:&nbsp; Review of MF85b]]<br />&nbsp;</td>
  <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Stability/BiPolytropes/Pt3|III:&nbsp; (5,1) Radial Oscillations]]<br />&nbsp;</td>
  <td align="center" bgcolor="lightblue"><br />[[SSC/Stability/BiPolytropes/Pt4|IV:&nbsp; Reconciliation]]<br />&nbsp;</td>
</tr>
<tr>
  <td align="left" width="100%" colspan="4">
&nbsp; &nbsp; These four chapters, labeled Parts I - IV, are segments of the much longer chapter titled, [[SSC/Stability/BiPolytropes/PlannedApproach|SSC/Stability/BiPolytropes/PlannedApproach]].  An [[SSC/Stability/BiPolytropes/Index|accompanying organizational index]] has helped us write this chapter succinctly.
  </td>
</tr>
</table>

=In Search of Marginally Unstable (n<sub>c</sub>,n<sub>e</sub>) = (5,1) Bipolytropes=

<!-- Our aim is to determine whether or not there is a relationship between (1) equilibrium models at turning points along bipolytrope sequences and (2) bipolytropic models that are marginally (dynamically) unstable toward collapse (or dynamical expansion).-->


==Overview==
<!--
<table border="0" cellpadding="8" align="right">
<tr>
  <th align="center">[[File:DataFileButton02.png|right|60px|file = Dropbox/WorkFolder/Wiki edits/EmbeddedPolytropes/CombinedSequences.xlsx --- worksheet = EqSeqCombined]]Figure 1: &nbsp;Equilibrium Sequences<br />of Pressure-Truncated Polytropes
  </th>
</tr>
<tr>
  <td align="center" colspan="1">
[[File:MassVsRadiusCombined02.png|300px|Equilibrium sequences of Pressure-Truncated Polytropes]]
  </td>
</tr>
</table>
-->
<table border="0" cellpadding="8" align="right">
<tr>
  <th align="center">[[File:DataFileButton02.png|right|60px|file = Dropbox/WorkFolder/Wiki edits/EmbeddedPolytropes/CombinedSequences.xlsx --- worksheet = EqSeqCombined2]]Figure 1: &nbsp;Equilibrium Sequences<br />of Pressure-Truncated Polytropes
  </th>
</tr>
<tr>
  <td align="center" colspan="1">
[[File:DFBsequenceB.png|300px|Equilibrium sequences of Pressure-Truncated Polytropes]]
  </td>
</tr>
</table>
We expect the  content of this chapter &#8212; which examines the relative stability of bipolytropes &#8212; to parallel in many ways the content of an [[SSC/Stability/InstabilityOnsetOverview#Marginally_Unstable_Pressure-Truncated_Gas_Clouds|accompanying chapter in which we have successfully analyzed the relative stability of pressure-truncated polytopes]].  Figure 1, shown here on the right, has been copied from [[SSC/Structure/PolytropesEmbedded#Additional.2C_Numerically_Constructed_Polytropic_Configurations|a closely related discussion]].  The curves show the mass-radius relationship for pressure-truncated model sequences having a variety of polytropic indexes, as labeled, over the range <math>1 \le n \le 6</math>.  ([[SSC/Stability/InstabilityOnsetOverview#Turning_Points_along_Sequences_of_Pressure-Truncated_Polytropes|Another version of this figure]] includes the isothermal sequence.)  On each sequence for which <math>n \ge 3</math>, the green filled circle identifies the model with the largest mass.  We have shown ''analytically'' that the oscillation frequency of the fundamental-mode of radial oscillation is precisely zero<sup>&dagger;</sup> for each one of these maximum-mass models.  As a consequence, we know that each green circular marker identifies the point along its associated sequence that separates dynamically stable (larger radii) from dynamically unstable (smaller radii) models.

<sup>&dagger;</sup>In each case, the fundamental-mode oscillation frequency is precisely zero if, and only if, the adiabatic index governing expansions/contractions is related to the underlying ''structural'' polytropic index via the relation, <math>\gamma_g = (n + 1)/n</math>, and if a constant surface-pressure boundary condition is imposed.

<font color="red">'''Key Realization:'''</font>  ''Along sequences of pressure-truncated polytropes, the maximum-mass models identify precisely where the onset of dynamical instability occurs.''

----


In another [[SSC/Structure/BiPolytropes/Analytic51#BiPolytrope_with_nc_.3D_5_and_ne_.3D_1|accompanying chapter]], we have used purely analytic techniques to construct equilibrium sequences of spherically symmetric bipolytropes that have, <math>(n_c,n_e) = (5,1)</math>.   For a given choice of <math>~\mu_e/\mu_c</math> &#8212; the ratio of the mean-molecular weight of envelope material to the mean-molecular weight of material in the core &#8212; a physically relevant sequence of models can be constructed by steadily increasing the value of the dimensionless radius at the core/envelope interface, <math>\xi_i</math>, from zero to infinity.  Figure 2, whose content is essentially the same as [[SSC/Structure/BiPolytropes/Analytic51#Model_Sequences|Figure 1 of this separate chapter]], 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 seven different values of <math>\mu_e/\mu_c</math>, as labeled: 1 (black), &frac12; (dark blue), 0.345 (brown), &#8531; (dark green), 0.316943 (purple), 0.309 (orange),  and &frac14; (light blue).  

When modeling bipolytropes, the default 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 Figure 2 sequences that have the largest mean-molecular weight ratios:  <math>\mu_e/\mu_c</math> = 1 and &frac12;.  But the behavior is different along the other five illustrated sequences.  For sufficiently large <math>\xi_i</math>, the relative radius of the core begins to decrease; along each sequence, a solid purple circular marker identifies the location of this ''turning point'' in radius.  Furthermore, along sequences for which <math>\mu_e/\mu_c < \tfrac{1}{3}</math>, eventually the fractional mass of the core reaches a maximum and, thereafter, decreases even as the value of <math>\xi_i</math> continues to increase.  In Figure 2, a solid green circular marker identifies the location of this ''maximum mass turning point'' along each of these sequences; the analytically determined values of <math>\xi_i, q </math> and <math>\nu</math> that are associated with each of these ''turning points'' are provided in the table adjacent to Figure 2.  (Additional properties of these equilibrium sequences are discussed in [[SSC/FreeEnergy/PolytropesEmbedded#Behavior_of_Equilibrium_Sequence|yet another accompanying chapter]].)  

<span id="PrincipalQ">
<font color="red">'''The principal question is:'''</font>  ''Along bipolytropic sequences, are maximum-mass models associated with the onset of dynamical instabilities?''</span>


<table border="1" align="center" cellpadding="3">
<tr>
  <td align="center" rowspan="1">
'''Figure 2:  Equilibrium Sequences of Bipolytropes''' <br /><p>
'''with <math>(n_c,n_e) = (5,1)</math> and Various <math>\mu_e/\mu_c</math>'''
  </td>
  <td align="center" colspan="4">
<b>Analytically Determined Parameters<sup>&dagger;</sup><br />for Models that have the Maximum Fractional Core Mass<br />(solid green circular markers)<br />Along Various Equilibrium Sequences
  </td>
</tr>
<tr>
  <td align="center" rowspan="8">
[[File:TurningPoints51Bipolytropes.png|450px|Extrema along Various Equilibrium Sequences]]
  </td>
  <td align="center">
<math>\frac{\mu_e}{\mu_c}</math>
  </td>
  <td align="center">
<math>\xi_i</math>
  </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>
</tr>

<tr>
  <td align="center">
<math>\frac{1}{3}</math>
  </td>
  <td align="center">
<math>\infty</math> </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.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.068652714  </td>
  <td align="right">
0.382383875  </td>
</tr>

<tr>
  <td align="center">
0.31
  </td>
  <td align="right">
9.014959766  </td>
  <td align="right">
0.0755022550  </td>
  <td align="right">
0.3372170064  </td>
</tr>

<tr>
  <td align="center">
0.3090
  </td>
  <td align="right">
8.8301772  </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.084824137  </td>
  <td align="right">
0.139370157  </td>
</tr>

<tr>
  <td align="left" colspan="4">
<sup>&dagger;</sup>Additional model parameters [[SSC/Stability/BiPolytropes/51Models#Structure|can be found here]].
  </td>
</tr>
</table>

==Planned Approach(es)==

In an effort to answer the '''<font color="red">principal question</font>''' posed above, we have pursued each stability-analysis approach described in the introductory section of {{ B-KB74full }}.

<table border="0" align="center" cellpadding="8" width="95%">
<tr><td align="left" rowspan="1">
<font color="darkgreen"><span id="BKB74pt1">"Three different approaches are used in the study of the hydrodynamical stability of stars</span> and other gravitating objects &hellip;" &nbsp; 
<ul><li>"The first approach is based on the use of the equations of small oscillations.  In that case the problem is reduced to a search for the solution of the boundary-value problem of the Stourme-Liuville type for the linearised system of equations of small oscillations.  The solutions consist of a set of eigenfrequencies and eigenfunctions."</font> </li>

<li>Second, one can derive <font color="darkgreen">"a variational principle from the equations of small oscillations &hellip;</font>  <!-- This principle replaces the straightforward solution of these equations:</font>  In the context of rotating Newtonian systems, see, for example, [http://adsabs.harvard.edu/abs/1964ApJ...140.1045C Clement (1964)], [http://adsabs.harvard.edu/abs/1968ApJ...152..267C Chandrasekhar &amp; Lebovitz (1968)], [http://adsabs.harvard.edu/abs/1967MNRAS.136..293L Lynden-Bell and Ostriker (1967)], or [http://adsabs.harvard.edu/abs/1972ApJS...24..319S Schutz (1972)]. -->  <font color="darkgreen">With the aid of the variational principle, the problem is reduced to the search of the best trial functions; this leads to approximate eigenvalues of oscillations.  In spite of the simplifications introduced by the use of the variational principle and by not solving the equations of motion exactly, the problem still remains complicated &hellip;"</font></li>

<li>The third approach is what we usually refer to as a free-energy &#8212; and associated virial theorem &#8212; analysis.  <font color="darkgreen">"When this method is used, it is not necessary to use the equations of small oscillations but, instead, the functional expression for the total energy of the momentarily stationary (but not necessarily in equilibrium) star is sufficient.  The condition that the first variation of the energy vanishes, determines the state of equilibrium of the star and the positiveness of a second variation indicates stability."</font> </li>
</ul>

<span id="BKB74pt2"><font color="darkgreen">"If one wants</span> to know from a stability analysis the answer to only one question &#8212; whether the model is stable or not &#8212; then the most straightforward procedure is to use the third, static method &hellip;  For the application of this method, one needs to construct only equilibrium, stationary models, with no further calculation."</font>

<font color="darkgreen">"Generally the static analysis gives no information about the shape of the modes of oscillation, but, in the vicinity of critical points, where instability sets in, this method makes it possible to find the eigenfunction of the mode which becomes unstable at the critical point."</font>
</td></tr>
<tr>
<td align="right">
— Drawn from pp. 391 - 392 of {{ B-KB74 }} 
</td>
</tr>
</table>

==Supplemental Chapters==
<ol>
  <li>Contains [[SSC/Stability/InstabilityOnsetOverview#Displacement_Functions_Summary|Displacement Functions Summary]]</li>
  <li>[[SSC/Stability/BiPolytropes/PlannedApproach|Earlier Planned Approach]]</li>
  <li>[[SSC/Stability/BiPolytropes/HeadScratching|Headscratching]]</li>
  <li>[[SSC/Stability/BiPolytropes/SuccinctDiscussion|Succinct Discussions]]</li>
  <li>[[SSC/Stability/BiPolytropes/51Models|51Models]]</li>
</ol>

=Equilibrium Models=

The original text contained in this section has been commented out.
<!--
The equilibrium sequences displayed in Figure 1 &#8212; for pressure-truncated polytropes &#8212; and in Figure 2 &#8212; for <math>(n_c, n_e) = (5, 1)</math> bipolytropes &#8212; were constructed by solving the second-order ODE that governs detailed force-balance throughout spherically symmetric configurations that obey polytropic equations of state.  In Figure 1, the sequences labeled <math>n = 1</math> and <math>n = 5</math> can be specified entirely through analytical expressions; these are well-known analytic solutions to the Lane-Emden equation.  In Figure 2, ''all seven equilibrium sequences'' are specifiable analytically; the location(s) of ''turning points'' along each sequence (when they exist) are also completely specifiable analytically.

In this chapter we examine the relative stability of various equilibrium <math>(n_c, n_e) = (5, 1)</math> bipolytropes.  A model is uniquely selected once we specify the radial location, <math>\xi_i</math>, of the core/envelope interface, and specify the ratio <math>\mu_e/\mu_c \le 1</math> of the mean-molecular weights of the gas in the envelope to the gas in the core.  Once this pair of parameter values has been specified, the model's radius <math>R_\mathrm{eq}</math> and total mass <math>M_\mathrm{tot}</math> are known analytically, as is the relative size of the core &#8212; both in terms of its fractional radius, <math>q \equiv r_i/R_\mathrm{eq}</math>, and its fractional mass, <math>\nu \equiv m_\mathrm{core}/M_\mathrm{tot}</math>.  Knowledge of the <math>(q, \nu)</math> parameter pair precisely identifies the selected model's location in the Figure 2 diagram.

In line with the '''<font color="red">principal question</font>''' stated above, we are particularly interested in examining the relative stability of models that are associated with the "maximum-mass turning points" along bipolytropic sequences for which <math>\mu_e/\mu_c \le \tfrac{1}{3}</math> &#8212; see the solid green circular markers in Figure 2.
-->

=Third Approach: Free-Energy Analysis=

==Old Material==
<table border="0" align="right" cellpadding="5"><tr><td align="center">[[File:VirialStability01.png|350px|Virial Stability01]]</td></tr></table>
The following set of menu tiles include links to chapters where this approach has been applied to: (a) uniform-density configurations, (b) pressure-truncated isothermal spheres, (c) an isolated n = 3 polytrope, (d) pressure-truncated n = 5 configurations, and (e) bipolytropes having <math>(n_c, n_e) = (1, 5)</math>.

One menu tile, below, links to a chapter in which an analytic (''exact'') demonstration of the variational principle's utility is provided in the context pressure-truncated n = 5 polytropes.

Ideally we would like to answer the just-stated "principal question" using purely analytic techniques.  But, to date, we have been unable to fully address the relevant issues analytically, even in what would be expected to be the simplest case: &nbsp; [[SSC/Stability/BiPolytrope00#Radial_Oscillations_of_a_Zero-Zero_Bipolytrope|bipolytropic models that have <math>(n_c,n_e) = (0, 0)</math>]].  Instead, we will streamline the investigation a bit and proceed &#8212; at least initially &#8212; using a blend of techniques.  We will investigate the relative stability of bipolytropic models having <math>(n_c,n_e) = (5,1) </math> whose ''equilibrium structures'' are completely defined analytically; then the eigenvectors describing radial modes of oscillation will be determined, one at a time, by solving the relevant LAWE(s) numerically.  We are optimistic that this can be successfully accomplished because we have had experience numerically integrating the LAWE that governs the oscillation of:
* [[SSC/Stability/n3PolytropeLAWE#Radial_Oscillations_of_n_.3D_3_Polytropic_Spheres|Isolated n = 3 polytropes]] &#8212; including a quantitative comparison against the published work of {{ Schwarzschild41full }};
* [[SSC/Stability/Isothermal#Radial_Oscillations_of_Pressure-Truncated_Isothermal_Spheres|Pressure-truncated isothermal spheres]] &#8212; including a quantitative comparison against the published analysis of {{ TVH74full }}; and
* [[SSC/Stability/n5PolytropeLAWE#Radial_Oscillations_of_n_.3D_5_Polytropic_Spheres|Pressure-truncated n = 5 polytropes]].

A key reference throughout this investigation will be the paper by {{ MF85bfull }}.  They studied ''Radial Pulsations and Vibrational Stability of a Sequence of Two Zone Polytropic Stellar Models.''  Specifically, their underlying equilibrium models were bipolytropes that have <math>(n_c,n_e) = (1, 5)</math>.  In an [[SSC/Structure/BiPolytropes/Analytic15#BiPolytrope_with_nc_.3D_1_and_ne_.3D_5|accompanying chapter]], we describe in detail how {{ MF85b }} obtained these equilibrium bipolytropic structures and detail some of their equilibrium properties.

Here are the steps we initially plan to take:

* Governing LAWEs:
** Identify the relevant LAWEs that govern the behavior of radial oscillations in the <math>~n_c = 5</math> core and, separately, in the <math>~n_e = 1</math> envelope.  Check these LAWE specifications against the published work of {{ MF85b }}.
** Determine the matching conditions that must be satisfied across the core/envelope interface.  Be sure to take into account the critical interface ''jump'' conditions spelled out by {{ LW58full }}, as we have already discussed in the context of an [[SSC/Stability/BiPolytrope00#Radial_Oscillations_of_a_Zero-Zero_Bipolytrope|analysis of radial oscillations in zero-zero bipolytropes]].
* Determine what surface boundary condition should be imposed on physically relevant LAWE solutions, i.e., on the physically relevant radial-oscillation eigenvectors.
* Initial Analysis:
** Choose a maximum-mass model along the bipolytropic sequence that has, for example, <math>\mu_e/\mu_c = 1/4</math>.  Hopefully, we will be able to identify precisely (analytically) where this maximum-mass model lies along the sequence.  <font color="red">'''Yes!'''</font>  Our [[SSC/Structure/BiPolytropes/Analytic51#Limiting_Mass|earlier analysis]] does provide an analytic prescription of the model that sits at the maximum-mass location along the chosen sequence.
** Solve the relevant eigenvalue problem for this specific model, initially for <math>(\gamma_c, \gamma_e) = (6/5, 2)</math> and initially for the fundamental mode of oscillation.

==New Material==
Drawing from our [[SSC/Structure/BiPolytropes/FreeEnergy51#Free_Energy_of_BiPolytrope_with|accompanying detailed discussion]] &#8212; see also an [[SSC/Stability/BiPolytropes/PlannedApproach#Free_Energy_and_Its_Derivatives|accompanying summary]] &#8212; the normalized free-energy associated with each of our spherically symmetric, bipolytropic configurations is given by the expression,

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

<tr>
  <td align="right">
<math>\mathfrak{g} \equiv \frac{\mathfrak{G}}{E_\mathrm{norm}}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
-\biggl(\frac{\chi_\mathrm{eq}}{\chi}\biggr)\mathfrak{w} 
+  \biggl[ \frac{2}{3(\gamma_c - 1)} \biggr] \biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{3 - 3\gamma_c}\mathfrak{s}_\mathrm{core}  
+  \biggl[ \frac{2}{3(\gamma_e - 1)} \biggr] \biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{3 - 3\gamma_e}\mathfrak{s}_\mathrm{env}  
\, ,
</math>
  </td>
</tr>
</table>

where: <math>\chi_\mathrm{eq}</math> is the dimensionless radius of the configuration when it is in equilibrium; <math>\mathfrak{s}_\mathrm{core}</math> and <math>\mathfrak{s}_\mathrm{env}</math> are the appropriately normalized thermal energy content of the <math>n = 5</math> core and of the <math>n = 1</math> envelope, respectively, in the configuration's equilibrium state; and <math>\mathfrak{w}</math> is the absolute value of the normalized total gravitational potential energy of the equilibrium configuration.  For every <math>(n_c, n_e) = (5, 1)</math> bipolytropic model, the values of these four terms can be obtained via [[SSC/Stability/BiPolytropes/PlannedApproach#Free_Energy_and_Its_Derivatives|our derived analytic expressions]]. And the value of <math>\mathfrak{g}</math> in an equilibrium state is obtained by setting the configuration's dimensionless radius, <math>\chi</math>, equal to <math>\chi_\mathrm{eq}</math>.  

This expression for <math>\mathfrak{g}</math> has been written in such a way that we can readily assess how the free energy varies while the configuration undergoes homologous (<math>\mathfrak{s}_\mathrm{core}, \mathfrak{s}_\mathrm{env}, \mathfrak{w}</math> all held fixed) expansion/contraction <math>\chi \gtrless \chi_\mathrm{eq}</math> about its equilibrium state.  Specifically the ''first variation'' is,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\frac{d\mathfrak{g}}{d\chi} </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
+\biggl(\frac{\chi_\mathrm{eq}}{\chi^2}\biggr)\mathfrak{w} 
~-~\frac{2}{\chi_\mathrm{eq}} \biggl[\biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{2 - 3\gamma_c}\mathfrak{s}_\mathrm{core}  
~+~\biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{2 - 3\gamma_e}\mathfrak{s}_\mathrm{env}  \biggr]  \, ;
</math>
  </td>
</tr>
</table>
and the ''second variation'' is,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\frac{d^2\mathfrak{g}}{d\chi^2} </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
-2\biggl(\frac{\chi_\mathrm{eq}}{\chi^3}\biggr)\mathfrak{w} 
~-~\frac{2}{\chi_\mathrm{eq}^2} \biggl[(2 - 3\gamma_c) \biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{1 - 3\gamma_c}\mathfrak{s}_\mathrm{core}  
~+~(2 - 3\gamma_e)\biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{1 - 3\gamma_e}\mathfrak{s}_\mathrm{env}  \biggr] \ .
</math>
  </td>
</tr>
</table>
The condition that the first variation of the energy vanishes, determines the state of equilibrium of the star and the positiveness of a second variation indicates stability.

=Virial Analysis=

In an [[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|accompanying chapter]] we have examined the viral stability of <math>~(n_c, n_e) = (5, 1)</math> bipolytropes &#8212; see also an [[SSC/Structure/BiPolytropes/FreeEnergy51#Free_Energy_of_BiPolytrope_with|associated discussion of the free-energy of these configurations]].  We present, here, an overview of this separate discussion in an effort to provide a broader perspective to the focus of this current chapter.

==Various Energy Expressions==
Drawing from our [[SSC/Structure/BiPolytropes/FreeEnergy51#Free_Energy_of_BiPolytrope_with|accompanying detailed discussion]], we have,
<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>~W_\mathrm{grav}\biggr|_\mathrm{core} + \mathfrak{S}_A\biggr|_\mathrm{core} + W_\mathrm{grav}\biggr|_\mathrm{env} + \mathfrak{S}_A\biggr|_\mathrm{env} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~W_\mathrm{grav}\biggr|_\mathrm{core} +  \biggl[ \frac{2}{3(\gamma_c - 1)} \biggr] S_\mathrm{core} + W_\mathrm{grav}\biggr|_\mathrm{env} + \biggl[ \frac{2}{3(\gamma_e - 1)} \biggr] S_\mathrm{env}  \, .
</math>
  </td>
</tr>
</table>
</div>

Recognizing that,
<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>~\frac{r_i}{R} = \frac{\eta_i}{\eta_s} \, ,</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>

<tr>
  <td align="right">
<math>~\ell</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{\xi}{\sqrt{3}}  \, ,</math>
  </td>
</tr>

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

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

<tr>
  <td align="right">
<math>~- \chi \biggl[ \frac{W_\mathrm{grav}}{E_\mathrm{norm}} \biggr]_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
 <td align="left">
<math>
\biggl[\biggl(\frac{3}{2^4}\biggr) \biggl( \frac{q}{\ell_i}\biggr)^{5}\biggl(\frac{\nu}{q^3} \biggr)^2 \biggl( 1 + \ell_i^2 \biggr)^{3}  \biggr]_\mathrm{eq}
\biggl[
\ell_i (\ell_i^4 - \frac{8}{3}\ell_i^2 - 1) (\ell_i^2 +1)^{-3} + \tan^{-1}\ell_i
\biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \biggl[ \frac{W_\mathrm{grav}}{E_\mathrm{norm}} \biggr]_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
 <td align="left">
<math>- \frac{\chi_\mathrm{eq}}{\chi} 
\biggl( \frac{3^8}{2^5 \pi}\biggr)^{1 / 2}
\biggl[
\ell_i (\ell_i^4 - \frac{8}{3}\ell_i^2 - 1) (\ell_i^2 +1)^{-3} + \tan^{-1}\ell_i
\biggr] \, .
</math>
  </td>
</tr>
</table>
Next, given that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~b_\eta = \eta_s</math>
  </td>
  <td align="center">
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
and
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
 </td>
  <td align="left">
<math>~b_\eta q = \eta_i = 3\biggl( \frac{\mu_e}{\mu_c} \biggr) \ell_i (1 + \ell_i^2)^{-1}  \, ,</math>
  </td>
</tr>
</table>
</div>
we have,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\biggl[ \frac{W_\mathrm{grav}}{E_\mathrm{norm}} \biggr]_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
 <td align="left">
<math>
- \frac{\chi_\mathrm{eq}}{\chi} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3}  \biggl( \frac{1}{2^3\pi} \biggr)^{1/2} A^2 
\biggl[6 b_\eta x  - 3\sin[2(b_\eta x-B)] - 4b_\eta x \sin^2(b_\eta x - B)  \biggr]_q^1 
</math>
  </td>
</tr>

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

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

<tr>
  <td align="right">
<math>~q a_\xi^{1 / 2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\xi_i}{\sqrt{3}} \, ,</math>
  </td>
</tr>
</table>
we have,
<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{2}{3({\gamma_c}-1)} \biggl( \frac{\chi}{\chi_\mathrm{eq}} \biggr)^{3-3\gamma_c} 
\biggl\{\biggl( \frac{3^8}{2^7\pi} \biggr)^{1/2} \biggl[ 
\tan^{-1}[a_\xi^{1/2}q] - a_\xi^{1/2}q ~\frac{(1 - a_\xi q^2)}{(1 + a_\xi q^2)^2} \biggr]
\biggr\} 
</math>
  </td>
</tr>

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

Finally, then, we can write,

<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{env}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
  <td align="left">
<math>
\frac{2}{3({\gamma_e}-1)} \biggl( \frac{\chi}{\chi_\mathrm{eq}} \biggr)^{3-3\gamma_e} \biggl\{
\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2  \biggl( \frac{3^2}{2^5\pi} \biggr)^{1/2}  \biggl[ 2b_\eta x -\sin[2(b_\eta x - B)] \biggr]_q^1  \biggr\} 
</math>
  </td>
</tr>

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

==Free Energy and Its Derivatives==

We can now rewrite the free-energy expression in the form,

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

<tr>
  <td align="right">
<math>~\mathfrak{g} \equiv \frac{\mathfrak{G}}{E_\mathrm{norm}}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~-\biggl(\frac{\chi_\mathrm{eq}}{\chi}\biggr)\mathfrak{w} 
+  \biggl[ \frac{2}{3(\gamma_c - 1)} \biggr] \biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{3 - 3\gamma_c}\mathfrak{s}_\mathrm{core}  
+  \biggl[ \frac{2}{3(\gamma_e - 1)} \biggr] \biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{3 - 3\gamma_e}\mathfrak{s}_\mathrm{env}  
\, ,
</math>
  </td>
</tr>
</table>
</div>
where,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{s}_\mathrm{core}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{3^8}{2^7\pi} \biggr)^{1/2} \biggl[ 
\tan^{-1}\biggl( \frac{\xi_i}{\sqrt{3}}\biggr) - \biggl( \frac{\xi_i}{\sqrt{3}}\biggr)\frac{(1 - \xi_i^2/3)}{(1 + \xi_i^2/3)^2} \biggr]\, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathfrak{s}_\mathrm{env}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\bigg( \frac{\mu_e}{\mu_c} \biggr)^{-3} A^2  \biggl( \frac{3^2}{2^5\pi} \biggr)^{1/2}  \biggl\{ 
\biggl[ 2\eta_s -\sin[2(\eta_s - B)] \biggr]  
-
\biggl[ 2\eta_i -\sin[2(\eta_i - B)] \biggr]  
\biggr\} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathfrak{w}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{3^8}{2^5 \pi}\biggr)^{1 / 2}
\biggl[
\ell_i (\ell_i^4 - \frac{8}{3}\ell_i^2 - 1) (\ell_i^2 +1)^{-3} + \tan^{-1}\ell_i
\biggr] 
</math>
  </td>
</tr>

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

The first derivative is,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{d\mathfrak{g}}{d\chi} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
+~\biggl(\frac{\chi_\mathrm{eq}}{\chi^2}\biggr)\mathfrak{w} 
~-~\frac{2}{\chi_\mathrm{eq}} \biggl[\biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{2 - 3\gamma_c}\mathfrak{s}_\mathrm{core}  
~+~\biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{2 - 3\gamma_e}\mathfrak{s}_\mathrm{env}  \biggr]  \, .
</math>
  </td>
</tr>
</table>
And the second derivative is,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{d^2\mathfrak{g}}{d\chi^2} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
-~2\biggl(\frac{\chi_\mathrm{eq}}{\chi^3}\biggr)\mathfrak{w} 
~-~\frac{2}{\chi_\mathrm{eq}^2} \biggl[(2 - 3\gamma_c) \biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{1 - 3\gamma_c}\mathfrak{s}_\mathrm{core}  
~+~(2 - 3\gamma_e)\biggl(\frac{\chi}{\chi_\mathrm{eq}}\biggr)^{1 - 3\gamma_e}\mathfrak{s}_\mathrm{env}  \biggr] \ .
</math>
  </td>
</tr>
</table>

==What to Expect for Equilibrium Configurations==
In equilibrium we should set <math>~d\mathfrak{g}/d\chi = 0</math> and <math>~\chi = \chi_\mathrm{eq}</math>.  In which case we should expect to find that,
<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>
\frac{1}{\chi_\mathrm{eq}}\biggl[\mathfrak{w}
~-~2(\mathfrak{s}_\mathrm{core}  
~+~\mathfrak{s}_\mathrm{env})  \biggr]  \, .
</math>
  </td>
</tr>
</table>
In addition, an evaluation of the second derivative should give,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{\chi_\mathrm{eq}^2}{2} \cdot \frac{d^2\mathfrak{g}}{d\chi^2} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
-~\mathfrak{w} 
~-~ \biggl[(2 - 3\gamma_c) \mathfrak{s}_\mathrm{core}  
~+~(2 - 3\gamma_e)\mathfrak{s}_\mathrm{env}  \biggr] \, .
</math>
  </td>
</tr>
</table>
Then, the transition from stable to unstable configurations occurs when <math>~d^2\mathfrak{g}/d\chi^2 = 0</math>, that is, when,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{w} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~-~ \biggl[(2 - 3\gamma_c) \mathfrak{s}_\mathrm{core}  
~+~(2 - 3\gamma_e)\mathfrak{s}_\mathrm{env}  \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~(4 - 3\gamma_c) \mathfrak{s}_\mathrm{core}  
~+~(4 - 3\gamma_e)\mathfrak{s}_\mathrm{env}  
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow ~~~\frac{\mathfrak{s}_\mathrm{core} }{\mathfrak{s}_\mathrm{env}} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~-~\frac{(4 - 3\gamma_e) }{(4 - 3\gamma_c) } \, .
</math>
  </td>
</tr>
</table>
For example, if we set <math>~\gamma_c = 6/5</math> and <math>~\gamma_e = 2</math>, this implies,

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

<tr>
  <td align="right">
<math>~\frac{\mathfrak{s}_\mathrm{core} }{\mathfrak{s}_\mathrm{env}} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\frac{2 }{(4 - 18/5) } = 5 \, .
</math>
  </td>
</tr>
</table>


==Free-Energy Stability Evaluation==

Here we pull together excerpts from several different [[H_BookTiledMenu#Tiled_Menu|H_Book Chapters]] in which we have presented, from several different perspectives, an analysis of the free-energy of bipolytropes.  
* In [[SSC/Structure/BiPolytropes/Analytic51#BiPolytrope_with_nc_.3D_5_and_ne_.3D_1|one chapter]], using purely analytic techniques, we have derived expressions that detail the structural properties of bipolytropes having <math>~(n_c, n_e) = (5, 1)</math>.  Among these are analytic expressions for various terms that make up the free-energy expression:  <math>~\mathfrak{s}_\mathrm{core}</math>, <math>~\mathfrak{s}_\mathrm{env}</math>, <math>~\mathfrak{w}_\mathrm{core}</math>, <math>~\mathfrak{w}_\mathrm{core}</math>, and <math>~P_iV_\mathrm{core}</math>. Equilibrium model ''sequences'' are defined by fixing the ratio, <math>~\mu_e/\mu_c</math>, then varying the radial location, <math>~r_i</math>, of the core-envelope interface; note that the volume of the core is, then, <math>~V_\mathrm{core} \equiv 4\pi r_i^3/3</math>.  
<div align="left">
* <font color="red">[Virial Equilibrium]</font> In a subsection of this same chapter titled, ''[[SSC/Structure/BiPolytropes/Analytic51#Global|Equilibrium Condition: Global]]'', we have shown that a statement of virial equilibrium &#8212; obtained by setting the first derivative of the free-energy expression to zero &#8212; is,<table border="0" cellpadding="4" align="center">
<tr>
  <td align="right">
<math>~( 2S + W )_\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>In another subsection of this same chapter titled, ''[[SSC/Structure/BiPolytropes/Analytic51#Global|Equilibrium Condition: In Parts]]'', we showed that for each bipolytropic equilibrium structure, 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>also hold separately. Therefore, for every equilibrium configuration we should expect the <b>CASE1</b> expression (see Table XXX) to precisely sum to unity.
</div>

<div align="left">
* <font color="red">[Marginally Unstable Model]</font> Near the bottom of this same chapter, in a subsection titled, ''[[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|Stability Condition]]'', we point out that the model along each sequence that is marginally (dynamically) unstable &#8212; obtained setting the second derivative of the free-energy expression to zero &#8212; is identified by the configuration for which,<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~
2(\gamma_e - \gamma_c) \mathfrak{s}_\mathrm{core} + ( \mathfrak{w}_\mathrm{core} +  \mathfrak{w}_\mathrm{env} )(\gamma_e - \tfrac{4}{3})
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0  \, .</math>
  </td>
</tr>
</table>Therefore, along each equilibrium sequence, the marginally unstable model can be identified by the configuration for which the <b>CASE2</b> expression (see Table XXX) precisely sums to zero.  Immediately above, in a subsection titled, ''[[SSC/Stability/BiPolytropes#What_to_Expect_for_Equilibrium_Configurations|What to Expect for Equilibrium Configurations]]'', we have shown that this same marginally unstable model can be identified by the configuration for which the <b>CASE3</b> energy ratio,  <math>~\mathfrak{s}_\mathrm{core}/\mathfrak{s}_\mathrm{env}</math>, has a value that is precisely 5. And, thirdly, as highlighted in our accompanying [[SSC/SynopsisStyleSheet#Bipolytropes|Tabular Overview]], this same marginally unstable model can be identified by the configuration for which the <b>CASE4</b> expression precisely sums to zero.
</div>


<div align="center">
<table border="1" cellpadding="8" align="center">
<tr>
  <th align="center" colspan="16">
'''Table XXX:'''  Properties of Marginally Unstable Bipolytropes Having<br /><br /><math>~(n_c, n_e) = (5, 1)</math> and <math>~(\gamma_c, \gamma_e) = (\tfrac{6}{5}, 2)</math><br /><br />Determined from Free-Energy Arguments
  </th>
</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>~R^*_\mathrm{surf}</math></td>
  <td align="center"><math>~q \equiv \frac{r_\mathrm{core}}{R_\mathrm{surf}}</math></td>
  <td align="center"><math>~\nu \equiv \frac{M_\mathrm{core}}{M_\mathrm{tot}}</math></td>
  <td align="center"><math>~\frac{\rho_c}{\bar\rho}</math></td>
  <td align="center"><math>~P_i</math></td>
  <td align="center"><math>~\mathfrak{s}_\mathrm{core}</math></td>
  <td align="center"><math>~\mathfrak{w}_\mathrm{core}</math></td>
  <td align="center"><math>~\mathfrak{s}_\mathrm{env}</math></td>
  <td align="center"><math>~\mathfrak{w}_\mathrm{env}</math></td>
  <td align="center"><math>~3P_i V_\mathrm{core}</math></td>
  <td align="center"><b>CASE1</b></td>
  <td align="center"><b>CASE2</b></td>
  <td align="center"><b>CASE3</b></td>
  <td align="center"><b>CASE4</b></td>
</tr>
<tr>
  <td align="center">1</td>
  <td align="left">2.41610822</td>
  <td align="right">2.8049</td>
  <td align="center">0.59520261</td>
  <td align="center">0.68306067</td>
  <td align="center">16.3788</td>
  <td align="center">0.039116848</td>
  <td align="center">4.446748782</td>
  <td align="center">- 6.606135366</td>
  <td align="center">0.889349762</td>
  <td align="center">- 4.066061722</td>
  <td align="center">2.287362198</td>
  <td align="center">1</td>
  <td align="center">0</td>
  <td align="center">5.0000000</td>
  <td align="center">2.1 &times; 10<sup>-8</sup></td>
</tr>
<tr>
  <td align="center"><math>~\tfrac{1}{2}</math></td>
  <td align="left">4.1853093</td>
  <td align="right">8.8058</td>
  <td align="center">0.328419479</td>
  <td align="center">0.70131896</td>
  <td align="center">354.089</td>
  <td align="center">0.003126324</td>
  <td align="center">5.76978580</td>
  <td align="center">- 10.58931853</td>
  <td align="center">1.153956968</td>
  <td align="center">- 3.258165567</td>
  <td align="center">0.95025163</td>
  <td align="center">1</td>
  <td align="center">6.3 &times; 10<sup>-8</sup></td>
  <td align="center">5.0000002</td>
  <td align="center">0</td>
</tr>
<tr>
  <td align="center">0.345</td>
  <td align="left">7.64325</td>
  <td align="right">44.116</td>
  <td align="center">0.119714454</td>
  <td align="center">0.52700045</td>
  <td align="center">2.85 &times; 10<sup>4</sup></td>
  <td align="center">0.000116533</td>
  <td align="center">6.230343527</td>
  <td align="center">- 12.24495934</td>
  <td align="center">1.1246068658</td>
  <td align="center">- 2.707865028</td>
  <td align="center">0.215727713</td>
  <td align="center">1</td>
  <td align="center">6.4 &times; 10<sup>-8</sup></td>
  <td align="center">5.0000002</td>
  <td align="center">0</td>
</tr>
<tr>
  <td align="center"><math>~\tfrac{1}{3}</math></td>
  <td align="left">8.548103</td>
  <td align="right">59.643</td>
  <td align="center">0.099032423</td>
  <td align="center">0.47901529</td>
  <td align="center">6.30 &times; 10<sup>4</sup></td>
  <td align="center">6.1337 &times; 10<sup>-5</sup></td>
  <td align="center">6.261548334</td>
  <td align="center">- 12.36425897</td>
  <td align="center">1.252309682</td>
  <td align="center">- 2.663457063</td>
  <td align="center">0.158837699</td>
  <td align="center">1</td>
  <td align="center">0</td>
  <td align="center">4.9999999</td>
  <td align="center">6.0 &times; 10<sup>-8</sup></td>
</tr>
<tr>
  <td align="center">0.316943</td>
  <td align="left">10.7441565</td>
  <td align="right">108.14</td>
  <td align="center">0.068655205</td>
  <td align="center">0.38238387</td>
  <td align="center">2.93 &times; 10<sup>5</sup></td>
  <td align="center">1.6252 &times; 10<sup>-5</sup></td>
  <td align="center">6.301810768</td>
  <td align="center">- 12.52005323</td>
  <td align="center">1.260362204</td>
  <td align="center">- 2.604292714</td>
  <td align="center">0.083568307</td>
  <td align="center">1</td>
  <td align="center">0</td>
  <td align="center">4.9999998</td>
  <td align="center">2.0 &times; 10<sup>-7</sup></td>
</tr>
<tr>
  <td align="center">0.309</td>
  <td align="left">12.77156</td>
  <td align="right">166.06</td>
  <td align="center">0.053145011</td>
  <td align="center">0.31696879</td>
  <td align="center">8.70 &times; 10<sup>5</sup></td>
  <td align="center">5.8905 &times; 10<sup>-6</sup></td>
  <td align="center">6.318902171</td>
  <td align="center">- 12.58692884</td>
  <td align="center">1.26378042</td>
  <td align="center">- 2.57843634</td>
  <td align="center">0.050875500</td>
  <td align="center">1</td>
  <td align="center">1.9 &times; 10<sup>-8</sup></td>
  <td align="center">5.0000001</td>
  <td align="center">0</td>
</tr>
<tr>
  <td align="center" colspan="16">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<b>CASE1</b>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\frac{1}{6P_i V_\mathrm{core}}\biggl[
(2 \mathfrak{s}_\mathrm{core} + \mathfrak{w}_\mathrm{core})  
-(2 \mathfrak{s}_\mathrm{core} + \mathfrak{w}_\mathrm{core})  
\biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
<b>CASE2</b>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
2(\gamma_e - \gamma_c) \mathfrak{s}_\mathrm{core} + ( \mathfrak{w}_\mathrm{core} +  \mathfrak{w}_\mathrm{env} )(\gamma_e - \tfrac{4}{3})
</math>
  </td>
</tr>

<tr>
  <td align="right">
<b>CASE3</b>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\frac{\mathfrak{s}_\mathrm{core}}{\mathfrak{s}_\mathrm{env}}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<b>CASE4</b>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
(4 - 3\gamma_c)\mathfrak{w}_\mathrm{core} + (4 - 3\gamma_e)\mathfrak{w}_\mathrm{env} + 3^2(\gamma_c - \gamma_e)P_i V_\mathrm{core}
</math>
  </td>
</tr>
</table>

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


The left-hand panel of Figure 3 is identical to Figure 2, above.  It displays in the <math>~q-\nu</math> parameter space, the behavior of <math>~(n_c, n_e) = (5, 1)</math> bipolytropic equilibrium sequences that have, as labeled, seven different values of the ratio of mean-molecular-weights, <math>~\mu_e/\mu_c</math>.  Using a numerical root-finding technique, we have determined where the virial stability condition, <math>~\mathfrak{s}_\mathrm{core}/\mathfrak{s}_\mathrm{env} = 5</math>, is satisfied along each of these sequences &#8212; as well as along a number of additional equilibrium sequences.  Key properties of each of these identified models have been recorded in [[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|Table 1 of an accompanying discussion]]; see also an [[SSC/Structure/BiPolytropes/FreeEnergy51#Free_Energy_of_BiPolytrope_with|associated discussion of the free-energy of these configurations]].  Pulling from this tabulated data, the solid-red circular markers that appear in the right-hand panel of Figure 3 identify where this virial stability condition is satisfied along the separate equilibrium sequences while the accompanying dashed red curve identifies more broadly how the <math>~q-\nu</math> parameter space is divided into stable (below and to the right) versus unstable (above and to the left) regions.  

<table border="0" align="center" cellpadding="8">
<tr>
  <th align="center">Figure 3</th>
</tr>
<tr>
  <td align="center">
[[File:CompositeAlabeled.png|800px|dynamical stability in qNu plane]]
  </td>
</tr>
</table>

In what follows we use a complementary &#8212; and more quantitatively rigorous &#8212; approach to evaluating the stability of equilibrium models, and contrast the results of that analysis with the virial-analysis results presented graphically here in Figure 3.

=See Also=

* [http://adsabs.harvard.edu/abs/2018Sci...362..201D K. De et al. (12 October 2018, Science, Vol. 362, No. 6411, pp. 201 - 206)], ''A Hot and Fast Ultra-stripped Supernova that likely formed a Compact Neutron Star Binary.''


{{ SGFfooter }}