Hayashi Limit
Total Page:16
File Type:pdf, Size:1020Kb
HAYASHI LIMIT Let us consider a simple ”model atmosphere” of a star. The equation of hydrostatic equilibrium and the definition of optical depth are dP dτ = −gρ, = −κρ, (s2.24) dr dr and may be combined to write dP g = . (s2.25) dτ κ Assuming κ = const we may integrate this equation to obtain 2 g GM Pτ=2/3 = , g ≡ 2 , (s2.26) 3 κτ=2/3 R where the subscript τ =2/3 indicates that we evaluate the particular quantity at the photosphere. Let us consider a cool star with the negative hydrogen ion H− dominating opacity in the atmosphere. When temperature is low we may neglect radiation pressure in the atmosphere. Adopting k 0.5 7.7 −25 0.5 P = ρT, κ = κ0ρ T , κ0 = 10 Z , (s2.27) µH we may write the equation (s2.26) as k 2 1 GM ρT = 0 5 7 7 2 , (s2.28) µH 3 κ0ρ . T . R which may be rearranged to have 1.5 8.7 2 µH G M ρ T = 2 . (s2.29) 3 k κ0 R We know that a star with the H− opacity in the atmosphere becomes convective below optical depth τ =0.775, i.e. very close to the photosphere. Let us suppose that the convection extends all the way to the stellar center, and let us ignore here all complications due to hydrogen and helium ionization. Convective star is adiabatic, and if it is made of a perfect gas with the equation of state (s2.27) then it is a polytrope with an index n =1.5. Therefore, we expect a polytropic relation all the way from the photosphere down to the center, and we have ρ ρc 3M µH GM 1.5 = 1.5 , ρc =5.99ρav =5.99 3 , Tc =0.539 . (s2.30) T Tc 4πR k R Combining equations (s2.29) and (s2.30) we obtain 3.25 10.95 −1.5 2.25 µHG M 0.10 µHG 1.75 0.25 T ≈ ρc Tc 2 ≈ M R . (s2.31) kκ0 R κ0 k Let us make an approximation that convection begins at the photosphere, i.e. at T = Teff , and let us replace stellar radius with the combination of effective temperature and luminosity according to 2 4 L =4πR σTeff : 3.25 11.45 ≈ 0.07 µHG 1.75 1/8 Teff 1 8 M L , (s2.32a) κ0σ / k 0.15 0.01 −0.04 3 M L Z Teff ≈ 2 × 10 , (s2.32b) M⊙ L⊙ 0.02 To the right of the Hayashi limit no stars in a hydrostatic equilibrium can exist. s2 — 1 Ignition Mass for Low-Mass Stars Low mass stars are supported by gas pressure, while radiation pressure is unimportant. Electron gas may be partly degenerate. Numerical models demonstrate that very low mass stars, with M < 0.3M⊙ are fully convective, and may be very well approximated with n =1.5 polytropes. We shall approximate the non.-rel. equation of state with the following formula 1/2 2 2 k 5/3 P ≈ ρT + K1ρ , (lms.1) " µH # 8 −1 −1 13 −5/3 −5/3 2 where k/H = 0.825 × 10 [erg g K ], K1 = 0.991 × 10 µe [erg g cm ]. Algebraic approximation to the stellar structure equations gives M GM 2 ρ ≈ , P ≈ . (lms.2) R3 R4 Combining the last two equations we obtain 2 4 2 2 2 10/3 G M k 2 10/3 kT M 2 M ≈ ρT + K1 ρ ≈ + K1 , (lms.3) R8 µH µH R6 R10 which may be written as 2 1/2 µH GM K1 T ≈ 1 − . (lms.4) k R GRM 1/3 " # Fully degenerate, the stellar radius satisfies a polytropic ( n =1.5 ) mass-radius relation: K1 5/3 k R = , K1ρ ≫ ρT, (lms.5) min 0.4242 GM 1/3 µH where Rmin is the minimum radius that a star with a mass M may have. This relation is recovered 2 from the equation (lms.4) when we replace the square bracket with [1 − (Rmin/R) ]. In the limit when degeneracy is negligible we should have µH GM 5/3 k T =0.539 , K1ρ ≪ ρT, (lms.6) c k R µH where Tc is the central temperature of an n =1.5 polytrope with pressure provided by non-degenerate gas. Combining the two limiting cases with the equation (lms.4) we may write it as 2 1/2 µH GM R K1 T =0.539 1 − min , R = . (lms.7) c k R R min 0.4242 GM 1/3 " # The central temperature reaches its maximum where dTc/dR =0 for R = RTmax , i.e. 1/2 RTmax =2 Rmin, (lms.8) and the corresponding maximum central temperature is 0.539 µH GM 0.539 µH GM Tc,max = 1/2 = = (lms.9) 2 k RTmax 2 k Rmin 2 4/3 4/3 µHG 4/3 × 8 5/3 M ≈ × 6 M 0.1143 M =1.56 10 µµe 6 10 [ K ], kK1 M⊙ 0.1 M⊙ lms — 2 Sch¨onberg-Chandrasekhar Limit After core hydrogen burning, the quasi-isothermal helium core grows by the accumulation of helium ashes from the hydrogen burning shell surrounding it. However, the core can maintain itself in hydrostatic equilbrium only if the mass of the helium core is below the so-called “Schonberg-Chandrasekhar"¨ limit. If it exceeds this limit, the core starts to contract on Kelvin-Helmholtz timescales, inaugurating the next phase of stellar evolution. We can derive the SC critical mass by using the Virial theorem with the core boundary pressure retained and assuming an ideal gas, finding the maximum surface pressure the core can sustain (calculated in terms of core quantities), calculating the surface pressure due to the weight of the overlying envelope, and then setting these two pressures equal. A rearrangement of the resulting equation yields the ratio of the critical core mass to the total stellar mass (as long as this ratio is small) in terms of the mean molecular weights in the envelope and core. The derivation proceeds as follows: We state the Virial theorem: 3 4πRcPc − 2Uc =Ωc , (sc.1) 3 where U = Mc kT and the subscript c stands for “core." We will assume that Ω ∼ c 2 µcmp c c 3 2 GMc . Then, we solve for P , and take dP /dM , and set it to zero. The result, written for 3Rc c c c Rc, is 2µ m GM R = c p c . c 5kT c Plugging this radius into eq. (sc.1) and solving for Pc, we have 4 375 kT P (crit)= c . c 64πG3M 2 µ m c c p The pressure at the surface of the core due to the envelope weight is M 2 2 GM G M − Mc Pc,env = − 4 dM ∼ 4 . (sc.2) M 4πr 8π < r > Z c sc—3 We assume that < r4 >∼ R4/2, where R is the outer radius of the star and M is the total mass of the star. The pressure on the core due to the envelope, Pc,env, can be written in terms of the mass density (ρc,env), temperature (Tc), and envelope mean molecular mass (µenv), and ∼ 3M if we assume that ρc,env 4πR3 , we find (solving for R) µ m GM R ∼ env p . 3kTc Plugging this into the eq. (sc.2) for the pressure, we find 4 81 kT P (crit)= c . c 4πG3M 2 µ m env p Setting the two pressures equal (note that Tc is in both expressions), we find 2 M µ c ∼ 0.54 env . M µ c Doingthisabitmorerigorouslyyieldsacoefficientof 0.37, not much different. If µenv ∼ 0.63 ∼ Mc ∼ ∼ and µc 1.34, we find that M 0.08 ( 8%). This is the Schonberg-Chandrasekhar¨ ratio. Note that the subsequent KH collapse of the core leads to a self-bound object with central pressures much larger than the pressure in the hydrogen-burning shell. As we will state again, “once an isothermal core becomes self-gravitating, it remains self-gravitating forever, even if a new nuclear fuel ignites in the core, and the core is no longer isothermal." sc—4.