View source for Appendix/Ramblings/MacSphCriticalPoints

__FORCETOC__ 
<!-- __NOTOC__ will force TOC off -->
=Critical Points along the Maclaurin Spheroid Sequence=
{| class="MaclaurinSequence" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
|- 
! style="height: 150px; width: 150px; background-color:maroon;" |[[H_BookTiledMenu#Spheroidal_.26_Spheroidal-Like|<font color="white">Maclaurin<br />Spheroid<br />Sequence</font>]]
|}

The thoughts and recommendations that are presented in this chapter stem from [[Appendix/Ramblings/ForCohlHoward|collaborative discussions with Howard Cohl]] and associated colleagues.


==Introduction==

As has been explained in, for example, our [[ThreeDimensionalConfigurations/RiemannStype#Fig2|accompanying discussion of Riemann S-type ellipsoids]], the "EFE Diagram"  refers to a two-dimensional parameter space defined by the pair of axis ratios (b/a, c/a), ''usually'' covering the ranges, 0 &le; b/a &le; 1 and 0 &le; c/a &le; 1.  The classic/original version of this diagram appears as Figure 2 on p. 902 of {{ Chandrasekhar65_XXV }}; a somewhat less cluttered version appears on p. 147 of Chandrasekhar's [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  

The version of the EFE Diagram shown in the left-hand panel of the following figure, highlights &hellip;
<ul>
  <li>The Maclaurin spheroid sequence &#8212; the vertical line that runs from <math>(b/a, c/a) = (1, 1)</math> (a nonrotating, sphere) to <math>(b/a, c/a) = (1, 0)</math> (an infinitesimally thin disk), aligning with the right-hand boundary of the Diagram;</li>
  <li>The Jacobi ellipsoid sequence &#8212; the curve running through the purple circular markers that connects the origin of the diagram <math>(b/a, c/a) = (0, 0)</math> to the point <math>(b/a, c/a) = (1, 0.582724)</math> where it intersects (and bifurcates from) the Maclaurin spheroid sequence.</li>
  <li>The pair of self-adjoint sequences (USA and LSA) as [[#Self-Adjoint_Sequences|discussed in detail above]] &#8212; the black dot-dashed curves that run from the origin of the diagram to points identified by yellow circular markers where they intersect the Maclaurin spheroid sequence: at, respectively, <math>(b/a, c/a) = (1, 1)</math> for the USA and <math>(b/a, c/a) = (1, 0.30333)</math> for the LSA.</li>
</ul>

To aid subsequent discussion, in this EFE diagram we have broken the Maclaurin spheroid sequence into three differently colored segments: &nbsp; <b>Blue</b> extends from the nonrotating sphere (point of USA intersection) to the point where the Jacobi ellipsoid sequence intersects; <b>Orange</b> extends from the Jacobi intersection point to the point of LSA intersection; and <b>Black</b> extends from the LSA intersection to the end of the Maclaurin sequence (the infinitesimally thin disk).

<table border="1" align="center" cellpadding="5">
<tr>
  <td align="center">
<b>EFE Diagram</b><br />

[[File:OurEFEannotated.png|400px|OurEFE]]
  </td>
  <td align="center">
<b>&Omega;<sup>2</sup> vs. j<sup>2</sup> Diagram</b><br />

[[File:OurHE84Fig1annotated.png|525px|OurHE84Fig1]]
  </td>
</tr>
</table>

Now, as we consider examining the stability of individual models &#8212; or the behavior of equilibrium model sequences &#8212; whose configurations are not constrained to have purely ellipsoidal shapes, the two-dimensional parameter space <math>(b/a, c/a)</math> associated with the EFE diagram proves to be of little use.  In Figures 1 through 6 of an [[Apps/MaclaurinSpheroidSequence|accompanying overview]], we show how the Maclaurin spheroid sequence behaves when displayed in six alternate 2D-parameter-space diagrams that have been proposed/used &#8212; to varying degrees of success &#8212; by various research groups over the years.  Following the [[Apps/EriguchiHachisu/Models|extensive set of related work published by Eriguchi, Hachisu and their collaborators over the decade of the 1980s]], here we adopt the &Omega;<sup>2</sup> versus j<sup>2</sup> diagram as displayed in the right-hand panel of the above figure, where, <math>\Omega^2 \equiv \omega_0^2/(4\pi G \rho)</math>.

The blue, orange, and black segments of the (vertical straight line) Maclaurin spheroid sequence that have been highlighted in our version of the EFE diagram map respectively to the blue, orange, and black segments of the (curved, double-valued) Maclaurin spheroid sequence that appears in our version of the &Omega;<sup>2</sup> versus j<sup>2</sup> diagram.  Also shown (purple circular markers) is the Jacobi ellipsoid sequence, along with its intersection with the Maclaurin spheroid sequence; and (see the yellow circular markers) the points where the USA and LSA sequences intersect the Maclaurin sequence.  As an additional reference point, in the diagram on the right, the small green square marker identifies the [[Apps/MaclaurinSpheroidSequence#OmegaMax|point along the Maclaurin spheroid sequence where &Omega;<sup>2</sup> attains its maximum value]] &#8212; <math>(\Omega^2, j^2)= (0.112333, 0.010105)</math>; in our EFE diagram, a green square marker identifies this same Maclaurin spheroid &#8212; <math>(b/a, c/a) = (1, 0.36769)</math>.

==Model01==

<table border="1" align="center" cellpadding="8" width="80%">
<tr>
  <td align="center">Model</td>
  <td align="center"><math>\frac{b}{a}</math></td>
  <td align="center"><math>\frac{c}{a}</math></td>
  <td align="center"><math>e</math></td>
  <td align="center"><math>\frac{\Omega^2}{\pi G\rho}</math></td>
  <td align="center"><math>\frac{L}{(GM^3\bar{a})^{1 / 2}}</math></td>
  <td align="center"><math>\frac{T_\mathrm{rot}}{|W_\mathrm{grav}|}</math></td>
</tr>
<tr>
  <td align="center">01</td>
  <td align="center"><math>1</math></td>
  <td align="center"><math>\frac{\sqrt{3}}{2} = 0.866025404</math></td>
  <td align="center"><math>\sin(\pi/6) = 0.5</math></td>
  <td align="center"><math>\biggl[\frac{10\pi}{ 3^{1 / 2}} - 18\biggr]</math><br />&nbsp;<br /><math>= 0.137993640</math></td>
  <td align="center"><math>\frac{1}{5}\biggl( \frac{2^{14}}{3}\biggr)^{1 / 12} \biggl[ 5\pi - 3^{5/2}  \biggr]^{1 / 2}</math><br />&nbsp;<br /><math> = 0.141633637</math></td>
  <td align="center"><math>\biggl[\frac{3}{\pi}\biggl( 2\pi - 3^{3/2}\biggr) - 1 \biggr]</math><br />&nbsp;<br /><math>=0.038039942</math></td>
</tr>

<tr>
  <td align="left" colspan="7">

Model01 is a Maclaurin spheroid for which <math>(a, b, c) = (1, 1, \sqrt{3}/2)</math>.  This specific model has been chosen because &hellip;
<ul>
<li>It is axisymmetric &#8212; hence, it can, in principle, be modeled with a two-dimensional (rather than 3D) code;</li>
<li>It is rotating rather slowly, so it is only mildly distorted from a sphere &#8212; in the EFE Diagram it lies between a spherically symmetric model (upper-right corner) and the point at which the Jacobi/Dedekind sequence bifurcates from the Maclaurin-spheroid sequence (see our discussion of Model03, below);</li>
<li>Its eccentricity, <math>e</math>, is precisely <math>1/2</math> and is among the models for which the numerical value of various physical parameters can be found in Chapter 5, Table I (p. 78) of EFE &#8212; this provides a cross-check for our model;</li>
<li>Because <math>\sin^{-1}e = \pi/6</math> precisely, the expression for <math>\Omega^2</math> can be written in algebraic form.</li>
</ul>
  </td>
</tr>
</table>

For this chosen value of <math>c/a</math>, we appreciate that the eccentricity is, 
<table border="0" align="center" cellpadding="8">
<tr>
 <td align="right"><math>e</math></td>
  <td align="center"><math>\equiv</math></td>
  <td align="left">
<math>
\biggl[1 - \frac{c^2}{a^2}\biggr]^{1 / 2}
=
\biggl[1 - \frac{3}{4}\biggr]^{1 / 2}
=
\frac{1}{2} \, ,
</math>
  </td>
</tr>
</table>
in which case,
<table border="0" align="center" cellpadding="8">
<tr>
 <td align="right"><math>\frac{\sin^{-1}e}{e}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{(\pi/6)}{(1/2)} = \frac{\pi}{3}\, .
</math>
  </td>
</tr>
</table>
Hence, from our [[Apps/MaclaurinSpheroidSequence#MaclaurinFrequency|accompanying discussion of the Maclaurin spheroid sequence]],

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

<tr>
  <td align="right">
<math>
\frac{\Omega^2}{\pi G \rho }
</math>
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\frac{2}{e^2}(3-2e^2)(1-e^2)^{1 / 2} \cdot \frac{\sin^{-1}e}{e} 
- \frac{6(1-e^2)}{e^2}
</math>
  </td>
</tr>

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

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
2^3\biggl[ \frac{5}{2} \biggr]\frac{\pi}{2\cdot 3^{1 / 2}} 
- 18
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\frac{10\pi}{ 3^{1 / 2}} - 18 \, .
</math>
  </td>
</tr>
</table>

The [[Apps/MaclaurinSpheroidSequence#Corresponding_Total_Angular_Momentum|corresponding total angular momentum]] is,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\frac{L}{(GM^3\bar{a})^{1 / 2}}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{6^{1 / 2}}{5} \biggl[  A_1 - A_3 (1-e^2) \biggr]^{1 / 2}(1 - e^2)^{-1 / 3} \, ,</math>
  </td>
</tr>
</table>
where,

<table align="center" border=0 cellpadding="3">
<tr>
  <td align="right">
<math>
A_1
</math>
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\frac{1}{e^2} \biggl[\frac{\sin^{-1}e}{e} - (1-e^2)^{1/2}  \biggr](1-e^2)^{1/2} 
= 0.627598728\, ,
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
A_3
</math>
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\frac{2}{e^2} \biggl[(1-e^2)^{-1/2} -\frac{\sin^{-1}e}{e} \biggr](1-e^2)^{1/2} 
= 0.744802542\, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\Rightarrow ~~~ A_1 - A_3 (1-e^2) 
</math>
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\frac{1}{e^2} \biggl[\frac{\sin^{-1}e}{e} - (1-e^2)^{1/2}  \biggr](1-e^2)^{1/2}
-
\frac{2}{e^2} \biggl[(1-e^2)^{-1/2} -\frac{\sin^{-1}e}{e} \biggr](1-e^2)^{3/2} 
</math>
  </td>
</tr>

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

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

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

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\frac{1}{\sqrt{3}} \biggl[ 5\pi - 3^{5/2}  \biggr]  = 0.068996821 \, .
</math>
  </td>
</tr>
</table>
Hence,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\frac{L}{(GM^3\bar{a})^{1 / 2}}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{6^{1 / 2}}{5} \biggl\{  3^{-1 / 2} \biggl[ 5\pi - 3^{5/2}  \biggr]\biggr\}^{1 / 2}\biggl(\frac{4}{3} \biggr)^{1 / 3}</math>
  </td>
</tr>

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

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>  
\frac{1}{5}\biggl( \frac{2^{14}}{3}\biggr)^{1 / 12} \biggl[ 5\pi - 3^{5/2}  \biggr]^{1 / 2} 
= 0.141633637\, .
</math>
  </td>
</tr>
</table>

Finally, the [[Apps/MaclaurinSpheroidSequence#Alternate_Sequence_Diagrams|ratio of rotational to gravitational potential energy]] is,
<table border="0" cellpadding="5" align="center">

<tr>
  <td align="right">
<math>\tau \equiv \frac{T_\mathrm{rot}}{|W_\mathrm{grav}|}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>~ \frac{3}{2e^2}\biggl[ 1 - \frac{e(1-e^2)^{1 / 2}}{\sin^{-1} e}\biggr] - 1</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{3}{\pi}\biggl[ 2\pi - 3^{3/2}\biggr] - 1
= 0.038039942 \, .</math>
  </td>
</tr>
</table>

==Model04==

<table border="1" align="center" cellpadding="8">
<tr>
  <td align="center">Model</td>
  <td align="center"><math>\frac{b}{a}</math></td>
  <td align="center"><math>\frac{c}{a}</math></td>
  <td align="center"><math>e</math></td>
  <td align="center"><math>\frac{\Omega^2}{\pi G\rho}</math></td>
  <td align="center"><math>\frac{T_\mathrm{rot}}{|W_\mathrm{grav}|}</math></td>
</tr>
<tr>
  <td align="center">04</td>
  <td align="center"><math>1</math></td>
  <td align="center"><math>0.258819045</math></td>
  <td align="center"><math>\sin(5\pi/12) =</math><br /><math>0.965925827</math></td>
  <td align="center"><math>0.421802734</math></td>
  <td align="center"><math>0.300647998</math></td>
</tr>
</table>

Here we examine a Maclaurin spheroid for which <math>(a, b) = (1, 1)</math>, and,
<table border="0" align="center" cellpadding="8">
<tr>
 <td align="right"><math>c</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\biggl[ 1 - \frac{1}{8} \biggl(\sqrt{3}+1\biggr)^2 \biggr]^{1 / 2} \, .
</math>
  </td>
</tr>
</table>

For this chosen value of <math>c/a</math>, we appreciate that the eccentricity is, 
<table border="0" align="center" cellpadding="8">
<tr>
 <td align="right"><math>e</math></td>
  <td align="center"><math>\equiv</math></td>
  <td align="left">
<math>
\biggl[1 - \frac{c^2}{a^2}\biggr]^{1 / 2}
=
\biggl\{
\frac{1}{8} \biggl(\sqrt{3}+1\biggr)^2 
\biggr\}^{1 / 2}
=
\frac{\sqrt{3}+1}{2\sqrt{2}} \, ;
</math>
  </td>
</tr>
</table>
and from our accompanying examples of [[Appendix/Ramblings/TrigFunctions#Integer-Degree_Angles|analytic functions for selected trigonometric functions]] we appreciate that, for this value of the eccentricity,
<table border="0" align="center" cellpadding="8">

<tr>
 <td align="right"><math>\sin^{-1}e</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{5\pi}{12} 
</math>
  </td>
</tr>

<tr>
 <td align="right"><math>\Rightarrow ~~~ \frac{\sin^{-1}e}{e}</math></td>
  <td align="center"><math>=</math></td>
  <td align="left">
<math>
\frac{5\pi}{6} \biggl[ \frac{\sqrt{2}}{\sqrt{3}+1}  \biggr]
\, .
</math>
  </td>
</tr>
</table>
Hence, from our [[Apps/MaclaurinSpheroidSequence#MaclaurinFrequency|accompanying discussion of the Maclaurin spheroid sequence]],

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

<tr>
  <td align="right">
<math>
\frac{\Omega^2}{\pi G \rho }
</math>
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\frac{2}{e^2}(3-2e^2)(1-e^2)^{1 / 2} \cdot \frac{\sin^{-1}e}{e} 
- \frac{6(1-e^2)}{e^2}
</math>
  </td>
</tr>

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\frac{2}{e^2}\biggl\{ 
(3-2e^2)(1-e^2)^{1 / 2} \cdot \frac{\sin^{-1}e}{e} 
+3e^2 - 3
\biggr\}
</math>
  </td>
</tr>

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

<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>
=
</math>
  </td>
  <td align="left">
<math>
\biggl[ \frac{2^4}{\ell^2} \biggr]\biggl\{ 
\biggl[3- \frac{\ell^2}{4}\biggr]
\biggl[ 1 - \frac{\ell^2}{8}  \biggr]^{1 / 2} \cdot \frac{5\pi}{6} \biggl[ \frac{\sqrt{2}}{\ell}  \biggr] 
+ \frac{3\ell^2}{8} - 3
\biggr\} \, ,
</math>
  </td>
</tr>
</table>
where,

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

<tr>
  <td align="right">
<math>
\ell
</math>
  </td>
  <td align="center">
<math>
\equiv
</math>
  </td>
  <td align="left">
<math>
\sqrt{3} + 1 \, .
</math>
  </td>
</tr>
</table>

==Potentially Interesting Models==

===Handy Maclaurin Spheroid Formulae===
Drawing principally from an [[Apps/MaclaurinSpheroidSequence#Maclaurin_Spheroid_Sequence|accompanying discussion of equilibrium models along the Maclaurin spheroid sequence]], once a configuration's eccentricity, <math>e</math>, is specified, we know the following:

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

<tr>
  <td align="right">
<math>
\frac{c}{a }
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
(1 - e^2)^{1 / 2}
\, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>
\frac{\omega_0^2}{4\pi G \rho }
</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\frac{1}{2e^2} \biggl[
(3-2e^2)(1-e^2)^{1 / 2} \cdot \frac{\sin^{-1}e}{e} 
- 3(1-e^2)
\biggr]
\, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>[L_*^2]_\mathrm{MPT77} \equiv \frac{L^2}{(GM^3\bar{a})}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\frac{6}{25 e^2} \biggl[ 
(3 - 2e^2) \frac{\sin^{-1}e}{e} - 3(1-e^2)^{1/2}  
\biggr] (1-e^2)^{- 1 / 6} \, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>j^2 \equiv \frac{L^2 \rho^{1 / 3}}{4\pi G M^{10/3} } </math>
  </td>
  <td align="center">
<math>= </math>
  </td>
  <td align="left">
<math>
 \biggl( \frac{3}{2^8 \pi^4}  \biggr)^{1/3}  \frac{L^2}{(GM^3\bar{a})} \, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>x_\mathrm{Wong74} </math>
  </td>
  <td align="center">
<math>= </math>
  </td>
  <td align="left">
<math>
5^2\biggl( \frac{2^2\pi^4}{3^4}\biggr)^{1 / 3} j^2 \, ;
</math>
  </td>
</tr>

<tr>
  <td align="right">
<math>E_\mathrm{Wong74} </math>
  </td>
  <td align="center">
<math>\equiv </math>
  </td>
  <td align="left">
<math>
\frac{3}{5} \biggl(\frac{4\pi}{3}\biggr)^{1 / 2} \frac{GM^{5 / 3}}{\rho^{1 / 3}} \, .
</math>
  </td>
</tr>
</table>

===First Ideas===

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

<tr>
  <td align="center" colspan="1">
<b>First Ideas</b><br />

[[File:FirstIdeasAnnotated.png|550px|First Ideas]]
  </td>
</tr>
</table>


<table border="1" align="center" cellpadding="8" width="80%">
<tr>
  <td align="center" colspan="7">
Selected Models Along the (Axisymmetric) Maclaurin Spheroid Sequence<br /><font size="01">Listed in Order of Increasing <math>e</math> and Increasing <math>j^2</math></font>
  </td>
</tr>
<tr>
  <td align="center"><math>\frac{b}{a}</math></td>
  <td align="center"><math>\frac{c}{a}</math></td>
  <td align="center"><math>e</math></td>
  <td align="center"><math>\frac{\Omega^2}{4\pi G\rho}</math></td>
  <td align="center"><math>j^2 = \biggl(\frac{3}{2^8\pi^4}\biggr)^{1 / 3}\frac{L^2}{(GM^3\bar{a})}</math></td>
  <td align="center"><math>\frac{T_\mathrm{rot}}{|W_\mathrm{grav}|}</math></td>
  <td align="center">NOTES</td>
</tr>
<tr>
  <td align="center"><math>1</math></td>
  <td align="center"><math>\frac{\sqrt{3}}{2} = 0.866025404</math></td>
  <td align="center"><math>\sin(\pi/6) = 0.5</math></td>
  <td align="center"><math>
\frac{1}{4}\biggl[\frac{10\pi}{ 3^{1 / 2}} - 18\biggr]</math><br />&nbsp;<br /><math>= 0.034498410</math></td>
  <td align="center"><math>\frac{1}{5^2}\biggl( \frac{3}{2^2\pi^8}\biggr)^{1 / 6} \biggl[ 5\pi - 3^{5/2}  \biggr]</math><br />&nbsp;<br /><math> = 9.902847\times 10^{-4}</math></td>
  <td align="center"><math>\biggl[\frac{3}{\pi}\biggl( 2\pi - 3^{3/2}\biggr) - 1 \biggr]</math><br />&nbsp;<br /><math>=0.038039942</math></td>
  <td align="center">see<br />[[#Model01|Model01]]<br />above</td>
</tr>

<tr>
  <td align="center"><math>1</math></td>
  <td align="center"><math>\cos(5\pi/12) =</math><br /><math>0.258819045</math></td>
  <td align="center"><math>\sin(5\pi/12) =</math><br /><math>0.965925827</math></td>
  <td align="center"><math>0.105450684</math></td>
  <td align="center"><math>0.0493659 \times 0.553964^2</math><br /><math>= 0.015149</math></td>
  <td align="center"><math>0.300647998</math></td>
  <td align="center">see<br />[[#Model04|Model04]]<br />above</td>
</tr>

</table>

===Drawn From CKST95d===

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

<tr>
  <td align="center" colspan="1">
<b>Drawn From CKST95d</b><br />

[[File:CKST95dAnnotated.png|550px|Drawn from CKST95d]]
  </td>
</tr>
</table>



<table border="1" align="center" cellpadding="8" width="80%">
<tr>
  <td align="center" colspan="14">
Selected Models Along the (Axisymmetric) Maclaurin Spheroid Sequence<br /><font size="01">Listed in Order of Increasing <math>e</math> and Increasing <math>j^2</math></font>

----

Similar to Table 1 (p. 502) in &hellip;<br />
{{ HE84figure }}
<br />and, as well,<br />
Similar to Table 1 (p. 513) in &hellip;<br />
{{ CKST95dfigure }}
  </td>
</tr>
<tr>
  <td align="center" rowspan="2">'''Bifurcating<br />Sequence'''</td>
  <td align="center" rowspan="1" colspan="2">[[Apps/EriguchiHachisu/Models#HE84|HE84]]<br />notation</td>
  <td align="center" rowspan="2"><math>\frac{b}{a}</math></td>
  <td align="center" rowspan="2"><math>\frac{c}{a}</math></td>
  <td align="center" rowspan="2"><math>e</math></td>
  <td align="center" rowspan="2"><math>\frac{\Omega^2}{4\pi G\rho}</math></td>
  <td align="center" rowspan="2"><math>j^2 = \biggl(\frac{3}{2^8\pi^4}\biggr)^{1 / 3}\frac{L^2}{(GM^3\bar{a})}</math></td>
<td align="center" rowspan="1" colspan="3">[<b>[[Appendix/References#T78|<font color="red">T78</font>]]</b>]'s Energy Normalization</td>
<td align="center" rowspan="1" colspan="1">{{ CKST95dhereafter }}'s<br />Normalization</td>
  <td align="center" rowspan="2"><math>\frac{T_\mathrm{rot}}{|W_\mathrm{grav}|}</math></td>
  <td align="center" rowspan="2">NOTES</td>
</tr>

<tr>
  <td align="center" rowspan="1" colspan="1"><math>(n, m)</math></td>
  <td align="center" rowspan="1" colspan="1"><math>\xi_0</math></td>
  <td align="center" rowspan="1" colspan="1"><math>T_\mathrm{rot}/E_\mathrm{T87}</math></td>
  <td align="center" rowspan="1" colspan="1"><math>W_\mathrm{grav}/E_\mathrm{T87}</math></td>
  <td align="center" rowspan="1" colspan="1"><math>E_\mathrm{tot}/E_\mathrm{T87}</math></td>
  <td align="center" rowspan="1" colspan="1"><math>E_\mathrm{tot}/E_\mathrm{CKST95d}</math></td>
</tr>

<tr>
  <td align="right">Jacobi</td>
  <td align="center"><math>(2, 2)</math></td>
  <td align="center"><math>0.717049</math></td>
  <td align="center"><math>1</math></td>
  <td align="center" bgcolor="yellow"><math>0.582724</math></td>
  <td align="center"><math>0.812670</math></td>
  <td align="center"><math>0.0935574</math></td>
  <td align="center"><math>0.00455473</math></td>
  <td align="center"><math>0.0804614</math></td>
  <td align="center"><math>-0.585054</math></td>
  <td align="center"><math>-0.504592</math></td>
  <td align="center"><math>-2.94820 \times 10^{-4}</math></td>
  <td align="center"><math>0.137528</math></td>
  <td align="center">(a)</td>
</tr>

<tr>
  <td align="right">Triangle</td>
  <td align="center"><math>(3, 3)</math></td>
  <td align="center"><math>0.48642</math></td>
  <td align="center"><math>1</math></td>
  <td align="center"><math>0.43741</math></td>
  <td align="center" bgcolor="yellow"><math>0.89926</math></td>
  <td align="center"><math>0.11004</math></td>
  <td align="center"><math>0.0078526</math></td>
  <td align="center"><math>0.11458</math></td>
  <td align="center"><math>-0.56628</math></td>
  <td align="center"><math>-0.45171</math></td>
  <td align="center"><math>-4.5502 \times 10^{-4}</math></td>
  <td align="center"><math>0.20233</math></td>
  <td align="center">(b)</td>
</tr>

<tr>
  <td align="right">Square</td>
  <td align="center"><math>(4, 4)</math></td>
  <td align="center"><math>0.38652</math></td>
  <td align="center"><math>1</math></td>
  <td align="center"><math>0.36052</math></td>
  <td align="center" bgcolor="yellow"><math>0.93275</math></td>
  <td align="center"><math>0.11231</math></td>
  <td align="center"><math>0.010371</math></td>
  <td align="center"><math>0.13303</math></td>
  <td align="center"><math>-0.55029</math></td>
  <td align="center"><math>-0.41726</math></td>
  <td align="center"><math>-5.5513 \times 10^{-4}</math></td>
  <td align="center"><math>0.24174</math></td>
  <td align="center">(c)</td>
</tr>

<tr>
  <td align="right">Dyn. Unstable<br /><math>2^\mathrm{nd}~\mathrm{Harmonic}</math></td>
  <td align="center">n/a</td>
  <td align="center"><math>0.31833</math></td>
  <td align="center"><math>1</math></td>
  <td align="center" bgcolor="yellow"><math>0.30333</math></td>
  <td align="center"><math>0.95289</math></td>
  <td align="center"><math>0.11006</math></td>
  <td align="center"><math>0.012796</math></td>
  <td align="center"><math>0.14627</math></td>
  <td align="center"><math>-0.53418</math></td>
  <td align="center"><math>-0.38791</math></td>
  <td align="center"><math>-6.3671 \times 10^{-4}</math></td>
  <td align="center"><math>0.27382</math></td>
  <td align="center">(d)</td>
</tr>

<tr>
  <td align="right">Dyn. Unstable<br /><math>3^\mathrm{rd}~\mathrm{Harmonic}</math></td>
  <td align="center">n/a</td>
  <td align="center"><math>0.26364</math></td>
  <td align="center"><math>1</math></td>
  <td align="center"><math>0.25493</math></td>
  <td align="center" bgcolor="yellow"><math>0.96696</math></td>
  <td align="center"><math>0.10491</math></td>
  <td align="center"><math>0.015380</math></td>
  <td align="center"><math>0.15657</math></td>
  <td align="center"><math>-0.51660</math></td>
  <td align="center"><math>-0.36003</math></td>
  <td align="center"><math>-7.1029 \times 10^{-4}</math></td>
  <td align="center"><math>0.30308</math></td>
  <td align="center">(e)</td>
</tr>

<tr>
  <td align="right">Ammonite</td>
  <td align="center"><math>(3, 1)</math></td>
  <td align="center"><math>0.25337</math></td>
  <td align="center"><math>1</math></td>
  <td align="center"><math>0.24560</math></td>
  <td align="center" bgcolor="yellow"><math>0.96937</math></td>
  <td align="center"><math>0.10353</math></td>
  <td align="center"><math>0.015950</math></td>
  <td align="center"><math>0.15840</math></td>
  <td align="center"><math>-0.51269</math></td>
  <td align="center"><math>-0.35429</math></td>
  <td align="center"><math>-7.2491 \times 10^{-4}</math></td>
  <td align="center"><math>0.30895</math></td>
  <td align="center">(f)</td>
</tr>

<tr>
  <td align="center">---</td>
  <td align="center"><math>(4, 2)</math></td>
  <td align="center"><math>0.19793</math></td>
  <td align="center"><math>1</math></td>
  <td align="center"><math>0.19416</math></td>
  <td align="center" bgcolor="yellow"><math>0.98097</math></td>
  <td align="center"><math>0.093348</math></td>
  <td align="center"><math>0.019674</math></td>
  <td align="center"><math>0.16704</math></td>
  <td align="center"><math>-0.48713</math></td>
  <td align="center"><math>-0.32009</math></td>
  <td align="center"><math>-8.0783 \times 10^{-4}</math></td>
  <td align="center"><math>0.34290</math></td>
  <td align="center">(g)</td>
</tr>

<tr>
  <td align="right">One-Ring</td>
  <td align="center"><math>(4, 0)</math></td>
  <td align="center"><math>0.17332</math></td>
  <td align="center"><math>1</math></td>
  <td align="center"><math>0.17078</math></td>
  <td align="center" bgcolor="yellow"><math>0.98531</math></td>
  <td align="center"><math>0.087121</math></td>
  <td align="center"><math>0.021788</math></td>
  <td align="center"><math>0.16982</math></td>
  <td align="center"><math>-0.47271</math></td>
  <td align="center"><math>-0.30289</math></td>
  <td align="center"><math>-8.4656 \times 10^{-4}</math></td>
  <td align="center"><math>0.35925</math></td>
  <td align="center">(h)</td>
</tr>

<tr>
  <td align="left" colspan="14">
NOTES:
<ol type="a">
<li>Primary data source (bgcolor = yellow): &nbsp;Table IV (Chapter 6, &sect;39, p. 103) of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  See also &hellip; the first line in <sup>&dagger;</sup>Table B1 (p. 446) of {{ Bardeen71 }}; and Appendices D.3 & D.4 (pp. 485-486) of [<b>[[Appendix/References#T78|<font color="red">T78</font>]]</b>].
</li>
<li>Primary data source (bgcolor = yellow): &nbsp;Eq. (63) (Chapter 6, &sect;41, p. 111) of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  
</li>
<li>Primary data source (bgcolor = yellow): &nbsp;bottom of p. 128 of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  
</li>
<li>Primary data source (bgcolor = yellow): &nbsp;Table VI (Chapter 7, &sect;48, p. 142) of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  From p. 141 of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>], we find <font color="darkgreen">"&hellip; the Maclaurin spheroid on the verge of dynamical instability is the first member of the</font> [Riemann S-type ellipsoid] <font color="darkgreen"> self-adjoint sequence <math>x = + 1</math>."</font>  In <sup>&dagger;</sup>Table B1 (p. 446) of {{ Bardeen71 }}, this is referred to as the "First nonaxisymmetric dynamical instability."</li>
<li>Primary data source (bgcolor = yellow): &nbsp;Eq. (105) (Chapter 6, &sect;43, p. 119) of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  
</li>
<li>Primary data source (bgcolor = yellow): &nbsp;Eq. (68) (Chapter 6, &sect;41, p. 112) of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  
</li>
<li>Primary data source (bgcolor = yellow): &nbsp;bottom of p. 128 of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  
</li>
<li>Primary data source (bgcolor = yellow): &nbsp;bottom of p. 128 of [<b>[[Appendix/References#EFE|<font color="red">EFE</font>]]</b>].  In <sup>&dagger;</sup>Table B1 (p. 446) of {{ Bardeen71 }}, this is referred to as the "First axisymmetric secular instability." 
</li>

</ol>

<sup>&dagger;</sup>Table B1 of {{ Bardeen71 }} lists numerically determined values of <math>c/a, \xi_0, e, T_\mathrm{rot}/W_\mathrm{grav}</math>, and a dimensionless measure of the squared angular momentum that is a factor of <math>4\pi</math> larger than our listed <math>j^2</math>.
  </td>
</tr>
</table>

==Axisymmetric Equilibrium Sequences that Display a Topological Change==

In the context of discussions of self-gravitating configurations that are rotating and have a uniform density, there are three well-known ''axisymmetric'' equilibrium sequences:
<ol>
  <li>The Maclaurin spheroid sequence &#8212; first constructed in 1742 by Maclaurin.</li>
  <li>A so-called one-ring (Dyson-Wong toroid) sequence &#8212; the ring-like segment was first mapped out in 1893 by Dyson</li>
  <li>A so-called Maclaurin toroid sequence &#8212; the toroidal segment was first constructed by {{ MPT77 }}.</li>
</ol>
Along the first of these sequences, every equilibrium configuration has a surface that is precisely a spheroid. The other two bifurcate from the Maclaurin spheroid sequence; contain a segment of concave, hamburger-shaped (spheroidal-like) configurations; and, upon further extension, blend into the (separately identified) sequence of toroidal-shaped configurations. In other words, moving along the respective model sequences, we encounter a spheroidal-like to ring-like topological change.  The overriding question that requires nonlinear dynamical modeling to answer is: "Is there a model (or a range of models) along the Maclaurin spheroid sequence that is unstable toward evolution away from that sequence and toward one of the ring-like sequences?"

'''TWO KEY CONCEPTS:'''
<ol type="A">
  <li>The model parameterization that distinguishes the two ring-like sequences from one another is the specified ''distribution of angular momentum.''  Dyson-Wong tori are uniformly rotating, whereas, each Maclaurin toroid has an angular momentum distribution that is the same as is present in all Maclaurin spheroids.  It is this single structural feature that drives the bifurcation points for these sequences to two quite different locations along the Maclaurin spheroid sequence.</li>
  <li>The bifurcation points of both ring-like sequences arise at positions on the Maclaurin spheroid sequence where the models are expected to be violently unstable toward the development of nonaxisymmetric distortions.  Hence, in order to study how unstable spheroids undergo a transition to a ring-like configuration (the topological change), a numerical code must be developed with the capability to suppress all nonaxisymmetric distortions.  This can be accomplished, for example, by adopting a cylindrical coordinate representation of the fluid equations then performing two-dimensional rather than three-dimensional simulations.</li>
</ol>


===One-Ring (Dyson-Wong) Sequence===

It is important to remember &#8212; as emphasized above &#8212; that ''all'' equilibrium models along the one-ring (Dyson-Wong) sequence are uniformly rotating.

====Background Storyline====

<ol>
  <li>Over 125 years ago via a pair of detailed publications &#8212; {{ Dyson1893full }} and {{ Dyson1893Part2full }} &#8212; Dyson demonstrated that a sequence of rapidly rotating, self-gravitating equilibrium models could be constructed that had a uniform density, were uniformly rotating, and had a toroidal (ring) shape. [Configurations that, in every respect except their ''shape'', were like Maclaurin  spheroids.]  See [[Apps/DysonPotential#Dyson_(1893)|our review and discussion]] of this work.</li>
  <li>{{ Wong74full }} tackled this same problem, improving on, and extending Dyson's work.  As a consequence, this sequence of models is often referred to as "Dyson-Wong tori."  See [[Apps/DysonWongTori#Self-Gravitating,_Incompressible_(Dyson-Wong)_Tori|our detailed review and discussion]].</li>
  <li>Motivated by the work of {{ FESB-K80full }}, {{ ES81full }} demonstrated that the Dyson-Wong toroidal sequence can be "smoothly connected" to the Maclaurin spheroid sequence via an intermediate branch of models having a <font color="darkgreen">&hellip; concave hamburger-like shape of equilibrium &hellip;"</font>
<ul>
  <li>Table I of {{ ES81 }} provides quantitative data describing the properties of eighteen models that lie along this combined "one-ring" sequence, such as:  <math>\Omega^2/(4\pi G \rho), j^2, T_\mathrm{rot}/|W_\mathrm{grav}|,</math> and <math>(T_\mathrm{rot} + W_\mathrm{grav})/E_0</math>.  We have copied the values of two of these parameters from their Table I into the first two (pink) columns of our table, immediately below.</li>
  <li>This data has been used to generate the pink-colored one-ring sequence shown in our plot, below; see especially the plot inset.  Over the years, the same set of data has been used to display the behavior of the one-ring sequence in numerous publications; see, for example, the reproduction of Figure 1 (p. 488) from {{ CKST95b }} that we have presented, below.</li>
  <li>Figure 2 of {{ ES81 }} displays meridional-plane cross-sections through five of their eighteen models in an effort to illustrate how the surface geometry smoothly changes along the complete sequence: from spheroid, to "hamburger" shape, to torus.</li>
</ul>
  </li>
  <li>{{ ES81 }} claim that the one-ring sequence bifurcates from the Maclaurin sequence precisely at the point where the spheroid has an eccentricity, <math>e = e_\mathrm{cr} = 0.98523</math> &#8212; in which case, also, <math>\Omega^2/(4\pi G \rho) = 0.08726</math> and <math>j^2 = 0.02174</math>.  In support of this conjecture, they point out that, Chandrasekhar (1967; publication XXX) and {{ Bardeen71 }} have shown that this is <font color="darkgreen">&hellip; a neutral point on the Maclaurin sequence against the perturbation of <math>P_4(\eta)</math> displacement at the surface where <math>\eta</math> is one of the spheroidal coordinates."</font>  This is also the "neutral point" on the Maclaurin sequence labeled "F" in Table I of {{ HE82 }}; and the "bifurcation point" along the Maclaurin sequence that is labeled by the quantum numbers, <math>(n, m) = (4, 0)</math> in Table 1 of {{ HE84 }}.</li>
</ol>

====Model-Sequence Details====

<table border="1" align="center" cellpadding="5">
<tr><td align="center" colspan="13"><b>One-Ring Sequence</b> (see figure inset)<br />as quantitatively described in three separate studies<sup>&dagger;</sup></td></tr>

<tr>
  <td align="center" colspan="3" rowspan="2">{{ ES81 }}<br /><font size="-1">Data extracted from their Table I</font></td>
<td align="center" width="2%" bgcolor="lightgrey" colspan="1" rowspan="21">&nbsp;</td>
  <td align="center" colspan="3">{{ HES82 }}</td></td>
<td align="center" width="2%" bgcolor="lightgrey" colspan="1" rowspan="21">&nbsp;</td>
  <td align="center" colspan="5">{{ Hachisu86a }}</td></td>
</tr>
<tr>
  <td align="center" colspan="2"><font size="-1">Data extracted from their Table I</font></td>
  <td align="center" colspan="1"><font size="-1">Implication</font></td>
  <td align="center" colspan="3"><font size="-1">Data extracted from his Table Ia</font></td>
  <td align="center" colspan="2"><font size="-1">Implication (assuming G = &rho; = 1)</font></td>
</tr>

<tr>
  <td align="center"><math>\frac{\Omega^2}{4\pi G \rho}</math></td>
  <td align="center"><math>j^2</math></td>
  <td align="center"><math>\frac{T_\mathrm{rot}}{|W_\mathrm{grav}|}</math></td>
  <td align="center"><math>\frac{\Omega^2}{4\pi G \rho}</math></td>
  <td align="center"><math>j</math></td>
  <td align="center"><math>j^2</math></td>
  <td align="center"><math>\Omega^2</math></td>
  <td align="center"><math>M</math></td>
  <td align="center"><math>J</math></td>
  <td align="center"><math>\frac{\Omega^2}{4\pi G \rho}</math></td>
  <td align="center"><math>j^2 = \frac{J^2 \rho^{1 / 3}}{4\pi G M^{10/3}}</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08506</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02243</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3648</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.09635</math></td>
  <td align="center"><math>0.13725</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.018837</math></td>

  <td align="center"><math>0.000</math></td>
  <td align="center"><math>4.17</math></td>
  <td align="center"><math>0.0</math></td>
  <td align="center"><math>0.00</math></td>
  <td align="center"><math>0.00</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08324</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02270</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3654</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.09104</math></td>
  <td align="center"><math>0.14360</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.020621</math></td>

  <td align="center"><math>1.31</math></td>
  <td align="center"><math>2.09</math></td>
  <td align="center"><math>0.960</math></td>
  <td align="center"><math>0.1042</math></td>
  <td align="center"><math>0.00628</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08206</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02269</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3637</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.08763</math></td>
  <td align="center"><math>0.14702</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.021619</math></td>

  <td align="center"><math>1.41</math></td>
  <td align="center"><math>1.40</math></td>
  <td align="center"><math>0.666</math></td>
  <td align="center"><math>0.1122</math></td>
  <td align="center"><math>0.0115</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08139</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02251</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3608</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.08646</math></td>
  <td align="center"><math>0.14799</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.021901</math></td>

  <td align="center"><math>1.31</math></td>
  <td align="center"><math>1.05</math></td>
  <td align="center"><math>0.483</math></td>
  <td align="center"><math>0.1042</math></td>
  <td align="center"><math>0.0158</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08113</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02224</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3576</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.08488</math></td>
  <td align="center"><math>0.14896</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.022189</math></td>

  <td align="center"><math>1.09</math></td>
  <td align="center"><math>0.723</math></td>
  <td align="center"><math>0.307</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0867</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0221</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08119</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02198</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3550</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.08379</math></td>
  <td align="center"><math>0.14927</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.022281</math></td>

  <td align="center"><math>1.01</math></td>
  <td align="center"><math>0.811</math></td>
  <td align="center"><math>0.377</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0804</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0227</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08139</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02181</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3534</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.08276</math></td>
  <td align="center"><math>0.14904</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.022213</math></td>

  <td align="center"><math>1.01</math></td>
  <td align="center"><math>0.856</math></td>
  <td align="center"><math>0.404</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0804</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0218</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08150</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02177</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3531</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.08207</math></td>
  <td align="center"><math>0.14822</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.021970</math></td>

  <td align="center"><math>1.01</math></td>
  <td align="center"><math>0.929</math></td>
  <td align="center"><math>0.447</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0804</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0203</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08183</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02163</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3519</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.08172</math></td>
  <td align="center"><math>0.14714</math></td>
  <td align="center" bgcolor="#7DCEA0"><math>0.021650</math></td>

  <td align="center"><math>0.938</math></td>
  <td align="center"><math>0.953</math></td>
  <td align="center"><math>0.461</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0746</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0199</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08236</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02104</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3461</math></td>
  <td align="center" colspan="3" rowspan="9">&nbsp;</td>

  <td align="center"><math>0.826</math></td>
  <td align="center"><math>0.924</math></td>
  <td align="center"><math>0.445</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0657</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0205</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.08174</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02037</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3378</math></td>

  <td align="center"><math>0.692</math></td>
  <td align="center"><math>0.845</math></td>
  <td align="center"><math>0.398</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0551</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0221</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.07944</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.01992</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3279</math></td>

  <td align="center"><math>0.548</math></td>
  <td align="center"><math>0.721</math></td>
  <td align="center"><math>0.326</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0436</math></td>
  <td align="center" bgcolor="#D7BDE2"><math>0.0252</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.07556</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.01980</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3228</math></td>

  <td align="center"><math>0.408</math></td>
  <td align="center"><math>0.571</math></td>
  <td align="center"><math>0241</math></td>
  <td align="center"><math>0.0325</math></td>
  <td align="center"><math>0.0299</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.07041</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02003</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3176</math></td>

  <td align="center"><math>0.278</math></td>
  <td align="center"><math>0.409</math></td>
  <td align="center"><math>00.154398</math></td>
  <td align="center"><math>0.0221</math></td>
  <td align="center"><math>0.0372</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.06473</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02063</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3140</math></td>

  <td align="center"><math>0.169</math></td>
  <td align="center"><math>0.255</math></td>
  <td align="center"><math>0.0815</math></td>
  <td align="center"><math>0.0134</math></td>
  <td align="center"><math>0.0503</math></td>
</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.05775</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02162</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3119</math></td>

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

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.05088</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02304</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3113</math></td>

</tr>

<tr>
  <td align="center" bgcolor="#FDE0DA"><math>0.04399</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.02497</math></td>
  <td align="center" bgcolor="#FDE0DA"><math>0.3122</math></td>

</tr>

<tr>
<td align="left" colspan="13">
<sup>&dagger;</sup>The data drawn from these three separate studies are displayed in the figure inset as follows:
<ul>
  <li>Pink circular markers and accompanying smooth curve: &nbsp; {{ ES81 }}</li>
  <li>Green square markers: &nbsp; {{ HES82 }}</li>
  <li>Light-purple triangular markers: &nbsp; {{ Hachisu86a }}</li>
</ul>
</td>
</tr>
</table>

<span id="PhaseTransition">&nbsp;</span>
<table border="1" align="center" cellpadding="5">
<tr>
  <td align="center" colspan="1" rowspan="2">
[[File:ES81OneRingWithInsetBox.png|800px|One-Ring Sequence]]
  </td>
  <td align="center" colspan="1" rowspan="1">
Figure 1 extracted from &sect;2.2, p. 488 of &hellip;<br />{{ CKST95bfigure }}
  </td>
</tr>
<tr>
  <td align="center" colspan="1">[[File:CKST95bFig1annotated.png|400px|CKST95b Figure 1]]</td>
</tr>
</table>

==Most Interesting Initial Configurations for Axisymmetric Simulations==
In what follows, keep in mind that,

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

<tr>
  <td align="right">
<math>j^2</math>
  </td>
  <td align="center">
<math>\equiv \frac{L^2}{4\pi G M^{10/3} \rho^{-1 / 3}} = </math>
  </td>
  <td align="left">
<math>
\biggl( \frac{3}{2^8 \pi^4}  \biggr)^{1/3}  \underbrace{\frac{L^2}{(GM^3\bar{a})}}_{L_*^2} 
= \frac{1}{3}\biggl(\frac{4\pi}{3}\biggr)^{-4 / 3} L_*^2
\, .
</math>
  </td>
</tr>
</table>

===Bardeen71===

{{ Bardeen71 }} highlights five "critical points" along the Newtonian Maclaurin spheroid sequence in his Table B1 (Appendix B, p. 446).

<table border="1" align="center" cellpadding="5">
<tr>
  <td align="center" colspan="5">
Data extracted from Table B1 (Appendix B, p. 446) of &hellip;<br />{{ Bardeen71figure }}
  </td>
  <td align="center" colspan="1">
Our Determination
  </td>
</tr>
<tr>
  <td align="center">
<math>\frac{T_\mathrm{rot}}{|W_\mathrm{grav}|}</math>
  </td>
  <td align="center">
<math>\frac{J^2\rho^{1 / 3}}{GM_0^{10/3}}</math>
  </td>
  <td align="center">
<math>\xi_s</math>
  </td>
  <td align="center">
<math>\frac{c}{a}</math>
  </td>
  <td align="center">
<math>e</math>
  </td>
  <td align="center">
<math>j^2 = \frac{1}{4\pi}\biggl[\frac{J^2\rho^{1 / 3}}{GM_0^{10/3}}\biggr]</math>
  </td>
</tr>

<tr>
  <td align="center"><math>0.13753^a</math></td>
  <td align="center"><math>0.057236</math></td>
  <td align="left"><math>0.71705</math></td>
  <td align="center"><math>0.58272</math></td>
  <td align="center"><math>0.81267</math></td>
  <td align="center"><math>4.555\times 10^{-3}</math></td>
</tr>

<tr>
  <td align="center"><math>0.23790^b</math></td>
  <td align="center"><math>0.126991</math></td>
  <td align="left"><math>0.39537</math></td>
  <td align="center"><math>0.36767</math></td>
  <td align="center"><math>0.92996</math></td>
  <td align="center"><math>1.0106\times 10^{-2}</math></td>
</tr>

<tr>
  <td align="center"><math>0.27383^c</math></td>
  <td align="center"><math>0.160802</math></td>
  <td align="left"><math>0.31831</math></td>
  <td align="center"><math>0.30332</math></td>
  <td align="center"><math>0.95289</math></td>
  <td align="center"><math>1.280\times 10^{-2}</math></td>
</tr>

<tr>
  <td align="center"><math>0.35890^d </math></td>
  <td align="center"><math>0.273205</math></td>
  <td align="left"><math>0.17383</math></td>
  <td align="center"><math>0.17126</math></td>
  <td align="center"><math>0.98523</math></td>
  <td align="center"><math>2.174\times 10^{-2}</math></td>
</tr>

<tr>
  <td align="center"><math>0.4512^e</math></td>
  <td align="center" colspan="5">n/a</td>
</tr>

<tr>
  <td align="center"><math>0.45742^f </math></td>
  <td align="center"><math>0.577894</math></td>
  <td align="left"><math>0.04657</math></td>
  <td align="center"><math>0.04657</math></td>
  <td align="center"><math>0.99892</math></td>
  <td align="center"><math>4.599\times 10^{-2}</math></td>
</tr>
</table>

'''NOTES:'''
<ol type="a">
  <li>''Secular'' bifurcation point to the (ellipsoidal) Jacobi sequence.</li>
  <li>Maximum value of <math>\omega_0^2/(4\pi G\rho)</math> along the Maclaurin spheroid sequence.</li>
  <li>First nonaxisymmetric (ellipsoidal) ''dynamical'' instability.</li>
  <li>''Secular'' bifurcation point to the Dyson-Wong toroid.</li>
  <li>n/a &#8212; Not identified as a critical point by {{ Bardeen71 }}.</li>
  <li>
''Dynamical'' identified by {{ Bardeen71 }} &hellip;
  <ul>
<li>Also identified at <math>(e, j^2) = (0.998917, 0.04599)</math> as a <math>P_4</math> bifurcation point in Table 2 (p. 292) of {{ EH85 }}.</li>
  </ul>
  </li>
</ol>

===HTE87===

{{ HTE87 }} highlight six "critical points" obtained by "local analysis" in their Table 4 (Appendix A, p. 611).

<table border="1" align="center" cellpadding="5">
<tr>
  <td align="center" colspan="3">
Data extracted from Table 4 (Appendix A, p. 611) of &hellip;<br />{{ HTE87figure }}
  </td>
  <td align="center" colspan="5">
Our Determination
  </td>
</tr>
<tr>
  <td align="center">
<math>\frac{T_\mathrm{rot}}{|W_\mathrm{grav}|}</math>
  </td>
  <td align="center">
<math>\log_{10} F</math>
  </td>
  <td align="center">
Type
  </td>
  <td align="center">
<math>\xi_s = \biggl[\frac{(1-e^2)}{e^2}\biggr]^{1 / 2}</math>
  </td>
  <td align="center">
<math>\frac{c}{a}</math>
  </td>
  <td align="center">
<math>e</math>
  </td>
  <td align="center">
<math>L_*^2 = \biggl(\frac{4\pi}{3}\biggr)^{1 / 3} F^{-1 / 3}</math>
  </td>
  <td align="center">
<math>j^2 = \frac{1}{3}\biggl(\frac{4\pi}{3}\biggr)^{-4 / 3} L_*^2</math>
  </td>
</tr>

<tr>
  <td align="center"><math>0.1375^a</math></td>
  <td align="center"><math>3.727</math></td>
  <td align="left"><math>P_2^2(\eta)\cos(2\phi)</math></td>
  <td align="center"><math>0.71705</math></td>
  <td align="center"><math>0.582724</math></td>
  <td align="center"><math>0.812670</math></td>
  <td align="center"><math>0.09226</math></td>
  <td align="center"><math>4.555\times 10^{-3}</math></td>
</tr>

<tr>
  <td align="center"><math>0.2379^b</math></td>
  <td align="center" colspan="7">n/a</td>
</tr>

<tr>
  <td align="center"><math>0.2738^c</math></td>
  <td align="center"><math>2.381</math></td>
  <td align="left"><math>P_2^2(\eta)\cos(2\phi)</math></td>
  <td align="center"><math>0.31831</math></td>
  <td align="center"><math>0.30333</math></td>
  <td align="center"><math>0.95289</math></td>
  <td align="center"><math>0.25924</math></td>
  <td align="center"><math>1.280\times 10^{-2}</math></td>
</tr>

<tr>
  <td align="center"><math>0.3589^d </math></td>
  <td align="center"><math>1.691</math></td>
  <td align="left"><math>P_4(\eta)</math></td>
  <td align="center"><math>0.17386</math></td>
  <td align="center"><math>0.17129</math></td>
  <td align="center"><math>0.98522</math></td>
  <td align="center"><math>0.44025</math></td>
  <td align="center"><math>2.173\times 10^{-2}</math></td>
</tr>

<tr>
  <td align="center"><math>0.4512^e </math></td>
  <td align="center"><math>0.801</math></td>
  <td align="left"><math>P_6(\eta)</math></td>
  <td align="center"><math>0.053724</math></td>
  <td align="center"><math>0.05365</math></td>
  <td align="center"><math>0.99856</math></td>
  <td align="center"><math>0.87169</math></td>
  <td align="center"><math>4.303\times 10^{-2}</math></td>
</tr>

<tr>
  <td align="center"><math>0.4574^f </math></td>
  <td align="center"><math>0.714</math></td>
  <td align="left"><math>P_4(\eta)</math></td>
  <td align="center"><math>0.046621</math></td>
  <td align="center"><math>0.046571</math></td>
  <td align="center"><math>0.998915</math></td>
  <td align="center"><math>0.93189</math></td>
  <td align="center"><math>4.600\times 10^{-2}</math></td>
</tr>
</table>

'''NOTES:'''
<ol type="a">
  <li>''Secular'' bifurcation point to the (ellipsoidal) Jacobi sequence.</li>
  <li>n/a &#8212; Not identified as a critical point by {{ HTE87 }}.</li>
  <li>First nonaxisymmetric (ellipsoidal) ''dynamical'' instability.</li>
  <li>''Secular'' bifurcation point to the Dyson-Wong toroid.</li>
  <li>''Dynamical''; first ring mode instability and bifurcation point to the Maclaurin toroid &hellip;
  <ul>
<li>Also identified at <math>(e, j^2) = (0.998556, 0.04305)</math> as a <math>P_6</math> bifurcation point in Table 2 (p. 292) of {{ EH85 }}.</li>
  </ul>
  </li>
  <li>
''Dynamical'' identified by {{ Bardeen71 }} &hellip;
  <ul>
<li>Also identified at <math>(e, j^2) = (0.998917, 0.04599)</math> as a <math>P_4</math> bifurcation point in Table 2 (p. 292) of {{ EH85 }}.</li>
  </ul>
  </li>
</ol>

===(Temporary)===
<table border="0" align="center" cellpadding="5"><tr><td align="center">
<table border="0" align="center" cellpadding="3">
<tr>
  <td align="center">'''Figure 7'''</td>
<td align="center" rowspan="3">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</td>
  <td align="center">'''Figure 8'''<br />Duplicate of Fig. 1 from [https://ui.adsabs.harvard.edu/abs/1984PASJ...36..497H/abstract Hachisu &amp; Eriguchi (1984)]</td>
</tr>
<tr>
<td align="center" width="60%">
[[File:JacobiWithDynamicalLimit02.png|center|350px|Jacobi Ellipsoid Sequence & Dynamical Limit]]
</td>
<td align="left" rowspan="2">
<div align="center">[[File:BifurcationPointsHE84.png|center|500px|Bifurcation Points Along Maclaurin Spheroid Sequence]]</div><br />
Bifurcation points on the Maclaurin sequence for the deformation type <math>P_n^m(\eta) \times \cos(m\phi)</math>, plotted in the <math>\omega^2-j^2</math> plane.  The numbers in the parentheses denote the deformation type of <math>(n, m)</math>.  Other computed sequences are also plotted [taken from Eriguchi and Hachisu (1982)].  One-ring sequence starts from the bifurcation point of <math>(4, 0)</math> and two-ring sequence bifurcates from the point of <math>(6, 0)</math>.
</td>
</tr>
<tr>
  <td align="center">
This solid black curve also appears in:
<div align="center">
Fig. 3 (p. 1134) of [https://ui.adsabs.harvard.edu/abs/1983PThPh..69.1131E/abstract Eriguchi &amp; Hachisu (1983)]<br />
Fig. 3 (p. 487) of [https://ui.adsabs.harvard.edu/abs/1986ApJS...61..479H/abstract Hachisu (1986)]<br />
Fig. 4 (p. 4507) of [https://ui.adsabs.harvard.edu/abs/2019MNRAS.487.4504B/abstract Basillais &amp; Hur&eacute; (2019)]
</div>
  </td>
</tr>
</table>
</td></tr></table>

=See Also=
* [[Apps/MaclaurinSpheroids#Maclaurin_Spheroids_.28axisymmetric_structure.29|Properties of Maclaurin Spheroids]]
* [[Apps/MaclaurinSpheroids/GoogleBooks#Excerpts_from_A_Treatise_of_Fluxions|Excerpts from Maclaurin's (1742) ''A Treatise of Fluxions'']]
* [[ThreeDimensionalConfigurations/EFE_Energies#Properties_of_Homogeneous_Ellipsoids_.282.29|Properties of Homogeneous Ellipsoids]]
* Equilibrium Configurations and Sequences Generated by [[Apps/EriguchiHachisu/Models|Eriguchi, Hachisu, and their various colleagues]]
* [[Appendix/Ramblings/ForCohlHoward#Recommended_Riemann_S-Type_Ellipsoids_for_Modeling|Ongoing Collaborative Discussions]] with Howard Cohl.

{{ SGFfooter }}