Lecture 1 Overview Time Scales, Temperature-density Scalings, Critical Masses
I. Preliminaries The life of any star is a continual struggle between the force of gravity, seeking to reduce the star to a point, and pressure, which holds it up. Stars are gravitationally confined thermonuclear reactors. So long as they remain non-degenerate and have not encountered the pair instability, overheating leads to expansion and cooling. Cooling, on the other hand, leads to contraction and heating. Hence stars are generally stable. The Virial Theorem works. But, since ideal gas pressure depends on temperature, stars must remain hot. By being hot, they are compelled to radiate. In order to replenish the energy lost to radiation, stars must either contract or obtain energy from nuclear reactions. Since nuclear reactions change their composition, stars must evolve.
The Virial Theorm implies that if a star is neither too degenerate nor too relativistic (radiation or pair dominated) GM 2 R MN AkT (for an ideal gas) T GM N A kr R So T 3M 4πρ GM 2/3 1/3 N A k ρ1/3 for constant density That is as stars of ideal gas contract, they get hotter and since a given fuel (H, He, C etc) burns at about the same temperature, more massive stars will burn their fuels at lower density, i.e., higher entropy.
Central Conditions at death the iron cores of massive stars are somewhat degenerate ρ 3 T
Four important time scales for a (non-rotating) star to adjust its structure can be noted. The shortest by far is the time required to approach and maintain hydrostatic equilibrium. Stars not in a state of dynamical implosion or explosion maintain a balance between pressure and gravity on a few sound crossing times. The sound crossing time is typically comparable to the free fall time scale. and dp dr = GM ρ r 2 P cent GM ρ R P c s = γ ρ v esc 1/2 GM R 1/2 So the escape speed and the sound speed in the deep interior are comparable GM R 1/2
The free fall time scale is ~R/v esc so τ esc R v esc = R 3 2GM 1/2 3 8πGρ 1/2 The number out front depends upon how the time scale is evaluated. The e-folding time for the density is 1/3 of this τ HD 1 24πGρ 1/2 446 ρ sec which is often used by the nuclear astrophysicist
The second relevant adjustment time, moving up in scale, is the thermal time scale given either by radiative diffusion, appropriately modified for convection. This is the time it takes for a star to come into and maintain thermal steady state, e.g., for the energy generated in the interior to balance that emitted in the form of radiation at the surface. τ therm = R2 D therm nb. These two time scales, HD and therm sometimes place limits on the accuracy with which jobs run in one code can be remapped into another (e.g., 1D -> 3D). KH can also be important
The thermal diffusion coefficient (also called the thermal diffusivity) is defined as D = conductivity heat capacity = K C P ρ where K appears in Fourier's equation heat flow = -K T For radiative diffusion the "conductivity", K, is given by K = 4acT 3 3κρ (see Clayton 3-12) where κ is the opacity (cm 2 g 1 ), thus D = 4a ct 3 3κC P ρ 2 where C P is the heat capacity (erg g 1 K 1 ) D here thus has units cm 2 s 1 Note that D c κρ at 4 ρc P T = c κρ radiation energy content total heat content
If radiation energy density is a substantial fraction of the internal energy (not true in the usual case), D ~ c/κρ with κ the opacity, and taking advantage of the fact that in massive stars electron scattering dominates so that κ = 0.2 to 0.4 cm 2 g 1, the thermal time scale then scales like τ rad R2 κρ c ~ 0.2R2 c 3M 4π R 3 M R however, massive stars have convective cores so the thermal time is generally governed by the diffusion time in their outer layers. Since the dimensions are still (several) solar radii while the densities are less and the opacity about the same, the radiative time scales are somewhat less than the sun ( 10 5 yr; Mitalas and Sills, ApJ, 401, 759 (1992)).
Possible point of confusion. The radiative diffusion time vs the thermal time. The former is the time for a photon to diffuse out τ rad ~ R 2 mfp c ~ R 2 κρ c The latter is the time for an appreciable fraction of the heat content to diffuse out. It is longer if most of the heat is in the ions, not the radiation τ thermal ~τ rad In the sun these differ by a factor of about 100 and the thermal time is of order the Kelvin Helmholtz time τ thermal ~ τ rad total thermal energy in star radiation content of star E matter E rad ~ E matter L rad αgm 2 = τ KH 2RL rad http://link.springer.com/article/10.1023/a%3a1022952621810
In general a diffusion coefficient is given by D ~ 1 3 v l where v is a typical speed (light or a convective element) and l is a characteristic length scale (mean free path or pressure scale height)
A closely related time scale is the Kelvin Helmholtz time scale τ KH GM 2 RL M 5/3 ρ 1/3 L i.e. R ~ ( M / ρ) 1/3 Except for very massive stars, L on the main sequance is proportional to M to roughly the power 2 to 4, and ρ decreases with M so the Kelvin Helmholtz time scale is shorter for more massive stars. Note that there are numerous Kelvin Helmholtz time scales for massive stars since they typically go through six stages of nuclear burning. During the stages after helium burning, L in the heavy element core is given by pair neutrino emission and the Kelvin Helmhotz time scale becomes quite short - e.g. a protoneutron star evolves in a few seconds.
e.g., the sun is in thermal steady state. A presupernova star is not.
3 ρ T
Final composition of a 25 M 0 star: up to Si burned up to O burned up to Ne-burned up to He burned up to H-burned unburned mass fraction interior mass (M 0 )
25 M Presupernova Star (typical for 9-130 M ) 1400 R 0.5 R H, He He 240,000 L C,O O,Mg,Ne Si, S, Ar, Ca Fe He, C 0.1 R 0.003 R
Examples of Time Scales In stellar explosions, the relevant time scale is the hydrodynamic one. Explosive nucleosynthesis happens when HD ~ nuc During the first stages of stellar evolution KH < nuc. The evolution occurs on a Kelvin Helmholtz time In a massive presupernova star in its outer layers, nuc < thermal. The outer layers are not in thermal equilibrium with the interior. It also happens in these same stars that nuc (core) < KH (envelope). The core evolves like a separate star. Rotation and accretion can add additional time scales. E.g. Eddington Sweet vs nuclear. Accretion vs nuclear. Convection also has its own time scale.
Central Conditions for Polytropes and hence m(r) = const*r 3 GR 2 ρ 2 since 1 R = 4πρ c 3M 1/3 for constant density
Eq. 2-313
= 24 most stars have 1.5 < n < 3
Actually the dimensions of Y are Mole/gm and N A has dimensions particles per Mole. P = Σ n i kt = N A k µ ρt
X i = mass fraction of species i ρn A Y e = n e = Z i n i 0.5 < µ < 2 = ρn A Z i Y i (The limit µ=2 is achieved as A goes to infinity and Z = A/2, i.e. electrons dominate)
For advanced stages of evolution where A > 1, most of the pressure is due to the electrons(and radiation) P c ρ c T c µ This would suggest that the ratio would increase as the star evolved and µ became greater. T c µρ c 1/3
The decline of T 3 / mostly reflects the fact that beyond H burning the star becomes a red giant and no longer is a single polytrope. M core is essentially reduced. This more than compensates for the increase in µ. A decreasing also decreases T 3 /. n the final stages degeneracy becomes important. ρ 3 T
= P ideal P total P tot P ideal β ρ c T c µβ Decrease in as star evolves acts to suppress T 3 /.
CRITICAL MASSES Ideal gas plus radiation
note cancellation of ρ c = 1.45 M if Y e = 0.50 neglecting Coulomb corrections and relativistic corrections. 1.39 M if they are included.
For lower densities and hence degenerate cores significantly less than the Chandrasekhar mass non-relativistic degeneracy pressure gives another solution ( = 5/3) This is the well known central density mass relation for (non-relativistic) white dwarfs
Specifying a core mass gives a maximum temperature achieved before degeneracy supports the star. If that maximum temperature exceeds a critical value, burning ignites M 1 > M 2 M i <M Ch
He core CO Core Ne-O Core 100% He 50%C- 50% O 75% O- 25% Ne Ne-O cores CO-cores He-cores
Critical Temperatures
These calculations give critical masses: C ~1.0 Ne ~1.25 O 1.39 Si 1.39 These are core masses. The corresponding main sequence masses are larger. All stars with main sequence mass above the Chandrasekhar mass could in principle go on to burn Si. In fact, that never happens below 9 solar masses. Stars develop a red giant structure with a low density surrounding a compact core. The convective envelope dredges up helium core material and causes it to shrink. Only for stars above about 7 or 8 solar masses does the He core stay greater than the Chandrasekhar mass after helium burning.
Main Sequence Critical Masses 0.08 M Lower limit for hydrogen ignition 0.45 M helium ignition 7.25 M carbon ignition 9.00 M neon, oxygen, silicon ignition (off center) 10.5 M ignite all stages at the stellar center ~ 80 M First encounter the pair instability (neglecting mass loss) ~35 M Lose envelope if solar metallicity star These are for models calculated with the KELER code including semiconvection and convective overshoot mixing but ignoring rotation. With rotation the numbers may be shifted to slightly lower values. Low metallicity may raise the numbers slightly since less initial He means a smaller helium core. Other codes give different results typically to within 1 solar mass.
Between 7.25 and 10.5 solar masses the evolution can be quite complicated and code dependent owing to the combined effects of degeneracy and neutrino losses. Off-center ignition is the norm for the post-carbon burning stages. Mass loss introduces additional uncertainty, especially with regard to final outcome. Does a 8.5 solar mass main sequence star produce a NeO white dwarf or an electron-capture supernova (TBD). Above 9 solar masses an iron core eventually forms on up to the pair instability limit around 150 solar masses.
Off-Center Ignition
The death of a star and how it may potentially explode is also very sensitive to: The density structure surrounding the iron core The rotation rate of the core and that material The density structure depends on the entropy of presupernova cores (TBD). Higher entropy cores occur for higher masses and are less degenerate and less centrally condensed.
Density Profiles of Supernova Progenitor Cores These make the heavy elements These should be easy to explode 2D SASI-aided, Neutrino-Driven Explosion?
Density Profiles of Supernova Progenitor Cores These make the heavy elements These should be easy to explode 2D SASI-aided, Neutrino-Driven Explosion?
Sukhbold and Woosley (2015) Hard to Explode Easy to Explode Hard to Explode Easy to Explode rapid rise above here Below 13 solar masses very easy to explode
Mass (solar masses) End point Remnant < 7 to ~8 planetary nebula CO white dwarf ~8 to ~11 ~11 - ~20 degenerate core low energy SN neutrino powered supernova; SN Ibc in binary. Islands of explosion at higher mass Ne-O WD below 9? neutron star above 9 neutron stars and black holes 20-85 85 150 150-260 > 260 without mass loss probably no SN (unless rotationally powered); with mass loss SN Ibc pulsational pair SN if low mass loss pair instability SN if low mass loss pair induced collapse if low mass loss black hole if enough mass loss neutron star black hole none black hole