SSC/FreeEnergy/PolytropesEmbedded/Pt3B

From jetwiki
Jump to navigation Jump to search

Background

Index to original, very long chapter

Free-Energy of Bipolytropes

In this case, the Gibbs-like free energy is given by the sum of four separate energies,

𝔊

=

[Wgrav+𝔖therm]core+[Wgrav+𝔖therm]env.

In addition to specifying (generally) separate polytropic indexes for the core, nc, and envelope, ne, and an envelope-to-core mean molecular weight ratio, μe/μc, we will assume that the system is fully defined via specification of the following five physical parameters:

  • Total mass, Mtot;
  • Total radius, R;
  • Interface radius, Ri, and associated dimensionless interface marker, qRi/R;
  • 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

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


From Detailed Force-Balance Models

Equilibrium Radius

First View

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

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

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

First View

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

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

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

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 |