Editing SSC/Virial/PolytropesEmbedded/SecondEffortAgain/Pt2


==Mass-Radius Relation==
Up to this point in our discussion, we have focused on an analysis of the pressure-radius relationship that defines the equilibrium configurations of pressure-truncated polytropes.  In effect, we have viewed the problem through the same lens as did {{ Horedt70full }} and, separately, {{ Whitworth81full }}, defining variable normalizations in terms of the polytropic constant, <math>K</math>, and the configuration mass, <math>M_\mathrm{tot}</math>, which were both assumed to be held fixed throughout the analysis.  Here we switch to the approach championed by {{ Stahler83full }}, defining variable normalizations in terms of <math>K</math> and <math>P_e</math>, and examining the ''mass-radius'' relationship of pressure-truncated polytropes.

===Detailed Force-Balanced Solution===
As has been summarized in our [[SSC/Structure/PolytropesEmbedded#Stahler.27s_Presentation|accompanying review]] of detailed force-balanced models of pressure-truncated polytropes,  {{ Stahler83 }} &#8212; hereafter, {{ Stahler83hereafter }} &#8212; found that a spherical configuration's equilibrium radius is related to its mass through the following pair of parametric equations:
<div align="center">
<table border="0" cellpadding="3">

<tr>
  <td align="right">
<math>
~\frac{M_\mathrm{limit}}{M_\mathrm{SWS} }
</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{n^3}{4\pi} \biggr)^{1/2} \tilde\theta^{(n-3)/2} (- \tilde\xi^2 \tilde\theta^') \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
~\frac{R_\mathrm{eq}}{R_\mathrm{SWS} }
</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{n}{4\pi} \biggr)^{1/2} \tilde\xi \tilde\theta^{(n-1)/2} \, ,
</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<math>M_\mathrm{SWS} \equiv 
\biggl( \frac{n+1}{nG} \biggr)^{3/2} K^{2n/(n+1)} P_\mathrm{e}^{(3-n)/[2(n+1)]} \, ,</math>
</div>
<div align="center">
<math>
R_\mathrm{SWS} \equiv \biggl( \frac{n+1}{nG} \biggr)^{1/2} K^{n/(n+1)} P_\mathrm{e}^{(1-n)/[2(n+1)]} \, .
</math>
</div>

===Mapping from Above Discussion===
====Deriving Concise Virial Theorem Mass-Radius Relation====
Looking back on the definitions of <math>~\Pi_\mathrm{ad}</math> and <math>~\Chi_\mathrm{ad}</math> that we introduced in connection with our initial [[#ConciseVirial|concise algebraic expression of the virial theorem]], we can write,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>~P_e </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~P_\mathrm{norm} \biggl( \frac{3}{4\pi} \biggr) \Pi_\mathrm{ad}  \biggl[ \frac{\mathcal{B}^{4n}}{\mathcal{A}^{3(n+1)}} \biggr]^{1/(n-3)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\biggl( \frac{3}{4\pi} \biggr) \Pi_\mathrm{ad}  \biggl[ \frac{\mathcal{B}^{4n}}{\mathcal{A}^{3(n+1)}} \biggr]^{1/(n-3)}
\biggl[ \frac{K^{4n}}{G^{3(n+1)} M_\mathrm{tot}^{2(n+1)}} \biggr]^{1/(n-3)} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~R_\mathrm{eq} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~R_\mathrm{norm} \Chi_\mathrm{ad} \biggl[ \frac{\mathcal{A}}{\mathcal{B}} \biggr]^{n/(n-3)} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\Chi_\mathrm{ad} \biggl[ \frac{\mathcal{A}}{\mathcal{B}} \biggr]^{n/(n-3)} 
\biggl[ \biggl( \frac{G}{K} \biggr)^n M_\mathrm{tot}^{n-1} \biggr]^{1/(n-3)} \, .
</math>
  </td>
</tr>

</table>
</div>
The first of these two expressions can be flipped around to give an expression for <math>~M_\mathrm{tot}</math> in terms of <math>~P_e</math> and, then, as normalized to <math>~M_\mathrm{SWS}</math>.  Specifically,
<div align="center">
<table border="0" cellpadding="5">

<tr>
  <td align="right">
<math>~ M_\mathrm{tot}^{2(n+1)}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{n-3}  \biggl[ \frac{\mathcal{B}^{4n}}{\mathcal{A}^{3(n+1)}} \biggr]
\biggl[ \frac{K^{4n}}{G^{3(n+1)}P_e^{n-3} } \biggr] 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~M_\mathrm{SWS}^{2(n+1)} \biggl( \frac{n}{n+1} \biggr)^{3(n+1)} 
\biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{n-3}  \biggl[ \frac{\mathcal{B}^{4n}}{\mathcal{A}^{3(n+1)}} \biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~ \Rightarrow~~~ \frac{M_\mathrm{tot}}{M_\mathrm{SWS}}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\biggl( \frac{n}{n+1} \biggr)^{3/2} 
\biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{(n-3)/[2(n+1)]}  \biggl[ \frac{\mathcal{B}^{2n/(n+1)}}{\mathcal{A}^{3/2}} \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>

This means, as well, that we can rewrite the equilibrium radius as,
<div align="center">
<table border="0" cellpadding="5">

<tr>
  <td align="right">
<math>~R_\mathrm{eq}^{n-3} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\Chi_\mathrm{ad}^{n-3} \biggl[ \frac{\mathcal{A}}{\mathcal{B}} \biggr]^{n} 
\biggl( \frac{G}{K} \biggr)^n M_\mathrm{tot}^{n-1} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\Chi_\mathrm{ad}^{n-3} \biggl[ \frac{\mathcal{A}}{\mathcal{B}} \biggr]^{n} 
\biggl( \frac{G}{K} \biggr)^n \biggl\{   
\biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{n-3}  \biggl[ \frac{\mathcal{B}^{4n}}{\mathcal{A}^{3(n+1)}} \biggr]
\biggl[ \frac{K^{4n}}{G^{3(n+1)}P_e^{n-3} } \biggr] 
\biggr\}^{(n-1)/[2(n+1)]} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\Chi_\mathrm{ad}^{n-3} \biggl[ \frac{\mathcal{A}}{\mathcal{B}} \biggr]^{n} 
\biggl[ \frac{\mathcal{B}^{4n}}{\mathcal{A}^{3(n+1)}} \biggr]^{(n-1)/[2(n+1)]}
\biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{(n-3)(n-1)/[2(n+1)]} \biggl( \frac{G}{K} \biggr)^n 
\biggl\{  
\biggl[ \frac{K^{4n}}{G^{3(n+1)}P_e^{n-3} } \biggr] 
\biggr\}^{(n-1)/[2(n+1)]} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\Chi_\mathrm{ad}^{n-3} \biggl\{ \biggl[ \frac{\mathcal{A}}{\mathcal{B}} \biggr]^{2n(n+1)} 
\biggl[ \frac{\mathcal{B}^{4n(n-1)}}{\mathcal{A}^{3(n+1)(n-1)}} \biggr]\biggr\}^{1/[2(n+1)]}
\biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{(n-3)(n-1)/[2(n+1)]} 
\biggl\{ \biggl( \frac{G}{K} \biggr)^{2n(n+1)} 
\biggl[ \frac{K^{4n(n-1)}}{G^{3(n+1)(n-1)}P_e^{(n-3)(n-1)} } \biggr] 
\biggr\}^{1/[2(n+1)]} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\Chi_\mathrm{ad}^{n-3} \biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{(n-3)(n-1)/[2(n+1)]} 
\biggl[ \mathcal{A}^{-(n+1)(n-3)} \mathcal{B}^{2n(n-3)} \biggr]^{1/[2(n+1)]}
\biggl[ G^{(3-n)(n+1)} K^{2n(n-3)} P_e^{(n-3)(1-n)} \biggr]^{1/[2(n+1)]}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~R_\mathrm{SWS}^{n-3} \biggl( \frac{n}{n+1} \biggr)^{(n-3)/2} 
\Chi_\mathrm{ad}^{n-3} \biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{(n-3)(n-1)/[2(n+1)]} 
\biggl[ \mathcal{A}^{-(n+1)(n-3)} \mathcal{B}^{2n(n-3)} \biggr]^{1/[2(n+1)]}
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Rightarrow~~~\frac{R_\mathrm{eq}}{R_\mathrm{SWS} } </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\biggl( \frac{n}{n+1} \biggr)^{1/2} 
\Chi_\mathrm{ad} \biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{(n-1)/[2(n+1)]} 
\biggl[ \frac{\mathcal{B}^{n/(n+1)}}{\mathcal{A}^{1/2}} \biggr] \, .
</math>
  </td>
</tr>

</table>
</div>
Flipping both of these expressions around, we see that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\Pi_\mathrm{ad} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\frac{4\pi}{3} \biggl\{ \biggl( \frac{n+1}{n} \biggr)^{3(n+1)}   \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{2(n+1)} 
\biggl[ \frac{\mathcal{A}^{3(n+1)}}{\mathcal{B}^{4n}} \biggr] \biggr\}^{1/(n-3)} \, ,
</math>
  </td>
</tr>
</table>
</div>
and,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\Chi_\mathrm{ad} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\frac{R_\mathrm{eq}}{R_\mathrm{SWS} }  \biggl( \frac{n+1}{n} \biggr)^{1/2} 
\biggl[ \frac{\mathcal{A}^{1/2}}{\mathcal{B}^{n/(n+1)}} \biggr] 
\biggl( \frac{3\Pi_\mathrm{ad} }{4\pi} \biggr)^{(1-n)/[2(n+1)]} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\frac{R_\mathrm{eq}}{R_\mathrm{SWS} }  \biggl( \frac{n+1}{n} \biggr)^{1/2} 
\biggl[ \frac{\mathcal{A}^{1/2}}{\mathcal{B}^{n/(n+1)}} \biggr] 
\biggl\{ \biggl( \frac{n+1}{n} \biggr)^{3(n+1)}   \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{2(n+1)} 
\biggl[ \frac{\mathcal{A}^{3(n+1)}}{\mathcal{B}^{4n}} \biggr] \biggr\}^{(1-n)/[2(n+1)(n-3)]} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\frac{R_\mathrm{eq}}{R_\mathrm{SWS} } \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{(1-n)/(n-3)}  \biggl( \frac{n}{n+1} \biggr)^{n/(n-3)} 
\biggl[ \frac{\mathcal{B}}{\mathcal{A}} \biggr]^{n/(n-3)} \, .
</math>
  </td>
</tr>
</table>
</div>
Hence, our earlier derived [[#ConciseVirial3|compact expression for the virial theorem]] becomes,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl\{ \frac{R_\mathrm{eq}}{R_\mathrm{SWS} } \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{(1-n)/(n-3)}  \biggl( \frac{n}{n+1} \biggr)^{n/(n-3)} 
\biggl[ \frac{\mathcal{B}}{\mathcal{A}} \biggr]^{n/(n-3)} \biggr\}^{(n-3)/n} 
</math> 
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>
-~ \frac{4\pi}{3} \biggl\{ \biggl( \frac{n+1}{n} \biggr)^{3(n+1)}   \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{2(n+1)} 
\biggl[ \frac{\mathcal{A}^{3(n+1)}}{\mathcal{B}^{4n}} \biggr] \biggr\}^{1/(n-3)} 
\biggl\{ \frac{R_\mathrm{eq}}{R_\mathrm{SWS} } \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{(1-n)/(n-3)}  \biggl( \frac{n}{n+1} \biggr)^{n/(n-3)} 
\biggl[ \frac{\mathcal{B}}{\mathcal{A}} \biggr]^{n/(n-3)} \biggr\}^4 </math> 
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS} } \biggr)^{(n-3)/n} \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{(1-n)/n}  \biggl( \frac{n}{n+1} \biggr)
\biggl[ \frac{\mathcal{B}}{\mathcal{A}} \biggr]
-~ \frac{4\pi}{3} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS} }  \biggr)^4 \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{-2} \biggl( \frac{n}{n+1} \biggr)
\frac{1}{\mathcal{A}} \, .
</math> 
  </td>
</tr>
</table>
</div>

Or, rearranged,

<div align="center" id="CompactStahlerVirial">
<table border="1" cellpadding="10" align="center">

<tr>
  <td align="right">
<math>\frac{4\pi}{3} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS} }  \biggr)^4  - 
\mathcal{B} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS} } \biggr)^{(n-3)/n} \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{(n+1)/n}  
+~ \mathcal{A} \biggl( \frac{n+1}{n} \biggr) \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{2} = 0 \, .
</math> 
  </td>
</tr>
</table>
</div>

After adopting the modified coefficient definitions,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>~\mathcal{A}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\mathcal{A} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}} \biggr)^{-2}  =
\frac{1}{5} \cdot \frac{\tilde\mathfrak{f}_W}{\tilde\mathfrak{f}_M^2}  \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathcal{B}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\mathcal{B} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{-(n+1)/n} =
\biggl( \frac{3}{4\pi}\biggr)^{1/n} 
\frac{\tilde\mathfrak{f}_A}{\tilde\mathfrak{f}_M^{(n+1)/n}} \, ,
</math>
  </td>
</tr>
</table>
</div>

<span id="modNormalizations">as well as the modified length- and mass-normalizations,</span> <math>~R_\mathrm{mod}</math> and <math>~M_\mathrm{mod}</math>, such that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{M_\mathrm{SWS}}{M_\mathrm{mod}}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{4\pi}{3} \biggr)^{2n/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{3/2} \frac{\mathcal{A}_{M_\ell}^{3/2}}{\mathcal{B}_{M_\ell}^{2n/(n+1)}} 
\, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\frac{R_\mathrm{SWS}}{R_\mathrm{mod}}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\biggl( \frac{4\pi}{3} \biggr)^{n/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{1/2} \frac{\mathcal{A}_{M_\ell}^{1/2}}{\mathcal{B}_{M_\ell}^{n/(n+1)}} \, ,</math>
  </td>
</tr>
</table>
</div>
we obtain the
<div align="center" id="ConciseVirialMR">
<font color="#770000">'''Virial Theorem in terms of Mass and Radius'''</font><br />

<math>
\biggl( \frac{R_\mathrm{eq}}{R_\mathrm{mod}} \biggr)^4 
- \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{mod}} \biggr)^{(n-3)/n} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{mod}} \biggr)^{(n+1)/n} 
+ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{mod}} \biggr)^2 = 0 
\, .
</math>
</div>


For later use we note as well that, with these modified coefficient definitions, we can write,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\Pi_\mathrm{ad}^{n-3} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\biggl[ \biggl( \frac{4\pi}{3} \biggr)^{n-3} \biggl( \frac{n+1}{n} \biggr)^{3(n+1)}    
\frac{\mathcal{A}_{M_\ell}^{3(n+1)}}{\mathcal{B}_{M_\ell}^{4n}} \biggr] \mathcal{Y}^{2(n+1)} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\Chi_\mathrm{ad}^{n-3} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~ 
\biggl[ \frac{n}{n+1} \biggl( \frac{\mathcal{B}_{M_\ell}}{\mathcal{A}_{M_\ell}} \biggr)\biggr]^n 
\mathcal{X}^{n-3} \mathcal{Y}^{1-n} \, ,
</math>
  </td>
</tr>
</table>
</div>
where <math>~\mathcal{X}</math> and <math>~\mathcal{Y}</math> are defined [[#Confirmation|immediately below]].

====Corresponding Concise Free-Energy Expression====
Let's also rewrite the [[#FreeEnergyExpression|algebraic free-energy function]] in terms of Stahler's normalized mass and radius variables.  Expressed in terms of the polytropic index, the free-energy function is,
<div align="center">
<math>
\mathfrak{G}^* = 
-3\mathcal{A} \chi^{-1} +~ n\mathcal{B} \chi^{-3/n} +~ \mathcal{D}\chi^3 \, .
</math>
</div>
First, we recognize that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\chi \equiv \frac{R}{R_\mathrm{norm}}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(\frac{R}{R_\mathrm{SWS}}\biggr) \frac{R_\mathrm{SWS}}{R_\mathrm{norm}} \, .</math>
  </td>
</tr>
</table>
</div>
From the definition of <math>~R_\mathrm{norm}</math> &#8212; reprinted, for example, [[#Detailed_Force-Balanced_Solution|here]] &#8212; we can write,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\biggl( \frac{R_\mathrm{SWS}}{R_\mathrm{norm}}\biggr)^{n-3}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~R_\mathrm{SWS}^{n-3} \biggl[ G^{-n} K^n M_\mathrm{tot}^{1-n} \biggr]
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~R_\mathrm{SWS}^{n-3} M_\mathrm{SWS}^{1-n} \biggl[ \biggl(\frac{K}{G}\biggr)^n \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}}\biggr)^{1-n} \biggr] \, ;
</math>
  </td>
</tr>
</table>
</div>
and from the definitions of <math>~R_\mathrm{SWS}</math> and <math>~M_\mathrm{SWS}</math> &#8212; reprinted, for example, [[#Detailed_Force-Balanced_Solution_2|here]] &#8212; we have,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\biggl( \frac{R_\mathrm{SWS}}{R_\mathrm{norm}}\biggr)^{n-3}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \biggl(\frac{K}{G}\biggr)^n \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}}\biggr)^{1-n} \biggr]
\biggl\{\biggl( \frac{n+1}{nG} \biggr)^{1/2} K^{n/(n+1)} P_\mathrm{e}^{(1-n)/[2(n+1)]} \biggr\}^{n-3}
\biggl\{ \biggl( \frac{n+1}{nG} \biggr)^{3/2} K^{2n/(n+1)} P_\mathrm{e}^{(3-n)/[2(n+1)]} \biggr\}^{1-n}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \biggl(\frac{K}{G}\biggr)^n \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}}\biggr)^{1-n} \biggr]
\biggl( \frac{n+1}{n} \biggr)^{[(n-3) +3(1-n)]/2} G^{[(3-n) + 3(n-1)]/2} K^{n[(n-3)+2(1-n)]/(n+1)} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}}\biggr)^{1-n} \biggl( \frac{n}{n+1} \biggr)^n \, .
</math>
  </td>
</tr>
</table>
</div>
Hence, in each term in the free-energy expression we can make the substitution,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\chi </math>
  </td>
  <td align="center">
<math>~~~\rightarrow~~~</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr) \biggl\{ \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}}\biggr)^{1-n} \biggl( \frac{n}{n+1} \biggr)^n \biggr\}^{1/(n-3)}
= 
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)\biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{(1-n)/(n-3)}  
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{n-1} \biggl( \frac{n}{n+1} \biggr)^n \biggr\}^{1/(n-3)} \, .
</math>
  </td>
</tr>
</table>
</div>

Next, drawing on the definition of <math>P_\mathrm{norm}</math> &#8212; reprinted, for example, [[#Detailed_Force-Balanced_Solution|here]] &#8212; along with the definition of <math>M_\mathrm{SWS}</math>, we recognize that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathcal{D} \equiv \frac{4\pi}{3} \cdot \frac{P_e}{P_\mathrm{norm}}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{4\pi}{3} \cdot P_e \biggl[ K^{-4n} G^{3(n+1)} M_\mathrm{tot}^{2(n+1)} \biggr]^{1/(n-3)} </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{4\pi}{3} \cdot \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \biggr)^{2(n+1)/(n-3)} P_e \biggl[ K^{-4n} G^{3(n+1)}\biggr]^{1/(n-3)} M_\mathrm{SWS}^{2(n+1)/(n-3)}  </math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{4\pi}{3} \cdot \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{2(n+1)/(n-3)} 
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}} \biggr)^{-2}\biggl( \frac{n+1}{n}\biggr)^3 \biggr\}^{(n+1)/(n-3)} \, .</math>
  </td>
</tr>
</table>
</div>

After making these substitutions into the free-energy function, as well as replacing <math>~\mathcal{A}</math> and <math>~\mathcal{B}</math> with <math>~\mathcal{A}_{M_\ell}</math> and <math>~\mathcal{B}_{M_\ell}</math>, respectively, 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>~
-3\mathcal{A}_{M_\ell}\biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}} \biggr)^{2} 
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-1} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{-(1-n)/(n-3)}  
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{n-1} \biggl( \frac{n}{n+1} \biggr)^n \biggr\}^{-1/(n-3)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+~ n\mathcal{B}_{M_\ell} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}} \biggr)^{(n+1)/n}
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-3/n} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{-3(1-n)/[n(n-3)]}  
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{n-1} \biggl( \frac{n}{n+1} \biggr)^n \biggr\}^{-3/[n(n-3)]}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+~ \frac{4\pi}{3} \cdot \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{2(n+1)/(n-3)} 
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}} \biggr)^{-2}\biggl( \frac{n+1}{n}\biggr)^3 \biggr\}^{(n+1)/(n-3)}
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^3 \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{3(1-n)/(n-3)}  
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{n-1} \biggl( \frac{n}{n+1} \biggr)^n \biggr\}^{3/(n-3)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-3\mathcal{A}_{M_\ell} 
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-1} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{-(1-n)/(n-3)}  
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{[(n-1)-2(n-3)]} \biggl( \frac{n}{n+1} \biggr)^n \biggr\}^{-1/(n-3)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+~  n\mathcal{B}_{M_\ell}  
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-3/n} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{-3(1-n)/[n(n-3)]}  
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{[(n+1)(n-3)-3(n-1)]} \biggl( \frac{n}{n+1} \biggr)^{-3n} \biggr\}^{1/[n(n-3)]}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+~ \frac{4\pi}{3} \cdot \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^3 \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{[2(n+1)+3(1-n)]/(n-3)}   
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}} \biggr)^{[(3(n-1) -2(n+1)]}\biggl( \frac{n}{n+1}\biggr)^{[3n-3(n+1)]} \biggr\}^{1/(n-3)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-3\mathcal{A}_{M_\ell} 
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-1} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{-(1-n)/(n-3)}  
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{(n-5)} \biggl( \frac{n+1}{n} \biggr)^n \biggr\}^{1/(n-3)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+~  n\mathcal{B}_{M_\ell}  
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-3/n} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{-3(1-n)/[n(n-3)]}  
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}}\biggr)^{(n-5)} \biggl( \frac{n+1}{n} \biggr)^{3} \biggr\}^{1/(n-3)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+~ \frac{4\pi}{3} \cdot \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^3 \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{(5-n)/(n-3)}   
\biggl\{ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}} \biggr)^{(n-5)}\biggl( \frac{n+1}{n}\biggr)^{3} \biggr\}^{1/(n-3)}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl[ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{tot}} \biggr)^{(n-5)}\biggl( \frac{n+1}{n}\biggr)^{3} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{(5-n)} \biggr]^{1/(n-3)} 
\biggl\{
-3\mathcal{A}_{M_\ell}\biggl( \frac{n+1}{n}\biggr) \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-1} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{2} 
+~  n\mathcal{B}_{M_\ell}  
\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-3/n} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}\biggr)^{(n+1)/n}  
+~ \frac{4\pi}{3} \cdot \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{3}   
\biggr\}\, .
</math>
  </td>
</tr>
</table>
</div>

Hence, after defining,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{G}^*_\mathrm{SWS}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~
\frac{\mathfrak{G}}{[G^{-3} K^n M_\mathrm{SWS}^{n-5}]^{1/(n-3)}} \biggl( \frac{n}{n+1}\biggr)^{3/(n-3)} 
=
\frac{\mathfrak{G}}{[K^{6n} P_e^{5-n}]^{1/[2(n+1)]}} \biggl( \frac{nG}{n+1}\biggr)^{3/2} \, ,
</math>
  </td>
</tr>
</table>
</div>

we can write,
<div align="center" id="ConciseFreeEnergyExpression">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathfrak{G}^*_\mathrm{SWS} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
-3\mathcal{A}_{M_\ell}\biggl( \frac{n+1}{n}\biggr) \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{2} \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-1} 
+~ n\mathcal{B}_{M_\ell}  \cdot \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{(n+1)/n}  \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-3/n} 
+~ \frac{4\pi}{3} \cdot \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{3}  \, . 
</math>
  </td>
</tr>
</table>
</div>
Setting the first derivative of this function equal to zero should produce the virial theorem expression.  Let's see &hellip;
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{\partial\mathfrak{G}^*_\mathrm{SWS}}{\partial \mathcal{X}}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
3\mathcal{A}_{M_\ell}\biggl( \frac{n+1}{n}\biggr) \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{2} \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-2} 
-~ 3\mathcal{B}_{M_\ell} \cdot \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{(n+1)/n}  \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-(3+n)/n} 
+~ 4\pi \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{2}  
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ 3\biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-2}
\biggl[ \mathcal{A}_{M_\ell}\biggl( \frac{n+1}{n}\biggr) \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{2} 
-~ \mathcal{B}_{M_\ell} \cdot \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{(n+1)/n}  \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{(n-3)/n} 
+~ \frac{4\pi}{3} \cdot \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{4} \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>
Replacing <math>\mathcal{A}_{M_\ell}</math> and <math>\mathcal{B}_{M_\ell}</math> with <math>\mathcal{A}</math> and <math>\mathcal{B}</math>, as prescribed by their defined relationships, and setting the expression inside the square brackets equal to zero does, indeed, produce the [[#CompactStahlerVirial|above, boxed-in ''viral theorem'' mass-radius relationship]].

====Plotting Concise Mass-Radius Relation====
Our derived, [[#ConciseVirialMR|concise analytic expression for the virial theorem]], namely,
<div align="center">
<math>
\biggl( \frac{R_\mathrm{eq}}{R_\mathrm{mod}} \biggr)^4 
- \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{mod}} \biggr)^{(n-3)/n} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{mod}} \biggr)^{(n+1)/n} 
+ \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{mod}} \biggr)^2 = 0 
\, ,
</math>
</div>

is plotted for seven different values of the polytropic index, <math>~n</math>, as indicated, in the lefthand diagram of the following composite figure.  For comparison, the ''schematic'' diagram displayed on the righthand side of the figure is a reproduction of Figure 17 from Appendix B of {{ Stahler83 }}.  It seems that our derived, analytically prescribable, mass-radius relationship &#8212; which is, in essence, a statement of the scalar virial theorem &#8212; embodies most of the attributes of the mass-radius relationship for pressure-truncated polytropes that were already understood, and conveyed schematically, by Stahler in 1983. 
<table border="1" width="100%" cellpadding="3" align="center">
<tr>
  <td align="center" colspan="2">
'''Virial Theorem Mass-Radius Relationships'''
  </td>
</tr>
<tr>
  <td align="center" rowspan="2">
[[File:MassRadiusVirialLabeled.png|350px|Virial Theorem Mass-Radius Relation]]
  </td>
  <td align="center">
Digital copy of Figure 17 from &hellip;<br />
{{ Stahler83figure }}
  </td>
</tr>
<tr>
  <td align="center">
[[File:Stahler_MRdiagram1.png|300px|center|Stahler (1983) Figure 17]]
  </td>
</tr>
</table>
<!--
<table border="1" align="center" cellpadding="3">
<tr>
  <td align="center" rowspan="2">
[[File:MassRadiusVirialLabeled.png|350px|Virial Theorem Mass-Radius Relation]]
  </td>
  <td align="center">
[[File:Stahler1983TitlePage0.png|300px|center|Stahler (1983) Title Page]]
  </td>
</tr>
<tr>
  <td align="center" bgcolor="white">
[[File:Stahler_MRdiagram1.png|300px|center|Stahler (1983) Figure 17 (edited)]]
  </td>
</tr>
</table>
-->

Let's do this again using the mass-radius relation as written explicitly in terms of the normalizations, <math>~M_\mathrm{SWS}</math> and <math>~R_\mathrm{SWS}</math>.  The relevant, generic nonlinear equation is,
<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>~
\biggl( \frac{R}{R_\mathrm{SWS}} \biggr)^4 
- \biggl( \frac{R}{R_\mathrm{SWS}} \biggr)^{(n-3)/n} \biggl[ \biggl(\frac{3}{4\pi} \biggr) \frac{M}{M_\mathrm{SWS}} \biggr]^{(n+1)/n} 
+ \frac{3}{20\pi} \biggl( \frac{n+1}{n}\biggr) \biggl( \frac{M}{M_\mathrm{SWS}} \biggr)^2 \, .
</math>
  </td>
</tr>
</table>
</div>

Analytically determined roots:
* <math>~n=1</math>
<div align="center">
<math>~\frac{M}{M_\mathrm{SWS}} = \biggl( \frac{10\pi}{3}\biggr)^{1 / 2} \biggl(\frac{R}{R_\mathrm{SWS}} \biggr)^3 \biggl[ \frac{3\cdot 5}{2^3\pi} - \biggl(\frac{R}{R_\mathrm{SWS}} \biggr)^2
\biggr]^{-1/2} </math> &nbsp; &nbsp; &nbsp; &nbsp; for, &nbsp; &nbsp; <math>~0 \le \frac{R}{R_\mathrm{SWS}} \le \biggl(\frac{3\cdot 5}{2^3\pi}\biggr)^{1/2} \, .</math>
</div>
* <math>~n=3</math>
<div align="center">
<math>~\frac{R}{R_\mathrm{SWS}} = \biggl\{ \biggl[ \biggl(\frac{3}{4\pi}\biggr)\frac{M}{M_\mathrm{SWS}}\biggr]^{4/3} 
- \biggl(\frac{1}{5\pi}\biggr) \biggl( \frac{M}{M_\mathrm{SWS}}\biggr)^2 \biggr\}^{1/4} </math> &nbsp; &nbsp; &nbsp; &nbsp; for, &nbsp; &nbsp; <math>~0 \le \frac{M}{M_\mathrm{SWS}} \le \biggl(\frac{3^4\cdot 5^3}{2^8\pi}\biggr)^{1/2} \, .</math>
</div>
* <span id="Isothermal">''Isothermal''</span> (explained [[#IsothermalExplained|immediately below]])
<div align="center">
<math>~\frac{M}{M_\mathrm{SWS}} = \frac{5}{2} \biggl( \frac{R}{R_\mathrm{SWS}} \biggr) \biggl\{
1 \pm \biggl[ 1 - \frac{16\pi}{15}\biggl( \frac{R}{R_\mathrm{SWS}} \biggr)^2  \biggr]^{1 / 2} \biggr\} </math>
&nbsp; &nbsp; &nbsp; &nbsp; for, &nbsp; &nbsp; <math>~0 \le \frac{R}{R_\mathrm{SWS}} \le \biggl(\frac{3\cdot 5}{2^4\pi}\biggr)^{1/2} \, .</math>
</div>


<span id="TabulatedValues">First, we'll create a table of the normalized coordinate values that satisfy this nonlinear expression.</span>
<div align="center">
<table border="1" align="center" cellpadding="5">
<tr>
  <td align="center" colspan="1"><math>~n =2</math></td>
  <td align="center" colspan="1"><math>~n =2.8</math></td>
  <td align="center" colspan="1"><math>~n =3.5</math></td>
  <td align="center" colspan="1"><math>~n = 4</math></td>
  <td align="center" colspan="1"><math>~n = 5</math></td>
</tr>
<tr>
<td align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="center"><math>~ \frac{R}{R_\mathrm{SWS}}</math></td>
  <td align="center">&nbsp;</td>
  <td align="center"><math>~ \frac{M}{M_\mathrm{SWS}}</math></td>
</tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr>
  <td align="right">0.3800</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.26562</td>
</tr>
<tr>
  <td align="right">0.4500</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.477153</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.70919</td>
</tr>
<tr>
  <td align="right">0.5500</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.063602</td>
</tr>
<tr>
  <td align="right">0.5800</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.39755</td>
</tr>
<tr>
  <td align="right">0.5950</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.64662</td>
</tr>
<tr>
  <td align="right">0.6050</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.893915</td>
</tr>
<tr>
  <td align="right">0.6120</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.22372</td>
</tr>
<tr>
  <td align="right"><font color="darkgreen">0.6131721</font></td>
  <td align="center">&nbsp;</td>
  <td align="left"><font color="darkgreen">2.433375</font></td>
</tr>
<tr>
  <td align="right">0.6120</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.64923</td>
</tr>
<tr>
  <td align="right">0.6050</td>
  <td align="center">&nbsp;</td>
  <td align="left">3.01688</td>
</tr>
<tr>
  <td align="right">0.5950</td>
  <td align="center">&nbsp;</td>
  <td align="left">3.32037</td>
</tr>
<tr>
  <td align="right">0.5800</td>
  <td align="center">&nbsp;</td>
  <td align="left">3.658702</td>
</tr>
<tr>
  <td align="right">0.5500</td>
  <td align="center">&nbsp;</td>
  <td align="left">4.19097</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">4.94599</td>
</tr>
<tr>
  <td align="right">0.4700</td>
  <td align="center">&nbsp;</td>
  <td align="left">5.38791</td>
</tr>
<tr>
  <td align="right">0.4500</td>
  <td align="center">&nbsp;</td>
  <td align="left">5.69164</td>
</tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
</table>
</td>
<td align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="center"><math>~ \frac{R}{R_\mathrm{SWS}}</math></td>
  <td align="center">&nbsp;</td>
  <td align="center"><math>~ \frac{M}{M_\mathrm{SWS}}</math></td>
</tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr>
  <td align="right">0.3800</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.266134</td>
</tr>
<tr>
  <td align="right">0.4500</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.47971</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.71765</td>
</tr>
<tr>
  <td align="right">0.5250</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.831825</td>
</tr>
<tr>
  <td align="right">0.5600</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.20977</td>
</tr>
<tr>
  <td align="right">0.5750</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.427183</td>
</tr>
<tr>
  <td align="right">0.5850</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.653232</td>
</tr>
<tr>
  <td align="right">0.5900</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.89304</td>
</tr>
<tr>
  <td align="right"><font color="darkgreen">0.5904492</font></td>
  <td align="center">&nbsp;</td>
  <td align="left"><font color="darkgreen">1.989927</font></td>
</tr>
<tr>
  <td align="right">0.5900</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.086584</td>
</tr>
<tr>
  <td align="right">0.5850</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.32394</td>
</tr>
<tr>
  <td align="right">0.5750</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.54527</td>
</tr>
<tr>
  <td align="right">0.5600</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.75612</td>
</tr>
<tr>
  <td align="right">0.5250</td>
  <td align="center">&nbsp;</td>
  <td align="left">3.07134</td>
</tr>
<tr>
  <td align="right">0.4500</td>
  <td align="center">&nbsp;</td>
  <td align="left">3.460304</td>
</tr>
<tr>
  <td align="right">0.3500</td>
  <td align="center">&nbsp;</td>
  <td align="left">3.75881</td>
</tr>
<tr>
  <td align="right">0.2500</td>
  <td align="center">&nbsp;</td>
  <td align="left">3.97835</td>
</tr>
<tr>
  <td align="right">0.2000</td>
  <td align="center">&nbsp;</td>
  <td align="left">4.09302</td>
</tr>
<tr>
  <td align="right">0.1500</td>
  <td align="center">&nbsp;</td>
  <td align="left">4.232786</td>
</tr>
<tr>
  <td align="right">0.1000</td>
  <td align="center">&nbsp;</td>
  <td align="left">4.430303</td>
</tr>
<tr>
  <td align="right">0.0700</td>
  <td align="center">&nbsp;</td>
  <td align="left">4.60984</td>
</tr>
<tr>
  <td align="right">0.0400</td>
  <td align="center">&nbsp;</td>
  <td align="left">4.9057</td>
</tr>
<tr>
  <td align="right">0.0150</td>
  <td align="center">&nbsp;</td>
  <td align="left">5.47056</td>
</tr>
</table>
</td>
<td align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="center"><math>~ \frac{R}{R_\mathrm{SWS}}</math></td>
  <td align="center">&nbsp;</td>
  <td align="center"><math>~ \frac{M}{M_\mathrm{SWS}}</math></td>
</tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr>
  <td align="right">0.3800</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.26639</td>
</tr>
<tr>
  <td align="right">0.4500</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.481072</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.722406</td>
</tr>
<tr>
  <td align="right">0.5250</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.89152</td>
</tr>
<tr>
  <td align="right">0.5600</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.246123</td>
</tr>
<tr>
  <td align="right">0.5650</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.32113</td>
</tr>
<tr>
  <td align="right">0.5750</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.52651</td>
</tr>
<tr>
  <td align="right">0.5800</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.745165</td>
</tr>
<tr>
  <td align="right"><font color="darkgreen">0.5803836</font></td>
  <td align="center">&nbsp;</td>
  <td align="left"><font color="darkgreen">1.823995</font></td>
</tr>
<tr>
  <td align="right">0.5800</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.90201</td>
</tr>
<tr>
  <td align="right">0.5780</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.01647</td>
</tr>
<tr>
  <td align="right">0.5750</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.11019</td>
</tr>
<tr>
  <td align="right">0.5600</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.35906</td>
</tr>
<tr>
  <td align="right">0.5400</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.543602</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.75557</td>
</tr>
<tr>
  <td align="right">0.4500</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.543602</td>
</tr>
<tr>
  <td align="right">0.3800</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.94829</td>
</tr>
<tr>
  <td align="right" bgcolor="yellow"><font color="black">0.3749583</font></td>
  <td align="center">&nbsp;</td>
  <td align="left" bgcolor="yellow"><font color="black">2.948526</font></td>
</tr>
<tr>
  <td align="right">0.3300</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.93161</td>
</tr>
<tr>
  <td align="right">0.2500</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.829401</td>
</tr>
<tr>
  <td align="right">0.1500</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.578605</td>
</tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
</table>
</td>
<td align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="center"><math>~ \frac{R}{R_\mathrm{SWS}}</math></td>
  <td align="center">&nbsp;</td>
  <td align="center"><math>~ \frac{M}{M_\mathrm{SWS}}</math></td>
</tr>
<tr>
  <td align="right">0.1000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.004224</td>
</tr>
<tr>
  <td align="right">0.2000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.034709</td>
</tr>
<tr>
  <td align="right">0.3000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.1230901</td>
</tr>
<tr>
  <td align="right">0.4000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.31735</td>
</tr>
<tr>
  <td align="right">0.4500</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.48177</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.72493</td>
</tr>
<tr>
  <td align="right">0.5250</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.89686</td>
</tr>
<tr>
  <td align="right">0.5400</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.028495</td>
</tr>
<tr>
  <td align="right">0.5500</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.13574</td>
</tr>
<tr>
  <td align="right">0.5600</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.26965</td>
</tr>
<tr>
  <td align="right">0.5730</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.55527</td>
</tr>
<tr>
  <td align="right"><font color="darkgreen">0.5756189</font></td>
  <td align="center">&nbsp;</td>
  <td align="left"><font color="darkgreen">1.750930</font></td>
</tr>
<tr>
  <td align="right">0.5730</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.93949</td>
</tr>
<tr>
  <td align="right">0.5600</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.18983</td>
</tr>
<tr>
  <td align="right">0.5400</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.376318</td>
</tr>
<tr>
  <td align="right">0.5250</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.46661</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.56895</td>
</tr>
<tr>
  <td align="right">0.4600</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.657809</td>
</tr>
<tr>
  <td align="right" bgcolor="yellow"><font color="black">0.41184646</font></td>
  <td align="center">&nbsp;</td>
  <td align="left" bgcolor="yellow"><font color="black">2.688999</font></td>
</tr>
<tr>
  <td align="right">0.4100</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.68895</td>
</tr>
<tr>
  <td align="right">0.3800</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.677703</td>
</tr>
<tr>
  <td align="right">0.3000</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.56612</td>
</tr>
<tr>
  <td align="right">0.2500</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.44565</td>
</tr>
<tr>
  <td align="right">0.2000</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.28789</td>
</tr>
<tr>
  <td align="right">0.1500</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.08747</td>
</tr>
<tr>
  <td align="right">0.1000</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.82708</td>
</tr>
</table>
</td>
<td align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="center"><math>~ \frac{R}{R_\mathrm{SWS}}</math></td>
  <td align="center">&nbsp;</td>
  <td align="center"><math>~ \frac{M}{M_\mathrm{SWS}}</math></td>
</tr>
<tr>
  <td align="right">0.1000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.004224</td>
</tr>
<tr>
  <td align="right">0.2000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.03471</td>
</tr>
<tr>
  <td align="right">0.3000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.123115</td>
</tr>
<tr>
  <td align="right">0.4000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.31766</td>
</tr>
<tr>
  <td align="right">0.4500</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.48278</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.72866</td>
</tr>
<tr>
  <td align="right">0.5250</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.905006</td>
</tr>
<tr>
  <td align="right">0.5400</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.042907</td>
</tr>
<tr>
  <td align="right">0.5500</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.15886</td>
</tr>
<tr>
  <td align="right">0.5600</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.313712</td>
</tr>
<tr>
  <td align="right">0.5675</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.511304</td>
</tr>
<tr>
  <td align="right"><font color="darkgreen">0.5692185</font></td>
  <td align="center">&nbsp;</td>
  <td align="left"><font color="darkgreen">1.657839</font></td>
</tr>
<tr>
  <td align="right">0.5675</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.798532</td>
</tr>
<tr>
  <td align="right">0.5600</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.97061</td>
</tr>
<tr>
  <td align="right">0.5400</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.17282</td>
</tr>
<tr>
  <td align="right">0.5250</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.25888</td>
</tr>
<tr>
  <td align="right">0.5000</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.34793</td>
</tr>
<tr>
  <td align="right">0.4600</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.410374</td>
</tr>
<tr>
  <td align="right" bgcolor="yellow"><font color="black">0.4391754</font></td>
  <td align="center">&nbsp;</td>
  <td align="left" bgcolor="yellow"><font color="black">2.417330</font></td>
</tr>
<tr>
  <td align="right">0.4000</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.396465</td>
</tr>
<tr>
  <td align="right">0.3000</td>
  <td align="center">&nbsp;</td>
  <td align="left">2.19848</td>
</tr>
<tr>
  <td align="right">0.2000</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.84195</td>
</tr>
<tr>
  <td align="right">0.1000</td>
  <td align="center">&nbsp;</td>
  <td align="left">1.31421</td>
</tr>
<tr>
  <td align="right">0.0500</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.930314</td>
</tr>
<tr>
  <td align="right">0.0200</td>
  <td align="center">&nbsp;</td>
  <td align="left">0.58847</td>
</tr>
<tr><td align="center" colspan="3">&nbsp;</td></tr>
</table>
</td>

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

<span id="IsothermalExplained">From a</span> [[SSC/Virial/Isothermal#Bonnor.27s_.281956.29_Equivalent_Relation|free-energy analysis of isothermal spheres]], we have demonstrated that, when the structural form factors are all set to unity, the statement of virial equilibrium is,
<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>~
\biggl( \frac{R}{R_\mathrm{SWS}} \biggr)^4 
- \frac{3}{4\pi} \biggr( \frac{M}{M_\mathrm{SWS}} \biggr)  \biggl( \frac{R}{R_\mathrm{SWS}} \biggr)
+ \frac{3}{20\pi} \biggl( \frac{M}{M_\mathrm{SWS}} \biggr)^2 \, ,
</math>
  </td>
</tr>
</table>
</div>
where, in order to be consistent with the above polytropic normalizations, we have adopted the ''isothermal'' normalizations,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~M_\mathrm{SWS}\biggr|_\mathrm{isothermal} \equiv \biggl( \frac{c_s^8}{G^3 P_e} \biggr)^{1 / 2} \, ,</math>
  </td>
  <td align="center">
&nbsp; &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; &nbsp; 
  </td>
  <td align="left">
<math>~R_\mathrm{SWS}\biggr|_\mathrm{isothermal} \equiv \biggl( \frac{c_s^4}{G P_e} \biggr)^{1 / 2} \, .</math>
  </td>
</tr>
</table>
</div>
This is a quadratic equation that can be readily solved to provide an analytic expression for the ''isothermal'' mass-radius relation; the relevant expression has already been [[#Isothermal|provided, above]].

====Confirmation====

Rewriting the [[#ConciseVirialMR|just-derived virial theorem expression]] in terms of Stahler's dimensionless radius and mass variables, written in the abbreviated form,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathcal{X}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathcal{Y}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \, ,</math>
  </td>
</tr>
</table>
</div>
we have,

<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>~\biggl[ \mathcal{X} \cdot \frac{R_\mathrm{SWS}}{R_\mathrm{mod}} \biggr]^4 
- \biggl[ \mathcal{X} \cdot \frac{R_\mathrm{SWS}}{R_\mathrm{mod}} \biggr]^{(n-3)/n} 
\biggl[ \mathcal{Y} \cdot \frac{M_\mathrm{SWS}}{M_\mathrm{mod}} \biggr]^{(n+1)/n} 
+ \biggl[ \mathcal{Y} \cdot \frac{M_\mathrm{SWS}}{M_\mathrm{mod}} \biggr]^2</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\mathcal{X}^4 \biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{n/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{1/2} \frac{\mathcal{A}_{M_\ell}^{1/2}}{\mathcal{B}_{M_\ell}^{n/(n+1)}} \biggr\}^4 
+ \mathcal{Y}^2 \biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{2n/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{3/2} \frac{\mathcal{A}_{M_\ell}^{3/2}}{\mathcal{B}_{M_\ell}^{2n/(n+1)}} \biggr\}^2</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
- \mathcal{X}^{(n-3)/n} \mathcal{Y}^{(n+1)/n} \biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{n/(n+1)} 
\biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{1/2} \frac{\mathcal{A}_{M_\ell}^{1/2}}{\mathcal{B}_{M_\ell}^{n/(n+1)}}\biggr\}^{(n-3)/n} 
\biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{2n/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{3/2} \frac{\mathcal{A}_{M_\ell}^{3/2}}{\mathcal{B}_{M_\ell}^{2n/(n+1)}}  \biggr\}^{(n+1)/n} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\mathcal{X}^4 \biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{4n/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{2} 
\frac{\mathcal{A}_{M_\ell}^{2}}{\mathcal{B}_{M_\ell}^{4n/(n+1)}} \biggr\} 
+ \mathcal{Y}^2 \biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{4n/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{3} \frac{\mathcal{A}_{M_\ell}^{3}}{\mathcal{B}_{M_\ell}^{4n/(n+1)}} \biggr\}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
- \mathcal{X}^{(n-3)/n} \mathcal{Y}^{(n+1)/n} \biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{[(n-3)+2(n+1)]/(n+1)} 
\biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{[(n-3)+3(n+1)]/2n} \frac{\mathcal{A}_{M_\ell}^{[(n-3)+3(n+1)]/2n}}{\mathcal{B}_{M_\ell}^{[(n-3)+2(n+1)]/(n+1)}}\biggr\} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl\{\mathcal{X}^4 + \biggl[\frac{3(n+1)}{4\pi n}\biggr] \mathcal{A}_{M_\ell} \mathcal{Y}^2 \biggr\} 
\biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{4n/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{2} \frac{\mathcal{A}_{M_\ell}^{2}}{\mathcal{B}_{M_\ell}^{4n/(n+1)}} \biggr\}</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
- \mathcal{X}^{(n-3)/n} \mathcal{Y}^{(n+1)/n} \biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{(3n-1)/(n+1)} 
\biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{2} \frac{\mathcal{A}_{M_\ell}^{2}}{\mathcal{B}_{M_\ell}^{(3n-1)/(n+1)}}\biggr\} 
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl[ \biggl(\frac{4\pi}{3}\biggr) \mathcal{X}^4 - \mathcal{B}_{M_\ell} \cdot \mathcal{X}^{(n-3)/n} \mathcal{Y}^{(n+1)/n} 
+ \biggl(\frac{n+1}{n}\biggr) \mathcal{A}_{M_\ell} \mathcal{Y}^2 \biggr] 
\biggl\{ \biggl( \frac{4\pi}{3} \biggr)^{(3n-1)/(n+1)} \biggl[ \frac{3(n+1)}{4\pi n} \biggr]^{2} \frac{\mathcal{A}_{M_\ell}^{2}}{\mathcal{B}_{M_\ell}^{4n/(n+1)}} \biggr\} \, .
</math>
  </td>
</tr>
</table>
</div>

Replacing <math>\mathcal{A}_{M_\ell}</math> and <math>\mathcal{B}_{M_\ell}</math> with <math>\mathcal{A}</math> and <math>\mathcal{B}</math>, as prescribed by their defined relationships, the expression inside the square brackets becomes the [[#CompactStahlerVirial|above, boxed-in mass-radius relationship]], namely,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~
\frac{4\pi}{3} \cdot \mathcal{X}^4  
- \mathcal{B} \cdot \mathcal{X}^{(n-3)/n} \biggl[ \mathcal{Y} \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{limit}} \biggr) \biggr]^{(n+1)/n}  
+~ \mathcal{A} \biggl( \frac{n+1}{n} \biggr) \biggl[ \mathcal{Y} \biggl( \frac{M_\mathrm{tot}}{M_\mathrm{limit}} \biggr) \biggr]^{2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0 \, .</math>
  </td>
</tr>
</table>
</div>

====In Terms of Structural Form-Factors====

Alternatively, replacing <math>~\mathcal{A}_{M_\ell}</math> and <math>~\mathcal{B}_{M_\ell}</math> by their expressions in terms of the structural form factors gives, 

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

<tr>
  <td align="right">
<math>~
\frac{4\pi}{3} \cdot \mathcal{X}^4  
- \mathcal{X}^{(n-3)/n} \mathcal{Y}^{(n+1)/n} \biggl( \frac{4\pi}{3} \biggr)^{-1/n}
\frac{\tilde\mathfrak{f}_A}{\tilde\mathfrak{f}_M^{(n+1)/n}} 
+ \mathcal{Y}^2\biggl( \frac{n+1}{5n} \biggr) \frac{\tilde\mathfrak{f}_W  }{\tilde\mathfrak{f}_M^2} 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0 \, .</math>
  </td>
</tr>
</table>
</div>

Finally, inserting into this relation the [[#PTtable|expressions presented above]] for the structural form-factors, <math>\tilde\mathfrak{f}_M</math> and <math>\tilde\mathfrak{f}_A</math>, namely,

<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\tilde\mathfrak{f}_M</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~ \biggl[ - \frac{3\Theta^'}{\xi} \biggr]_{\tilde\xi} </math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~
\tilde\mathfrak{f}_A 
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\tilde\theta^{n+1} + \biggl[ \frac{(n+1)}{3\cdot 5} \biggr] \tilde\xi^2 \cdot \tilde\mathfrak{f}_W
</math>
  </td>
</tr>
</table>
</div>

gives us the desired,
<div align="center" id="ConciseVirialXY">
<table border="1" cellpadding="8" align="center">
<tr><td align="center">
<font color="#770000">'''Virial Theorem written in terms of <math>~\mathcal{X}</math>, <math>~\mathcal{Y}</math>, and <math>~\tilde\mathfrak{f}_W</math>'''</font><br />

<math>~
4\pi \cdot \mathcal{X}^4   ~- ~
\mathcal{X}^{(n-3)/n} \mathcal{Y}^{(n+1)/n} ( 4\pi)^{-1/n}
\biggl[\frac{\tilde\xi}{(-\tilde\theta^')}\biggr]^{(n+1)/n} \biggl[\tilde\theta^{n+1} +  \frac{(n+1)\tilde\xi^2}{3\cdot 5} \cdot \tilde\mathfrak{f}_W \biggr]
~+ ~ \mathcal{Y}^2\biggl( \frac{n+1}{3\cdot 5n} \biggr) \frac{\tilde\xi^2}{(- \tilde\theta^')^2} \cdot \tilde\mathfrak{f}_W 
= 0 \, .
</math>
</td></tr>
</table>
</div>

===Relating and Reconciling Two Mass-Radius Relationships for n = 5 Polytropes===

Now, let's examine the case of pressure-truncated, <math>n=5</math> polytropes.  As we have discussed in the context of [[SSC/Structure/PolytropesEmbedded#Tabular_Summary_.28n.3D5.29|detailed force-balanced models]], [http://adsabs.harvard.edu/abs/1983ApJ...268..165S Stahler (1983)] has deduced that all <math>n=5</math> equilibrium configurations obey the mass-radius relationship,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^2 - 5 \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)\biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)
+ \frac{2^2 \cdot 5 \pi}{3} \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \biggr)^4 </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0 \, ,</math>
  </td>
</tr>
</table>
</div>
where, as [[SSC/Virial/PolytropesSummary#Detailed_Force-Balanced_Solution_2|reviewed above]], the mass and radius normalizations, <math>M_\mathrm{SWS}</math> and <math>~R_\mathrm{SWS}</math>, may be treated as constants once the parameters <math>K</math> and <math>P_e</math> are specified.  In contrast to this, the mass-radius relationship that we have just derived ''from the virial theorem'' for pressure-truncated, <math>~n=5</math> polytropes is,
<div align="center">
<math>
\biggl( \frac{M_\mathrm{limit}}{M_\mathrm{mod}} \biggr)^2
- \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{mod}} \biggr)^{2/5} \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{mod}} \biggr)^{6/5} 
+ \biggl( \frac{R_\mathrm{eq}}{R_\mathrm{mod}} \biggr)^4 = 0 \, ,
</math>
</div>
where the mass and radius normalizations,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~M_\mathrm{mod}\biggr|_{n=5}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~M_\mathrm{SWS} \biggl( \frac{3\mathcal{B}_{M_\ell}}{4\pi} \biggr)^{5/3} \biggl[ \frac{2\cdot 5\pi}{3^2 \mathcal{A}_{M_\ell}} \biggr]^{3/2} 
\, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~R_\mathrm{mod}\biggr|_{n=5}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~R_\mathrm{SWS} \biggl( \frac{3\mathcal{B}_{M_\ell}}{4\pi}\biggr)^{5/6} \biggl[ \frac{2\cdot 5\pi}{3^2\mathcal{A}_{M_\ell}} \biggr]^{1/2} \, ,</math>
  </td>
</tr>
</table>

depend, not only on <math>~K</math> and <math>~P_e</math> via the definitions of <math>~M_\mathrm{SWS}</math> and <math>~R_\mathrm{SWS}</math>, but also on the structural form factors via the free-energy coefficients, <math>~\mathcal{A}_{M_\ell}</math> and <math>~\mathcal{B}_{M_\ell}</math>.  While these two separate mass-radius relationships are similar, they are not identical.  In particular, the middle term involving the cross-product of the mass and radius contains different exponents in the two expressions.  It is not immediately obvious how the two different polynomial expressions can be used to describe the same physical sequence.

This apparent discrepancy is reconciled as follows:  The structural form factors &#8212; and, hence, the free-energy coefficients &#8212; vary from equilibrium configuration to equilibrium configuration.  So it does not make sense to discuss ''evolution along the sequence'' that is defined by the second of the two polynomial expressions.  If you want to know how a given system's equilibrium radius will change ''as its mass changes'', the first of the two polynomials will do the trick.  However, the equilibrium radius of ''a given system'' can be found by looking for extrema in the free-energy function while holding the free-energy coefficients, <math>~\mathcal{A}_{M_\ell}</math> and <math>~\mathcal{B}_{M_\ell}</math>, constant; more importantly, the relative stability ''of a given equilibrium system'' can be determined by analyzing the behavior of the system's free energy ''while holding the free-energy coefficients constant''.  Dynamically stable versus dynamically unstable configurations can be readily distinguished from one another along the sequence that is defined by the second polynomial expression; they cannot be readily distinguished from one another along the sequence that is defined by the first polynomial expression.  It is useful, therefore, to determine how to map a configuration's position on one of the sequences to the other. 

====Plotting Stahler's Relation====

<span id="MRplot">[[File:CorrectedStahlerN5.png|thumb|300px|Mass (Y) vs. radius (X) plot]]Switching, again, to the shorthand notation,</span>
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathcal{X}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{R_\mathrm{eq}}{R_\mathrm{SWS}} \, ,</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathcal{Y}</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{M_\mathrm{tot}}{M_\mathrm{SWS}} \, ,</math>
  </td>
</tr>
</table>
</div>
the equilibrium mass-radius relation defined by the first of the two polynomial expressions can be plotted straightforwardly in either of two ways.  

=====Quadratic Equation=====
One way is to recognize that the polynomial is a quadratic equation whose solution is,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathcal{Y}_\pm</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{5}{2} \mathcal{X} \biggl\{ 1 \pm \biggl[ 1 - \biggl( \frac{2^4\cdot \pi}{3\cdot 5} \biggr) \mathcal{X}^2 \biggr]^{1/2} \biggr\} \, .</math>
  </td>
</tr>
</table>
</div>
In the figure shown here on the right &#8212; see also the bottom panel of [[SSC/Structure/PolytropesEmbedded#Stahler1983Fig17|Figure 2 in our accompanying discussion of detailed force-balance models]] &#8212; Stahler's mass-radius relation has been plotted using the solution to this quadratic equation; the green segment of the displayed curve was derived from the ''positive'' root while the segment derived from the ''negative'' root is shown in orange.  The two curve segments meet at the maximum value of the normalized equilibrium radius, namely, at
<div align="center">
<math>\mathcal{X}_\mathrm{max} \equiv \biggl[ \frac{3\cdot 5}{2^4 \pi} \biggr]^{1/2} \approx 0.54627 \, .</math>
</div>
We note that, when <math>~\mathcal{X} = \mathcal{X}_\mathrm{max}</math>, <math>~\mathcal{Y} = (5\mathcal{X}_\mathrm{max}/2) \approx 1.36569</math>.  Along the entire sequence, the maximum value of <math>~\mathcal{Y}</math> occurs at the location where <math>~d\mathcal{Y}/d\mathcal{X} = 0</math> along the segment of the curve corresponding to the ''positive'' root.  This occurs along the upper segment of the curve where <math>~\mathcal{X}/\mathcal{X}_\mathrm{max} = \sqrt{3}/2</math>, at the location,
<div align="center">
<math>\mathcal{Y}_\mathrm{max} \equiv \biggl[ \frac{3^3 \cdot 5^2}{2^6 } \biggr]^{1/2} \mathcal{X}_\mathrm{max} 
= \biggl[ \frac{3^4 \cdot 5^3}{2^{10} \pi } \biggr]^{1/2}  \approx 1.77408 \, .</math>
</div>

=====Parametric Relations=====
The other way is to determine the normalized mass and normalized radius individually through Stahler's pair of parametric relations.  Drawing partly from our [[SSC/Virial/PolytropesSummary#Detailed_Force-Balanced_Solution_2|above discussion]] and partly from a separate discussion where we provide a [[SSC/Structure/PolytropesEmbedded#Tabular_Summary_.28n.3D5.29|tabular summary of the properties of pressure-truncated <math>n=5</math> polytropes]], these are,
<div align="center">
<table border="0" cellpadding="3">

<tr>
  <td align="right">
<math>
~\mathcal{X}\biggr|_{n=5}
</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{5}{4\pi} \biggr)^{1/2} \tilde\xi \tilde\theta^{2} =
\biggl\{ \frac{3\cdot 5}{2^2 \pi} \biggl[ \frac{\tilde\xi^2/3}{(1+\tilde\xi^2/3)^{2}} \biggr] \biggr\}^{1/2} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
~\mathcal{Y}\biggr|_{n=5}
</math>
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{5^3}{4\pi} \biggr)^{1/2} \tilde\theta (- \tilde\xi^2 \tilde\theta^') =
\biggl[  \biggl( \frac{3 \cdot 5^3}{2^2\pi} \biggr) \frac{(\tilde\xi^2/3)^3}{(1+\tilde\xi^2/3)^{4}} \biggr]^{1/2} \, .
</math>
  </td>
</tr>
</table>
</div>
The entire sequence will be traversed by varying the Lane-Emden parameter, <math>~\tilde\xi</math>, from zero to infinity.  Using the first of these two expressions, we have determined, for example, that the point along the sequence corresponding to the maximum normalized equilibrium radius, <math>~\mathcal{X}_\mathrm{max}</math>, is associated with an embedded <math>~n=5</math> polytrope whose truncated, dimensionless Lane-Emden radius is,
<div align="center">
<math>
~\tilde\xi \biggr|_{\mathcal{X}_\mathrm{max}} = 3^{1/2} \, .
</math>
</div>
Similarly, we have determined that the point along the sequence that corresponds to the maximum dimensionless mass, <math>~\mathcal{Y}_\mathrm{max}</math>, is associated with an embedded <math>~n=5</math> polytrope whose truncated, dimensionless Lane-Emden radius is, precisely,
<div align="center">
<math>
~\tilde\xi \biggr|_{\mathcal{Y}_\mathrm{max}} = 3 \, .
</math>
</div>

<div id="KimuraApplication">
<table border="1" width="90%" align="center" cellpadding="8">
<tr><td align="left">
Referring back to [[SSC/Structure/PolytropesEmbedded#Turning_Points|our review of turning points]] along equilibrium sequences and, especially, the work of {{ Kimura81bfull }}, we appreciate that the point that corresponds to the maximum mass, <math>~\mathcal{Y}_\mathrm{max}</math>, is the turning point that Kimura refers to as the "extremum in M<sub>1</sub>" along a p<sub>1</sub> sequence.  [[SSC/Structure/PolytropesEmbedded#Location_of_Pressure_Limit|As we have highlighted]], according to Kimura, this point should occur along the sequence where <math>~h_G=0</math>, that is, where the following condition applies:
<div align="center">
<math>~\frac{\tilde\theta^{n+1}}{(\tilde\theta^')^2} = \frac{(n-3)}{2}  \, .</math>
</div>
For the specific case being studied here, namely, <math>~n = 5</math> polytropic configurations, we therefore expect from Kimura's work that <math>~[\tilde\theta^6/(\tilde\theta^')^2] = 1</math> at the "maximum mass" turning point.  Given that,
<div align="center">
<math>~\tilde\xi \biggr|_{\mathcal{Y}_\mathrm{max}} = 3</math>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>~\Rightarrow</math>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>\tilde\theta_{n=5} = \frac{1}{2}</math>&nbsp;&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;&nbsp; <math>\tilde\theta_{n=5}^' = -\frac{1}{8} \, ,</math>
</div>
we see that Kimura's condition holds and, hence, that our identification of the location along the sequence of the maximum mass matches Kimura's identification of the location of that turning point.

We appreciate, as well, that the point corresponding to the maximum normalized equilibrium radius, <math>\mathcal{X}_\mathrm{max}</math>, is the turning point that Kimura would reference as the "extremum in r<sub>1</sub>" along a p<sub>1</sub> sequence.  Following Kimura's analysis [[SSC/Structure/PolytropesEmbedded#TurningPointXmax|we have shown that this point occurs along the sequence where the following condition applies]]:
<div align="center">
<math>~\frac{\xi (-\theta^')}{\tilde\theta} = \frac{2}{(n-1)} \, ,</math>
</div>
that is, for the specific case being studied here, we should expect <math>~[\tilde\xi (-\tilde\theta^')/\tilde\theta] = 1/2</math> at the "maximum radius" turning point.  Given that,
<div align="center">
<math>~\tilde\xi \biggr|_{\mathcal{X}_\mathrm{max}} = 3^{1/2}</math>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>~\Rightarrow</math>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<math>\tilde\theta_{n=5} = 2^{-1/2}</math>&nbsp;&nbsp;&nbsp;&nbsp; and &nbsp;&nbsp;&nbsp;&nbsp; 
<math>\tilde\theta_{n=5}^' = -(2^3 \cdot 3)^{-1/2} \, ,</math>
</div>
we see that Kimura's condition holds and, hence, that our identification of the location along the sequence of the maximum radius matches Kimura's identification of the location of that turning point.
</td></tr>
</table>
</div>

=====Discussion=====
Notice that if the quadratic equation is used to map out the mass-radius relationship, the parameter, <math>~\tilde\xi</math>, never explicitly enters the discussion.  Instead, a radius <math>~0 \le \mathcal{X} \le \mathcal{X}_\mathrm{max}</math> is specified and the ''two'' equilibrium masses associated with <math>~\mathcal{X}</math> &#8212; call them, <math>~\mathcal{Y}_+</math> and <math>~\mathcal{Y}_-</math> &#8212; are determined.  (The values of the two masses are degenerate at both limiting values of <math>~\mathcal{X}</math>.)  If the pair of parametric relations is used, instead, only ''one'' value of the mass is obtained for each specified value of <math>~\tilde\xi</math>.  As <math>~\tilde\xi</math> is increased from <math>~0</math> to <math>~\sqrt{3}</math>, <math>~\mathcal{X}</math> increases monotonically from <math>~0</math> to <math>~\mathcal{X}_\mathrm{max}</math> and the corresponding mass is (only) <math>~\mathcal{Y}_-</math>; that is, as <math>~\tilde\xi</math> is increased from <math>~0</math> to <math>~\sqrt{3}</math>, we move away from the origin in a counter-clockwise direction along the lower segment (colored orange in the [[#MRplot|above figure]]) of the plotted equlibrium sequence.  Then, as <math>~\tilde\xi</math> is increased from <math>~\sqrt{3}</math> to <math>~\infty</math>, we continue to move in a counter-clockwise direction along the equilibrium sequence, but now along the upper segment (colored green in the [[#MRplot|above figure]]) of the sequence, back to the origin; that is to say, <math>~\mathcal{X}</math> steadily decreases from <math>~\mathcal{X}_\mathrm{max}</math> back to <math>~0</math> and this time the relevant associated mass is the positive root of the quadratic relation, <math>~\mathcal{Y}_+</math>. 

Clearly, then, each value of <math>~\mathcal{X}</math> is associated with two different values of the parametric parameter, <math>~\tilde\xi</math>.  By inverting the <math>~\mathcal{X}(\tilde\xi)</math> parametric expression we see that, the two values of <math>~\tilde\xi</math> associated with a given equilibrium radius are,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\tilde\xi_\pm</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl\{\frac{3}{\alpha} \biggl[ 1 \pm \sqrt{1 - \alpha^2} \biggr]  \biggr\}^{1/2} \, ,</math>
  </td>
</tr>
</table>
</div>

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

<tr>
  <td align="right">
<math>~\alpha</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{(\mathcal{X}/\mathcal{X}_\mathrm{xmax})^2}{2-(\mathcal{X}/\mathcal{X}_\mathrm{xmax})^2} \, .</math>
  </td>
</tr>
</table>
</div>
We note as well that, for a given equilibrium radius, <math>~\mathcal{X}</math>, the ''ratio'' of the two mass solutions is given by a very simple expression, namely,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\frac{\mathcal{Y}_-}{\mathcal{Y}_+} = \frac{\tilde\xi_-^2}{3}</math>
  </td>
  <td align="center">
&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; or &nbsp; &nbsp; &nbsp; &nbsp; 
  </td>
  <td align="left">
<math>~\frac{\mathcal{Y}_+}{\mathcal{Y}_-} = \frac{\tilde\xi_+^2}{3} \, .</math>
  </td>
</tr>
</table>
</div>
This implies, as well, that,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\tilde\xi_+ \cdot \tilde\xi_-</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~3\, .</math>
  </td>
</tr>
</table>
</div>

====Plotting the Virial Theorem Relation====

[[#ConciseFreeEnergyExpression|Drawing from our above derivations, the concise free-energy expression]] that reflects the properties of pressure-truncated <math>n = 5</math> polytropic configurations is,
<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>~
- \frac{18}{5}\cdot \mathcal{A}_{M_\ell}  \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{2} \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-1} 
+~ 5\mathcal{B}_{M_\ell}  \cdot \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{6/5}  \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-3/5} 
+~ \frac{4\pi}{3} \cdot \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{3}  \, , 
</math>
  </td>
</tr>
</table>
</div>
where,
<div align="center">
<table border="0" cellpadding="5">
<tr>
  <td align="right">
<math>~\mathcal{A}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{1}{5} \cdot \frac{\tilde\mathfrak{f}_W}{\tilde\mathfrak{f}_M^2}  
=\frac{1}{3^2\cdot 5} \biggl( \frac{\tilde\xi}{-\tilde\theta^'} \biggr)^2 \tilde\mathfrak{f}_W \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathcal{B}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( \frac{3}{4\pi}\biggr)^{1/5} \frac{\tilde\mathfrak{f}_A}{\tilde\mathfrak{f}_M^{6/5}} 
= \frac{1}{3}( 4\pi)^{-1/5} \biggl( \frac{\tilde\xi}{-\tilde\theta^'} \biggr)^{6/5}
\biggl[ \tilde\theta^6 + \frac{2}{5} \cdot \tilde\xi^2 \tilde\mathfrak{f}_W \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>

The virial theorem which is derived from this free-energy expression provides a mass-radius relationship to be compared with the detailed force-balance relationship presented by Stahler.
Because our intent is to make this comparison, we begin with the virial theorem as written in terms of the variables, <math>~\mathcal{X}</math> and <math>~\mathcal{Y}</math>, and specialized for the case of <math>~n = 5</math> polytropic configurations.  Written in terms of the (constant) coefficients in the free-energy expression, we have
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~
\mathcal{X}^4  
- \frac{3\mathcal{B}_{M_\ell}}{4\pi} \cdot ( \mathcal{X} \mathcal{Y}^3 )^{2/5}  
+~ \frac{9 \mathcal{A}_{M_\ell}}{10\pi}\cdot \mathcal{Y}^{2}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~0 \, ;</math>
  </td>
</tr>
</table>
</div>
or, [[#ConciseVirialXY|from above]], the
<div align="center">
<table border="0" cellpadding="8" align="center">
<tr><td align="center">
<font color="#770000">'''Virial Theorem written in terms of <math>~\mathcal{X}</math>, <math>~\mathcal{Y}</math>, and <math>~\tilde\mathfrak{f}_W</math>'''</font><br />

<math>
\mathcal{X}^4   ~- ~
(\mathcal{X} \mathcal{Y}^3)^{2/5} 
\biggl[\frac{\tilde\xi}{4\pi (-\tilde\theta^')}\biggr]^{6/5} \biggl[\tilde\theta^{6} +  \frac{2}{5} \cdot \tilde\xi^2 \tilde\mathfrak{f}_W \biggr]
~+ ~ \mathcal{Y}^2 \biggl[\frac{\tilde\xi}{(- \tilde\theta^')} \biggr]^2 \frac{\tilde\mathfrak{f}_W}{2\cdot 5^2 \pi} 
= 0 \, ,
</math>
</td></tr>
</table>
</div>
where, specifically for <math>~n = 5</math> polytropic configurations &#8212; see our [[SSC/Structure/PolytropesEmbedded#Tabular_Summary_.28n.3D5.29|summary of the radial profiles of physical variables]] and our [[SSCpt1/Virial/FormFactors#Summary_.28n.3D5.29|determination of expressions for the structural form-factors]],
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\tilde\theta</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~  
( 1+ \ell^2 )^{-1/2} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\tilde\theta^'</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~  
- 3^{-1/2} \ell ( 1+ \ell^2 )^{-3/2} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathfrak{f}_W</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~  
\frac{5}{2^4} \cdot \ell^{-5}  
\biggl[ \ell \biggl( \ell^4 - \frac{8}{3}\ell^2 - 1 \biggr)(1 + \ell^2)^{-3} + \tan^{-1}(\ell ) \biggr] \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\ell^2</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~  
\frac{\tilde\xi^2}{3} \, .
</math>
  </td>
</tr>
</table>
</div>

Once numerical values have been assigned to the free-energy coefficients, <math>~\mathcal{A}_{M_\ell}</math> and <math>~\mathcal{B}_{M_\ell}</math>, the mass-radius relationship given by the scalar virial theorem can be compared quantitatively with Stahler's (detailed force-balance) mass-radius relationship.  The simplest, physically reasonable approximation would be to assume uniform-density structures, in which case, <math>~\tilde\mathfrak{f}_M = \tilde\mathfrak{f}_W = \tilde\mathfrak{f}_A = 1</math>, and accordingly, <math>~\mathcal{A}_{M_\ell} = 5^{-1}</math> and <math>~\mathcal{B}_{M_\ell} = (4\pi/3)^{-1/5}</math>.  But a better approximation would be to assign values to the structural form-factors that properly represent the properties of at least one detailed force-balanced model.  By way of illustration, the following table details what the proper values are for the two free-energy coefficients, and other relevant parameters, specifically for the model along Stahler's sequence that sits at <math>~\mathcal{Y}_\mathrm{max}</math> &#8212; that is, the model whose truncation radius is <math>~\tilde\xi = 3</math>.  As is recorded in the table, in this case the precise values of the free-energy coefficients are,
<div align="center">
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>~\mathcal{A}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~  
\biggl( \frac{2^4 \pi^2}{3^7} \biggr)^{1/2} \, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>~\mathcal{B}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~  
\biggl( \frac{3}{2^{14}\pi} \biggr)^{1/5} \biggl[  1 +
\biggl( \frac{2^6\pi^2}{3^3} \biggr)^{1/2}\biggr] \, .
</math>
  </td>
</tr>
</table>
</div>
Notice that, by choosing <math>~\tilde\xi = 3</math>, the evaluation of <math>~\tilde\mathfrak{f}_W</math> is particularly simple, in part, because <math>~\tan^{-1}(\ell) = \tan^{-1} \sqrt{3} = \pi/3</math>, but also because the term <math>~(\ell^4 - 8\ell^2/3 - 1)</math> equals zero.


<div align="center" id="ExampleXi3">
<table border="1" align="center" cellpadding="8">
<tr>
  <th align="center" colspan="3">
Determination of Coefficient Values in the Specific Case of <math>~\tilde\xi = 3</math>
  </th>
</tr>
<tr>
  <td align="center">
Quantity
  </td>
  <td align="center">
Analytic Evaluation
  </td>
  <td align="center">
Numerical
  </td>
</tr>
<tr>
  <td align="center">
<math>~\ell</math>
  </td>
  <td align="center">
<math>~3^{1/2}</math>
  </td>
  <td align="left">
<math>~1.732051</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\tilde\theta</math>
  </td>
  <td align="center">
<math>~2^{-1}</math>
  </td>
  <td align="left">
<math>~0.5</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\tilde\theta^'</math>
  </td>
  <td align="center">
<math>~2^{-3}</math>
  </td>
  <td align="left">
<math>~0.125</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\mathcal{X}</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{3^2\cdot 5}{2^6\pi} \biggr)^{1/2}</math>
  </td>
  <td align="left">
<math>~0.473087</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\mathcal{Y}</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{3^4\cdot 5^3}{2^{10}\pi} \biggr)^{1/2}</math>
  </td>
  <td align="left">
<math>~1.774078</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\tilde\mathfrak{f}_W</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{5^2 \pi^2}{2^8\cdot 3^7} \biggr)^{1/2}</math>
  </td>
  <td align="left">
<math>~0.020993</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\mathcal{A}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{2^4 \pi^2}{3^7} \biggr)^{1/2}</math>
  </td>
  <td align="left">
<math>~0.268711</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~\mathcal{B}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{3}{2^{14}\pi} \biggr)^{1/5} \biggl[  1 +
\biggl( \frac{2^6\pi^2}{3^3} \biggr)^{1/2}\biggr]</math>
  </td>
  <td align="left">
<math>~0.830395</math>
  </td>
</tr>

<tr>
  <td align="center">
<math>~G^*</math>
  </td>
  <td align="center">
<math>~\biggl( \frac{3\cdot 5^3}{2^{12}\pi} \biggr)^{1/2} [  2^3\pi + 3^{5/2} ]</math>
  </td>
  <td align="left">
<math>~6.951544</math>
  </td>
</tr>

<tr>
  <td align="left">
Virial: 

<math>~\mathcal{X}^4 - \frac{3\mathcal{B}_{M_\ell}}{4\pi} \cdot (\mathcal{X}\mathcal{Y}^3 )^{2/5}</math>

&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>~~~~~~~~~~+\frac{9 \mathcal{A}_{M_\ell}}{10\pi} \cdot \mathcal{Y}^2</math>
  </td>
  <td align="left">
<math>~ \frac{3^4\cdot 5^2}{2^{12}\pi} - 
\frac{3^4\cdot 5^2}{2^{12}\pi} \biggl[  1 + \biggl( \frac{2^6\pi^2}{3^3} \biggr)^{1/2}\biggr]</math>

&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; <math>~~~~~~~~~~+ \biggl( \frac{3^5\cdot 5^4}{2^{18} \pi^2} \biggr)^{1/2}
</math>
  </td>
  <td align="left">
Sums to zero, 

exactly!
  </td>
</tr>
</table>
</div>


The curve traced out by the light-blue diamonds in each panel of the following comparison figure displays Stahler's analytically prescribed mass-radius relation; this curve is identical in all six panels and is the same as the [[#Plotting_Stahler.27s_Relation|curve displayed above]] in connection with our description of Stahler's mass-radius relation.  Each point along this "Stahler" curve identifies a model having a different truncation radius, <math>~\tilde\xi</math>; as plotted here, starting near the origin and moving counter-clockwise around the curve, <math>~\tilde\xi</math> is varied from 0.05 to 42.5.  As foreshadowed by the above discussion, the model having the greatest mass <math>~(\mathcal{Y}_\mathrm{max})</math> along the Stahler sequence &#8212; highlighted by the red filled circle in most of the figure panels &#8212;  is defined by <math>~\tilde\xi = 3</math>.  

In each figure panel, the curve traced out by the orange triangles &#8212; or, in one case, the orange triangles &amp; light purple diamonds &#8212; displays the mass-radius relation defined by the virial theorem.  These "Virial" curves are all defined by the same virial theorem polynomial expression, as just presented, but the coefficient of the <math>~\mathcal{Y}^2</math> term and the coefficient of the <math>~(\mathcal{X}\mathcal{Y}^3)^{2/5}</math> cross term &#8212; essentially, the value of <math>~\mathcal{A}_{M_\ell}</math> and the value of <math>~\mathcal{B}_{M_\ell}</math>, respectively &#8212; have different values in the six separate figure panels.  In each case, a value has been specified for the parameter, <math>~\tilde\xi</math> (as identified in the title of each figure panel), and this, in turn, has determined the values of the two (constant) free-energy coefficients.  For example, in the top-right figure panel whose title indicates <math>~\tilde\xi = 3</math>, the "Virial" curve traces the mass-radius relation prescribed by the virial theorem after the values of the free-energy coefficients have been set to values that correspond to a detailed force-balanced model with this specified truncation radius, that is (see the above table),  <math>~\mathcal{A}_{M_\ell} = 0.268711</math> and <math>~\mathcal{B}_{M_\ell} = 0.830395</math>.  Columns 2 and 3, respectively, of the table affixed to the bottom of the following comparison figure list the values of <math>~\mathcal{A}_{M_\ell}</math> and <math>~\mathcal{B}_{M_\ell}</math> that have been used to define the "Virial" curve in each of the six figure panels, in accordance with the value of <math>~\tilde\xi</math> listed in column 1 of the table.


<div align="center" id="OverlapPlots">
<table border="1" cellpadding="8" align="center">
<tr>
  <th align="center" colspan="7">
Comparing Two Separate Mass-Radius Relations for Pressure-Truncated ''n = 5'' Polytropes
  </th>
</tr>
<tr>
<td align="center" colspan="3">
[[File:Xi1.26419.png|225px|Comparison of Two Mass-Radius Relations]]
</td>
<td align="center" colspan="2">
[[File:Xi2.0.png|225px|Comparison of Two Mass-Radius Relations]]
</td>
<td align="center" colspan="2">
[[File:Xi3.0.png|225px|Comparison of Two Mass-Radius Relations]]
</td>
</tr>
<tr>
<td align="center" colspan="3">
[[File:Xi3.5.png|225px|Comparison of Two Mass-Radius Relations]]
</td>
<td align="center" colspan="2">
[[File:Xi3.850652.png|225px|Comparison of Two Mass-Radius Relations]]
</td>
<td align="center" colspan="2">
[[File:Xi9.8461.png|225px|Comparison of Two Mass-Radius Relations]]
</td>
</tr>

<tr>
  <td align="center" rowspan="2">
<math>~\tilde\xi</math>
  </td>
  <td align="center" colspan="2">
Free-Energy Coefficients
  </td>
  <td align="center" colspan="2">
Primary Overlap
  </td>
  <td align="center" colspan="2">
Secondary Overlap

(see [[#SecondaryOverlap|further elaboration below]])
  </td>
</tr>

<tr>
  <td align="center">
<math>~\mathcal{A}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~\mathcal{B}_{M_\ell}</math>
  </td>
  <td align="center">
<math>~\mathcal{X} \equiv \frac{R_\mathrm{eq}}{R_\mathrm{SWS}}</math>
  </td>
  <td align="center">
<math>~\mathcal{Y} \equiv \frac{M_\mathrm{limit}}{M_\mathrm{SWS}}</math>
  </td>
  <td align="center">
<math>~\mathcal{X}</math>
  </td>
  <td align="center">
<math>~\mathcal{Y}</math>
  </td>
</tr>

<tr>
<td align="center">
1.26419
</td>
<td align="center">
0.214429
</td>
<td align="center">
0.758099
</td>
<td align="center" bgcolor="yellow">
0.520269
</td>
<td align="center" bgcolor="yellow">
0.904142
</td>
<td align="center">
0.388938
</td>
<td align="center">
0.289568
</td>
</tr>

<tr>
<td align="center">
2
</td>
<td align="center">
0.233842
</td>
<td align="center">
0.779836
</td>
<td align="center">
0.540671
</td>
<td align="center">
1.544775
</td>
<td align="center">
0.468918
</td>
<td align="center">
0.570916
</td>
</tr>

<tr>
<td align="center">
[[#ExampleXi3|3]]
</td>
<td align="center">
0.268711
</td>
<td align="center">
0.830395
</td>
<td align="center">
0.473087
</td>
<td align="center">
1.774078
</td>
<td align="center">
0.507387
</td>
<td align="center">
0.798441
</td>
</tr>

<tr>
<td align="center">
3.5
</td>
<td align="center">
0.288708
</td>
<td align="center">
0.861503
</td>
<td align="center">
0.434310
</td>
<td align="center">
1.744359
</td>
<td align="center">
0.515168
</td>
<td align="center">
0.859518
</td>
</tr>

<tr>
<td align="center">
3.850652
</td>
<td align="center">
0.303490
</td>
<td align="center">
0.884746
</td>
<td align="center">
0.408738
</td>
<td align="center">
1.699778
</td>
<td align="center">
0.518588
</td>
<td align="center">
0.888969
</td>
</tr>

<tr>
<td align="center">
9.8461
</td>
<td align="center">
0.601012
</td>
<td align="center">
1.313904
</td>
<td align="center">
0.186424
</td>
<td align="center">
0.904141
</td>
<td align="center" bgcolor="yellow">
0.520269
</td>
<td align="center" bgcolor="yellow">
0.904143
</td>
</tr>
</table>
</div>


As is [[#Quantitative_Study|discussed more fully, below]], in each of the six panels of the above comparison figure, the "Virial" curve intersects the "Stahler" curve at two locations.  These points of intersection are identified by the black filled circles in each figure panel.  In each case, the intersection point that is farthest along on the Stahler sequence &#8212; as determined by starting at the origin and moving counter-clockwise along the sequence &#8212; identifies the "Primary Overlap" between the two curves.  That is to say, the <math>~(\mathcal{X}, \mathcal{Y})</math> coordinates of this point (see columns 4 and 5 of the table affixed to the bottom of the figure) are the coordinate values that are obtained by plugging the specified value of <math>\tilde\xi</math> (see the title of the figure panel or column 1 of the affixed table) into [[#Parametric_Relations|Stahler's pair of parametric relations]].  The second point of intersection in each panel &#8212; which we will refer to as the "Secondary Overlap" points and whose coordinates are provided in columns 6 and 7 of the affixed table &#8212; appears to be fortuitous and of no particularly significant astrophysical interest.


The mass-radius diagram displayed in the top-right panel of the [[#OverlapPlots|above comparison figure]] has been reproduced in the upper-left panel of the [[#GraphicalDepictionXi3|following figure]] &#8212; in this case, with a coordinate aspect ratio that is closer to 1:1  &#8212; along with color images of the corresponding free-energy surface, viewed from two different perspectives, and a three-column table listing the 3D coordinates, <math>~(X,Y,Z) = (R_\mathrm{eq}/R_\mathrm{SWS}, M_\mathrm{limit}/M_\mathrm{SWS}, \mathfrak{G}^*)</math>, of the seventeen points that have been used to define the displayed "Virial" curve.  To be more explicit, the rainbow-colored free-energy surface, <math>~\mathfrak{G}^*(R,M_\mathrm{limit})</math>, has been defined by the free-energy function appropriate to pressure-truncated <math>~n=5</math> polytropic configurations [[#Plotting_the_Virial_Theorem_Relation|as defined above]], that is,
<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>~
- \frac{18}{5}\cdot \mathcal{A}_{M_\ell}  \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{2} \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-1} 
+~ 5\mathcal{B}_{M_\ell}  \cdot \biggl( \frac{M_\mathrm{limit}}{M_\mathrm{SWS}} \biggr)^{6/5}  \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{-3/5} 
+~ \frac{4\pi}{3} \cdot \biggl(\frac{R}{R_\mathrm{SWS}}\biggr)^{3}  \, , 
</math>
  </td>
</tr>
</table>
</div>

with the values of the two free-energy coefficients set to the values that correspond to a <math>~\tilde\xi = 3</math> virial curve as discussed above, namely, <math>~\mathcal{A}_{M_\ell} = 0.268711</math> and <math>~\mathcal{B}_{M_\ell} = 0.830395</math>.


<div align="center" id="GraphicalDepictionXi3">
<table border="1" cellpadding="8" align="center">
<tr>
  <th align="center" colspan="2">
Free-Energy Surface for ''n = 5'' Polytropic Configurations Truncated at <math>~\tilde\xi = 3</math>
  </th>
  <th align="center" colspan="1">
Radius
  </th>
  <th align="center" colspan="1">
Mass
  </th>
  <th align="center" colspan="1">
Free

Energy
  </th>
</tr>
<tr>
<td align="center" rowspan="8">
[[File:VirialEquilibria.png|150px|Virial Mass-Radius Relation]]
</td>
<td align="center" rowspan="8" bgcolor="CCFFFF">
[[File:FreeEnergyFaceOn.png|325px|Virial Mass-Radius Relation on top of Free-Energy Surface]]
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
0.3152
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
0.1345
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
0.8231
</td>
</tr>

<tr>
<td align="right" rowspan="1">
0.4159
</td>
<td align="right" rowspan="1">
0.3470
</td>
<td align="right" rowspan="1">
1.9947
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.4731
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
0.5735
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
3.1095
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="white">
0.5035
</td>
<td align="right" rowspan="1" bgcolor="white">
0.7661
</td>
<td align="right" rowspan="1" bgcolor="white">
3.9591
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="white">
0.5114
</td>
<td align="right" rowspan="1" bgcolor="white">
0.8347
</td>
<td align="right" rowspan="1" bgcolor="white">
4.2409
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.5297
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
1.0759
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
5.1453
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.5310
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
1.4036
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
6.1560
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="white">
0.5250
</td>
<td align="right" rowspan="1" bgcolor="white">
1.5010
</td>
<td align="right" rowspan="1" bgcolor="white">
6.4046
</td>
</tr>

<tr>
<td align="center" colspan="2" rowspan="9" bgcolor="#CCFFFF">
[[File:FreeEnergySurface.png|450px|Free-Energy Surface]]
</td>
<td align="right" rowspan="1" bgcolor="white">
0.5185
</td>
<td align="right" rowspan="1" bgcolor="white">
1.5680
</td>
<td align="right" rowspan="1" bgcolor="white">
6.5606
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.5114
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
1.6206
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
6.6741
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="black"><font color="white">
0.4731</font>
</td>
<td align="right" rowspan="1" bgcolor="black"><font color="white">
1.7741</font>
</td>
<td align="right" rowspan="1" bgcolor="black"><font color="white">
6.9515</font>
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.4343
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
1.8287
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
7.0247
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="red"><font color="white">
0.3984</font>
</td>
<td align="right" rowspan="1" bgcolor="red"><font color="white">
1.8333</font>
</td>
<td align="right" rowspan="1" bgcolor="red"><font color="white">
7.0301</font>
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.3379
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
1.7764
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
6.9923
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.2911
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
1.6909
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
6.9541
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.2260
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
1.5229
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
6.9091
</td>
</tr>

<tr>
<td align="right" rowspan="1" bgcolor="lightblue">
0.1545
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
1.2749
</td>
<td align="right" rowspan="1" bgcolor="lightblue">
6.8786
</td>
</tr>
</table>
</div>

In the bottom panel of this figure, the undulating free-energy surface is drawn in three dimensions and viewed from a vantage point that illustrates its "valley of stability" and "ridge of instability;" the surface color correlates with the value of the free energy.  Twelve small colored dots identify extrema &#8212; either the bottom of a valley or the top of a ridge &#8212; in the free-energy function and therefore trace the mass-radius relation defined by the scalar virial theorem.  The 3D coordinates of these twelve points are provided in the three-column table that is affixed to the righthand edge of the figure:  The coordinates of the (only) red dot are provided in row 13 of the table (red has also been assigned as the "bgcolor" of this table row); the equilibrium configuration having the greatest mass along the intersecting Stahler sequence is identified by the (only) black dot (coordinates are provided in row 11 of the table and, correspondingly, bgcolor="black" for that row); bgcolor="lightblue" has been assigned to the other rows of the affixed table that provide coordinates of the other 10, blue dots.

The upper-right panel of this figure presents the two-dimensional projection that results from viewing the identical free-energy surface "from above," along a line of sight that is parallel to the free-energy <math>~(Z)</math> axis and looking directly down onto the radius-mass <math>~(X-Y)</math> plane. From this vantage point, the twelve small colored dots cleanly trace out the <math>~M_\mathrm{limit}(R_\mathrm{eq})</math> equilibrium sequence that is defined by the scalar virial theorem, exactly reproducing the "Virial" curve that is depicted in the mass-radius diagram shown in the upper-left panel of the figure.