User:Tohline/SSC/Stability/BiPolytropes
From VisTrailsWiki
Marginally Unstable (n_{c},n_{e}) = (5,1) Bipolytropes
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
Overview
Figure 1: Equilibrium Sequences of PressureTruncated Polytropes 

We expect the content of this chapter — which examines the relative stability of bipolytropes — to parallel in many ways the content of an accompanying chapter in which we have successfully analyzed the relative stability of pressuretruncated polytopes. Figure 1, shown here on the right, has been copied from a closely related discussion. The curves show the massradius relationship for pressuretruncated model sequences having a variety of polytropic indexes, as labeled, over the range . (Another version of this figure includes the isothermal sequence.) On each sequence for which , the green filled circle identifies the model with the largest mass. We have shown analytically that the oscillation frequency of the fundamentalmode of radial oscillation is precisely zero^{†} for each one of these maximummass models. As a consequence, we know that each green circular marker identifies the point along its associated sequence that separates dynamically stable (larger radii) from dynamically unstable (smaller radii) models.
^{†}In each case, the fundamentalmode oscillation frequency is precisely zero if, and only if, the adiabatic index governing expansions/contractions is related to the underlying structural polytropic index via the relation, , and if a constant surfacepressure boundary condition is imposed.
In another accompanying chapter, we have used purely analytic techniques to construct equilibrium sequences of spherically symmetric bipolytropes that have, . For a given choice of — the ratio of the meanmolecular weight of envelope material to the meanmolecular weight of material in the core — a physically relevant sequence of models can be constructed by steadily increasing the value of the dimensionless radius at the core/envelope interface, , from zero to infinity. Figure 2, whose content is essentially the same as Figure 1 of this separate chapter, shows how the fractional core mass, , varies with the fractional core radius, , along sequences having seven different values of , as labeled: 1 (black), ½ (dark blue), 0.345 (brown), ⅓ (dark green), 0.316943 (purple), 0.309 (orange), and ¼ (light blue).
When modeling bipolytropes, the default expectation is that an increase in along a given sequence will correspond to an increase in the relative size — both the radius and the mass — of the core. This expectation is realized along the Figure 2 sequences that have the largest meanmolecular weight ratios: = 1 and ½. But the behavior is different along the other five illustrated sequences. For sufficiently large , the relative radius of the core begins to decrease; along each sequence, a solid purple circular marker identifies the location of this turning point in radius. Furthermore, along sequences for which , eventually the fractional mass of the core reaches a maximum and, thereafter, decreases even as the value of continues to increase; a solid green circular marker identifies the location of this maximum mass turning point along each of these sequences. (Additional properties of these equilibrium sequences are discussed in yet another accompanying chapter.)
The principal question is: Along bipolytropic sequences, are maximummass models associated with the onset of dynamical instabilities?
Planned Approach
Figure 2: Equilibrium Sequences of Bipolytropes with and Various 


Ideally we would like to answer the juststated "principal question" using purely analytic techniques. But, to date, we have been unable to fully address the relevant issues analytically, even in what would be expected to be the simplest case: bipolytropic models that have . Instead, we will streamline the investigation a bit and proceed — at least initially — using a blend of techniques. We will investigate the relative stability of bipolytropic models having whose equilibrium structures are completely defined analytically; then the eigenvectors describing radial modes of oscillation will be determined, one at a time, by solving the relevant LAWE(s) numerically. We are optimistic that this can be successfully accomplished because we have had experience numerically integrating the LAWE that governs the oscillation of:
 Isolated n = 3 polytropes — including a quantitative comparison against Schwarzschild's (1941) published work;
 Pressuretruncated isothermal spheres — including a quantitative comparison against the published analysis of Taff & Van Horn (1974); and
 Pressuretruncated n = 5 polytropes.
A key reference throughout this investigation will be the paper by J. O. Murphy & R. Fiedler (1985b, Proc. Astr. Soc. of Australia, 6, 222). They studied Radial Pulsations and Vibrational Stability of a Sequence of Two Zone Polytropic Stellar Models. Specifically, their underlying equilibrium models were bipolytropes that have . In an accompanying chapter, we describe in detail how Murphy & Fiedler obtained these equilibrium bipolytropic structures and detail some of their equilibrium properties.
Here are the steps we initially plan to take:
 Governing LAWEs:
 Identify the relevant LAWEs that govern the behavior of radial oscillations in the core and, separately, in the envelope. Check these LAWE specifications against the published work of Murphy & Fiedler (1985b).
 Determine the matching conditions that must be satisfied across the core/envelope interface. Be sure to take into account the critical interface jump conditions spelled out by P. Ledoux & Th. Walraven (1958), as we have already discussed in the context of an analysis of radial oscillations in zerozero bipolytropes.
 Determine what surface boundary condition should be imposed on physically relevant LAWE solutions, i.e., on the physically relevant radialoscillation eigenvectors.
 Initial Analysis:
 Choose a maximummass model along the bipolytropic sequence that has, for example, . Hopefully, we will be able to identify precisely (analytically) where this maximummass model lies along the sequence. Yes! Our earlier analysis does provide an analytic prescription of the model that sits at the maximummass location along the chosen sequence.
 Solve the relevant eigenvalue problem for this specific model, initially for and initially for the fundamental mode of oscillation.
Virial Analysis
In an accompanying chapter we have examined the viral stability of bipolytropes — see also an associated discussion of the freeenergy of these configurations. We present, here, an overview of this separate discussion in an effort to provide a broader perspective to the focus of this current chapter.
Various Energy Expressions
Drawing from our accompanying detailed discussion, we have,






Recognizing that,












we have,






Next, given that,

and 

we have,






Also, given that,



we have,






Finally, then, we can write,







Free Energy and Its Derivatives
We can now rewrite the freeenergy expression in the form,



where,












The first derivative is,



And the second derivative is,



What to Expect for Equilibrium Configurations
In equilibrium we should set and . In which case we should expect to find that,



In addition, an evaluation of the second derivative should give,



Then, the transition from stable to unstable configurations occurs when , that is, when,









For example, if we set and , this implies,



FreeEnergy Stability Evaluation
Here we pull together excerpts from several different H_Book Chapters in which we have presented, from several different perspectives, an analysis of the freeenergy of bipolytropes.
 In one chapter, using purely analytic techniques, we have derived expressions that detail the structural properties of bipolytropes having . Among these are analytic expressions for various terms that make up the freeenergy expression: , , , , and . Equilibrium model sequences are defined by fixing the ratio, , then varying the radial location, , of the coreenvelope interface; note that the volume of the core is, then, .
 [Virial Equilibrium] In a subsection of this same chapter titled, Equilibrium Condition: Global, we have shown that a statement of virial equilibrium — obtained by setting the first derivative of the freeenergy expression to zero — is,
and
 [Marginally Unstable Model] Near the bottom of this same chapter, in a subsection titled, Stability Condition, we point out that the model along each sequence that is marginally (dynamically) unstable — obtained setting the second derivative of the freeenergy expression to zero — is identified by the configuration for which,
Table XXX: Properties of Marginally Unstable Bipolytropes Having 


CASE1  CASE2  CASE3  CASE4  
1  2.41610822  2.8049  0.59520261  0.68306067  16.3788  0.039116848  4.446748782   6.606135366  0.889349762   4.066061722  2.287362198  1  0  5.0000000  2.1 × 10^{8}  
4.1853093  8.8058  0.328419479  0.70131896  354.089  0.003126324  5.76978580   10.58931853  1.153956968   3.258165567  0.95025163  1  6.3 × 10^{8}  5.0000002  0  
0.345  7.64325  44.116  0.119714454  0.52700045  2.85 × 10^{4}  0.000116533  6.230343527   12.24495934  1.1246068658   2.707865028  0.215727713  1  6.4 × 10^{8}  5.0000002  0  
8.548103  59.643  0.099032423  0.47901529  6.30 × 10^{4}  6.1337 × 10^{5}  6.261548334   12.36425897  1.252309682   2.663457063  0.158837699  1  0  4.9999999  6.0 × 10^{8}  
0.316943  10.7441565  108.14  0.068655205  0.38238387  2.93 × 10^{5}  1.6252 × 10^{5}  6.301810768   12.52005323  1.260362204   2.604292714  0.083568307  1  0  4.9999998  2.0 × 10^{7}  
0.309  12.77156  166.06  0.053145011  0.31696879  8.70 × 10^{5}  5.8905 × 10^{6}  6.318902171   12.58692884  1.26378042   2.57843634  0.050875500  1  1.9 × 10^{8}  5.0000001  0  

The lefthand panel of Figure 3 is identical to Figure 2, above. It displays in the parameter space, the behavior of bipolytropic equilibrium sequences that have, as labeled, seven different values of the ratio of meanmolecularweights, . Using a numerical rootfinding technique, we have determined where the virial stability condition, , is satisfied along each of these sequences — as well as along a number of additional equilibrium sequences. Key properties of each of these identified models have been recorded in Table 1 of an accompanying discussion; see also an associated discussion of the freeenergy of these configurations. Pulling from this tabulated data, the solidred circular markers that appear in the righthand panel of Figure 3 identify where this virial stability condition is satisfied along the separate equilibrium sequences while the accompanying dashed red curve identifies more broadly how the parameter space is divided into stable (below and to the right) versus unstable (above and to the left) regions.
Figure 3 

In what follows we use a complementary — and more quantitatively rigorous — approach to evaluating the stability of equilibrium models, and contrast the results of that analysis with the virialanalysis results presented graphically here in Figure 3.
Review of the Analysis by Murphy & Fiedler (1985b)
As we have detailed separately, the boundary condition at the center of a polytropic configuration is,
and the boundary condition at the surface of an isolated polytropic configuration is,


at 
But this surface condition is not applicable to bipolytropes. Instead, let's return to the original, more general expression of the surface boundary condition:



Utilizing an accompanying discussion, let's examine the frequency normalization used by Murphy & Fiedler (1985b) (see the top of the lefthand column on p. 223):
For a given radial quantum number, , the factor inside the square brackets in this last expression is what Murphy & Fiedler (1985b) refer to as . Keep in mind, as well, that, in the notation we are using,
This also means that the surface boundary condition may be rewritten as,

Let's apply these relations to the core and envelope, separately.
Interface Conditions
Here, we will simply copy the discussion already provided in the context of our attempt to analyze the stability of bipolytropes; specifically, we will draw from STEP 4: in the Piecing Together subsection. Following the discussion in §§57 & 58 of P. Ledoux & Th. Walraven (1958), the proper treatment is to ensure that fractional perturbation in the gas pressure (see their equation 57.31),



is continuous across the interface. That is to say, at the interface , we need to enforce the relation,









In the context of this interfacematching constraint (see their equation 62.1), P. Ledoux & Th. Walraven (1958) state the following: In the static (i.e., unperturbed equilibrium) model … discontinuities in or in might occur at some [radius]. In the first case — that is, a discontinuity only in density, while — the interface conditions imply the continuity of at that [radius]. In the second case — that is, a discontinuity in the adiabatic exponent — the dynamical condition may be written as above. This implies a discontinuity of the first derivative at any discontinuity of .
The algorithm that Murphy & Fiedler (1985b) used to "… [integrate] through each zone …" was designed "… with continuity in and being imposed at the interface …" Given that they set , their interface matching condition is consistent with the one prescribed by P. Ledoux & Th. Walraven (1958).
Radial Oscillations of (n_{c}, n_{e}) = (5, 1) Models
Foundation
In an accompanying discussion, we derived the socalled,
whose solution gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, selfgravitating fluid configurations. Assuming that the underlying equilibrium structure is that of a bipolytrope having , it makes sense to adopt the normalizations used when defining the equilibrium structure, namely,



; 






; 



We note as well that,









Hence, multiplying the LAWE through by gives,


















Profile
Now, referencing the derived bipolytropic model profile, we should incorporate the following relations:
Table 1: Radial Profile of Various Physical Variables
Variable 
Throughout the Core 
Throughout the Envelope^{†} 
Plotted Profiles 






















^{†}In order to obtain the various envelope profiles, it is necessary to evaluate φ(η) and its first derivative using the information presented in Step 6, above. 
Therefore, throughout the core we have,






and, throughout the envelope we have,






NOTE on 15 May 2019: Prior to this date the last RHS expression had an incorrect exponent on . It previously (incorrectly) read,

Numerical Integration
General Approach
Here, we begin by recognizing that the 2^{nd}order ODE that must be integrated to obtain the desired eigenvectors has the generic form,



where,



and 



Adopting the same approach as before when we integrated the LAWE for pressuretruncated polytropes, we will enlist the finitedifference approximations,



and 



The finitedifference representation of the LAWE is, therefore,









In what follows we will also find it useful to rewrite in the form,
The relevant coefficient expressions for all regions of the configuration are,



, 



and 



Special Handling at the Center
In order to kickstart the integration, we set the displacement function value to at the center of the configuration , then draw on the derived powerseries expression to determine the value of the displacement function at the first radial grid line, , away from the center. Specifically, we set,



Special Handling at the Interface
Integrating outward from the center, the general approach will work up through the determination of when "j+1" refers to the interface location. In order to properly transition from the core to the envelope, we need to determine the value of the slope at this interface location. Let's do this by setting j = i, then projecting forward to what would be — that is, to what the amplitude just beyond the interface would be — if the core were to be extended one more zone. Then, the slope at the interface (as viewed from the perspective of the core) will be,












Conversely, as viewed from the envelope, if we assume that we know and , we can determine the amplitude, , at the first zone beyond the interface as follows:












Splitting Analysis Into Separate Core and Envelope Components
Core:
Given that, , lets multiply the LAWE through by . This gives,



Specifically for the core, therefore, the finitedifference representation of the LAWE is,









This also means that, as viewed from the perspective of the core, the slope at the interface is



Envelope:
Given that,
let's multiply the LAWE through by . This gives,



Specifically for the envelope, therefore, the finitedifference representation of the LAWE is,









This also means that, once we know the slope at the interface (see immediately below), the amplitude at the first zone outside of the interface will be given by the expression,



Interface
If we consider only cases where , then at the interface we expect,









Switching at the interface from to therefore means that,



If, however, we want to consider values for the adiabatic index that are different in the two regions, we have to follow the aboveoutlined guidelines, that is,









Eigenvectors for Marginally Unstable Models with (γ_{c}, γ_{e}) = (6/5, 2)
We now have the tools in hand to identify the eigenvectors — that is, various radial eigenfunctions and the corresponding eigenfrequency for each — associated with various modes of oscillation in bipolytropes. Which models should we examine?
In our accompanying review of the bipolytrope stability analysis presented by Murphy & Fiedler (1983b), our primary objective was to show that we were able to match their results quantitatively. We therefore set = 1 — the only ratio that they considered — and picked values of the coreenvelope interface radius, , that were listed among their set of chosen models. For a fixed value of , we integrated the relevant LAWE from the center toward the surface for many different eigenfrequency guesses until an eigenfunction was found whose behavior at the surface matched with high precision the physically justified surface boundary condition.
With the above virial stability analysis in mind (see especially Figure 3), here we have chosen to focus on models that reside along six of the equilibrium sequences that have already been analytically identified, above — specifically, the sequences for which = 1, Â½, 0.345, â…“, 0.309, and Â¼ — and to examine oscillation modes under the assumption that,
and
Fundamental Modes
We decided to examine, first, whether any model along each sequence marks a transition from dynamically stable to dynamically unstable configurations. We accomplished this by setting = 0, then integrating the relevant LAWE from the center toward the surface for many different guesses of the coreenvelope interface radius until an eigenfunction with no radial nodes — i.e., an eigenfunction associated with the fundamental mode of radial oscillation — was found whose behavior at the surface matched with high precision the physically desired surface boundary condition. We were successful in this endeavor. A marginally unstable model was identified on each of the six separate equilibrium sequences.
Equilibrium Properties of Marginally Unstable Models
Table 2 summarizes some of the equilibrium properties of these six models. For example, the second column of the table gives the value of the coreenvelope interface radius, , associated with each marginally unstable model. The table also lists: the value of the model's dimensionless radius, , the key structural parameters, & , and the centraltomean density associated with each model; and in each case the dimensionless thermal energy and dimensionless gravitational potential energy associated, separately, with the core and the envelope. Note that, once the pair of parameters, , has been specified, we can legitimately assign highprecision values to all of the other model parameters because they are analytically prescribed.
Table 2: Properties of Marginally Unstable Bipolytropes Having 


1  1.6686460157  2.139737  0.53885819  0.497747626  8.51704656  3.021916335  3.356583022  1.47780476  5.642859167 
2.27925811317  5.146499  0.306021732  0.401776274  63.29514949  4.241287117  6.074241035  4.284931508  10.97819621  
0.345  2.560146865247  9.554041  0.185160563  0.234302525  209.7739052  4.639705843  7.125754184  11.72861751  25.61089252 
2.582007485476  10.120558  0.176288391  0.218241608  230.4125398  4.667042505  7.200966267  13.15887139  28.45086152  
0.309  2.6274239687695  11.464303  0.158362807  0.184796947  279.0788798  4.722277318  7.354156963  17.1374434  36.36528446 
2.7357711469398  15.895632  0.118924863  0.11071211  430.0444648  4.84592201  7.70305421  37.84289623  77.67458196 
As was expected from our above discussion of virial equilibrium conditions, we found that to high precision for each of these equilibrium models,



However, contrary to expectations, in no case did we find that . That is to say, we found that none of the models lies on the (reddashed) curve in the parameter space that separates stable from unstable models as defined by our above freeenergybased stability analysis. The lefthand panel of Figure 4 shows this (reddashed) demarcation curve; for all intents and purposes, it is a reproduction of the righthand panel of Figure 3, above — turningpoint markers have been removed to minimize clutter, the equilibrium sequences have been labeled, and the horizontal axis has been extended to unity in order to include a longer portion of the sequence. The orange triangular markers that appear in the righthand panel of Figure 4 pinpoint where each of the Table 2 "marginally unstable" models resides in this plane. Clearly, all six of the orange triangles lie well off of — and to the stable side of — the reddashed demarcation curve. This discrepancy, which has resulted from our use of two separate approaches to stability analysis, will be discussed further and gratifyingly resolved, below.
Figure 4  

Eigenfunction Details
Here we examine some of the properties of the fundamentalmode eigenfunctions that we have found are associated with marginally unstable, bipolytropes.
Figure 5 

Consider the model on the sequence for which ; key properties of this specific equilibrium model are enumerated in the first row of numbers provided in Table 2, above. Figure 5 shows how our numerically derived, fundamentalmode eigenfunction, , varies with the fractional radius over the entire range, . By prescription, the eigenfunction has a value of unity and a slope of zero at the center . Integrating the LAWE outward from the center, through the model's core (blue curve segment), drops smoothly to the value at the interface . Our numerical integration of the LAWE showed that, at the interface, the logarithmic slope of the core (blue) segment of the eigenfunction is,



Next, following the above discussion of matching conditions at the interface, we determined that, from the perspective of the envelope, the slope of the eigenfunction at the interface must therefore be,



Adopting this "env" slope along with the amplitude, , as the appropriate interface boundary conditions, we integrated the LAWE from the interface to the surface, obtaining the greencolored segment of the eigenfunction that is shown in Figure 5. The amplitude continued to steadily decrease, reaching a value of , at the model's surface . At the surface, this envelope (green) segment of the eigenfunction exhibits a logarithmic slope that matches to eight significant digits the value that is expected from astrophysical arguments for this marginally unstable model, namely,
Key Reminder: We were able to find an eigenfunction whose surface boundary condition matched the desired value — in this particular case, a logarithmic slope of negative one — to this high level of precision only by iterating many times and, at each step, finetuning our choice of the equilibrium model's radial interface location, before performing a numerical integration of the LAWE. 
The discontinuous jump that occurs in the slope of the eigenfunction at the interface results from our assumption that the effective adiabatic index of material in the core is different from the effective adiabatic index of the envelope material . In an effort to emphasize and more clearly illustrate the behavior of this fundamentalmode eigenfunction as it crosses the core/envelope interface, we have added a pair of dashed line segments to the Figure 5 plot. The reddashed line segment touches, and is tangent to, the blue segment of the eigenfunction at the location of the core/envelope interface; it has a slope,



On the other hand, the purpledashed line segment touches, and is tangent to, the green segment of the eigenfunction at the location of the core/envelope interface; it has a slope,



For comparison purposes, the eigenfunction shown in Figure 5 has been presented again in Figure 6, along with several other of our numerically derived eigenfunctions, but in Figure 6 the plotted amplitude has been renormalized to give a surface value — rather than a central value — of unity.
In Figure 6 we show the behavior of the fundamentalmode eigenfunction for each of the marginally unstable models identified in Table 2. In the top figure panel, each curve shows — on a linearlinear plot — how the amplitude varies with radius; in the bottom figure panel, the amplitude is plotted on a logarithmic scale. On each curve, the black plus sign marks the radial location of the coreenvelope interface; in the bottom panel, these markers are accompanied by the values of that are associated with each corresponding model (see also the second column of Table 2). Each eigenfunction has been normalized such that the surface amplitude is unity. In the top panel, the value of the central amplitude of the eigenfunction that results from this normalization is recorded near the point where each eigenfunction touches the vertical axis. (In each case, the value provided on the plot is simply the inverse of the value of given in Table 3, below.)
Figure 6: Eigenfunctions Associated with the FundamentalMode of Radial Oscillation 

Notice that, especially as they approach the surface, the "envelope" segments of these six marginally unstable eigenfunction appear to merge into the same curve, irrespective of their value of the ratio of mean molecular weights. Note as well that the discontinuous jump that occurs in the slope of each eigenfunction at the radial location of the core/envelope interface — resulting from our choice to adopt a different adiabatic index, , in the core from the one in the envelope — becomes less and less noticeable for smaller and smaller values of the ratio of mean molecular weights.
Is There an Analytic Expression for the Eigenfunction?
After noticing that, in Figure 6, the envelope segments of all of the marginally unstable eigenfunctions merge into the same curve, we began to wonder whether a single expression — and, even better, an analytically defined expression — would perfectly describe the eigenfunction. We had reason to believe that this might actually be possible because, in pressuretruncated polytropic configurations, we have derived analytic expressions for the marginally unstable, fundamentalmode eigenfunctions of both and systems.
Very quickly, we convinced ourselves that a parabolic function does indeed perfectly match the "core" segment of each displayed eigenfunction. Specifically, throughout the core ,









The envelope segment posed a much greater challenge. In the context of our discussion of Radial Oscillations of n = 1 Polytropic Spheres, and in an accompanying Ramblings Appendix chapter we have detailed some trial derivations that are mostly blind alleyways. Twice — once in January, 2019 and again (independently) in April 2019 — we have analytically demonstrated that the following appears to work for the envelope:
Given that the Structural Properties of the envelope are described by the LaneEmden function,






the relevant LAWE is satisfied by the fractional displacement function,



where, is an arbitrary scale factor.
Note that,

But, as far as we have been able to determine (as of 16 April 2019), this analytic displacement function does not match the displacement function that has been generated through numerical integration of the LAWE (see the lightgreen segment of the eigenfunction displayed above in Figure 5). It remains unclear whether (a) the numerical integration is at fault, (b) we are imposing an incorrect slope at the coreenvelope interface, or ( c) we are misinterpreting how to compare the two separately derived (one, numerical, and the other, analytic) envelope eigenfunctions.
Other Modes
Our Determinations for Marginally Unstable Model Having 

Mode  
core  env  expected  measured  
1 (Fundamental) 
0.00  0.00  +0.81437470  0.455872  1.473523  +0.3820  1  0.999999992  n/a  n/a  n/a  n/a  n/a  n/a 
2  2.51513333  10.7107538  0.20482050  7.09124  5.4547441   0.9962  4.355376917  4.35537692  0.64133  0.3502  n/a  n/a  n/a  n/a 
3  5.72371888  24.3745901  0.14269277  +8.046019  +3.627611  +0.9308  11.18729505  11.18729506  0.4837  0.5864  0.842  0.0854  n/a  n/a 
4  10.3458476  44.0622916  0.20845197  0.6949966  1.61699793  1.1443  21.03114578  21.03114577  0.3939  0.7154  0.6902  0.2777  0.9115  0.0284 
Reconciliation
Figure 7: Conflicting Instability Regions  

Figure 7, shown here on the right, is identical to the righthand panel of Figure 4, as displayed above. From a standard global, freeenergy analysis — such as the one summarized above — we have determined that the reddashed curve shown in the right panel of Figure 6 divides the plane into dynamically stable (below and to the right) and unstable (above and to the left) regions.
Variational Principle
Setup
Let's follow the guidelines of the variational principle. Instead of starting with the form of the LAWE given above, namely,
we will start with a form that is more amenable to the variational principle, namely,



ASIDE: Let's show that these two expressions are equivalent. Remembering that,
the second expression becomes,
Hence, we must multiply the first expression through by in order to obtain the second expression. 
From above, we realize that multiplying the second expression through by gives,









That is, multiplying the second expression through by, , should give a desirable, totally dimensionless version of the LAWE. Remembering that,



; 






; 






; 



let's try it.









Now, guided by the accompanying summary, if we multiply through by and integrate over the entire volume, we obtain the governing variational relation, namely,


















Energy Normalization:





















Hence, the dimensionless governing variational relation becomes,
























where,









Or, for inclusion in our accompanying Tabular Overview,









Implementation
Table 3 


Core  Envelope  Virial  
Integrals over  Integrals over  Integrals over  Integrals over  Numerical  
Analytic 
Numerical 
TERM1  Analytic 
Numerical 
TERM2  Analytic 
Numerical 
TERM3  Analytic 
Numerical 
TERM4  
1  1.6686460157  3.021916335  3.021921  0.116389175  3.356583022  3.35666  2.649752079  1.47780476  1.47791  1.0720821  5.642859167  5.642820  1.91142893  0.000020  
2.27925811317  4.241287117  4.241410819  0.440878529  6.074241035  6.074317546  4.150731169  4.284931508  4.28547195  1.44651932  10.97819621  10.97847622  0.92598634  0.000057  
0.345  2.560146865247  4.639705843  4.6399114  0.6794857  7.125754184  7.125854025  4.5487829  11.72861751  11.730381  1.51410084  25.61089252  25.6115597  0.4496513  0.000097  
2.582007485476  4.667042505  4.667254935  0.700414598  7.200966267  7.201068684  4.57274936  13.15887139  13.1608467  1.51408246  28.45086152  28.45153761  0.4170461  0.000101  
0.309  2.6274239687695  4.722277318  4.722504339  0.744964507  7.354156963  7.3542507  4.61961058  17.1374434  17.1399773  1.51055838  36.36528446  36.36591543  0.3524855  0.000110  
2.7357711469398  4.84592201  4.846185027  0.857001395  7.70305421  7.703178009  4.7163542  37.84289623  37.8479208  1.47966673  77.67458196  77.67408155  0.2194152  0.000128  








































Table 4 


Analytic 
Numerical 
TERM5  
1  1.6686460157  4.818155928  4.818145  13.14874521  0.814374698  1.153014872  0.139506172  +0.455871977  +1.47352  0.0002973  
2.27925811317  9.020985415  9.021084268  23.0612702  0.653665497  1.574940686  0.049058481  +1.059668912  +1.835801347  0.0000705  
0.345  2.560146865247  17.41399388  17.4141672  24.63990825  0.563043202  1.769031527  0.030957085  +1.552125296  +2.131275177  0.0000695  
2.582007485476  18.8449906  18.8451614  24.52624951  0.555549156  1.78413696  0.029889634  +1.600041467  +2.16002488  0.0004312  
0.309  2.6274239687695  22.61541791  22.6155686  24.13633675  0.539776219  1.815519219  0.027798189  +1.705239188  +2.223143513  0.0001077  
2.7357711469398  39.12088278  39.1208075  22.36902623  0.501037082  1.890385851  0.023427588  +1.991720516  +2.39503231  0.0002679  

Revised FreeEnergy Analysis
If we set = constant in the variational principle relation, we have,






Notice the similarity between this last expression and the pair of expressions — numbered ⑥ and ⑦ — that arise in the context of pressuretruncated polytropes.
See Also
 K. De et al. (12 October 2018, Science, Vol. 362, No. 6411, pp. 201  206), A Hot and Fast Ultrastripped Supernova that likely formed a Compact Neutron Star Binary.
© 2014  2021 by Joel E. Tohline 