SSC/FreeEnergy/PolytropesEmbedded/Pt3B

From JETohlineWiki
Jump to navigation Jump to search

Free Energy of Embedded Polytropes[edit]


Part I:   Synopsis
 

Part II:  Truncated Polytropes
 

Part III:  Free-Energy of Bipolytropes
 

IIIA:  Focus on (5, 1) Bipolytropes
 

IIIB:  Focus on (0, 0) Bipolytropes
 

IIIC:  Overview
 


  • Core mass, Mc, and associated dimensionless mass fraction, νMc/Mtot;
  • Polytropic constant in the core, Kc.

In general, the warped free-energy surface drapes across a five-dimensional parameter "plane" such that,

𝔊

=

𝔊(R,Kc,Mtot,q,ν).

Focus on Zero-Zero Free-Energy Expression[edit]

Here, we will draw heavily from the following accompanying chapters:


From Detailed Force-Balance Models[edit]

Equilibrium Radius[edit]

First View[edit]

In an accompanying chapter we find,

P0Req4GMtot2

=

(323π)(νq3)2[q2+25q5(f1𝔉)]

where,

f

1+52(ρeρc)(1q21)+(ρeρc)2[1q51+52(11q2)],

𝔉

52(ρeρc)1q5[(2q2+3q3q5)+35(ρeρc)(1+5q25q3+q5)],

ρeρc

=

q3(1ν)ν(1q3).

Here, we prefer to normalize the equilibrium radius to Rnorm. So, let's replace the central pressure with its expression in terms of Kc. Specifically,

P0

=

Kcρcγc=Kc[3Mcore4πRi3]γc=Kc[3νMtot4πq3Req3](nc+1)/ncP0Pnorm=[34π(νq3)1χeq3](nc+1)/nc

Kc[3νMtot4πq3Req3](nc+1)/ncReq4GMtot2

=

(323π)(νq3)2[q2+25q5(f1𝔉)]

Req(nc3)/nc

=

(GKc)Mtot(nc1)/nc[3ν4πq3](nc+1)/nc(323π)(νq3)2[q2+25q5(f1𝔉)]

χeq(nc3)/nc[ReqRnorm](nc3)/nc

=

12(4π3)1/nc(νq3)(nc1)/nc[q2+25q5(f1𝔉)].

Or, in terms of γc,

χeq43γc

=

12(34π)1γc(νq3)2γc[q2+25q5(f1𝔉)].

Second View[edit]

Alternatively, from our derivation and discussion of analytic detailed force-balance models,

[R4GMtot2]P0

  = 

(323π)ν2g2q4,

where,

[g(ν,q)]2

1+(ρeρ0)[2(1ρeρ0)(1q)+ρeρ0(1q21)].

In order to show that this expression is the same as the other one, above, we need to show that,

(323π)(νq3)2[q2+25q5(f1𝔉)]

=

(323π)ν2g2q4

f1𝔉

=

52q3[g21]

 

=

52q3(ρeρ0)[2(1ρeρ0)(1q)+ρeρ0(1q21)]

 

=

52q5(ρeρ0){2(q2q3)+ρeρ0[13q2+2q3]}.

Let's see …

f1𝔉

=

52(ρeρc)(1q21)+(ρeρc)2[1q51+52(11q2)]52(ρeρc)1q5[(2q2+3q3q5)+35(ρeρc)(1+5q25q3+q5)]

 

=

52(ρeρc)(1q21)52(ρeρc)1q5[(2q2+3q3q5)]

 

 

52(ρeρc)1q5[35(ρeρc)(1+5q25q3+q5)]+(ρeρc)2[1q51+52(11q2)]

 

=

52(ρeρc)1q5{(q3q5)+(2q23q3+q5)}

 

 

+12(ρeρc)21q5[3(15q2+5q3q5)]+12(ρeρc)21q5[22q5+5(q5q3)]

 

=

52(ρeρc)1q5[(q3q5)+(2q23q3+q5)]+12(ρeρc)21q5[3(15q2+5q3q5)+22q5+5(q5q3)]

 

=

52(ρeρc)1q5[2q22q3]+52q5(ρeρc)2[13q2+2q3].

Q.E.D.

Hence, the equilibrium radius can also be written as,

χeq43γc

=

12(34π)1γc(νq3)2γcq2g2;

or, in terms of the polytropic index,

χeqnc3

=

12nc(4π3)(νq3)nc1(qg)2nc.

Gravitational Potential Energy[edit]

Also from our accompanying discussion, we have,

WgravEnorm

=

X1(35)(νq3)2q5[12nc(4π3)(νq3)nc1(qg)2nc]1/(nc3)f(ν,q)

 

=

X1(65)q5f[2nc(nc3)(34π)(νq3)(1nc)+2(nc3)bξnc]1/(nc3)

 

=

X1(65)q5f[(6π)(νq3)nc5bξnc]1/(nc3).

Internal Energy Components[edit]

First View[edit]

Before writing out the expressions for the internal energy of the core and of the envelope, we note from our separate detailed derivation that, in either case,

[Piχ3γPnorm]eqχ33γ

=

[(PiP0)(P0Pnorm)χ3]eq[χχeq]33γ

 

=

{(PiP0)[34π(νq3)]γcχ33γc}eqX33γ,

where, in equilibrium,

(PiP0)eq

=

1bξq2

bξq2

=

{25q3f+[125q3(1+𝔉)]}1

 

=

[1+25q3(f1𝔉)]1

So, copying from our accompanying detailed derivation, we have,

(𝔖AEnorm)core

=

4π/3(γc1){(PiP0)[34π(νq3)]γcχ33γc}eqX33γc{(P0Pic)[q3(3bξ5)q5]}

 

=

1(γc1)[(4π3)1γc(νq3)γcχeq33γc]X33γcq3[1(35)bξq2],

(𝔖AEnorm)env

=

4π/3(γe1){(PiP0)[34π(νq3)]γcχ33γc}eqX33γe{(1q3)+bξ(P0Pie)[25q5𝔉]}

 

=

1(γe1)[(4π3)1γc(νq3)γcχeq33γc]X33γe(PiP0){(1q3)+bξ(P0Pie)[25q5𝔉]}

 

=

1(γe1)[(4π3)1γc(νq3)γcχeq33γc]X33γe{(1bξq2)(1q3)+bξ[25q5𝔉]}

 

=

1(γe1)[(4π3)1γc(νq3)γcχeq33γc]X33γe(1q3){1[125(q31q3)𝔉]bξq2}.

Furthermore,

[(4π3)1γc(νq3)γcχeq33γc]

=

(34π)γc1(νq3)γc{χeq43γc}(33γc)/(43γc)

 

=

(34π)γc1(νq3)γc{12(34π)1γc(νq3)2γc[q2+25q5(f1𝔉)]}(33γc)/(43γc)

 

=

(34π)(γc1)/(43γc)(νq3)(65γc)(43γc){q22[1+25q3(f1𝔉)]}(33γc)/(43γc)

 

=

(34π)1/(nc3)(νq3)(nc5)(nc3){q22[1+25q3(f1𝔉)]}3/(nc3)

 

=

[(23π)(νq3)(nc5)bξ3]1/(nc3).

Hence, we have,

(𝔖AEnorm)core

=

nc[(4π3)1γc(νq3)γcχeq33γc]X3/ncq3[1(35)bξq2]

 

=

nc[(23π)(νq3)(nc5)bξ3]1/(nc3)X3/ncq3[1(35)bξq2],

(𝔖AEnorm)env

=

ne[(23π)(νq3)(nc5)bξ3]1/(nc3)X3/ne(1q3){1[125(q31q3)𝔉]bξq2}.


Second View[edit]

In our accompanying discussion of energies associated with detailed force balance models, we used the notation,

Π

(323π)GMtot2R4(νq3)2=Pnormχ4(323π)(νq3)2,

which allows us to rewrite the above quoted relationship between the central pressure and the radius of the bipolytrope as,

P0=Π(qg)2.

We also showed that, in equilibrium, the relationship between the central pressure and the interface pressure is,

P0=Pi+Πeqq2.

This means that, in equilibrium, the ratio of the interface pressure to the central pressure is,

(PiP0)eq

=

1Πeqq2P0=11g2,

or given that (see above),

52q3[g21]

=

f1𝔉

g2

=

1+25q3(f1𝔉),

we have,

(PiP0)eq

=

1Πeqq2P0=1[1+25q3(f1𝔉)]1.

This is exactly the pressure-ratio expression presented in our "first view" and unveils the notation association,

bξq2

1g2.

From our separate derivation, we have, in equilibrium,

𝔊core=(2nc3)Score

=

(2nc3)(4π5)Req3q5(5Pi2q2+Π)eq

 

=

(q5nc5)Req3(23π3)Πeq[52q2(PiΠ)eq+1]

 

=

(nc5)[Rnorm3Pnorm]χeq1(ν2q)[52q2(PiP0)eq(P0Π)eq+1]

[𝔊coreEnorm]eq

=

(nc5)(ν2q)[52q2(11g2)(q2g2)+1]χeq1

 

=

(nc2)(ν2q)[g235]{12nc(4π3)(νq3)nc1(qg)2nc}1/(nc3)

 

=

nc[1(35)1g2](12)(ν2q)g2{2nc(34π)(νq3)1nc(qg)2nc}1/(nc3)

 

=

nc[1(35)1g2]{2nc2(3nc)(34π)(νq3)1nc(νq3)2(nc3)q5(nc3)q2ncg2ncg2(nc3)}1/(nc3)

 

=

nc[1(35)1g2]{(23π)(νq3)nc5q3nc15g6}1/(nc3).

Finally, switching from the g notation to the bξ notation gives,

[𝔊coreEnorm]eq

=

nc[1(35)bξq2]{(23π)(νq3)nc5q3nc15bξ3q6}1/(nc3)

 

=

ncq3[1(35)bξq2]{(23π)(νq3)nc5bξ3}1/(nc3),

which, after setting X=1, precisely matches the above, "first view" expression. Also from our previous derivation, we can write,

𝔊env=(2ne3)Senv

=

2π(2ne3)Req3Πeq{(1q3)(PiΠ)eq+(ρeρ0)[(2q2+3q3q5)+35(ρeρ0)(1+5q25q3+q5)]}

 

=

2π(2ne3)Req3[Pnormχ4(323π)(νq3)2]eq{(1q3)q2(g21)+(25)q5𝔉}

 

=

[PnormRnorm3]ne2(ν2q4)(1q3){(g21)+25(q31q3)𝔉}χeq1

[𝔊envEnorm]eq

=

ne(1q3){(g21)+25(q31q3)𝔉}q22(νq3)2[12nc(4π3)(νq3)nc1(qg)2nc]1/(nc3)

 

=

ne(1q3){(g21)+25(q31q3)𝔉}[2[nc(nc3)](34π)(νq3)(1nc)+2(nc3)q2(nc3)2ncg2nc]1/(nc3)

 

=

ne(1q3){(g21)+25(q31q3)𝔉}[(23π)(νq3)nc5q6g2nc]1/(nc3).

And, finally, switching from the g notation to the bξ notation gives,

[𝔊envEnorm]eq

=

ne(1q3)(bξq2)1{1[125(q31q3)𝔉]bξq2}[(23π)(νq3)nc5q6(bξq2)nc]1/(nc3)

 

=

ne(1q3){1[125(q31q3)𝔉]bξq2}[(23π)(νq3)nc5q62(nc3)+2ncbξ3nc+nc]1/(nc3)

 

=

ne[(23π)(νq3)nc5bξ3]1/(nc3)(1q3){1[125(q31q3)𝔉]bξq2},

which, after setting X=1, precisely matches the above, "first view" expression.


Summary00[edit]

In summary, the desired out of equilibrium free-energy expression is,

𝔊Enorm

=

A0X3/nc+B0X3/neC0X1

where,

A0(𝔖coreEnorm)eq

=

ncbξ[(23π)(νq3)(nc5)bξnc]1/(nc3)q3[1(35)bξq2],

B0(𝔖envEnorm)eq

=

nebξ[(23π)(νq3)(nc5)bξnc]1/(nc3)(1q3){1[125(q31q3)𝔉]bξq2},

C0(WgravEnorm)eq

=

(65)q5f[(23π)(νq3)nc5bξnc]1/(nc3).

Or, in a more compact form,

𝔊*[(23π)(νq3)(nc5)bξnc]1/(nc3)[𝔊Enorm]

=

ncA1X3/nc+neB1X3/ne3C1X1

where,

A1

1bξ(q3)[1(35)bξq2],

B1

1bξ(1q3){1[125(q31q3)𝔉]bξq2},

C1

(25)q5f.

Let's examine the behavior of the first radial derivative.

𝔊*X

=

3X[A1X3/ncB1X3/ne+C1X1].

Let's see whether the sum of terms inside the square brackets is zero at the derived equilibrium radius, that is, when X=1 and, hence, when

χ=χeq

=

[12nc(4π3)(νq3)nc1(qg)2nc]1/(nc3)

 

=

[12nc(4π3)(νq3)nc1bξnc]1/(nc3).

C1A1B1

=

(25)q5f1bξ(q3)[1(35)bξq2]1bξ(1q3){1[125(q31q3)𝔉]bξq2}

 

=

(25)q5f1bξ{1[125(q31q3)𝔉]bξq2}+q3bξ{1[125(q31q3)𝔉]bξq2}q3bξ[1(35)bξq2]

 

=

(25)q5f1bξ+[125(q31q3)𝔉]q2+q3bξ[125(q31q3)𝔉]q5q3bξ+(35)q5

 

=

q2{(25)q3f1bξq2+[125(q31q3)𝔉](1q3)+(35)q3}

 

=

q2{(25)q3f[1+25q3(f1𝔉)]+[(1q3)25q3𝔉]+(35)q3}

 

=

q2{0}.

Q.E.D.

Even slightly better:

1q2[(π23)(νq3)(5nc)bξnc]1/(nc3)[𝔊Enorm]

=

ncA2X3/nc+neB2X3/ne3C2X1,

or, better yet,

Out-of-Equilibrium, Free-Energy Expression for BiPolytropes with Structural (nc,ne)=(0,0)

2(q2ν)2χeq[𝔊Enorm]

=

ncA2X3/nc+neB2X3/ne3C2X1

where, keeping in mind that,

1(bξq2)

=

[1+25q3(f1𝔉)],

we have,

A2A1q2

q3(bξq2)[1(35)bξq2]

 

=

q3{[1+25q3(f1𝔉)](35)}

 

=

25q3[1+q3(f1𝔉)],

B2B1q2

1(bξq2)(1q3){1[125(q31q3)𝔉]bξq2}

 

=

(1q3){1(bξq2)1+25(q31q3)𝔉}

 

=

(1q3){[1+25q3(f1𝔉)]1+25(q31q3)𝔉}

 

=

25q3{(1q3)(f1𝔉)+𝔉}

 

=

25q3{f[1+q3(f1𝔉)]}

 

=

25q3fA2,

C2C1q2

25q3f.

As before, the equilibrium system is dynamically unstable if 2𝔊/X2<0. We have deduced that the system is unstable if,

ne3[3nencne]

<

A2C2=1f[1+q3(f1𝔉)].

See Also[edit]

In October 2023, this very long chapter was subdivided in order to more effectively accommodate edits. Here is a list of the resulting set of shorter chapters:

  1. Free-Energy Synopsis
  2. Free-Energy of Truncated Polytropes
  3. Free-Energy of BiPolytropes


Tiled Menu

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