Appendix/Ramblings/BdHN

From jetwiki
Jump to navigation Jump to search


Binary-driven Hpernovae

The material presented here builds on our separate discussion of close binary stars.

Setup

Initial Carbon-Oxygen Core

The total angular momentum associated with a Maclaurin spheroid of eccentricity, e, is given by the expression,

L*2L2(GM3a¯)

=

652[(32e2)(1e2)1/2sin1ee33(1e2)e2](1e2)2/3.

📚 Marcus, Press, & Teukolsky (1977), §IVa, p. 591, Eq. (4.2)

According to an accompanying table

Jacobi Sequence Bifurcation:

  • [e,ω02/(4πGρ),j2]=[0.812670,0.0935574,0.00455473][e,ω02/(πGρ),L*]=[0.812670,0.374230,0.303751]

Dynamically Unstable 2nd Harmonic:

  • [e,ω02/(4πGρ),j2]=[0.95289,0.11006,0.012796][e,ω02/(πGρ),L*]=[0.95289,0.44024,0.25921]

Now, given that, ω02=(2π/P)2, if we specify P, then we know ρCO from the expression,

ρCO

=

(ω02πGρ)1(2π/P)2πG;

and if we specify MCO, we know a¯CO via the expression,

a¯CO

=

(3MCO4πρCO)1/3.

Initial CO Core at the Jacobi Bifurcation Point
with MCO=10M and various spin periods, P
P
(min)
P
(s)
ρCO=factor1/P2 a¯CO=factor2ρCO1/3 LCO=factor3a¯1/2
factor1 (g cm-3) factor2 (cm) factor3 (g cm2 s-1)
100 6600 1.155 × 101 7.435 × 1010 6.001 × 1052
5 300 5.592 × 103 9.469 × 109 2.142 × 1052
  • G=6.672×108 cm3 g-1 s-2
  • M=1.989×1033 g
  • factor1=(0.37423)14πG=5.033×108 cgs
  • factor2=1.681×1011 cgs
  • factor3=(0.303751)(GMCO3)1/2=2.201×1047 cgs

Related Binary

Consider the simple model of two spherical stars in circular orbit about one another, as depicted here on the right. In addition to the physical parameters explicitly labeled in this diagram, we adopt the following variable notation:

  • The total system mass is,

MtotM+M';

  • The ratio of the primary to secondary mass is,

λMM';

  • And the separation between the two centers is,

drcm+rcm'.

Simple Binary Model

For a circular orbit, the angular velocity is related to the the system mass and separation via the Kepler relation,

ω2d3=GMtot,

and the distances, rcm and rcm', between the center of each star and the center of mass (cm) of the system must be related to one another via the expression,

rcm'rcm

=

MM'=λ.

Note that the following relations also hold:

M=Mtot(λ1+λ)

      and      

M'=Mtot(11+λ)

rcm=d(11+λ);

      and      

rcm'=d(λ1+λ).

Hence, the orbital angular momentum is,

Lorb

=

[Mrcm2+M'(rcm')2]ω

 

=

Mtotd2[(λ1+λ)(11+λ)2+(11+λ)(λ1+λ)2][GMtotd3]1/2

 

=

(GMtot3d)1/2[λ(1+λ)2].

Assuming that both stars are rotating synchronously with the orbit, their respective spin angular momenta are,

LM=IMω

=

25MR2ω

 

=

25Mtot(λ1+λ)R2[GMtotd3]1/2

 

=

25(GMtot3d)1/2(λ1+λ)(Rd)2,

LM'=IM'ω

=

25M'(R')2ω

 

=

25Mtot(11+λ)(R')2[GMtotd3]1/2

 

=

25(GMtot3d)1/2(11+λ)(R'd)2.

Hence, the total angular momentum of the system is,

Ltot=Lorb+LM+LM'

=

(GMtot3d)1/2[λ(1+λ)2]+25(GMtot3d)1/2(λ1+λ)(Rd)2

 

 

+25(GMtot3d)1/2(11+λ)(R'd)2

 

=

(GMtot3R)1/2{[λ(1+λ)2](dR)1/2+25(11+λ)[λ+(R'R)2](dR)3/2}

Solutions

Here we assume that the unprimed parameters (primary star) refer to the Maclaurin spheroid component while the primed parameters (secondary star) refer to the Jacobi ellipsoid component. Also assume that the mass ratio is,

λMM'

=

MCOfMJTE=8.51.5=173=5.6667.

Assume that the total angular momentum and total mass are conserved, and that the three frequencies that appear in the binary system — orbital, and two stellar spins — are identical (co-rotation). This means that the orbital separation, d, can be expressed in terms of ωf via the Keperian relation,

d3

=

GMtotωf2.

And the radii of the binary components are given, respectively, by the relations,

R3

=

(c2)primary[λGMtot(1+λ)ωf2];

(R')3

=

(c2)secondary[GMtot(1+λ)ωf2].

Hence, the length-ratios that appear in the expression for Ltot may be re-expressed as follows:

(R'R)3

=

(c2)secondary[GMtot(1+λ)ωf2]{(c2)primary[λGMtot(1+λ)ωf2]}1

 

=

(c2)secondaryλ(c2)primary;

(dR)3

=

GMtotωf2{(c2)primary[λGMtot(1+λ)ωf2]}1

 

=

(1+λ)λ(c2)primary.

Also,

Ltot6(GMtot3R)3=L*6(GMtot3a¯CO3)(GMtot3R)3

=

L*6a¯CO3{(c2)primary[λGMtot(1+λ)ωf2]}1

 

=

(1+λ)L*6λ(c2)primary[a¯CO3ωf2GMtot].

As a result, we find,

(1+λ)L*6λ(c2)primary[a¯CO3ωf2GMtot]

=

{[λ(1+λ)2](dR)1/2+25(11+λ)[λ+(R'R)2](dR)3/2}6

[a¯CO3ωf2GMtot]

=

λ(c2)primary(1+λ)L*6{[λ(1+λ)2](dR)1/2+25(11+λ)[λ+(R'R)2](dR)3/2}6

Critique

  1. Ill-advised to refer to the new NS as "νNS" because, in this context, readers might reasonably associate the greek letter, ν, with neutrinos.

See Also

Tiled Menu

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