SSC/FreeEnergy/Powerpoint: Difference between revisions

From jetwiki
Jump to navigation Jump to search
Created page with "__FORCETOC__ =Supporting Derivations for Free-Energy PowerPoint Presentation= The derivations presented here are an extension of our User:Tohline/SSC/FreeEnergy/Polytrope..."
 
 
Line 4: Line 4:


The derivations presented here are an extension of our  
The derivations presented here are an extension of our  
[[User:Tohline/SSC/FreeEnergy/PolytropesEmbedded#Free-Energy_Synopsis|accompanying free-energy synopsis]].  These additional details proved to be helpful while developing an overarching PowerPoint presentation.
[[SSC/FreeEnergy/PolytropesEmbedded#Free-Energy_Synopsis|accompanying free-energy synopsis]].  These additional details proved to be helpful while developing an overarching PowerPoint presentation.


==General Free-Energy Expression==
==General Free-Energy Expression==
Line 29: Line 29:
</div>
</div>


As [[User:Tohline/SSC/FreeEnergy/PolytropesEmbedded#Equilibrium_Radii_and_Critical_Radii|we have shown]], setting,
As [[SSC/FreeEnergy/PolytropesEmbedded#Equilibrium_Radii_and_Critical_Radii|we have shown]], setting,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 127: Line 127:


====Case M====
====Case M====
Set <math>~K</math> and <math>~M_\mathrm{tot}</math> constant and examine how the free-energy behaves as a function of the coordinates, <math>~(R,P_e)</math>.  In this case (see, for example, [[User:Tohline/SSC/FreeEnergy/PolytropesEmbedded#Case_M|here]]),
Set <math>~K</math> and <math>~M_\mathrm{tot}</math> constant and examine how the free-energy behaves as a function of the coordinates, <math>~(R,P_e)</math>.  In this case (see, for example, [[SSC/FreeEnergy/PolytropesEmbedded#Case_M|here]]),
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 171: Line 171:
</div>
</div>


where the structural form factors for pressure-truncated polytropes are precisely defined [[User:Tohline/SSC/Virial/FormFactors#PTtable|here]].  And (see, for example, [[User:Tohline/SSC/FreeEnergy/PolytropesEmbedded#Case_M_Free-Energy_Surface|here]]),
where the structural form factors for pressure-truncated polytropes are precisely defined [[SSCpt1/Virial/FormFactors#PTtable|here]].  And (see, for example, [[SSC/FreeEnergy/PolytropesEmbedded#Case_M_Free-Energy_Surface|here]]),
<div align="center">
<div align="center">
<table border="0" cellpadding="3">
<table border="0" cellpadding="3">
Line 622: Line 622:


====Case P====
====Case P====
Set <math>~K</math> and <math>~P_e</math> constant and examine how the free-energy behaves as a function of the coordinates, <math>~(R,M_\mathrm{tot})</math>.  In this case (see, for example, [[User:Tohline/SSC/FreeEnergy/PolytropesEmbedded#Case_P|here]]),
Set <math>~K</math> and <math>~P_e</math> constant and examine how the free-energy behaves as a function of the coordinates, <math>~(R,M_\mathrm{tot})</math>.  In this case (see, for example, [[SSC/FreeEnergy/PolytropesEmbedded#Case_P|here]]),
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 694: Line 694:
</div>
</div>


where the structural form factors for pressure-truncated polytropes are precisely defined [[User:Tohline/SSC/Virial/FormFactors#PTtable|here]]. If we set all three structural form-factors to unity, we have,
where the structural form factors for pressure-truncated polytropes are precisely defined [[SSCpt1/Virial/FormFactors#PTtable|here]]. If we set all three structural form-factors to unity, we have,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 1,535: Line 1,535:
=====Structural Form Factors=====
=====Structural Form Factors=====


The following table of structural form factors has been drawn from [[User:Tohline/SSC/Virial/FormFactors#PTtable|here]],
The following table of structural form factors has been drawn from [[SSCpt1/Virial/FormFactors#PTtable|here]],


<div align="center" id="PTtable">
<div align="center" id="PTtable">
Line 1,643: Line 1,643:
</div>
</div>


and [[User:Tohline/SSC/Virial/FormFactors#Summary_.28n.3D5.29|here]],
and [[SSCpt1/Virial/FormFactors#Summary_.28n.3D5.29|here]],


<div align="center">
<div align="center">
Line 1,708: Line 1,708:


=====Case M Equilibrium Conditions=====
=====Case M Equilibrium Conditions=====
Employing the [[User:Tohline/SphericallySymmetricConfigurations/Virial#Choices_Made_by_Other_Researchers|renormalization factors]],
Employing the [[SSCpt1/Virial#Choices_Made_by_Other_Researchers|renormalization factors]],


<div align="center">
<div align="center">
Line 1,739: Line 1,739:
</div>
</div>


we find from detailed force-balance analyses that the [[User:Tohline/SSC/Structure/PolytropesEmbedded#General_Properties|equilibrium radius and corresponding external pressure]] for "Case M" configurations are,
we find from detailed force-balance analyses that the [[SSC/Structure/PolytropesEmbedded#General_Properties|equilibrium radius and corresponding external pressure]] for "Case M" configurations are,
<div align="center">
<div align="center">
<table border="0" cellpadding="3">
<table border="0" cellpadding="3">
Line 1,760: Line 1,760:
</table>
</table>
</div>
</div>
which matches the expression derived in an [[User:Tohline/SSC/Structure/Polytropes#Lane-Emden_Equation|ASIDE box found with our introduction of the Lane-Emden equation]], and
which matches the expression derived in an [[SSC/Structure/Polytropes#Lane-Emden_Equation|ASIDE box found with our introduction of the Lane-Emden equation]], and
<div align="center">
<div align="center">
<table border="0" cellpadding="3">
<table border="0" cellpadding="3">
Line 1,815: Line 1,815:
</div>
</div>


There are two turning points:  One associated with a maximum in <math>~P_e</math> and one associated with a maximum in <math>~R_\mathrm{eq}</math>.  According to [[User:Tohline/SSC/Structure/PolytropesEmbedded#Turning_Points|Kimura's discussion]], the first of these occurs in the configuration for which,
There are two turning points:  One associated with a maximum in <math>~P_e</math> and one associated with a maximum in <math>~R_\mathrm{eq}</math>.  According to [[SSC/Structure/PolytropesEmbedded#Turning_Points|Kimura's discussion]], the first of these occurs in the configuration for which,


<div align="center">
<div align="center">
Line 2,075: Line 2,075:
</div>
</div>
-->
-->
This agrees with the [[User:Tohline/SSC/FreeEnergy/PolytropesEmbedded#Case_M|expression derived in a separate <font color="red">ASIDE</font>]]; as was pointed out in that context, the root of this equation is:  <math>~\ell \approx 2.223175</math>, that is, <math>~\tilde\xi \approx 3.85065</math>.  This point along the equilibrium sequence is identified by the red circular dot in Figure 3, below.  It is almost &#8212; but definitely not &#8212; coincident with the configuration along the sequence (marked by the yellow circular dot) that is associated with the minimum-radius turning point.
This agrees with the [[SSC/FreeEnergy/PolytropesEmbedded#Case_M|expression derived in a separate <font color="red">ASIDE</font>]]; as was pointed out in that context, the root of this equation is:  <math>~\ell \approx 2.223175</math>, that is, <math>~\tilde\xi \approx 3.85065</math>.  This point along the equilibrium sequence is identified by the red circular dot in Figure 3, below.  It is almost &#8212; but definitely not &#8212; coincident with the configuration along the sequence (marked by the yellow circular dot) that is associated with the minimum-radius turning point.


<div align="center" id="Figure3">
<div align="center" id="Figure3">
Line 2,101: Line 2,101:
=====Case P Equilibrium Conditions=====
=====Case P Equilibrium Conditions=====


The [[User:Tohline/SSC/Structure/PolytropesEmbedded#Stahler.27s_Presentation|equilibrium radius and corresponding configuration mass]] from a "Case P" analysis are,
The [[SSC/Structure/PolytropesEmbedded#Stahler.27s_Presentation|equilibrium radius and corresponding configuration mass]] from a "Case P" analysis are,
<div align="center">
<div align="center">
<table border="0" cellpadding="3">
<table border="0" cellpadding="3">
Line 2,140: Line 2,140:
</div>
</div>


According to our [[User:Tohline/SSC/Structure/PolytropesEmbedded#Other_Limits|review of, especially, Kimura's work]], the turning point associated with <math>~M_\mathrm{max}</math> occurs where,
According to our [[SSC/Structure/PolytropesEmbedded#Other_Limits|review of, especially, Kimura's work]], the turning point associated with <math>~M_\mathrm{max}</math> occurs where,


<div align="center">
<div align="center">
Line 2,386: Line 2,386:
-->
-->


Comparing the general free-energy expression at the beginning of this chapter with the free-energy expression provided via our [[User:Tohline/SSC/FreeEnergy/PolytropesEmbedded#BiPolytropeFreeEnergy|accompanying summary discussion of five-one bipolytropes]], we find that,
Comparing the general free-energy expression at the beginning of this chapter with the free-energy expression provided via our [[SSC/FreeEnergy/PolytropesEmbedded#BiPolytropeFreeEnergy|accompanying summary discussion of five-one bipolytropes]], we find that,


<div align="center">
<div align="center">
Line 2,606: Line 2,606:
</div>
</div>


Consistent with our [[User:Tohline/SSC/BipolytropeGeneralization#Free_Energy_and_Its_Derivatives|generic discussion of the stability of bipolytropes]] and the ''specific'' discussion of [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#Stability_Condition|the stability of bipolytropes having]] <math>~(n_c, n_e) = (5, 1)</math>, it can straightforwardly be shown that <math>~\partial \mathfrak{G}/\partial \chi = 0</math> is satisfied by setting <math>~\Chi = 1</math>; that is, the equilibrium condition is,
Consistent with our [[SSC/BipolytropeGeneralization#Free_Energy_and_Its_Derivatives|generic discussion of the stability of bipolytropes]] and the ''specific'' discussion of [[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|the stability of bipolytropes having]] <math>~(n_c, n_e) = (5, 1)</math>, it can straightforwardly be shown that <math>~\partial \mathfrak{G}/\partial \chi = 0</math> is satisfied by setting <math>~\Chi = 1</math>; that is, the equilibrium condition is,


Furthermore, the equilibrium configuration is unstable whenever <math>~\partial^2 \mathfrak{G}/\partial \chi^2 < 0</math>, that is, the transition from stable to unstable configurations whenever,
Furthermore, the equilibrium configuration is unstable whenever <math>~\partial^2 \mathfrak{G}/\partial \chi^2 < 0</math>, that is, the transition from stable to unstable configurations whenever,
Line 2,626: Line 2,626:
</table>
</table>
</div>
</div>
[[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#Stability_Condition|Table 1 of an accompanying chapter]] &#8212; and the red-dashed curve in the figure adjacent to that table &#8212; identifies some key properties of the model that marks the transition from stable to unstable configurations along equilibrium sequences that have various values of the mean-molecular weight ratio, <math>~\mu_e/\mu_c</math>.
[[SSC/Structure/BiPolytropes/Analytic51#Stability_Condition|Table 1 of an accompanying chapter]] &#8212; and the red-dashed curve in the figure adjacent to that table &#8212; identifies some key properties of the model that marks the transition from stable to unstable configurations along equilibrium sequences that have various values of the mean-molecular weight ratio, <math>~\mu_e/\mu_c</math>.


====Coincidence Between Points of Secular and Dynamical Instability====
====Coincidence Between Points of Secular and Dynamical Instability====
From the [[User:Tohline/SSC/Structure/BiPolytropes/Analytic5_1#Model_Sequences|accompanying graphical display of equilibrium sequences]], it seems that a <math>~\nu_\mathrm{max}</math> turning point will only exist in five-one bipolytropes for <math>~\mu_e/\mu_c</math> less than some value &#8212; call it, <math>~(\mu_e/\mu_c)_\mathrm{begin}</math> &#8212; which is less than but approximately equal to <math>~\tfrac{1}{3}</math>.  As we move along any sequence for which <math>~\mu_e/\mu_c < (\mu_e/\mu_c)_\mathrm{begin}</math>, in the direction of increasing <math>~\ell_i</math>, it is fair to ask whether the system becomes dynamically unstable (at <math>~[x_\mathrm{eq}]_\mathrm{crit}</math>) before or after it encounters the point of secular instability marked by <math>~\nu_\mathrm{max}</math>.
From the [[SSC/Structure/BiPolytropes/Analytic51#Model_Sequences|accompanying graphical display of equilibrium sequences]], it seems that a <math>~\nu_\mathrm{max}</math> turning point will only exist in five-one bipolytropes for <math>~\mu_e/\mu_c</math> less than some value &#8212; call it, <math>~(\mu_e/\mu_c)_\mathrm{begin}</math> &#8212; which is less than but approximately equal to <math>~\tfrac{1}{3}</math>.  As we move along any sequence for which <math>~\mu_e/\mu_c < (\mu_e/\mu_c)_\mathrm{begin}</math>, in the direction of increasing <math>~\ell_i</math>, it is fair to ask whether the system becomes dynamically unstable (at <math>~[x_\mathrm{eq}]_\mathrm{crit}</math>) before or after it encounters the point of secular instability marked by <math>~\nu_\mathrm{max}</math>.


{{ SGFfooter }}
{{ SGFfooter }}

Latest revision as of 19:40, 31 July 2021


Supporting Derivations for Free-Energy PowerPoint Presentation

The derivations presented here are an extension of our accompanying free-energy synopsis. These additional details proved to be helpful while developing an overarching PowerPoint presentation.

General Free-Energy Expression

We're considering a free-energy function of the following form:

𝔊type*

=

ax1+bx3/n+cx3/j+𝔊0,

where,

xRR0.

As we have shown, setting,

𝔊type*x

=

0,

generates a mathematical statement of virial equilibrium, namely,

bncxeq(n3)/na3c+1jxeq(j3)/j

=

0.

And equilibrium configurations for which the second (as well as first) derivative of the free energy is zero are found at "critical" radii given by the expression,

[xeq(j3)/j]crit

=

a32c[j2(n3)nj].

Pressure-Truncated Polytropes

For pressure-truncated polytropes, set j=1 and let n be the chosen polytropic index. In this case, the statement of virial equilibrium is,

bncxeq(n3)/na3cxeq4

=

0;

and the critical equilibrium configuration has,

[xeq]crit

=

[a(n3)32c(n+1)]1/4.

Case M

Set K and Mtot constant and examine how the free-energy behaves as a function of the coordinates, (R,Pe). In this case (see, for example, here),

a

35𝔣~W𝔣~M2,

b

n(4π3)1/n𝔣~A𝔣~M(n+1)/n,

c

4π3(PePnorm),

where the structural form factors for pressure-truncated polytropes are precisely defined here. And (see, for example, here),

R0=Rnorm

[(GK)nMtotn1]1/(n3),

Pnorm

[K4nG3(n+1)Mtot2(n+1)]1/(n3).

If we set all three structural form-factors to unity, we have,

a3c

=

3225π(PePnorm)1,

bnc

=

(34π)(n+1)/n(PePnorm)1.

Virial Equilibrium
Figure 1

Simplest Case M sequence

So the statement of virial equilibrium becomes,

xeq4

=

[(34π)(n+1)/nxeq(n3)/n3225π](PePnorm)1

 

=

3225π[5(34π)1/nxeq(n3)/n1](PePnorm)1

PePnorm

=

3225πxeq4[5(34π)1/nxeq(n3)/n1].

The light-blue dots in Figure 1 trace the equilibrium sequence that is defined by this virial equilibrium function in the case of n=5.

Dynamical Instability

Along the "Case M" equilibrium sequence, the transition from stable to unstable configurations occurs at,

[xeq]crit4

=

[(n3)3(n+1)]3225π(PePnorm)1

225π3(PePnorm)[xeq]crit4

=

[(n3)3(n+1)]

which, in combination with the virial equilibrium condition gives,

5(34π)1/n[xeq]crit(n3)/n1

=

[(n3)3(n+1)]

[xeq]crit

=

[4n35(n+1)(4π3)1/n]n/(n3).

The location of this critical configuration along the n=5 equilibrium sequence is marked by the red circular dot in Figure 1.

Turning Point

Let's examine the curvature of the equilibrium sequence.

ddx(PePnorm)

=

35πx5[5(34π)1/nx(n3)/n1]+3(n3)22nπx4(34π)1/nx3/n

 

=

35πx5+34π(34π)1/n[(n3)n4]x(n3)/nx5

 

=

35πx53(34π)(n+1)/n[n+1n]x(n3)/nx5.

Setting this derivative to zero let's us identify the location of the turning point that identifies Pmax.

[xeq(n3)/n]turn

=

15π[nn+1](4π3)(n+1)/n

[xeq]turn

=

[4n15(n+1)(4π3)1/n]n/(n3).

And, returning to the virial equilibrium expression, we find that, associated with this equilibrium radius,

PmaxPnorm

=

3225πxturn4[5(34π)1/nxturn(n3)/n1]

225πxturn4(PmaxPnorm)

=

15(34π)1/n[4n15(n+1)(4π3)1/n]3

 

=

(n3n+1)

PmaxPnorm

=

120π(n3n+1)[15(n+1)4n(34π)1/n]4n/(n3).

Notice that, under the assumption that all three structural filling-factors are unity, [xeq]turn=[xeq]crit, that is, the location of the turning point coincides precisely with the point along the equilibrium sequence where the transition from stable to unstable equilibrium configurations occurs (marked by the red circular dot in Figure 1).

Case M Summary

Case M

Order-of-Magnitude Analysis: Assume 𝔣~M=𝔣~W=𝔣~A=1

Virial Equilibrium:

PePnorm

=

3225π(ReqRnorm)4[5(34π)1/n(ReqRnorm)eq(n3)/n1]

Dynamical Instability:

(n>3)

ReqRnorm|crit

=

[4n35(n+1)(4π3)1/n]n/(n3)

Turning Point (Pmax):

(n>3)

PmaxPnorm

=

120π(n3n+1)[15(n+1)4n(34π)1/n]4n/(n3)

ReqRnorm|turn

=

[4n35(n+1)(4π3)1/n]n/(n3)

Case P

Set K and Pe constant and examine how the free-energy behaves as a function of the coordinates, (R,Mtot). In this case (see, for example, here),

a

=

35𝔣~W𝔣~M2(n+1n)n/(n3)(MtotMSWS)(n1)/(n3),

b

=

n(4π3)1/n𝔣~A𝔣~M(n+1)/n(n+1n)3/(n3)(MtotMSWS)3(n1)/[n(n3)],

c

=

4π3(n+1n)3/(n3)(MtotMSWS)(5n)/(n3),

R0=RSWS

(n+1n)1/2G1/2Knn/(n+1)Pe(1n)/[2(n+1)],

MSWS

(n+1n)3/2G3/2Kn2n/(n+1)Pe(3n)/[2(n+1)].

where the structural form factors for pressure-truncated polytropes are precisely defined here. If we set all three structural form-factors to unity, we have,

a3c

=

320π(n+1n)(MtotMSWS)2

bnc

=

[(34π)(MtotMSWS)](n+1)/n

Virial Equilibrium

So, the statement of virial equilibrium becomes,

[(34π)(MtotMSWS)](n+1)/nxeq(n3)/n320π(n+1n)(MtotMSWS)2xeq4

=

0.


ξ1


Known Analytic Lane-Emden Functions

n

θn(ξ)

θn'

ξ1

0

1ξ26

ξ3

6

1

sinξξ

sinξξ2cosξξ

π

5

[1+ξ23]1/2

ξ3[1+ξ23]3/2

Dynamical Instability

Along the "Case P" equilibrium sequence, the transition from stable to unstable configurations occurs at,

[xeq]crit4

=

(n3)3(n+1)[320π(n+1n)(MtotMSWS)2]

 

=

(n3)20πn(MtotMSWS)2=22π(n3)325n[34π(MtotMSWS)]2,

which, in combination with the "Case P" virial equilibrium expression gives,

0

=

{[325n22π(n3)]1/2xcrit2}(n+1)/nxcrit(n3)/n320π(n+1n){[20πn(n3)]xcrit4}xcrit4

 

=

[325n22π(n3)](n+1)/(2n)xcrit(3n1)/nxcrit4{[3(n+1)(n3)]+1}

 

=

[325n22π(n3)](n+1)/(2n)xcrit(3n1)/n[4n(n3)]xcrit4

xcrit(n+1)/n

=

[(n3)4n][325n22π(n3)](n+1)/(2n)

Turning Points

Let's simplify the notation, defining,

m34π(MtotMSWS).

The statement of virial equilibrium becomes,

m(n+1)/nx(n3)/nc0m2x4

=

0,

where,

c0[4π(n+1)15n].

Differentiating gives,

0

=

(n+1n)m1/nx(n3)/ndm+(n3n)m(n+1)/nx3/ndx2c0mdm4x3dx

 

=

[(n+1)m1/nx(n3)/n2c0nm]dm+[(n3)m(n+1)/nx3/n4nx3]dx

dmdx

=

4nx3(n3)m(n+1)/nx3/n(n+1)m1/nx(n3)/n2c0nm.


One turning point occurs where the numerator is zero, that is,

4nx3

=

(n3)m(n+1)/nx3/n

4nx3(n+1)/n

=

(n3)m(n+1)/n

mx3

=

[4n(n3)]n/(n+1).

Plugging this into the virial equilibrium expression gives,

0

=

{[4n(n3)]n/(n+1)x3}(n+1)/nx(n3)/nc0{[4n(n3)]n/(n+1)x3}2x4

 

=

[4n(n3)]x4c0[4n(n3)]2n/(n+1)x6x4

c0[4n(n3)]2n/(n+1)x2

=

[4n(n3)]1

x2

=

[15n4π(n+1)][3(n+1)(n3)][(n3)4n]2n/(n+1)

 

=

[325n4π(n3)][(n3)4n]2n/(n+1).

The associated mass is,

(MmaxMSWS)=(4π3)m

=

4π3[n34n]n/(n+1)[325n4π(n3)]3/2[(n3)4n]3n/(n+1)

 

=

4π3[n34n]2n/(n+1)[325n4π(n3)]3/2.

Notice that, for n=3,

(MmaxMSWS)n=3

=

[345328π]1/2.


Another turning point occurs where the denominator is zero, that is,

(n+1)m1/nx(n3)/n

=

2c0nm

(n+1)x(n3)/n

=

2c0nm(n1)/n

xn3mn1

=

[2c0n(n+1)]n

 

=

[8π15]n

m

=

[x(n3)/(n1)(158π)n/(n1)].

Plugging this into the virial equilibrium expression gives,

x4

=

[x(n3)/(n1)(158π)n/(n1)](n+1)/nx(n3)/nc0[x(n3)/(n1)(158π)n/(n1)]2

 

=

(158π)(n+1)/(n1)x(n3)(n+1)/[n(n1)]x(n3)/nc0(158π)2n/(n1)x2(n3)/(n1)

 

=

(158π)(n+1)/(n1)x2(n3)/(n1)[(n+1)2n](158π)(n+1)/(n1)x2(n3)/(n1)

 

=

(158π)(n+1)/(n1)x2(n3)/(n1)[1(n+1)2n]

x2(n+1)/(n1)

=

(n1)2n(158π)(n+1)/(n1)

RmaxRSWS

=

[(n1)2n](n1)/[2(n+1)](158π)1/2.

And the associated mass is,

MtotMSWS|turn

=

4π3(158π)n/(n1){[(n1)2n](n1)/[2(n+1)](158π)1/2}(n3)/(n1)

 

=

4π3(158π)3/2[(n1)2n](n3)/[2(n+1)].

Case P Summary

Case P

Order-of-Magnitude Analysis: Assume 𝔣~M=𝔣~W=𝔣~A=1

Virial Equilibrium:

[(34π)(MtotMSWS)](n+1)/n(ReqRSWS)(n3)/n320π(n+1n)(MtotMSWS)2(ReqRSWS)4

=

0

Dynamical Instability:

(n>3)

ReqRSWS|crit

=

[(n3)4n]n/(n+1)[325n22π(n3)]1/2

Turning Point (Mmax):

(n>3)

MmaxMSWS

=

4π3[n34n]2n/(n+1)[325n4π(n3)]3/2

ReqRSWS|turn

=

[(n3)4n]n/(n+1)[325n22π(n3)]1/2

Turning Point (Rmax):

(n>1)

RmaxRSWS

=

[(n1)2n](n1)/[2(n+1)](158π)1/2

MtotMSWS|turn

=

4π3(158π)3/2[(n1)2n](n3)/[2(n+1)]


Figure 2
Stahler (1983) Title Page
Stahler (1983) Title Page
Case P equilibrium sequences with key configurations highlighted
Case P equilibrium sequences with key configurations highlighted
Stahler (1983) Figure 17 (edited)
Stahler (1983) Figure 17 (edited)

The right-hand panel of Figure 2 presents substantial segments of Case P virial equilibrium sequences for a range of polytropic indexes (n = 1, 2, 2.8, 3, 3.5, 4, 5). For each sequence, the location of the Rmax and Mmax turning points — if they exist — are denoted by a yellow or red circular dot, respectively. The point along each (n3) sequence at which the transition from dynamically stable to dynamically unstable structures occurs coincides with the location of Mmax (i.e., with the red circular dot).

For display purposes, all normalized masses (Mtot/MSWS) have been further normalized to the maximum mass on the n = 3 sequence.

Detailed Force-Balance Models

Structural Form Factors

The following table of structural form factors has been drawn from here,

Structural Form Factors for Isolated Polytropes

Structural Form Factors for Pressure-Truncated Polytropes

𝔣M

=

[3θ'ξ]ξ1

𝔣W

=

3255n[θ'ξ]ξ12

𝔣A

=

3(n+1)(5n)[θ']ξ12

𝔣~M

=

(3θ~'ξ~)

𝔣~W

=

35(5n)ξ~2[θ~n+1+3(θ~')2𝔣~Mθ~]

𝔣~A

=

1(5n){6θ~n+1+(n+1)[3(θ~')2𝔣~Mθ~]}

and here,

Structural Form Factors for Pressure-Truncated n = 5 Polytropes

𝔣~M

=

(1+2)3/2

𝔣~W

=

5245[(48321)(1+2)3+tan1()]

𝔣~A

=

3233[tan1()+(21)(1+2)2]

where,       ξ~3

Case M Equilibrium Conditions

Employing the renormalization factors,

RHoredtRnorm

=

[4π(n+1)n]1/(n3),

PHoredtPnorm

=

[(n+1)34π](n+1)/(n3),

we find from detailed force-balance analyses that the equilibrium radius and corresponding external pressure for "Case M" configurations are,

ReqRnorm=ReqRHoredtRHoredtRnorm

=

[4π(n+1)n]1/(n3)ξ~(ξ~2θ~)(1n)/(n3),

which matches the expression derived in an ASIDE box found with our introduction of the Lane-Emden equation, and

PePnorm=PePHoredtPHoredtPnorm

=

[(n+1)34π](n+1)/(n3)θ~nn+1(ξ~2θ~)2(n+1)/(n3).

There are two turning points: One associated with a maximum in Pe and one associated with a maximum in Req. According to Kimura's discussion, the first of these occurs in the configuration for which,

θ~n+1(θ~')2

=

(n3)2

      …      

For n = 5, this occurs when ξ~=3

This point along the equilibrium sequence is identified by the dark green circular dot in Figure 3, below. The second occurs in the configuration for which,

ξ~θ~n(θ~')

=

(n3)(n1)

      …      

For n = 5, this occurs when ξ~=153.87298

This point along the equilibrium sequence is identified by the yellow circular dot in Figure 3, below. In addition, we have identified the point of dynamical instability.

(ReqRnorm)4

=

(n3)32(n+1)ac=133322π(PePnorm)135𝔣~W𝔣~M2

 

=

12235π(PePnorm)1524(1+2)35[(48321)(1+2)3+tan1]

(PePnorm)(ReqRnorm)4

=

1263π5[(48321)+(1+2)3tan1]

But the equilibrium condition for n = 5 configurations is,

(PePnorm)(ReqRnorm)4

=

[(n+1)34π](n+1)/(n3)θ~nn+1(ξ~2θ~)2(n+1)/(n3){[4π(n+1)n]1/(n3)ξ~(ξ~2θ~)(1n)/(n3)}4

 

=

[233π]3θ~n6ξ~12(θ~)6{16[π233]1/2ξ~3(θ~)2}4

 

=

12434[233π]θ~n6(θ~)2

 

=

1233π(1+2)3[31/2(1+2)3/2]2

 

=

123π2

Putting the two expressions together gives,

123π2

=

1263π5[(48321)+(1+2)3tan1]

2333

=

5833+(1+2)3tan1

[(1+2)3]tan1

=

1+[2453]24

This agrees with the expression derived in a separate ASIDE; as was pointed out in that context, the root of this equation is: 2.223175, that is, ξ~3.85065. This point along the equilibrium sequence is identified by the red circular dot in Figure 3, below. It is almost — but definitely not — coincident with the configuration along the sequence (marked by the yellow circular dot) that is associated with the minimum-radius turning point.

Figure 3:   Case M Equilibrium Sequence

Case M pressure-truncated n=5 polytrope
Case M pressure-truncated n=5 polytrope


Now for a movie!

Case M movie of free-energy structures

Case P Equilibrium Conditions

The equilibrium radius and corresponding configuration mass from a "Case P" analysis are,

ReqRSWS

=

(n4π)1/2ξ~θ~n(n1)/2,

MtotMSWS

=

(n34π)1/2θ~n(n3)/2(ξ~2θ~').

According to our review of, especially, Kimura's work, the turning point associated with Mmax occurs where,

θ~n+1(θ~')2

=

(n3)2

      …      

For n = 5, this occurs when ξ~=3.

And a turning point associated with Rmax occurs where,

ξ~(θ~')θ~

=

2(n1).

For n=5 configurations, this means,

12

=

ξ~(1+ξ~23)1/2ξ~3(1+ξ~23)3/2

 

=

2(1+2)

(1+2)

=

22

2

=

1

ξ~

=

3.

Five-One Bipolytropes

Basic Properties

For bipolytropes, in general, let n=nc and j=ne. The statement of virial equilibrium is, then,

bnccxeq(nc3)/nca3c+1nexeq(ne3)/ne

=

0.

And the critical equilibrium configuration has,

[xeq(ne3)/ne]crit

=

a32c[ne2(nc3)ncne].

Here we choose to set nc=5 and ne=1. Hence, these two conditions become, respectively,

b5cxeq2/5+xeq2

=

a3c,

and,

[xeq]crit

=

[232ca]1/2.


Comparing the general free-energy expression at the beginning of this chapter with the free-energy expression provided via our accompanying summary discussion of five-one bipolytropes, we find that,

a

=

(3𝔏i+12𝔎i)xeq[(3625π)1/2i3(1+i2)3]

b

=

(5𝔏i)xeq3/5[(3625π)1/2i3(1+i2)3]

c

=

(4𝔎i)xeq3[(3625π)1/2i3(1+i2)3]

where,

𝔏i

(i41)i2+(1+i2)3i3tan1i,

𝔎i

Λiηi+(1+Λi2)ηi[π2+tan1Λi],

Λi

1ηii,

ηi

=

3(μeμc)[i(1+i2)],

and,

xeq=ReqRnorm

=

(π23)1/2ν2q(1+i2)333i5

 

=

(π23)1/2(1+i2)333i3[1(1+Λi2)]{1+1ηi[π2+tan1Λi]}1


In rewriting this last expression, we have made use of the two relations,

qriR

=

{1+1ηi[π2+tan1Λi]}1,

νMcoreMtot

=

iq(1+Λi2)1/2.

Consistent with our generic discussion of the stability of bipolytropes and the specific discussion of the stability of bipolytropes having (nc,ne)=(5,1), it can straightforwardly be shown that 𝔊/χ=0 is satisfied by setting X=1; that is, the equilibrium condition is,

Furthermore, the equilibrium configuration is unstable whenever 2𝔊/χ2<0, that is, the transition from stable to unstable configurations whenever,

𝔏i𝔎i

=

20.

Table 1 of an accompanying chapter — and the red-dashed curve in the figure adjacent to that table — identifies some key properties of the model that marks the transition from stable to unstable configurations along equilibrium sequences that have various values of the mean-molecular weight ratio, μe/μc.

Coincidence Between Points of Secular and Dynamical Instability

From the accompanying graphical display of equilibrium sequences, it seems that a νmax turning point will only exist in five-one bipolytropes for μe/μc less than some value — call it, (μe/μc)begin — which is less than but approximately equal to 13. As we move along any sequence for which μe/μc<(μe/μc)begin, in the direction of increasing i, it is fair to ask whether the system becomes dynamically unstable (at [xeq]crit) before or after it encounters the point of secular instability marked by νmax.

Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |