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(f−1−𝔉)]

where,

f

≡

1+52(ρeρc)(1q2−1)+(ρeρc)2[1q5−1+52(1−1q2)],

𝔉

≡

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

ρeρc

=

q3(1−ν)ν(1−q3).

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)/nc⇒P0Pnorm=[34π(νq3)1χeq3](nc+1)/nc

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

=

(323π)(νq3)2[q2+25q5(f−1−𝔉)]

⇒Req(nc−3)/nc

=

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

⇒χeq(nc−3)/nc≡[ReqRnorm](nc−3)/nc

=

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

Or, in terms of γc,

χeq4−3γc

=

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

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)(1−q)+ρeρ0(1q2−1)].

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(f−1−𝔉)]

=

(323π)ν2g2q4

⇒f−1−𝔉

=

52q3[g2−1]

 

=

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

 

=

52q5(ρeρ0){2(q2−q3)+ρeρ0[1−3q2+2q3]}.

Let's see …

f−1−𝔉

=

52(ρeρc)(1q2−1)+(ρeρc)2[1q5−1+52(1−1q2)]−52(ρeρc)1q5[(−2q2+3q3−q5)+35(ρeρc)(−1+5q2−5q3+q5)]

 

=

52(ρeρc)(1q2−1)−52(ρeρc)1q5[(−2q2+3q3−q5)]

 

 

−52(ρeρc)1q5[35(ρeρc)(−1+5q2−5q3+q5)]+(ρeρc)2[1q5−1+52(1−1q2)]

 

=

52(ρeρc)1q5{(q3−q5)+(2q2−3q3+q5)}

 

 

+12(ρeρc)21q5[3(1−5q2+5q3−q5)]+12(ρeρc)21q5[2−2q5+5(q5−q3)]

 

=

52(ρeρc)1q5[(q3−q5)+(2q2−3q3+q5)]+12(ρeρc)21q5[3(1−5q2+5q3−q5)+2−2q5+5(q5−q3)]

 

=

52(ρeρc)1q5[2q2−2q3]+52q5(ρeρc)2[1−3q2+2q3].

Q.E.D.

Hence, the equilibrium radius can also be written as,

χeq4−3γc

=

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

or, in terms of the polytropic index,

χeqnc−3

=

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

Gravitational Potential Energy[edit]

Also from our accompanying discussion, we have,

WgravEnorm

=

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

 

=

−X−1(65)q5f[2nc−(nc−3)(34π)(νq3)(1−nc)+2(nc−3)bξnc]1/(nc−3)

 

=

−X−1(65)q5f[(6π)(νq3)nc−5bξnc]1/(nc−3).

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χ3−3γ

=

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

 

=

{(PiP0)[34π(νq3)]γcχ3−3γc}eqX3−3γ,

where, in equilibrium,

(PiP0)eq

=

1−bξq2

bξq2

=

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

 

=

[1+25q3(f−1−𝔉)]−1

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

(𝔖AEnorm)core

=

4π/3(γc−1){(PiP0)[34π(νq3)]γcχ3−3γc}eqX3−3γc{(P0Pic)[q3−(3bξ5)q5]}

 

=

1(γc−1)[(4π3)1−γc(νq3)γcχeq3−3γc]X3−3γcq3[1−(35)bξq2],

(𝔖AEnorm)env

=

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

 

=

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

 

=

1(γe−1)[(4π3)1−γc(νq3)γcχeq3−3γc]X3−3γe{(1−bξq2)(1−q3)+bξ[25q5𝔉]}

 

=

1(γe−1)[(4π3)1−γc(νq3)γcχeq3−3γc]X3−3γe(1−q3){1−[1−25(q31−q3)𝔉]bξq2}.

Furthermore,

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

=

(34π)γc−1(νq3)γc{χeq4−3γc}(3−3γc)/(4−3γc)

 

=

(34π)γc−1(νq3)γc{12(34π)1−γc(νq3)2−γc[q2+25q5(f−1−𝔉)]}(3−3γc)/(4−3γc)

 

=

(34π)(γc−1)/(4−3γc)(νq3)(6−5γc)(4−3γc){q22[1+25q3(f−1−𝔉)]}(3−3γc)/(4−3γc)

 

=

(34π)1/(nc−3)(νq3)(nc−5)(nc−3){q22[1+25q3(f−1−𝔉)]}−3/(nc−3)

 

=

[(2⋅3π)(νq3)(nc−5)bξ3]1/(nc−3).

Hence, we have,

(𝔖AEnorm)core

=

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

 

=

nc[(2⋅3π)(νq3)(nc−5)bξ3]1/(nc−3)X−3/ncq3[1−(35)bξq2],

(𝔖AEnorm)env

=

ne[(2⋅3π)(νq3)(nc−5)bξ3]1/(nc−3)X−3/ne(1−q3){1−[1−25(q31−q3)𝔉]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=1−1g2,

or given that (see above),

52q3[g2−1]

=

f−1−𝔉

⇒g2

=

1+25q3(f−1−𝔉),

we have,

(PiP0)eq

=

1−Πeqq2P0=1−[1+25q3(f−1−𝔉)]−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]χeq−1(ν2q)[52q2(PiP0)eq(P0Π)eq+1]

⇒[𝔊coreEnorm]eq

=

(nc5)(ν2q)[52q2(1−1g2)(q2g2)+1]χeq−1

 

=

(nc2)(ν2q)[g2−35]{12nc(4π3)(νq3)nc−1(qg)2nc}−1/(nc−3)

 

=

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

 

=

nc[1−(35)1g2]{2nc⋅2(3−nc)(34π)(νq3)1−nc(νq3)2(nc−3)q5(nc−3)q−2ncg−2ncg2(nc−3)}1/(nc−3)

 

=

nc[1−(35)1g2]{(2⋅3π)(νq3)nc−5q3nc−15g−6}1/(nc−3).

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

[𝔊coreEnorm]eq

=

nc[1−(35)bξq2]{(2⋅3π)(νq3)nc−5q3nc−15bξ3q6}1/(nc−3)

 

=

ncq3[1−(35)bξq2]{(2⋅3π)(νq3)nc−5bξ3}1/(nc−3),

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{(1−q3)(PiΠ)eq+(ρeρ0)[(−2q2+3q3−q5)+35(ρeρ0)(−1+5q2−5q3+q5)]}

 

=

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

 

=

[PnormRnorm3]ne2(ν2q4)(1−q3){(g2−1)+25(q31−q3)𝔉}χeq−1

⇒[𝔊envEnorm]eq

=

ne(1−q3){(g2−1)+25(q31−q3)𝔉}q22(νq3)2[12nc(4π3)(νq3)nc−1(qg)2nc]−1/(nc−3)

 

=

ne(1−q3){(g2−1)+25(q31−q3)𝔉}[2[nc−(nc−3)](34π)(νq3)(1−nc)+2(nc−3)q2(nc−3)−2ncg−2nc]1/(nc−3)

 

=

ne(1−q3){(g2−1)+25(q31−q3)𝔉}[(2⋅3π)(νq3)nc−5q−6g−2nc]1/(nc−3).

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

[𝔊envEnorm]eq

=

ne(1−q3)(bξq2)−1{1−[1−25(q31−q3)𝔉]bξq2}[(2⋅3π)(νq3)nc−5q−6(bξq2)nc]1/(nc−3)

 

=

ne(1−q3){1−[1−25(q31−q3)𝔉]bξq2}[(2⋅3π)(νq3)nc−5q−6−2(nc−3)+2ncbξ3−nc+nc]1/(nc−3)

 

=

ne[(2⋅3π)(νq3)nc−5bξ3]1/(nc−3)(1−q3){1−[1−25(q31−q3)𝔉]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

=

A0X−3/nc+B0X−3/ne−C0X−1

where,

A0≡(𝔖coreEnorm)eq

=

ncbξ[(2⋅3π)(νq3)(nc−5)bξnc]1/(nc−3)q3[1−(35)bξq2],

B0≡(𝔖envEnorm)eq

=

nebξ[(2⋅3π)(νq3)(nc−5)bξnc]1/(nc−3)(1−q3){1−[1−25(q31−q3)𝔉]bξq2},

C0≡(WgravEnorm)eq

=

(65)q5f[(2⋅3π)(νq3)nc−5bξnc]1/(nc−3).

Or, in a more compact form,

𝔊*≡[(2⋅3π)(νq3)(nc−5)bξnc]−1/(nc−3)[𝔊Enorm]

=

ncA1X−3/nc+neB1X−3/ne−3C1X−1

where,

A1

≡

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

B1

≡

1bξ(1−q3){1−[1−25(q31−q3)𝔉]bξq2},

C1

≡

(25)q5f.

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

∂𝔊*∂X

=

3X[−A1X−3/nc−B1X−3/ne+C1X−1].

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)nc−1(qg)2nc]1/(nc−3)

 

=

[12nc(4π3)(νq3)nc−1bξ−nc]1/(nc−3).

C1−A1−B1

=

(25)q5f−1bξ(q3)[1−(35)bξq2]−1bξ(1−q3){1−[1−25(q31−q3)𝔉]bξq2}

 

=

(25)q5f−1bξ{1−[1−25(q31−q3)𝔉]bξq2}+q3bξ{1−[1−25(q31−q3)𝔉]bξq2}−q3bξ[1−(35)bξq2]

 

=

(25)q5f−1bξ+[1−25(q31−q3)𝔉]q2+q3bξ−[1−25(q31−q3)𝔉]q5−q3bξ+(35)q5

 

=

q2{(25)q3f−1bξq2+[1−25(q31−q3)𝔉](1−q3)+(35)q3}

 

=

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

 

=

q2{0}.

Q.E.D.

Even slightly better:

1q2[(π2⋅3)(νq3)(5−nc)bξ−nc]1/(nc−3)[𝔊Enorm]

=

ncA2X−3/nc+neB2X−3/ne−3C2X−1,

or, better yet,

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

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

=

ncA2X−3/nc+neB2X−3/ne−3C2X−1

where, keeping in mind that,

1(bξq2)

=

[1+25q3(f−1−𝔉)],

we have,

A2≡A1q2

≡

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

 

=

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

 

=

25q3[1+q3(f−1−𝔉)],

B2≡B1q2

≡

1(bξq2)(1−q3){1−[1−25(q31−q3)𝔉]bξq2}

 

=

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

 

=

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

 

=

25q3{(1−q3)(f−1−𝔉)+𝔉}

 

=

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

 

=

25q3f−A2,

C2≡C1q2

≡

25q3f.

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

ne3[3−nenc−ne]

<

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

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 |