User:Tohline/SSC/Structure/BonnorEbert
From VisTrailsWiki
Contents 
PressureBounded Isothermal Sphere
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
As has been derived and discussed elsewhere, an isolated isothermal sphere has a density profile that extends to infinity and, correspondingly, an unbounded total mass. In an astrophysical context, neither of these properties is desirable. A more realistic isothermal configuration can be constructed by embedding the structure in a low density, but hot external medium whose pressure, , confines the isothermal configuration to a finite size. In a mathematical model, this can be accomplished by ripping off an outer layer of the isolated isothermal configuration down to the radius — label it — at which the configuration's original (internal) pressure equals ; the interior of the configuration that remains — containing mass — should be unaltered and in equilibrium. (This will work only for spherically symmetric configurations, as the gravitational acceleration at any location only depends on the mass contained inside that radius.) Ebert (1955) and Bonnor (1956) are credited with constructing the first such models and, most significantly, discovering that, for any specified sound speed and applied external pressure, there is a mass above which no equilibrium configuration exists. We present, here, the salient elements of these (essentially equivalent) derivations.
Prior to studying this discussion of pressurebounded isothermal spheres, we recommend studying our related discussion of pressurebounded = 5 polytropes. As with isolated isothermal spheres, isolated = 5 polytropes extend to infinity. But, unlike their isothermal counterparts, the structure of = 5 polytropes is describable analytically. Hence, an analysis of their structure and its extension to pressurebounded configurations avoids the clutter introduced by a model — such as the isothermal sphere — that can only be described numerically. As it turns out, the pressurebounded = 5 polytrope exhibits a BonnorEbert type limiting mass that is analytically prescribable. Its derivation is mathematically quite clean and provides a firm foundation for understanding the better known — but only numerically prescribable — BonnorEbert limiting mass.
Governing Relation
The equilibrium structure of an isolated isothermal sphere, as derived by Emden (1907), has been discussed elsewhere. From this separate discussion we appreciate that the governing ODE is,
where,
is the square of the isothermal sound speed. In their studies of pressurebounded isothermal spheres, Ebert (1955, ZA, 37, 217) and Bonnor (1956, MNRAS, 116, 351) both started with this governing ODE, but developed its solution in different ways. Here we present both developments while highlighting transformations between the two.
Derivation by Bonnor (1956) (edited)  translation  Derivation by Ebert (1955) (edited)  
Let us now transform (2.3) by making the following substitutions:
where is an arbitrary constant [we choose ], and
Then (2.3) becomes

Wir setzen unter Verwendung von (4) und mit . Für ergibt sich die Differentialgleichung der isothermen Gaskugel:
Der Strich bezeichnet die Differentiation nach . 

Both of these dimensionless governing ODEs — Bonnor's Eq. (2.8) and Ebert's Eq. (17) — are identical to the dimensionless expression derived by Emden (see the presentation elsewhere), namely,
The translation from EmdentoBonnortoEbert is straightforward:
In much of what follows, we will use Bonnor's notation rather than Emden's notation. This means that we will be referring to the isothermal LaneEmden function, , which provides a solution to the governing,
PV Diagram
Given a value for the isothermal sound speed, , and the central density, , our accompanying summary of the properties of an isothermal sphere, as derived by Emden (1907), provides the following structural relations:
 Radial Coordinate Position & Associated Volume:
;
;
 Pressure Profile:
;
 Enclosed Mass:
Bonnor and Ebert both asked how the surface pressure of a truncated isothermal sphere of a given sound speed — or, equivalently, how the external pressure, , required to confine that truncated configuration — will vary with volume if the mass of the configuration, rather than the central density, is held fixed. The above expression for the enclosed mass can be inverted to give in terms of ; specifically,
where the subscript denotes a truncated configuration of radius confined by an external pressure. This allows us to replace in favor of in the expressions for both the pressure and volume, obtaining,
 Pressure:
 or, adopting Ebert's (1955) normalization pressure, — see his equation (22),
 Volume:
 or, adopting a normalization volume, , defined in terms of Ebert's (1955) normalization radius, — again, see his equation (22),
 the corresponding normalized radius of the truncated configuration is,
The expressions shown here for and are identical to those presented by Bonnor (1956) for "" and "," immediately following his equation (2.17). They are also the set of parametric equations that lead to Bonnor's Fig. 1 PV diagram, which is reproduced here as the lefthand panel of our Figure 1 (with permission, as our caption to this figure documents).
Figure 1: Bonnor's PV Diagram 

Pertaining to Reprinted Figure (immediately above)

Pertaining to Our Original Figure (immediately above)
Caption: In an effort to demonstrate more explicitly how Bonnor's 1956 work ties back to the original work on isolated isothermal spheres that was published by Emden in 1907, the blue curve in the above panel displays a plot of versus that we have constructed using data drawn directly from Table 14 of Emden (1907). 
Ebert chose to illustrate the behavior of the sequence of pressurebounded isothermal spheres by plotting versus , rather than versus . The result — Figure 2 in Ebert (1955) — is reproduced here as the lefthand panel of our Figure 2. It should be clear that the information contained in Ebert's plot is identical to the information contained in the curve that is displayed in Bonner's PV diagram.
Figure 2: Ebert's PR Diagrams 

The reader may readily view this referenced image/text by opening (preferably in a separate browser window) the relevant article and scrolling to the relevant figure/page, as detailed immediately below.


Pertaining to Reprinted Figure (immediately above)

Pertaining to Our Original Figure (immediately above)
Caption: In an effort to demonstrate more explicitly how Ebert's 1955 work ties back to the original work on isolated isothermal spheres that was published by Emden in 1907, the blue curve in the above panel displays a plot of versus that we have constructed using data drawn directly from Table 14 of Emden (1907). 
In an effort to demonstrate more explicitly how Bonnor's 1956 work and Ebert's 1955 work tie back to the original work on isolated isothermal spheres that was published by Emden in 1907, the righthand panel of our Figure 1 shows a plot of versus , and the righthand panel of our Figure 2 shows a plot of versus that we have constructed using data drawn directly from Table 14 of Emden (1907).
Whitworth's Normalization
Whitworth (1981, MNRAS, 195, 967) examined the behavior of model sequences that result from embedding polytropes with various effective adiabatic indexes in an external medium of pressure . Following Ebert's lead, in his Figure 1b, Whitworth showed how the size of the configuration varies with along each sequence. A reproduction of Whitworth's Figure 1b has been displayed in our separate discussion of pressurebounded polytropic configurations that have analytic prescriptions, but we have again displayed the figure here, in the rightmost panel of our Figure 2, because Whitworth included pressurebounded isothermal configurations among his models. The red arrow identifies the isothermal sequence. Whitworth's Figure 1b effectively serves as an extension of Ebert's Figure 2, shown above in the leftmost panel of our Figure 2. In order to assist comparison with Ebert's figure, this time we have rotated and flipped Whitworth's figure relative to its original orientation because Whitworth chose to plot versus whereas Ebert plotted versus .
Figure 3: Whitworth's PR Diagrams 

Pertaining to Reprinted Figure


It should be noted that the scaling adopted along both axes in Whitworth's Figure 1b is different from the scaling used by Ebert. As is detailed in the accompanying ASIDE, Whitworth "referenced" P_{e} and R to, respectively,
In generating the righthand panel of our Figure 1, we have adopted Whitworth's scaling. That is to say, we have plotted versus , where,



and,



Limiting Pressure and Maximum Mass
No matter how you look at the bounded isothermal sphere sequence — whether plotted as a curve in Bonnor's PV diagram or as a curve in Whitworth's RP diagram — it is clear that, for a given choice of the sound speed and the mass, there is a value of the pressure above which no equilibrium configuration exists. (Alternatively, for a given sound speed and external pressure, there is a limiting mass above which no equilibrium configuration exists; see below.) The configuration identifying this limiting pressure resides at the position along Bonnor's PV diagram sequence where — or, in Whitworth's discussion, — first goes to zero.
Following Bonnor's Presentation
As is shown in his equation (3.1), Bonnor (1956) determined that, at a fixed mass and sound speed, goes to zero when,
In the following figure, the function defined by the lefthandside of this expression is plotted versus using Emden's (1907) tabulated data. As Bonnor noted, the function first crosses zero when .
An isothermal sphere that is truncated at this location will have a radius,
which matches equation (3.3) of Bonnor (1956); and, drawing on function values provided in Emden's (1907) Table 14, it will have a total mass,
Dividing this expression for by the expression for gives,
which matches equation (3.6) of Bonnor (1956). In order to maintain an equilibrium structure while truncating the isothermal model at this radius requires applying an external pressure of the following magnitude:
where, again, numerical values have been drawn from Emden's (1907) Table 14.
Finally, using this relation to eliminate from the expression for gives,

= 






This is the mass upper limit for a stable, pressurebounded isothermal sphere — the socalled BonnorEbert mass; see, for example, equation (1) of Shu (1977). In a separate discussion, we compare this result to the determinations of other related mass upperlimits.
Ebert's Corresponding Presentation
The expression derived by Bonnor for the ratio in the limiting configuration can be inverted to give,
This matches the expression for the critical radius, , that appears as equation (23) in Ebert's (1955) published derivation. Similarly, the above relation that expresses the BonnorEbert limiting mass in terms of , , and can be inverted to give,
which agrees with the limiting pressure that was derived by Ebert (1955) and that is also presented in his equation (23).
Mass versus Radius Equilibrium Sequence
Old Derivation
For use in a related discussion, let's examine how the mass varies with radius — at fixed and fixed — in the limit as . From above, we have,






and,






Hence,












Also,









So the derivative of interest is,














Now, we can also write,



and,



Hence,





New Derivation
For use in a related discussion, let's examine how the mass, , varies with radius, , at fixed and . From above, we have,






and,
;
Given that, , we can also write,






Hence,









Also,



So the derivative of interest is,











Now, we can also write,



and,



Hence,





Related Discussions
 Polytropes emdeded in an external medium
 Constructing BiPolytropes
 Analytic description of BiPolytrope with (n_{c},n_{e}) = (5,1)
 BonnorEbert spheres
 BonnorEbert Mass according to Wikipedia
 A MATLAB script to determine the BonnorEbert Mass coefficient developed by CheYu Chen as a graduate student in the University of Maryland Department of Astronomy
 SchönbergChandrasekhar limiting mass
 Relationship between BonnorEbert and SchönbergChandrasekhar limiting masses
© 2014  2019 by Joel E. Tohline 