SSC/Stability/NeutralMode: Difference between revisions

From jetwiki
Jump to navigation Jump to search
Line 2,164: Line 2,164:
</tr>
</tr>
</table>
</table>
We conclude, therefore, that the radial displacement function (''i.e.,'' the eigenfunction) for the neutral mode of all polytropic configurations is,
We conclude, therefore, that the radial displacement function (''i.e.,'' the eigenfunction) for the neutral <math>(\sigma_c^2 = 0)</math> mode of all polytropic configurations is,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">



Revision as of 17:36, 8 September 2021

LAWE

Most General Form

In an accompanying discussion, we derived the so-called,

Adiabatic Wave (or Radial Pulsation) Equation

d2xdr02+[4r0−(g0ρ0P0)]dxdr0+(ρ0γgP0)[ω2+(4−3γg)g0r0]x=0

where the gravitational acceleration,

g0 ≡

GMrr02=−1ρ0dP0dr0⇒g0ρ0r0P0=−dln⁡P0dln⁡r0.

The solution to this equation gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, self-gravitating fluid configurations. The boundary condition conventionally used in connection with the adiabatic wave equation is,

r0dln⁡xdr0

=

1γg(4−3γg+ω2R3GMtot)        at         r0=R.

Polytropic Configurations

Part 1

If the initial, unperturbed equilibrium configuration is a polytropic sphere whose internal structure is defined by the function, θ(ξ), that provides a solution to the,

Lane-Emden Equation

1ξ2ddξ(ξ2dΘHdξ)=−ΘHn

then,

r0

=

anξ,

ρ0

=

ρcθn,

P0

=

Kρ0(n+1)/n=Kρc(n+1)/nθn+1,

g0

=

GM(r0)r02=Gr02[4πan3ρc(−ξ2dθdξ)],

where,

an

=

[(n+1)K4πG⋅ρc(1−n)/n]1/2.

Hence, after multiplying through by an2, the above adiabatic wave equation can be rewritten in the form,

d2xdξ2+[4ξ−g0an(an2ρ0P0)]dxdξ+(an2ρ0γgP0)[ω2+(4−3γg)g0anξ]x

=

0.

In addition, given that,

g0an

=

4πGρc(−dθdξ),

and,

an2ρ0P0

=

(n+1)(4πGρc)θ=an2ρcPc⋅θcθ,

we can write,

0

=

d2xdξ2+[4−(n+1)V(ξ)ξ]dxdξ+[ω2(an2ρcγgPc)θcθ−(3−4γg)⋅(n+1)V(x)ξ2]x

 

=

d2xdξ2+[4−(n+1)V(ξ)ξ]dxdξ+(n+1)[ω2(an2ρcγgPc)ξ2θc(n+1)θ−(3−4γg)⋅V(x)]xξ2

where we have adopted the function notation,

V(ξ)

≡

−ξθdθdξ.

Part 2

Drawing from an accompanying discussion, we have the following:

Polytropic LAWE (linear adiabatic wave equation)

0=d2xdξ2+[4−(n+1)Q]1ξ⋅dxdξ+(n+1)[(σc26γg)ξ2θ−αQ]xξ2

where:    Q(ξ)≡−dln⁡θdln⁡ξ,    σc2≡3ω22πGρc,     and,     α≡(3−4γg)

In order to reconcile with the "Part 1" expression, we note first that V(ξ)↔Q(ξ). We note as well that since,

(an2ρcPc)θc

=

(n+1)4πGρc,

we have,

ω2(an2ρcγgPc)ξ2θc(n+1)θ

↔

ω2γg[(n+1)4πGρc]ξ2(n+1)θ=16γg[3ω22πGρc]ξ2θ=(σc26γg)ξ2θ.

All physically reasonable solutions are subject to the inner boundary condition,

dxdξ=0         at         ξ=0,

but the relevant outer boundary condition depends on whether the underlying equilibrium configuration is isolated (surface pressure is zero), or whether it is a "pressure-truncated" configuration. As is the case with the pressure-truncated isothermal spheres, discussed above, if the polytropic configuration is truncated by the pressure, Pe, of a hot, tenuous external medium, then the solution to the LAWE is subject to the outer boundary condition,

−dln⁡xdln⁡ξ=3         at         ξ=ξ~.

But, for isolated polytropes, the sought-after solution is subject to the more conventional boundary condition,

−dln⁡xdln⁡ξ=(3−nn+1)+nσc26(n+1)[ξθ']         at         ξ=ξsurf.

Radial Pulsation Neutral Mode

Background

The integro-differential version of the statement of hydrostatic balance is

dPdr=−GMrρr2

From our separate discussion, we have found that,

Exact Solution to the (3≤n<∞) Polytropic LAWE

σc2=0

      and      

xP≡3(n−1)2n[1+(n−3n−1)(1ξθn)dθdξ].

Let's rewrite the significant functional term in this expressions in terms of basic variables. That is,

(1ξθn)dθdξ

=

−(anρcr0ρ0)g04πGρcan

 

=

−M(r0)4πr03ρ0.

Trial Eigenfunction & Its Derivatives

Let's adopt the following trial solution:

xt

=

a−bMr4πr03ρ0=a−bg04πGr0ρ0.

Then we have,

−(1b)dxtdr0

=

ddr0[Mr4πr03ρ0]

 

=

[14πr03ρ0]dMrdr0−[Mr4πr03ρ02]dρ0dr0−[3Mr4πr04ρ0]

 

=

1r0−[3Mr4πr04ρ0]−[Mr4πr03ρ02]dρ0dr0

−(1b)d2xtdr02

=

ddr0{1r0−[3Mr4πr04ρ0]−[Mr4πr03ρ02]dρ0dr0}

 

=

−1r02−[34πr04ρ0]dMrdr0+[3Mr4πr04ρ02]dρ0dr0+4[3Mr4πr05ρ0]

 

 

−[Mr4πr03ρ02]d2ρ0dr02−[14πr03ρ02]dρ0dr0⋅dMrdr0+[3Mr4πr04ρ02]dρ0dr0+[2Mr4πr03ρ03](dρ0dr0)2

 

=

−4r02+3Mr4πr05ρ0[4+dln⁡ρ0dln⁡r0]+1r02[3Mr4πr03ρ0−1]dln⁡ρ0dln⁡r0−[Mr4πr03ρ02]d2ρ0dr02+[2Mr4πr05ρ0](dln⁡ρ0dln⁡r0)2.

Given that,

Δ≡Mr4πr03ρ0

=

14πG(g0r0ρ0)=−[P04πGr02ρ02⋅dln⁡P0dln⁡r0],

these expression can be rewritten as,

−(r02b)dxtdr0

=

r0{1−3Δ−Δ⋅dln⁡ρ0dln⁡r0},

and,

−(r02b)d2xtdr02

=

−4+3Δ[4+dln⁡ρ0dln⁡r0]+[3Δ−1]dln⁡ρ0dln⁡r0−Δ(r02ρ0)d2ρ0dr02+2Δ(dln⁡ρ0dln⁡r0)2.

Plug Trial Eigenfunction Into LAWE

 

LAWE

d2xdr02+[4r0−(g0ρ0P0)]dxdr0+(ρ0γgP0)[ω2+(4−3γg)g0r0]x=0


Plugging our trial radial displacement function, xt, into the LAWE gives,

LAWE

=

−(r02b)d2xtdr02−(r0b)[4+dln⁡P0dln⁡r0]dxtdr0−(r02b)(ρ0γgP0)[(4−3γg)g0r0+σc2(2πGρc3)]xt

 

=

−(r02b)d2xtdr02−(r0b)[4+dln⁡P0dln⁡r0]dxtdr0+(1b)1γg(dln⁡P0dln⁡r0)(4−3γg)xt−(r02b)(ρ0γgP0)[σc2(2πGρc3)]xt

 

=

−4+3Δ[4+dln⁡ρ0dln⁡r0]+[3Δ−1]dln⁡ρ0dln⁡r0−Δ(r02ρ0)d2ρ0dr02+2Δ(dln⁡ρ0dln⁡r0)2

 

 

+[4+dln⁡P0dln⁡r0]{1−3Δ−Δ⋅dln⁡ρ0dln⁡r0}+(1b)1γg(dln⁡P0dln⁡r0)(4−3γg)(a−bΔ)−(1b)(ρ0r02γgP0)[σc2(2πGρc3)](a−bΔ).

Now, if we set σc2=0 and dln⁡P0/dln⁡r0=γg(dln⁡ρ0/dln⁡r0), this expression becomes,

LAWE

=

−4+3Δ[4+dln⁡ρ0dln⁡r0]+[3Δ−1]dln⁡ρ0dln⁡r0−Δ(r02ρ0)d2ρ0dr02+2Δ(dln⁡ρ0dln⁡r0)2

 

 

+[4+dln⁡P0dln⁡r0]{1−3Δ−Δ⋅dln⁡ρ0dln⁡r0}+(1b)1γg(dln⁡P0dln⁡r0)(4−3γg)(a−bΔ)−(1b)(ρ0r02γgP0)[σc2(2πGρc3)](a−bΔ).

Notice that the key components of this last term may be rewritten as,

(ρ0r02γgP0)[σc2(2πGρc3)]

=

(4πGρ02r02P0)[σc26γg(ρcρ0)]

 

=

−(1Δ)dln⁡P0dln⁡r0[σc26γg(ρcρ0)].

So, for our trial eigenfunction, we have …

LAWE

=

[2Δ−1]dln⁡ρ0dln⁡r0−Δ(r02ρ0)d2ρ0dr02+2Δ(dln⁡ρ0dln⁡r0)2+dln⁡P0dln⁡r0⋅{[1−3Δ−Δ⋅dln⁡ρ0dln⁡r0]+(4−3γg)γg[ab−Δ]+σc26γg(ρcρ0)[abΔ−1]}.

Consider Polytropic Structures

Referring back to, for example, a separate review of polytropic structures, we recognize that,

Δ=14πG(g0r0ρ0)

=

1ξ3[(−ξ2dθdξ)]θ−n=1ξ(−dθdξ)θ−n=−θ'ξθn,

dln⁡ρ0dln⁡r0

=

n,

dln⁡P0dln⁡r0

=

(n+1).

Also,

(r02ρ0)d2ρ0dr02

=

(ξ2ρcθn)ddξ[nρcθn−1θ']

 

=

(nξ2θn)[(n−1)θn−2(θ')2+θn−1θ'′]

 

=

(nξ2θ2)[(n−1)(θ')2+θ⋅θ'′]

 

=

(nξ2θ2)[(n−1)(θ')2−(θn+1+2θθ'ξ)]

 

=

(nξ2θ2)[(n−1)(ξθnΔ)2+θn+1]

 

=

n(n−1)(ξn+1θn−1Δ)2+nξ2θn−1.

Hence,

LAWE

=

n(2Δ−1)−Δ(r02ρ0)d2ρ0dr02+2n2Δ+(n+1)[1−3Δ−nΔ]+(n+1){(4−3γg)γg[ab−Δ]+σc26γg(ρcρ0)[abΔ−1]}

 

=

2nΔ−n+2n2Δ+n−3nΔ−n2Δ+1−3Δ−nΔ−Δ(r02ρ0)d2ρ0dr02+(n+1){(4−3γg)γg[ab−Δ]+σc26γg(ρcρ0)[abΔ−1]}

 

=

1+n2Δ−(2n+3)Δ−Δ(r02ρ0)d2ρ0dr02+(n+1)[4γg−3][ab−Δ]+(n+1){σc26γg(ρcρ0)[abΔ−1]}.

If, γg=(n+1)/n, we can further simplify and obtain,

LAWE

=

1+n2Δ−(2n+3)Δ−Δ(r02ρ0)d2ρ0dr02+[n−3][ab−Δ]+(n+1){σc26γg(ρcρ0)[abΔ−1]}

 

=

1+(n−3)ab+n2Δ+(3−n)Δ−(2n+3)Δ−Δ(r02ρ0)d2ρ0dr02+(n+1){σc26γg(ρcρ0)[abΔ−1]}

Try Again

General Form of Wave Equation

 

LAWE

d2xdr02+[4r0−(g0ρ0P0)]dxdr0+(ρ0γgP0)[ω2+(4−3γg)g0r0]x=0

Employing the substitutions,

σc2 ≡

3ω22πGρc,

α ≡

3−4γg=3−nn+1,

g0 ≡

GMrr02=−1ρ0dP0dr0⇒g0ρ0r0P0=−dln⁡P0dln⁡r0,

Δ

≡

Mr4πr03ρ0=14πG(g0r0ρ0)=−[P04πGr02ρ02⋅dln⁡P0dln⁡r0],

we have,

LAWE =

d2xdr02+1r0[4−g0ρ0r0P0]dxdr0+[(4γg−3)g0ρ0r0P0]xr02+(ρ0P0)[4πGρc(σc26γg)]x

  =

d2xdr02+1r0[4+dln⁡P0dln⁡r0]dxdr0+[α⋅dln⁡P0dln⁡r0]xr02−1Δ[dln⁡P0dln⁡r0⋅ρcρ0(σc26γg)]xr02.


In the context of polytropic configurations (see more below), we appreciate that,

ρ0ρc =

θn,

dln⁡P0dln⁡r0 =

(n+1)dln⁡θdln⁡ξ=−(n+1)Q,      and,

1Δ⋅dln⁡P0dln⁡r0 =

−(n+1)ξ2θn−1.

Inserting these into the LAWE expression and multiplying through by the square of the polytropic length scale, an2, we obtain,

LAWE =

d2xdξ2+1ξ[4−(n+1)Q]dxdξ−[α(n+1)Q]xξ2+[(n+1)ξ2θ(σc26γg)]xξ2.


This is identical to what has been referred to in a separate discussion, as the

Polytropic LAWE

0=d2xdξ2+[4−(n+1)Q]1ξ⋅dxdξ+(n+1)[(σc26γg)ξ2θ−αQ]xξ2

where:    Q(ξ)≡−dln⁡θdln⁡ξ,    σc2≡3ω22πGρc,     and,     α≡(3−4γg)

Derivatives of Δ

Here we evaluate the first derivative of Δ with respect to r0,

dΔdr0 =

ddr0{Mr4πr03ρ0}

  =

14πr03ρ0⋅dMrdr0−3Mr4πr04ρ0−Mr4πr03ρ02⋅dρ0dr0

  =

1r0−3Mr4πr04ρ0−Mr4πr03ρ02⋅dρ0dr0

  =

1r0{1−Δ[3+dln⁡ρ0dln⁡r0]}

⇒r0⋅dΔdr0 =

[1−3Δ−Δ⋅dln⁡ρ0dln⁡r0];

and the second derivative of Δ with respect to r0,

d2Δdr02 =

ddr0{1r0}−ddr0{3Mr4πr04ρ0}−ddr0{Mr4πr03ρ02⋅dρ0dr0}

  =

−1r02−34π{1r04ρ0⋅dMrdr0−4Mrr05ρ0−Mrr04ρ02⋅dρ0dr0}

   

−14π{1r03ρ02⋅dρ0dr0⋅dMrdr0−3Mrr04ρ02⋅dρ0dr0−2Mrr03ρ03⋅[dρ0dr0]2+Mrr03ρ02⋅d2ρ0dr02}

  =

−1r02−3r02{1−Δ[4+dln⁡ρ0dln⁡r0]}−1r02{[1−3Δ]dln⁡ρ0dln⁡r0−2Δ⋅[dln⁡ρ0dln⁡r0]2+Δ(r02ρ0)d2ρ0dr02}

⇒−r02⋅d2Δdr02 =

1+3−3Δ[4+dln⁡ρ0dln⁡r0]+[1−3Δ]dln⁡ρ0dln⁡r0−2Δ⋅[dln⁡ρ0dln⁡r0]2+Δ(r02ρ0)d2ρ0dr02

  =

4−12Δ+[1−6Δ]dln⁡ρ0dln⁡r0−2Δ⋅[dln⁡ρ0dln⁡r0]2+Δ(r02ρ0)d2ρ0dr02.

Trial Eigenfunction

As above, let's adopt a trial eigenfunction of the form,

xt

=

a−bMr4πr03ρ0=a−bΔ.

Then we have,

1b[r02×LAWE]trial =

−r02⋅d2Δdr02−[4+dln⁡P0dln⁡r0]r0⋅dΔdr0+[α⋅dln⁡P0dln⁡r0](ab−Δ)−1Δ[dln⁡P0dln⁡r0⋅ρcρ0(σc26γg)](ab−Δ)

  =

{4−12Δ+[1−6Δ]dln⁡ρ0dln⁡r0−2Δ⋅[dln⁡ρ0dln⁡r0]2+Δ(r02ρ0)d2ρ0dr02}−[4+dln⁡P0dln⁡r0][1−3Δ−Δ⋅dln⁡ρ0dln⁡r0]

   

+[α⋅dln⁡P0dln⁡r0](ab−Δ)−1Δ[dln⁡P0dln⁡r0⋅ρcρ0(σc26γg)](ab−Δ).

Assume Polytropic Relations

If we assume that the equilibrium models are polytropes, then we know that,

ρ0∝θn

      ⇒     

dln⁡ρ0dln⁡r0=n⋅dln⁡θdln⁡ξ;

P0∝θn+1

      ⇒     

dln⁡P0dln⁡r0=(n+1)⋅dln⁡θdln⁡ξ.

We also deduce that,

(r02ρ0)d2ρ0dr02

=

(ξ2θn)ddξ[dθndξ]

 

=

(ξ2θn)ddξ[nθn−1θ']

 

=

(nξ2θn)[(n−1)θn−2(θ')2+θn−1θ'′]

 

=

[n(n−1)ξ2θ2](θ')2−(nξ2θ)[θn+2ξθ']

 

=

n(n−1)⋅[dln⁡θdln⁡ξ]2−nξ2θn−1−2n⋅dln⁡θdln⁡ξ

 

=

n(n−1)⋅Ξ2−nξ2θn−1−2n⋅Ξ,

where we have introduced the shorthand notation,

Ξ

≡

dln⁡θdln⁡ξ.


Drawing from our accompanying discussion, for example, we note as well that,

Δ

=

−[P04πGr02ρ02⋅dln⁡P0dln⁡r0]

⇒−1Δ⋅dln⁡P0dln⁡r0

=

4πGr02ρ02P0

 

=

4πG(an2ξ2)(ρcθn)2[K−1ρc−(n+1)/nθ−(n+1)]

 

=

4πG[(n+1)K4πG⋅ρc(1−n)/n](ρc)2[K−1ρc−(n+1)/n]ξ2θ(n−1)

 

=

(n+1)ξ2θ(n−1)

⇒Δ⋅ξ2θn−1

=

−1(n+1)⋅dln⁡P0dln⁡r0=−Ξ.

Hence,

1b[r02×LAWE]trial =

4−12Δ+[1−6Δ]nΞ−2Δ⋅[nΞ]2+Δ[n(n−1)⋅Ξ2−nξ2θn−1−2n⋅Ξ]−[4+(n+1)Ξ][1−3Δ−Δ⋅nΞ]

   

+[α⋅(n+1)Ξ](ab−Δ)−1Δ[(n+1)Ξ⋅ρcρ0(σc26γg)](ab−Δ)

  =

nΞ−6ΔnΞ−2Δ⋅[nΞ]2+Δ[n(n−1)⋅Ξ2−nξ2θn−1+2n⋅Ξ]+2(n+1)Ξ+n(n+1)Ξ2⋅Δ

   

+[α⋅(n+1)Ξ](ab)−[α⋅(n+1)Ξ]Δ−1Δ[(n+1)Ξ⋅ρcρ0(σc26γg)](ab−Δ)

  =

Ξ[1+n+2(n+1)+α⋅(n+1)(ab)]+Δ[−[4n+α⋅(n+1)]Ξ+n(n−1)⋅Ξ2+n(n+1)Ξ2−2n2Ξ2]

   

−1Δ[(n+1)Ξ⋅ρcρ0(σc26γg)](ab−Δ)

  =

Ξ[3(n+1)+α⋅(n+1)(ab)]−Δ[4n+α⋅(n+1)]Ξ−1Δ[(n+1)Ξ⋅ρcρ0(σc26γg)](ab−Δ)

Third Time

General Relations

Various general relations taken from above derivations:

LAWE =

d2xdr02+1r0[4+dln⁡P0dln⁡r0]dxdr0+[α⋅dln⁡P0dln⁡r0]xr02−1Δ[dln⁡P0dln⁡r0⋅ρcρ0(σc26γg)]xr02,

where,

Δ ≡

Mr4πr03ρ0

⇒1Δ⋅dln⁡P0dln⁡r0 =

−4πGr02ρ02P0;

r0⋅dΔdr0 =

[1−3Δ−Δ⋅dln⁡ρ0dln⁡r0];

−r02⋅d2Δdr02 =

4−12Δ+[1−6Δ]dln⁡ρ0dln⁡r0−2Δ⋅[dln⁡ρ0dln⁡r0]2+Δ(r02ρ0)d2ρ0dr02.

Polytropes

If polytropic relations are adopted:

Δ =

Qξ2θn−1;

1Δ⋅dln⁡P0dln⁡r0 =

−(n+1)ξ2θn−1;

r02×LAWE =

r02d2xdr02+[4−(n+1)Q]r0⋅dxdr0+[(n−3)Q]x+1Δ[(n+1)Q⋅1θn(σc26γg)]x;

r0⋅dΔdr0 =

[1−3Δ+nΔQ]=1+(nQ−3)Δ;

−r02⋅d2Δdr02 =

4−12Δ−n[1−6Δ]Q−2n2ΔQ2+Δ(r02ρ0)d2ρ0dr02

  =

4−12Δ−n[1−6Δ]Q−2n2ΔQ2+Δ[n(n−1)Q2+2nQ−nQΔ]

  =

4−nQ−12Δ+6nQΔ−2n2ΔQ2+Δ[n(n−1)Q2+2nQ]−nQ

  =

4−2nQ+Δ[6nQ−2n2Q2+n(n−1)Q2+2nQ−12].

  =

4−2nQ+Δ[8nQ−n2Q2−nQ2−12].

Eigenfunction Choice

Again, let's try the trial eigenfunction,

xt =

a−bΔ,

in which case,

1b[r02×LAWE] =

−r02d2Δdr02−[4−(n+1)Q]r0⋅dΔdr0+[(n−3)Q](ab−Δ)+1Δ[(n+1)Q⋅1θn(σc26γg)](ab−Δ)

  =

4−2nQ+Δ[8nQ−n2Q2−nQ2−12]−[4−(n+1)Q][1+(nQ−3)Δ]+[(n−3)Q](ab−Δ)

   

+1Δ[(n+1)Q⋅1θn(σc26γg)](ab−Δ)

  =

4−2nQ+Δ[8nQ−n2Q2−nQ2−12]−4−(4nQ−12)Δ+(n+1)Q+(n+1)(nQ2−3Q)Δ+(n−3)Q(ab)−(n−3)QΔ

   

+1Δ[(n+1)Q⋅1θn(σc26γg)](ab−Δ)

  =

[(n+1)−2n+(n−3)(ab)]Q+Δ[8nQ−n2Q2−nQ2−12−4nQ+12+(n+1)(nQ2)−3Q(n+1)+(3−n)Q]

   

+1Δ[(n+1)Q⋅1θn(σc26γg)](ab−Δ)

  =

[(1−n)+(n−3)(ab)]Q+Δ[0]+1Δ[(n+1)Q⋅1θn(σc26γg)](ab−Δ).

Hence, we are left with only the σc2 term if we set,

0 =

(1−n)+(n−3)(ab)

⇒(ab) =

n−1n−3.

We conclude, therefore, that the radial displacement function (i.e., the eigenfunction) for the neutral (σc2=0) mode of all polytropic configurations is,

xneutral =

1b[n−1n−3−Δ]=1b[n−1n−3−Qξ2θn−1]

  =

1b(n−1n−3)[1+(n−3n−1)1ξ2θn−1⋅dln⁡θdln⁡ξ]

  =

1b(n−1n−3)[1+(n−3n−1)1ξθn⋅dθdξ].

This last expression exactly matches our earlier result found for polytropic configurations if we choose an overall amplitude coefficient of the form,

b =

2n3(n−3).

Hooray!

See Also

Tiled Menu

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