Appendix/Ramblings/TurningPoints

Turning Points

Turning-Points (Broader Context)

The material presented below has been gathered together and expanded in preparation for developing a primary chapter that discusses instabilities associated with equilibrium-sequence turning points.

Introduction

Spherically symmetric, self-gravitating, equilibrium configurations can be constructed from gases exhibiting a wide variety of degrees of compressibility. When examining how the internal structure of such configurations varies with compressibility, or when examining the relative stability of such structures, it can be instructive to construct models using a polytropic equation of state,

$P=K_{\mathrm{n}}{\rho }^{1+1/n}$

because the degree of compressibility can be adjusted by simply changing the value of the polytropic index, $n$, across the range, $0\le n\le \mathrm{\infty }$. (Alternatively, one can vary the effective adiabatic exponent of the gas, $\mathrm{\Gamma }=1+1/n$.) In particular, $n=0(\mathrm{\Gamma }=\mathrm{\infty })$ represents a hard equation of state and describes an incompressible configuration, while $n=\mathrm{\infty }(\mathrm{\Gamma }=1)$ represents an isothermal and extremely soft equation of state.

As has been detailed in an accompanying discussion, the structural properties of spherical polytropes can be described entirely in terms of a dimensionless radial coordinate, $\xi$, and by the radial dependence of the dimensionless enthalpy function, ${\theta }_n(\xi )$, and its first radial derivative, ${\theta }_n^{\text{'}}(\xi )$. At the center of each configuration $(\xi =0)$, ${\theta }_n=1$ and ${\theta }_n^{\text{'}}=0$. The surface of each isolated polytrope is identified by the radial coordinate, ${\xi }_1$, at which ${\theta }_n$ first drops to zero. As a class, isolated polytropes exhibit three attributes that are especially key in the context of our present discussion:

1. The equilibrium structure is dynamically stable if $n<3$.
2. The equilibrium structure has a finite radius if $n<5$.
3. The equilibrium structure can be described in terms of closed-form analytic expressions for $n=0$, $n=1$, and $n=5$.

Isothermal spheres are discussed in a wide variety of astrophysical contexts because it is not uncommon for physical conditions to conspire to create an extended volume throughout which a configuration exhibits uniform temperature. But, as can be surmised from our list of three key polytrope attributes and recognition that equilibrium isothermal configurations are polytropes with index $n=\mathrm{\infty }$, mathematical models of isolated isothermal spheres are relatively cumbersome to analyze because they extend to infinity, they are dynamically unstable, and they are not describable in terms of analytic functions. In such astrophysical contexts, we have sometimes found it advantageous to employ an $n=5$ polytrope instead of an isothermal sphere. An isolated $n=5$ polytrope can serve as an effective surrogate for an isothermal sphere because it is both infinite in extent and dynamically unstable, but it is less cumbersome to analyze because its structure can be described by closed-form analytic expressions.

Bonnor-Ebert Sphere

In the mid-1950s, 📚 R. Ebert (1955, Zeitschrift für Astrophysik, Vol. 37, pp. 217 - 232) and 📚 W. B. Bonnor (1956, MNRAS, Vol. 116, pp. 351 - 359) independently realized that an equilibrium isothermal gas cloud can be constructed with a finite radius by embedding it in a hot, tenuous external medium. As has been described in an accompanying chapter, the relevant mathematical model is constructed by chopping off the isolated isothermal sphere at some finite radius — call it, ${\xi }_e$ — and imposing an externally applied pressure, $P_e$, that is equal to the pressure of the isothermal gas at the specified edge of the truncated sphere. By varying ${\xi }_e$, a sequence of equilibrium models can be constructed, as illustrated in Figure 1. (See also our more extensive discussion of isothermal equilibrium sequences.)

Reproduction of Figure 2 from R. Ebert (1955) Über die Verdichtung von HI-Gebieten Zeitschrift für Astrophysik, Vol. 37, pp. 217 - 232	Reproduction of Figure 1 from W. B. Bonnor (1956) Boyle's Law and gravitational instability Monthly Notices of the Royal Astronomical Society Vol. 116, pp. 351 - 359	Our Construction
Ebert (1955) Figure 2	Bonnor (1956, MNRAS, 116, 351)	Bonnor (1956, MNRAS, 116, 351)

Reproduction of Figure 2 from R. Ebert (1955) Über die Verdichtung von HI-Gebieten Zeitschrift für Astrophysik, Vol. 37, pp. 217 - 232	Reproduction of Figure 1 from W. B. Bonnor (1956) Boyle's Law and gravitational instability Monthly Notices of the Royal Astronomical Society Vol. 116, pp. 351 - 359
Ebert (1955) Figure 2	Bonnor (1956, MNRAS, 116, 351)
Our Analytic Analysis of n = 5 Polytropic Sequence]
P vs. R n = 5 sequence

Mass-Radius Relations for Pressure-Truncated Polytropes
Reproduction of Fig. 17 from p. 184 of S. W. Stahler (1983) The Equilibria of Rotating Isothermal Clouds. II. Structure and Dynamical Stability The Astrophysical Journal, Vol. 268, pp. 165 - 184	From Our Detailed Analyses
Stahler (1983)	M vs. R sequences for pressure-truncated polytropes	M vs. R sequences for pressure-truncated polytropes
Schematic	Virial-equilibrium configurations, assuming uniform density.	Exact configurations.

Schönberg-Chandrasekhar Mass

Reproduction of Figure 1 from L. R. Henrich & S. Chandrasekhar (1941) Stellar Models with Isothermal Cores The Astrophysical Journal, Vol. 94, pp. 525 - 536	Reproduction of Figure 1 from M. Schönberg & S. Chandrasekhar (1942) On the Evolution of the Main-Sequence Stars The Astrophysical Journal, Vol. 96, pp. 161 - 172	Bipolytropes with $(n_c,n_e)=(5,1)$

Reproduction of Figure 1 from L. R. Henrich & S. Chandrasekhar (1941) Stellar Models with Isothermal Cores The Astrophysical Journal, Vol. 94, pp. 525 - 536	Reproduction of Figure 1 from M. Schönberg & S. Chandrasekhar (1942) On the Evolution of the Main-Sequence Stars The Astrophysical Journal, Vol. 96, pp. 161 - 172	Bipolytropes with $(n_c,n_e)=(5,1)$

In the early 1940s, Chandrasekhar and his colleagues — see 📚 L. R. Henrich & S. Chandrasekhar (1941, ApJ, Vol. 94, pp. 525 - 536) and 📚 M. Schönberg & S. Chandrasekhar (1942, ApJ, Vol. 96, pp. 161 - 172)) — discovered that a star with an isothermal core will become unstable if the fractional mass of the core is above some limiting value. They discovered this by constructing models that are now commonly referred to as composite polytropes or bipolytropes, that is, models in which the star's core is described by a polytropic equation of state having one index — say, $n_c$ — and the star's envelope is described by a polytropic equation of state of a different index — say, $n_e$. In an accompanying discussion we explain in detail how the two structural components with different polytropic indexes are pieced together mathematically to build equilibrium bipolytropes. For a given choice of the two indexes, $n_c$ and $n_e$, a sequence of models can be generated by varying the radial location at which the interface between the core and envelope occurs. As the interface location is varied, the relative amount of mass enclosed inside the core, $\nu \equiv M_{\mathrm{c}\mathrm{o}\mathrm{r}\mathrm{e}}/M_{\mathrm{t}\mathrm{o}\mathrm{t}}$, quite naturally varies as well.

📚 Henrich & Chandrasekhar (1941) built structures of uniform composition having an isothermal core ($n_c=\mathrm{\infty }$) and an $n_e=3/2$ polytropic envelope and found that equilibrium models exist only for values of $\nu \le {\nu }_{\mathrm{m}\mathrm{a}\mathrm{x}}\approx 0.35$. 📚 Schönberg & Chandrasekhar (1942) extended this analysis to include structures in which the mean molecular weight of the gas changes discontinuously across the interface. Specifically, they used the same values of $n_c$ and $n_e$ as Henrich & Chandrasekhar, but they constructed models in which the ratio of the molecular weight in the core to the molecular weight in the envelope is ${\mu }_c/{\mu }_e=2$. This was done to more realistically represent stars as they evolve off the main sequence; they have inert, isothermal helium cores and envelopes that are rich in hydrogen. Note that introducing a discontinuous drop in the mean molecular weight at the core-envelope interface also introduces a discontinuous drop in the gas density across the interface. As the following excerpt from their article summarizes, in these models, 📚 Schönberg & Chandrasekhar (1942) found that ${\nu }_{\mathrm{m}\mathrm{a}\mathrm{x}}\approx 0.101$. This is commonly referred to as the Schönberg-Chandrasekhar mass limit, although it was Henrich & Chandrasekhar who were the first to identify the instability.

Text excerpt from p. 168 of M. Schönberg & S. Chandrasekhar (1942) On the Evolution of the Main-Sequence Stars The Astrophysical Journal, Vol. 96, pp. 161 - 172

"For the models with isothermal cores we have the following properties:  (a) There are no equilibrium configurations with cores containing … more than 0.101 of the stellar mass … [This] "upper limit is a decreasing function of ${\mu }_c/{\mu }_e$, since it is $\sim 0.35$ for the case of equal molecular weights, as was shown by" 📚 Henrich & Chandrasekhar (1941).

In an effort to develop a more complete appreciation of the onset of the instability associated with the Schönberg-Chandrasekhar mass limit, 📚 M. Beech (1988, Astrophysics and Space Science, Vol. 147, issue 2, pp. 219 - 227) matched an analytically prescribable, $n_e=1$ polytropic envelope to an isothermal core and, like Schönberg & Chandrasekhar, allowed for a discontinuous change in the molecular weight at the interface. [For an even more comprehensive generalization and discussion, see 📚 W. H. Ball, C. A. Tout, & A. N. Żytkow (2012, MNRAS, Vol. 421, issue 3, pp. 2713 - 2721)]. Beech's results were not significantly different from those reported by 📚 Schönberg & Chandrasekhar (1942); in particular, the value of ${\nu }_{\mathrm{m}\mathrm{a}\mathrm{x}}$ was still only definable numerically because an isothermal core cannot be described in terms of analytic functions.

In an accompanying derivation [see, also, 📚 P. P. Eggleton, J. Faulkner, & R. C. Cannon (1998, MNRAS, Vol. 298, issue 3, pp. 831 - 834)] we have gone one step farther, matching an analytically prescribable, $n_e=1$ polytropic envelope to an analytically prescribable, $n_c=5$ polytropic core. For this bipolytrope, we show that there is a limiting mass-fraction, ${\nu }_{\mathrm{m}\mathrm{a}\mathrm{x}}$, for any choice of the molecular weight ratio ${\mu }_c/{\mu }_e>3$ and that the interface location, ${\xi }_i$, associated with this critical configuration is given by the positive, real root of the following relation:

$(\frac{\pi }{2}+{\tan }^{-1}{\mathrm{\Lambda }}_i)(1+{\ell }_i^2)[3+(1-m_3{)}^2(2-{\ell }_i^2){\ell }_i^2]-m_3{\ell }_i[(1-m_3){\ell }_i^4-(m_3^2-m_3+2){\ell }_i^2-3]=0,$

where,

${\ell }_i\equiv \frac{{\xi }_i}{\sqrt{3}};$             $m_3\equiv 3(\frac{{\mu }_c}{{\mu }_e}{)}^{-1};$       and       ${\mathrm{\Lambda }}_i\equiv \frac{1}{m_3{\ell }_i}[1+(1-m_3){\ell }_i^2].$

Close Binary Stars

Preface

Text extracted from p. 162 of G. H. Darwin (1906) IV.  On the figure and stability of a liquid satellite Philosophical Transactions of the Royal Society A, Vol. 206, pp. 161 - 248

"More than half a century ago Édouard Roche wrote his celebrated paper on the form assumed by a liquid satellite when revolving, without relative motion, about a solid planet.* In consequence of the singular modesty of Roche's style, and also because the publication was made at Montpellier, this paper seems to have remained almost unnoticed for many years, but it has ultimately attained its due position as a classical memoir." "The laborious computations necessary for obtaining numerical results were carried out, partly at least, by graphical methods. Verification of the calculations, which as far as I know have never been repeated, forms part of the work of the present paper. The distance from a spherical planet which has been called 'Roche's limit' is expressed by the number of planetary radii in the radius vector of the nearest possible infinitesimal liquid satellite, of the same density as the planet, revolving so as always to present the same aspect to the planet. Our moon, if it were homogeneous, would have the form of one of Roche's ellipsoids; but its present radius vector is of course far greater than the limit. Roche assigned to the limit in question the numerical value 2.44; it the present paper I show that the true value is 2.455, and the closeness of the agreement with the previously accepted value affords a remarkable testimony to the accuracy with which he must have drawn his figures." *La figure d'une masse fluide soumise à l'attraction d'un point éloigné, 'Acad. des Sci. de Montpellier,' vol. 1, 1847-50, p. 243. Roche References The following articles by Édouard Roche have been found via "Google Book" searches: Tome Premier (1847 - 1850): (1848, pp. 113 - 124) — Calcul de l'Inégalité parallactique (1849, pp. 245 - 262) — Mémoire sur la figure d'une masse fluide, soumise à 'attraction d'un point éloigné. (1re Partie.) (1850, pp. 333 - 348) — Mémoire sur la figure d'une masse fluide, soumise à 'attraction d'un point éloigné. (2e Partie.) Tome Second (1851 - 1854): (1851, pp. 21 - 32) — Mémoire sur la figure d'une masse fluide, soumise à 'attraction d'un point éloigné.

Simplistic Illustrative Model

Using Our Adopted Notation

Consider the simple model of two spherical stars in circular orbit about one another, as depicted here on the right. In addition to the physical parameters explicitly labeled in this diagram, we adopt the following variable notation:

The total system mass is, $M_{\mathrm{t}\mathrm{o}\mathrm{t}}\equiv M+M^{\text{'}};$ The ratio of the primary to secondary mass is, $\lambda \equiv \frac{M}{M^{\text{'}}};$ And the separation between the two centers is, $d\equiv r_{\mathrm{c}\mathrm{m}}+r_{\mathrm{c}\mathrm{m}}^{\text{'}}.$

For a circular orbit, the angular velocity is related to the the system mass and separation via the Kepler relation,

${\omega }^2{d}^3=G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}},$

and the distances, $r_{\mathrm{c}\mathrm{m}}$ and $r_{\mathrm{c}\mathrm{m}}^{\text{'}}$, between the center of each star and the center of mass (cm) of the system must be related to one another via the expression,

$\frac{r_{\mathrm{c}\mathrm{m}}^{\text{'}}}{r_{\mathrm{c}\mathrm{m}}}$

$=$

$\frac{M}{M^{\text{'}}}=\lambda .$

Note that the following relations also hold:

$M=M_{\mathrm{t}\mathrm{o}\mathrm{t}}\left(\frac{\lambda }{1+\lambda }\right)$	and	$M^{\text{'}}=M_{\mathrm{t}\mathrm{o}\mathrm{t}}\left(\frac{1}{1+\lambda }\right)$
$r_{\mathrm{c}\mathrm{m}}=d\left(\frac{1}{1+\lambda }\right);$	and	$r_{\mathrm{c}\mathrm{m}}^{\text{'}}=d\left(\frac{\lambda }{1+\lambda }\right).$

Hence, the orbital angular momentum is,

$L_{\mathrm{o}\mathrm{r}\mathrm{b}}$	$=$	$[M{r}_{\mathrm{c}\mathrm{m}}^2+M^{\text{'}}(r_{\mathrm{c}\mathrm{m}}^{\text{'}}{)}^2]\omega$
	$=$	$M_{\mathrm{t}\mathrm{o}\mathrm{t}}{d}^2\left[\right(\frac{\lambda }{1+\lambda }\left)\right(\frac{1}{1+\lambda }{)}^2+\left(\frac{1}{1+\lambda }\right)\left(\frac{\lambda }{1+\lambda }{)}^2\right][\frac{G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}}{d^3}{]}^{1/2}$
	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}d{)}^{1/2}[\frac{\lambda }{(1+\lambda {)}^2}].$

Assuming that both stars are rotating synchronously with the orbit, their respective spin angular momenta are,

$L_M=I_M\omega$	$=$	$\frac{2}{5}M{R}^2\omega$
	$=$	$\frac{2}{5}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}\left(\frac{\lambda }{1+\lambda }\right)R^2[\frac{G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}}{d^3}{]}^{1/2}$
	$=$	$\frac{2}{5}(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}d{)}^{1/2}(\frac{\lambda }{1+\lambda }\left)\right(\frac{R}{d}{)}^2,$
$L_{M^{\text{'}}}=I_{M^{\text{'}}}\omega$	$=$	$\frac{2}{5}{M}^{\text{'}}(R^{\text{'}}{)}^2\omega$
	$=$	$\frac{2}{5}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}\left(\frac{1}{1+\lambda }\right)(R^{\text{'}}{)}^2[\frac{G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}}{d^3}{]}^{1/2}$
	$=$	$\frac{2}{5}(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}d{)}^{1/2}(\frac{1}{1+\lambda }\left)\right(\frac{R^{\text{'}}}{d}{)}^2.$

Hence, the total angular momentum of the system is,

$L_{\mathrm{t}\mathrm{o}\mathrm{t}}=L_{\mathrm{o}\mathrm{r}\mathrm{b}}+L_M+L_{M^{\text{'}}}$	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}d{)}^{1/2}[\frac{\lambda }{(1+\lambda {)}^2}]+\frac{2}{5}(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}d{)}^{1/2}\left(\frac{\lambda }{1+\lambda }\right)(\frac{R}{d}{)}^2$
		$+\frac{2}{5}(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}d{)}^{1/2}(\frac{1}{1+\lambda }\left)\right(\frac{R^{\text{'}}}{d}{)}^2$
	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}R{)}^{1/2}\{\left[\frac{\lambda }{(1+\lambda {)}^2}\right](\frac{d}{R}{)}^{1/2}+\frac{2}{5}(\frac{1}{1+\lambda }\left)\right[\lambda +\left(\frac{R^{\text{'}}}{R}{)}^2\right]\left(\frac{d}{R}{)}^{-3/2}\right\}$

If we assume that the two stars have the same (uniform) densities, $\rho$, then, following 📚 Darwin (1906) — see immediately below — the two stellar radii can be related to the mass ratio, $\lambda$ via the expressions,

$R=a\left(\frac{\lambda }{1+\lambda }\right),$

and

$R^{\text{'}}=a\left(\frac{1}{1+\lambda }\right),$

where, the characteristic length scale is,

$a\equiv (\frac{3M_{\mathrm{t}\mathrm{o}\mathrm{t}}}{4\pi \rho }{)}^{1/3}.$

Replacing $R$ and $R^{\text{'}}$ by these expressions in our equation for $L_{\mathrm{t}\mathrm{o}\mathrm{t}}$ results in the simplistic/illustrative expression for $L_1$, derived by 📚 Darwin (1906) and presented below. Darwin's expression for $L_2$ is obtained by using the same expression for $R$ but treating the secondary as a point mass, that is, setting $R^{\text{'}}=0$.

Darwin's (1906) Equivalent Illustration

From Pt. I, §1 (p. 164) of 📚 Darwin (1906) — verbatum text in green: It will be useful to make a rough preliminary investigation of the regions in which we shall have to look for cases of limiting stability in the two problems. For this purpose I consider [1] the case of two spheres as the analogue of [Darwin's] problem of the figure of equilibrium, and [2] the case of a sphere and a particle as the analogue of Roche's problem.

… let the mass of the whole system be $M_{\mathrm{t}\mathrm{o}\mathrm{t}}=\frac{4}{3}\pi \rho a^3$; let the masses of the two spheres be

$M=M_{\mathrm{t}\mathrm{o}\mathrm{t}}\left[\frac{\lambda }{1+\lambda }\right]$

and

$M^{\text{'}}=M_{\mathrm{t}\mathrm{o}\mathrm{t}}\left[\frac{1}{1+\lambda }\right]$

$\Rightarrow$

$\lambda =\frac{M}{M^{\text{'}}}$

or for Roche's problem let the latter $(M^{\text{'}})$ be the mass of the particle. Assuming that both spheres have the same characteristic density, $\rho$, that has been used to specify the total mass, we furthermore know that the radius of the first sphere is,

$a_M$

$=$

$(\frac{M}{\frac{4}{3}\pi \rho }{)}^{1/3}=a(\frac{\lambda }{1+\lambda }{)}^{1/3},$

and (for the Darwin problem) the radius of the second sphere is,

$a_{M^{\text{'}}}$

$=$

$(\frac{M^{\text{'}}}{\frac{4}{3}\pi \rho }{)}^{1/3}=a(\frac{1}{1+\lambda }{)}^{1/3}.$

Let $r$ be the distance from the centre of one sphere to that of the other, or to the particle, as the case may be; and $\omega$ the orbital angular velocity, where

${\omega }^2{r}^3$

$=$

$G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}.$

The centre of inertia of the two masses is distant $r/(1+\lambda )$ and $\lambda r/(1+\lambda )$ from their respective centres, and we easily find the orbital momentum to be

$L_{\mathrm{o}\mathrm{r}\mathrm{b}}$	$=$	$M(\frac{r}{1+\lambda }{)}^2\omega +M^{\text{'}}(\frac{\lambda r}{1+\lambda }{)}^2\omega$
	$=$	$M_{\mathrm{t}\mathrm{o}\mathrm{t}}\left[\right(\frac{\lambda }{1+\lambda }\left)\right(\frac{1}{1+\lambda }{)}^2\omega +\left(\frac{1}{1+\lambda }\right)\left(\frac{\lambda }{1+\lambda }{)}^2\right]r^2\omega$
	$=$	$M_{\mathrm{t}\mathrm{o}\mathrm{t}}\left[\frac{\lambda }{(1+\lambda {)}^2}\right]r^2\omega .$

In both problems the rotational momentum of the first sphere is

$L_{\mathrm{M}}=\frac{2}{5}M{a}_M^2\omega$

$=$

$\frac{2}{5}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}\left(\frac{\lambda }{1+\lambda }\right)(\frac{\lambda }{1+\lambda }{)}^{2/3}{a}^2\omega =\frac{2}{5}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}(\frac{\lambda }{1+\lambda }{)}^{5/3}{a}^2\omega .$

In the [Darwin] problem the rotational momentum of the second sphere is

$L_{{\mathrm{M}}^{\mathrm{\text{'}}}}=\frac{2}{5}{M}^{\text{'}}{a}_{M^{\text{'}}}^2\omega$

$=$

$\frac{2}{5}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}\left(\frac{1}{1+\lambda }\right)(\frac{1}{1+\lambda }{)}^{2/3}{a}^2\omega =\frac{2}{5}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}(\frac{1}{1+\lambda }{)}^{5/3}{a}^2\omega ,$

and in the [Roche] problem it is nil.

If, then, we write $L_1$ for the total angular momentum of the two spheres, and $L_2$ for that of the sphere and particle, we have

$L_1$	$=$	$L_{\mathrm{o}\mathrm{r}\mathrm{b}}+L_M+L_{M^{\text{'}}}$
	$=$	$M_{\mathrm{t}\mathrm{o}\mathrm{t}}\left[\frac{\lambda }{(1+\lambda {)}^2}\right]r^2\omega +\frac{2}{5}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}(\frac{\lambda }{1+\lambda }{)}^{5/3}{a}^2\omega +\frac{2}{5}{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}(\frac{1}{1+\lambda }{)}^{5/3}{a}^2\omega$
	$=$	$M_{\mathrm{t}\mathrm{o}\mathrm{t}}{a}^2\omega [\frac{\lambda r^2}{(1+\lambda {)}^2{a}^2}+\frac{2}{5}\frac{1+{\lambda }^{5/3}}{(1+\lambda {)}^{5/3}}],$
$L_2$	$=$	$L_{\mathrm{o}\mathrm{r}\mathrm{b}}+L_M$
	$=$	$M_{\mathrm{t}\mathrm{o}\mathrm{t}}{a}^2\omega [\frac{\lambda r^2}{(1+\lambda {)}^2{a}^2}+\frac{2}{5}\frac{{\lambda }^{5/3}}{(1+\lambda {)}^{5/3}}].$

For comparison, in the context of the LRS93b discussion of Compressible Roche Ellipsoids, the total angular momentum is, $J$ $=$ $[\frac{M{r}^2}{(1+\lambda )}+I]\mathrm{\Omega },$ where (see their equation 4.8), $I$ $=$ $\frac{1}{5}{\kappa }_n(a_1^2+a_2^2)M,$ and, for incompressible configurations, ${\kappa }_{n=0}=1$. This gives, $J$ $=$ $M_{\mathrm{t}\mathrm{o}\mathrm{t}}\left(\frac{\lambda }{1+\lambda }\right)a_1^2\mathrm{\Omega }[\frac{r^2}{(1+\lambda )a_1^2}+\frac{1}{5}(1+\frac{a_2^2}{a_1^2}\left)\right].$

On substituting for $\omega$ its value in terms of $r$, these expressions become

$L_1$	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}a{)}^{1/2}(\frac{a}{r}{)}^{3/2}[\frac{2}{5}\frac{1+{\lambda }^{5/3}}{(1+\lambda {)}^{5/3}}+\frac{\lambda r^2}{(1+\lambda {)}^2{a}^2}]$
	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}a{)}^{1/2}\frac{1}{(1+\lambda {)}^2}[\frac{2}{5}(1+{\lambda }^{5/3})(1+\lambda {)}^{1/3}(\frac{a}{r}{)}^{3/2}+\lambda \left(\frac{r}{a}{)}^{1/2}\right],$
$L_2$	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}a{)}^{1/2}(\frac{a}{r}{)}^{3/2}[\frac{2}{5}\frac{{\lambda }^{5/3}}{(1+\lambda {)}^{5/3}}+\frac{\lambda r^2}{(1+\lambda {)}^2{a}^2}]$
	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}a{)}^{1/2}\frac{1}{(1+\lambda {)}^2}[\frac{2}{5}{\lambda }^{5/3}(1+\lambda {)}^{1/3}(\frac{a}{r}{)}^{3/2}+\lambda \left(\frac{r}{a}{)}^{1/2}\right].$

After setting $G=1$, this matches the pair of equations that appears immediately following equation (1) in 📚 Darwin (1906).

In the following figure we have plotted, for both problems, how the total angular momentum varies with orbital separation. In order to facilitate a direct comparison with Figure 1 from LRS, in place of the dimensionless separation $r/a$ we plot along the abscissa the quantity, $r_{\mathrm{L}\mathrm{R}\mathrm{S}}/R$, where, $r_{\mathrm{L}\mathrm{R}\mathrm{S}}$ is the radius of the circular orbit and $R$ is the radius of the primary star; that is,

$r_{\mathrm{L}\mathrm{R}\mathrm{S}}=\frac{1}{2}r$	and	$R=a_M$	$\Rightarrow$	$\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}=\frac{1}{2}(\frac{1+\lambda }{\lambda }{)}^{1/3}\frac{r}{a}.$
$r_{\mathrm{L}\mathrm{R}\mathrm{S}}=\frac{1}{2}r$	and	$R=a_M$	$\Rightarrow$	$\frac{r}{a}=2(\frac{\lambda }{1+\lambda }{)}^{1/3}\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}.$

Rewriting Darwin's pair of angular momentum expressions in terms of this preferred dimensionless separation, we have,

$L_1$	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}a{)}^{1/2}\frac{1}{(1+\lambda {)}^2}\{\frac{2}{5}(1+{\lambda }^{5/3})(1+\lambda {)}^{1/3}[2(\frac{\lambda }{1+\lambda }{)}^{1/3}\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}{]}^{-3/2}+\lambda [2\left(\frac{\lambda }{1+\lambda }{)}^{1/3}\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}{]}^{1/2}\right\}$
	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}a{)}^{1/2}\frac{1}{(1+\lambda {)}^2}\{(\frac{1}{2\cdot 5^2}{)}^{1/2}(1+{\lambda }^{5/3})(1+\lambda {)}^{1/3}\left(\frac{1+\lambda }{\lambda }{)}^{1/2}\right(\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}{)}^{-3/2}+2^{1/2}\lambda \left(\frac{\lambda }{1+\lambda }{)}^{1/6}\right(\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}{)}^{1/2}\},$
$L_2$	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}a{)}^{1/2}\frac{1}{(1+\lambda {)}^2}\{\frac{2}{5}{\lambda }^{5/3}(1+\lambda {)}^{1/3}[2(\frac{\lambda }{1+\lambda }{)}^{1/3}\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}{]}^{-3/2}+\lambda [2\left(\frac{\lambda }{1+\lambda }{)}^{1/3}\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}{]}^{1/2}\right\}$
	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}a{)}^{1/2}\frac{1}{(1+\lambda {)}^2}\{(\frac{1}{2\cdot 5^2}{)}^{1/2}{\lambda }^{5/3}(1+\lambda {)}^{1/3}\left(\frac{1+\lambda }{\lambda }{)}^{1/2}\right(\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}{)}^{-3/2}+2^{1/2}\lambda \left(\frac{\lambda }{1+\lambda }{)}^{1/6}\right(\frac{r_{\mathrm{L}\mathrm{R}\mathrm{S}}}{R}{)}^{1/2}\}.$

Note that the two binary components come into contact when, for the Darwin problem,

$a_M+a_{M^{\text{'}}}=r$

$\Rightarrow$

$\frac{r}{a}=\frac{1+{\lambda }^{1/3}}{(1+\lambda {)}^{1/3}};$

and, for the Roche problem,

$a_M=r$

$\Rightarrow$

$\frac{r}{a}=\frac{{\lambda }^{1/3}}{(1+\lambda {)}^{1/3}}.$

Setup

Jeans (1919)

From § 50 (p. 46) of J. H. Jeans (1919) — verbatum text in green: Let the two bodies be spoken of as primary and secondary, and let their masses be $M$, $M^{\text{'}}$ respectively; let the distance apart of their centres of gravity be $R$, and let the angular velocity of rotation of the line joining them be $\omega$. It will be sufficient to fix our attention on the conditions of equilibrium of one of the two masses, say the primary. Let its centre of gravity be taken as origin, let the line joining it to the centre of the secondary be axis of $x$, and let the plane in which the rotation takes place be that of $xy$. Then the equation of the axis of rotation is

$x=\frac{M^{\text{'}}}{M+M^{\text{'}}}R$

and

$y=0.$

The problem may be reduced to a statical one (cf. § 31) by supposing the masses acted on by a field of force of [the centrifugal] potential

$\frac{1}{2}{\omega }^2\left[\right(x-\frac{M^{\text{'}}}{M+M^{\text{'}}}R{)}^2+y^2].$

Chandrasekhar (1969)

From pp. 189-190 of [EFE] — verbatum text in green: Let the masses of the primary and the secondary be $M$ and $M^{\text{'}}$, respectively; let the distance between their centers of mass be $R$; and let the constant angular velocity of rotation about their common center of mass be $\mathrm{\Omega }$. Choose a coordinate system in which the origin is at the center of mass of the primary, the $x_1-$axis points to the center of mass of the secondary, and the $x_3-$axis is parallel to the direction of $\overrightarrow{\mathrm{\Omega }}$. In this coordinate system, the equation of motion governing fluid elements of $M$ includes (see EFE's equation 1) a gradient of the centrifugal potential,

$\frac{1}{2}{\mathrm{\Omega }}^2\left[\right(x_1-\frac{M^{\text{'}}R}{M+M^{\text{'}}}{)}^2+x_2^2].$

Tassoul (1978)

From p. 449 of [T78] — verbatum text in green: Let the masses of the primary and the secondary be $M$ and $M^{\text{'}}$, respectively; let the distance between their centers of mass be $d$; and let the angular velocity of rotation about their common center of mass be $\mathrm{\Omega }$. Next choose a system of reference in which the origin is at the center of mass of the primary; for convenience, the $x_1-$axis points toward the center of mass of the secondary, and the $x_3-$axis is parallel to the direction of $\overrightarrow{\mathrm{\Omega }}$. Then, the equation of the rotation axis, which of course passes through the center of mass of the two bodies, is

$x_1=\frac{M^{\text{'}}}{M+M^{\text{'}}}d$

and

$x_2=0.$

Accordingly, the centrifugal force acting on the mass $M$ may be derived from the potential

$-\frac{1}{2}{\mathrm{\Omega }}^2\left[\right(x_1-\frac{M^{\text{'}}}{M+M^{\text{'}}}d{)}^2+x_2^2].$

Roche Ellipsoids

Jeans (1919)

From § 51 (p. 47) of J. H. Jeans (1919) — verbatum text in green: The simplest problem occurs when the secondary may be treated as a rigid sphere; this is the special problem dealt with by Roche. As in § 47 the tide-generating potential acting on the primary may be supposed to be

$\frac{M^{\text{'}}}{R}+\frac{M^{\text{'}}}{R^2}x+\frac{M^{\text{'}}}{R^3}(x^2-\frac{1}{2}{y}^2-\frac{1}{2}{z}^2)+\cdots$

We shall for the present be content to omit all terms beyond those written down. The correction required by the neglect of these terms will be discussed later, and will be found to be so small that the results now to be obtained are hardly affected.

On omitting these terms, and combining the two potentials … it appears that the primary may be supposed influenced by a statical field of potential

$\frac{M^{\text{'}}}{R}x(1-\frac{{\omega }^2{R}^3}{M+M^{\text{'}}})+\frac{M^{\text{'}}}{R^3}(x^2-\frac{1}{2}{y}^2-\frac{1}{2}{z}^2)+\frac{1}{2}{\omega }^2(x^2+y^2).$

The terms in $x$ may immediately be removed by supposing $\omega$ to have the appropriate value given by

${\omega }^2$

$=$

$\frac{M+M^{\text{'}}}{R^3}$

and the condition for equilibrium is now seen to be that we shall have, at every point of the surface,

$V_b+\mu (x^2-\frac{1}{2}{y}^2-\frac{1}{2}{z}^2)+\frac{1}{2}{\omega }^2(x^2+y^2)$

$=$

constant

where $\mu$ … stands for $M^{\text{'}}/R^3$ .

Chandrasekhar (1969)

From p. 190 of [EFE] — verbatum text in green: In Roche's particular problem, the secondary is treated as a rigid sphere. Then, over the primary, the tide-generating potential, ${\mathfrak{𝔅}}^{\text{'}}$ can be expanded in the form

${\mathfrak{𝔅}}^{\text{'}}$

$=$

$\frac{G{M}^{\text{'}}}{R}(1+\frac{x_1}{R}+\frac{x_1^2-\frac{1}{2}{x}_2^2-\frac{1}{2}{x}_3^2}{R^2}+\cdots );$

and the approximation which underlies this theory is to retain, in this expansion for ${\mathfrak{𝔅}}^{\text{'}}$, only the terms which have been explicitly written down and ignore all the terms which are of higher order. On this assumption, the equation of motion becomes

$\frac{d{u}_i}{dt}+\frac{1}{\rho }\frac{\partial p}{\partial x_i}$

$=$

$\frac{\partial }{\partial x_i}[\mathfrak{𝔅}+\frac{1}{2}{\mathrm{\Omega }}^2(x_1^2+x_2^2)+\mu (x_1^2-\frac{1}{2}{x}_2^2-\frac{1}{2}{x}_3^2)+(\frac{G{M}^{\text{'}}}{R^2}-\frac{M^{\text{'}}R}{M+M^{\text{'}}}{\mathrm{\Omega }}^2\left)x_1\right]+2\mathrm{\Omega }{ϵ}_{i\ell 3}{u}_{\ell }$

where we have introduced the abbreviation

$\mu =\frac{G{M}^{\text{'}}}{R^3}.$

So far, we have left ${\mathrm{\Omega }}^2$ unspecified. If we now let ${\mathrm{\Omega }}^2$ have the "Keplerian value"

${\mathrm{\Omega }}^2=\frac{G(M+M^{\text{'}})}{R^3}=\mu (1+\frac{M}{M^{\text{'}}}),$

the "unwanted" term in $x_1$, on the right-hand side of [this equation,] vanishes and we are left with

$\frac{d{u}_i}{dt}+\frac{1}{\rho }\frac{\partial p}{\partial x_i}$

$=$

$\frac{\partial }{\partial x_i}[\mathfrak{𝔅}+\frac{1}{2}{\mathrm{\Omega }}^2(x_1^2+x_2^2)+\mu (x_1^2-\frac{1}{2}{x}_2^2-\frac{1}{2}{x}_3^2)]+2\mathrm{\Omega }{ϵ}_{i\ell 3}{u}_{\ell }.$

This is the basic equation of this theory; and Roche's problem is concerned with the equilibrium and the stability of homogeneous masses governed by [this relation].

Tassoul (1978)

From pp. 449-450 of [T78] — verbatum text in green: In Roche's particular problem, the secondary is treated as a rigid sphere; hence, over the primary, the tide-generating potential can be expanded in the form

$-\frac{G{M}^{\text{'}}}{d}(1+\frac{x_1}{d}+\frac{x_1^2-\frac{1}{2}{x}_2^2-\frac{1}{2}{x}_3^2}{d^2}+\cdots ).$

The approximation that underlies the theory is to omit all terms beyond those written down. On this assumption, we find that, apart from its own gravitation, the primary may be supposed to be acted upon by a total field of force derived from the potential

$-\frac{1}{2}{\mathrm{\Omega }}^2(x_1^2+x_x^2)-\mu (x_1^2-\frac{1}{2}{x}_2^2-\frac{1}{2}{x}_3^2)-(\mu -\frac{M^{\text{'}}}{M+M^{\text{'}}}{\mathrm{\Omega }}^2)d{x}_1,$

where

$\mu$

$=$

$\frac{G{M}^{\text{'}}}{d^3}.$

Further letting ${\mathrm{\Omega }}^2$ have its "Keplerian value"

${\mathrm{\Omega }}^2$

$=$

$\frac{G(M+M^{\text{'}})}{d^3},$

we can thus write the conditions of relative equilibrium for the primary in the form

$\frac{1}{\rho }\mathrm{\nabla }p$

$=$

$-\mathrm{\nabla }[V-\frac{1}{2}{\mathrm{\Omega }}^2(x_1^2+x_x^2)-\mu (x_1^2-\frac{1}{2}{x}_2^2-\frac{1}{2}{x}_3^2)],$

where $V$ is the self-gravitating potential of the primary.

Incompressible Roche Ellipsoids (λ ≠ 0)

Let's see if we can understand the relationship between tabulated data presented by Lai, Rasio, & Shapiro (1993b, ApJS, 88, 205) — hereafter, LRS93S — for the case of incompressible Roche ellipsoids, when $\lambda =p=1$. After setting ${\kappa }_{n=0}=1$ in their equation (4.8), we have,

$I$

$=$

$\frac{1}{5}{a}_1^{2}M(1+\frac{a_2^2}{a_1^2}).$

And, from their equation (7.12),

$J_{\mathrm{t}\mathrm{o}\mathrm{t}}$

$=$

$[\frac{M{r}^2}{(1+p)}+I]\mathrm{\Omega }.$

If we adopt the Keplerian orbital frequency, the expression for the total angular momentum is,

$J_{\mathrm{K}\mathrm{e}\mathrm{p}}$	$=$	$[\frac{M{r}^2}{(1+p)}+\frac{1}{5}{a}_1^{2}M(1+\frac{a_2^2}{a_1^2}\left)\right][\frac{G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}}{r^3}{]}^{1/2}$
	$=$	$(G{M}^3{)}^{1/2}{a}_1^2[\frac{1}{(1+p)}(\frac{r}{a_1}{)}^2+\frac{1}{5}(1+\frac{a_2^2}{a_1^2}\left)\right]\left[\frac{1}{r^3}\right(\frac{1+p}{p}){]}^{1/2}$
	$=$	$(G{M}^{3}R{)}^{1/2}(\frac{a_1^4}{r^{3}R}{)}^{1/2}\left[\frac{1}{(1+p)}\right(\frac{r}{a_1}{)}^2+\frac{1}{5}(1+\frac{a_2^2}{a_1^2})\left]\right[\left(\frac{1+p}{p}\right){]}^{1/2}.$

For the specific case of an equal-mass binary sequence — that is, $\lambda =p=1$ — as considered in the following table, we have,

${\overline{J}}_{\mathrm{K}\mathrm{e}\mathrm{p}}{|}_{p=1}\equiv (G{M}^{3}R{)}^{-1/2}{J}_{\mathrm{K}\mathrm{e}\mathrm{p}}{|}_{p=1}$	$=$	$2^{-1/2}\left(\frac{a_1}{r}{)}^2\right(\frac{r^3}{a_1^3}\cdot \frac{a_1^3}{R^3}{)}^{1/6}\left[\right(\frac{r}{a_1}{)}^2+\frac{2}{5}(1+\frac{a_2^2}{a_1^2})]$
	$=$	$2^{-1/2}\left(\frac{a_1}{r}{)}^{3/2}\right(\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}{)}^{-1/6}\left[\right(\frac{r}{a_1}{)}^2+\frac{2}{5}(1+\frac{a_2^2}{a_1^2})],$

where we also have involved the expression for an equivalent spherical radius given just before their equation (7.21), namely,

$R^3$

$=$

$a_1{a}_2{a}_3=a_1^3(\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}).$

Table 1:  Incompressible $(n=0)$ Roche Ellipsoids with $\lambda =p=1$
Extracted from Table 1 of Chandrasekhar (1963) same as [EFE] Table XVI				EFE Check
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)
${\cos }^{-1}(a_3/a_1)$	$a_2/a_1$	$a_3/a_1$	${\mathrm{\Omega }}^2$	$r/a_1$	${\overline{J}}_{\mathrm{K}\mathrm{e}\mathrm{p}}=L_{\mathrm{t}\mathrm{o}\mathrm{t}}/(G{M}^{3}R{)}^{1/2}$	$r/R$	$\mathfrak{𝔍}$	$L_{\mathrm{t}\mathrm{o}\mathrm{t}}/(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}R{)}^{1/2}$
12°	0.98660	0.97815	0.009293	6.5181	1.8498	6.5959	1.0104	0.6540
24°	0.94376	0.91355	0.036152	3.9916	1.5168	4.1938	1.0436	0.5363
36°	0.86345	0.80902	0.076342	2.9005	1.3846	3.2689	1.1086	0.4895
48°	0.73454	0.66913	0.118726	2.2266	1.3353	2.8215	1.2360	0.4721
54°	0.64956	0.58779	0.134284	1.9645	1.3351	2.7080	1.3510	0.4720
59°	0.56892	0.51504	0.140854	1.7702	1.3494	2.6652	1.5003	0.4771
60°	0.55186	0.50000	0.141250	1.7335	1.3542	2.6627	1.5390	0.4788
61°	0.53451	0.48481	0.141298	1.6974	1.3597	2.6624	1.5816	0.4807
66°	0.44429	0.40674	0.135785	1.5253	1.4006	2.6980	1.8732	0.4952
69°	0.38813	0.35837	0.127424	1.4388	1.4278	2.7557	2.1105	0.5073
71°	0.35022	0.32557	0.119625	1.3647	1.4723	2.8144	2.3873	0.5205
72°	0.33119	0.30902	0.115054	1.3337	1.4919	2.8512	2.5357	0.5275
73°	0.31213	0.29237	0.110044	1.3028	1.5140	2.8938	2.7072	0.5353
75°	0.27405	0.25882	0.098753	1.2419	1.5663	3.0001	3.1371	0.5538
78°	0.21726	0.20791	0.078934	1.1513	1.6731	3.2327	4.1282	0.5915
81°	0.16126	0.15643	0.056499	1.0599	1.8353	3.6139	5.9641	0.6489

Table 2:  Incompressible $(n=0)$ Roche Ellipsoids with $\lambda =p=1$
LRS93 Supplements							LRS93 Check
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)	(9)	(10)
$r/a_1$	$r/R$	$a_2/a_1$	$a_3/a_1$	$\overline{\mathrm{\Omega }}$	$\overline{J}$	$\overline{E}$	$\left(\frac{r{p}^{1/3}}{a_1}{)}^3\right(\frac{R}{r{p}^{1/3}}{)}^3(\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}{)}^{-1}$	${\overline{\mathrm{\Omega }}}_{\mathrm{K}\mathrm{e}\mathrm{p}}$	${\overline{J}}_{\mathrm{K}\mathrm{e}\mathrm{p}}$
5.0	5.131	0.9707	0.9533	0.1406	1.653	-0.6943	1.0000	0.1405	1.6515
4.0	4.202	0.9441	0.9139	0.1901	1.522	-0.7128	0.9998	0.1896	1.5180
3.0	3.348	0.8750	0.8222	0.2690	1.408	-0.7349	1.0001	0.2666	1.3954
2.7	3.124	0.8345	0.7738	0.3000	1.386	-0.7404	0.9998	0.2958	1.3661
2.5	2.989	0.7981	0.7330	0.3222	1.377	-0.7427	1.0002	0.3160	1.3506
2.380	2.916	0.7715	0.7044	0.3358	1.375	-0.7432	1.0005	0.3279	1.3436
2.2	2.821	0.7236	0.6553	0.3556	1.380	-0.7418	1.0003	0.3446	1.3373
2.112	2.783	0.6960	0.6281	0.3648	1.386	-0.7399	0.9998	0.3518	1.3366
2.0	2.744	0.6561	0.5901	0.3753	1.399	-0.7358	1.0001	0.3592	1.3389
1.801	2.713	0.5708	0.5123	0.3886	1.441	-0.7218	1.0004	0.3654	1.3552
1.697	2.724	0.5184	0.4664	0.3908	1.477	-0.7097	1.0000	0.3632	1.3727
1.6	2.759	0.4644	0.4198	0.3884	1.524	-0.6939	1.0004	0.3563	1.3977
1.5	2.831	0.4040	0.3682	0.3798	1.590	-0.6717	1.0000	0.3428	1.4358
1.0	5.312	0.0823	0.0811	0.1685	2.888	-0.3772	0.9996	0.1334	2.2859

Digesting LRS93S Results

In column (12) of the above table we have combined the data from columns (5), (6), (7), & (8) to demonstrate that LRS93S indeed used the definition of $R$, given above for their incompressible, Roche ellipsoid configurations; these terms should zombie to give unity, and they appear to do so, within the accuracy presented by the data from LRS93S. In column 14 of the table, we have listed the value of ${\overline{J}}_{\mathrm{K}\mathrm{e}\mathrm{p}}$, as given by the above expression. Column 13 lists ${\overline{\mathrm{\Omega }}}_{\mathrm{K}\mathrm{e}\mathrm{p}}$, as follows:

${\mathrm{\Omega }}_K$	$=$	$[\frac{G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}}{r^3}{]}^{1/2}$
	$=$	$\left[\frac{GM}{R^3}\right(\frac{1+p}{p}\left)\right(\frac{R}{r}{)}^3{]}^{1/2}$
	$=$	$\left[\frac{4\pi G\rho }{3}{]}^{1/2}\right(\frac{1+p}{p}{)}^{1/2}\left[\right(\frac{a_1}{r}{)}^3\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}{]}^{1/2}$
$\Rightarrow {\overline{\mathrm{\Omega }}}_{\mathrm{K}\mathrm{e}\mathrm{p}}{|}_{p=1}\equiv \frac{{\mathrm{\Omega }}_{\mathrm{K}\mathrm{e}\mathrm{p}}}{(\pi G\rho {)}^{1/2}}{|}_{p=1}$	$=$	$\left[\frac{8}{3}\right(\frac{a_1}{r}{)}^3\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}{]}^{1/2}.$

The value of the LRS93S correction factor, $(1+\delta )$, can be obtained either from the ratio, $\overline{J}/{\overline{J}}_{\mathrm{K}\mathrm{e}\mathrm{p}}$, or from the ratio, $\overline{\mathrm{\Omega }}/{\overline{\mathrm{\Omega }}}_{\mathrm{K}\mathrm{e}\mathrm{p}}$.

Digesting the EFE Results

The EFE table lists values of ${\overline{\mathrm{\Omega }}}_{\mathrm{K}\mathrm{e}\mathrm{p}}$, but not values of $r/a_1$. Inverting the expression just provided gives,

$(\frac{a_1}{r}{)}^3\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}$	$=$	$\frac{3}{8}{\overline{\mathrm{\Omega }}}_{\mathrm{K}\mathrm{e}\mathrm{p}}^2$
$\Rightarrow \frac{r}{a_1}$	$=$	$\left[\frac{8}{3{\overline{\mathrm{\Omega }}}_{\mathrm{K}\mathrm{e}\mathrm{p}}^2}\right(\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}){]}^{1/3}.$

In the above "EFE Check" column, we've listed these inferred values of $r/a_1$. Or, we might prefer the ratio, $r/R$, which is obtained from the Keplerian frequency via the expression,

${\mathrm{\Omega }}_K^2$	$=$	$\frac{4\pi G\rho }{3}\left(\frac{1+p}{p}\right)[\frac{R}{r}{]}^3$
$\Rightarrow [\frac{R}{r}{]}^3$	$=$	$\left[\frac{3}{4}\right]\left(\frac{p}{1+p}\right)\left(\frac{{\mathrm{\Omega }}_K^2}{\pi G\rho }\right)$
$\Rightarrow (\frac{r}{R}{)}_{p=1}$	$=$	$\left[\frac{3}{8}\right(\frac{{\mathrm{\Omega }}_K^2}{\pi G\rho }{)}_{p=1}{]}^{-1/3}.$

Comparison

In the context of our simplistic spherical model, above, we derived the following expression for the total angular momentum:

$L_{\mathrm{t}\mathrm{o}\mathrm{t}}$

$=$

$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}R{)}^{1/2}\{\frac{1}{(1+\lambda )}(\frac{d}{R}{)}^{1/2}+\frac{2}{5}[1+{\xcancel{\frac{1}{\lambda }(\frac{R^{\text{'}}}{R}{)}^2}}^0\left]\right(\frac{d}{R}{)}^{-3/2}\left\}\right(\frac{\lambda }{1+\lambda }).$

Rewriting our just-derived "Keplerian" expression to emphasize the ratio $r/R$ instead of $r/a_1$, and to highlight the system's total mass in the leading dimensional coefficient, allows us to more readily recognize the overlap with this simpler expression.

$J_{\mathrm{K}\mathrm{e}\mathrm{p}}$	$=$	$(G{M}^{3}R{)}^{1/2}[\frac{1}{(1+p)}(\frac{r}{R}{)}^{1/2}+\frac{1}{5}(1+\frac{a_2^2}{a_1^2}\left)\right(\frac{a_1}{R}\cdot \frac{R}{r}{)}^2\left(\frac{r}{R}{)}^{1/2}\right](\frac{1+p}{p}{)}^{1/2}$
	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}R{)}^{1/2}(\frac{p}{1+p}{)}^{3/2}\left[\frac{1}{(1+p)}\right(\frac{r}{R}{)}^{1/2}+\frac{1}{5}(1+\frac{a_2^2}{a_1^2})\left(\frac{a_1}{R}{)}^2\right(\frac{r}{R}{)}^{-3/2}\left]\right(\frac{1+p}{p}{)}^{1/2}$
	$=$	$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}R{)}^{1/2}[\frac{1}{(1+p)}(\frac{r}{R}{)}^{1/2}+\frac{1}{5}(1+\frac{a_2^2}{a_1^2}\left)\right(\frac{a_1}{R}{)}^2\left(\frac{r}{R}{)}^{-3/2}\right]\left(\frac{p}{1+p}\right)$

It makes sense, then, to write the total angular momentum as,

$L_{\mathrm{t}\mathrm{o}\mathrm{t}}$

$=$

$(G{M}_{\mathrm{t}\mathrm{o}\mathrm{t}}^{3}R{)}^{1/2}\{\frac{1}{(1+\lambda )}(\frac{d}{R}{)}^{1/2}+\frac{2}{5}\cdot \mathfrak{𝔍}(\frac{d}{R}{)}^{-3/2}\left\}\right(\frac{\lambda }{1+\lambda }),$

where, $\mathfrak{𝔍}=1$ when one assumes that the primary star is spherical, but when tidal distortions are taken into account,

$\mathfrak{𝔍}=\frac{1}{2}(1+\frac{a_2^2}{a_1^2})(\frac{a_1}{R}{)}^2.$

Figure 1: "Roche" Binary Sequences with Point-Mass Secondary and $M/M^{\text{'}}=1$
Our Constructed Diagram	Extracted from Fig. 10 of LRS93Supplement
Left:   Curves showing how the total system angular momentum varies with binary separation when $n=0$ and the secondary star $(M^{\text{'}})$ is treated as a point mass. (Blue dashed curve) Primary star assumed to be a sphere and, hence, $\mathfrak{𝔍}=1$; (Green filled circular markers) Primary star is an (EFE) ellipsoidal configuration with axis ratios specified by columns 2 and 3 of our Table 1, normalized angular momentum specified by column 6 of our Table 1, and binary separation specified by column 7 of our Table 1; (Solid red curve connecting red filled circular markers) Primary star is an (LRS93S) ellipsoidal configuration with axis ratios specified by columns 3 and 4 of our Table 2, normalized angular momentum specified by column 6 of our Table 2, and binary separation specified by column 2 of our Table 2. The green filled circular markers define the same (EFE) sequence that is presented as a dot-dashed curve in the right-hand panel; the red filled circular markers and associated smoothed curve define the same (LRS93S) sequence that is presented as a solid curve in the right-hand panel. The purple filled circlular marker identifies the turning point along the (LRS93S) sequence associated with the minimum system angular momentum; the yellow filled circular marker identifies the turning point along the same sequence that is associated with the minimum separation — the so-called "Roche" limit. Right:   (The following text is largely taken from the Fig. 10 caption of LRS93S) Equilibrium curves generated by LRS93S showing total angular momentum as a function of binary separation along two incompressible, and three compressible Roche sequences with $M/M^{\text{'}}=1$. The various curves display results from polytropic configurations having $n=0$ (solid line), $n=1$ (dotted line), $n=1.5$ (short-dashed line), and $n=2.5$ (long-dashed line). For comparison, the sequence obtained by EFE for $n=0$ is also drawn (dotted-dashed line).

Incompressible Roche Ellipsoids (λ = 0)

Extracted from p. 229 of G. H. Darwin (1906) IV.  On the figure and stability of a liquid satellite Philosophical Transactions of the Royal Society A, Vol. 206, pp. 161 - 248

"Finally the solution for Roche's limit and for the ratio of the axes of the ellipsoid in limiting stability may be taken to be as follows … $\gamma$ ${\sin }^{-1}\kappa$ $\cos \gamma$ $\cos \beta$ $r/a$ 61° 8½' 78° 52' 0.4827 0.5114 2.4553 with uncertainty of unity in the las place of decimals in $r$ and of half a minute of arc in ${\sin }^{-1}\kappa$.

Extracted from p. 242 of G. H. Darwin (1906) IV.  On the figure and stability of a liquid satellite Philosophical Transactions of the Royal Society A, Vol. 206, pp. 161 - 248

"As stated in the Preface, the radius vector of limiting stability, which has been called 'Roche's limit,' is found to be 2.4553, and the axes of the critical ellipsoid are proportional to the numbers 10000, 5114, 4827. These may be compared with the 2.44 and 1000, 496, 469 determined by Roche himself. When we consider the methods which he employed, we must be struck with the closeness to accuracy to which he attained."

Here we examine the results presented by Roche, by Darwin, and by EFE for the case of a point-mass secondary ($(M^{\text{'}})$ and a primary whose mass $(M)$ is formally zero. In this case, we must use a different scheme for normalizing physical quantities. Because the secondary is not spinning and it has no orbital motion, only the primary contributes to the system's "angular momentum"; but because the primary has no mass, we need to examine its (and, hence, the system's) specific angular momentum. Specifically,

$\frac{I}{M}$

$=$

$\frac{1}{5}{a}_1^2(1+\frac{a_2^2}{a_1^2}),$

and,

$j\equiv \frac{J_{\mathrm{t}\mathrm{o}\mathrm{t}}}{M}$	$=$	$[\frac{r^2}{(1+{\xcancel{p}}^0)}+\frac{I}{M}]{\mathrm{\Omega }}_{\mathrm{K}\mathrm{e}\mathrm{p}}$
	$=$	$R^2\left[\right(\frac{r}{R}{)}^2+\frac{2}{5}\cdot \mathfrak{𝔍}\left]\right(\frac{G{M}^{\text{'}}}{r^3}{)}^{1/2}$
	$=$	$(G{M}^{\text{'}}R{)}^{1/2}[(\frac{r}{R}{)}^{1/2}+\frac{2}{5}\cdot \mathfrak{𝔍}(\frac{r}{R}{)}^{-3/2}],$

where, in order to ensure that the density of the primary remains constant along an equilibrium sequece, the adopted normalizing length scale is customarily,

$R^3\equiv a_1{a}_2{a}_3\Rightarrow \frac{R}{a_1}=(\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}{)}^{1/3},$

in which case,

$\mathfrak{𝔍}$	$=$	$\frac{1}{2}(1+\frac{a_2^2}{a_1^2})(\frac{a_1}{R}{)}^2$
	$=$	$\frac{1}{2}(1+\frac{a_2^2}{a_1^2})(\frac{a_2}{a_1}\cdot \frac{a_3}{a_1}{)}^{-2/3}$

Table 3:  Incompressible $(n=0)$ Roche Ellipsoids with $\lambda =p=0$
Extracted from Table 1 of Chandrasekhar (1963) same as [EFE] Table XVI				EFE Check
(1)	(2)	(3)	(4)	(5)	(6)	(7)	(8)
${\cos }^{-1}(a_3/a_1)$	$a_2/a_1$	$a_3/a_1$	${\mathrm{\Omega }}^2$	$r/R$	$R/a_1$	$\mathfrak{𝔍}$	$j/(G{M}^{\text{'}}R{)}^{1/2}$
24°	0.93188	0.91355	0.022624	3.8916	0.9478	1.0400	2.0269
36°	0.84112	0.80902	0.047871	3.0312	0.8796	1.1035	1.8247
48°	0.70687	0.66913	0.074799	2.6122	0.7791	1.2352	1.7333
57°	0.57787	0.54464	0.088267	2.4720	0.68022	1.4415	1.7206
60°	0.53013	0.50000	0.089946	2.4565	0.6424	1.5523	1.7286
61°	0.51373	0.48481	0.090068	2.4554	0.6292	1.5964	1.7329
62°	0.49714	0.46947	0.089977	2.4562	0.6157	1.6450	1.7382
63°	0.48040	0.45399	0.089689	2.4589	0.6019	1.6984	1.7443
66°	0.42898	0.40674	0.087201	2.48202	0.5588	1.8959	1.7694
71°	0.34052	0.32557	0.077474	2.5818	0.4804	2.4178	1.8399
72°	0.32254	0.30902	0.074648	2.6140	0.4636	2.5679	1.8598
75°	0.26827	0.25882	0.064426	2.7455	0.4110	3.1728	1.9359
79°	0.19569	0.19081	0.047111	3.0475	0.3342	4.6471	2.0951

References

• 📚 G. H. Darwin (1906, Philosophical Transactions of the Royal Society A, Vol. 206, pp. 402 - 412), On the Figure and Stability of a Liquid Satellite.
• J. H. Jeans (1919) Problems of Cosmogony and Stellar Dynamics -- see especially beginning in Ch. III (p. 46)
• Z. Kopal (1959), Close Binary Systems -- see especially, Ch. III (p. 125); but not particularly useful.
• Chapter 8 (pp. 189-241) of [EFE]
• Chapter 16 (pp. 448-469) of [T78]
• 📚 D. Lai, F. A. Rasio, & S. L. Shapiro (1993a, ApJLetters, Vol. 406, pp. L63 - L66), Hydrodynamic Instability and Coalescence of Close Binary Systems.
• 📚 D. Lai, F. A. Rasio, & S. L. Shapiro (1993b, ApJ Suppl., Vol. 88, pp. 205 - 252), Ellipsoidal Figures of Equilibrium: Compressible Models.
• Other LRS papers:
  • LRS (1994, ApJ, Vol. 420, p. 811 — Hydrodynamic Instability and Coalescence of Binary Neutron Stars
  • LRS (1994, ApJ, Vol. 423, p. 344 — Equilibrium, Stability, and Orbital Evolution of Close Binary Systems
  • LRS (1994, ApJ, Vol. 437, p. 742 — Hydrodynamics of Rotating Stars and Close Binary Interactions: Compressible Ellipsoid Models
• Y. Eriguchi & K. Uryu (1999, Prog. Theor. Phys. Suppl., 136, 199-215) Darwin-Riemann Problems in Newtonian Gravity; see also ADS reference

Reproduction of Fig. 1 (p. L64) from D. Lai, F. A. Rasio, & S. L. Shapiro (1993) Hydrodynamic Instability and Coalescence of Close Binary Systems Astrophysical Journal Letters, Vol. 406, pp. L63 - L66	Extracted from Fig. 1 (p. 315) of K. C. B. New & J. E. Tohline (1997) The Relative Stability Against Merger of Close, Compact Binaries The Astrophysical Journal, Vol. 490, pp. 311 - 327
Lai, Rasio & Shapiro (1993a)	New & Tohline (1997)

See Also

The discussion presented here is supported by detailed reviews and new derivations presented in the following associated chapters:

• A free-energy synopsis.
• Material associated with an overarching PowerPoint presentation.
• Limiting Masses
• An accompanying chapter on the Bonnor-Ebert sphere.

Overlapping discussions of this topic may also be found in the following key references:

• §6.8 of [ST83].

Other publications that are closely aligned with the analysis found in 📚 G. H. Darwin (1906, Philosophical Transactions of the Royal Society A, Vol. 206, pp. 402 - 412):

• G. H. Darwin (1901, Philosophical Transactions of the Royal Society A, Vol. 197, pp. 461 - 557), Ellipsoidal Harmonic Analysis.
• G. H. Darwin (1902, Philosophical Transactions of the Royal Society A, Vol. 198, pp. 301 - 331), The Pear-Shaped Figure of Equilibrium of a Rotating Mass of Liquid.
• G. H. Darwin (1903, Philosophical Transactions of the Royal Society A, Vol. 200, pp. 251 - 314), The Stability of the Pear-Shaped Figure of Equilibrium, &c.
• G. H. Darwin (1904, Philosophical Transactions of the Royal Society A, Vol. 203, pp. 111 - 137), Integrals of the Squares of Ellipsoidal Surface Harmonic Functions.
• G. H. Darwin (1908, Philosophical Transactions of the Royal Society A, Vol. 208, pp. 1 - 19), Further Consideration of the Stability of the Pear-Shaped Figure of a Rotating Mass of Liquid.   $\leftarrow$   Darwin points out that Liapounoff's (1905) work concludes that the "Pear-Shaped" figure is unstable, whereas in his own 1903 work, Darwin "… had arrived at an opposite conclusion." In this 1908 publication, Darwin has changed his conclusion and is now in agreement with Liapounoff.
• D. M. Christodoulou, D. Kazanas, I. Shlosman, and J. E. Tohline (1995, ApJ, Vol. 446, p. 500 - 509), Phase-Transition Theory of Instabilities.   III.   The Third-Harmonic Bifurcation on the Jacobi Sequence and the Fission Problem.