Ay 1 – Lecture 8

Stellar Structure and the 8.1 Basics How Work • : and radiation balance the

• Thermal Equilibrium: Energy generated = Energy radiated

• Thermonuclear Reactions: The source of the energy

• Energy Transport: How does it get from the core to the surface Equations of the Stellar Structure

Mass vs. radius: 2 dm = 4πr ρdr r dr dm m dm = 4πr2ρ ρ dr R M € vs. radius: dL = 4πr2ρdr × q € dL = 4πr2ρq dr q = rate of energy generation per unit mass

€ Hydrostatic Equilibrium At a given radius, the gravitational force on a shell is: dr GmΔm r m dm Fg = − 2 r ρ R M The of that mass shell over the area has to be the difference between the on an inner and an outer surface€ of the shell. dm = 4πr2ρdr dP Gm = − ρ dr r2 €

€ Pressure of What? Total pressure: P = Pgas + Pradiation

The for an is: Pgas = nkT n = the number of particles per unit volume (ions and )

Or: ρ Pgas = × kT € ρ = mass € µmH mH = the mass of atom μ = average particle mass in units of mH k The ideal gas constant: R ≡ m Depends on the € H chemical composition. R Typically μ ~ 0.8 Thus: P = ρT gas µ €

€

Momentum exchange between and electrons/ions results in a radiation pressure

1 For blackbody radiation: P = aT 4 r 3 …where a is the radiation constant: 8π 5k 4 4σ a = = = 7.565 ×10−15 erg cm-3 K-4 15c 3h 3 c Radiation pressure dominates €over the gas pressure inside very massive stars, which are hotter € Equations of Stellar Structure At radius r in a static, spherically symmetric and the density ρ: dm = 4πr2ρ Mass conservation dr dP Gm = − ρ Hydrostatic equilibrium dr r2 dT 3 κρ L Energy transport due to = − dr 4ac T 3 4πr2 radiation (only) dL = 4πr2ρq dr Energy generation

4 equations with 4 unknowns - enough for a solution once we know P(ρ ,T), opacity κ, and q € Some Order-of- Estimates

Let’s see if we can estimate roughly the conditions in the . Pressure P = F / A: 2 2 F ≈ G M / R  2 A ≈ 4 π R  2 4 P ≈ G M / 4 π R  33 10 - 8 (M ≈ 2  10 g, R  ≈ 7  10 cm, G ≈ 6.7  10 cgs) 15 2 Thus: Pest ~ 10 dyn / cm -- and surely an underestimate 17 2 True value: Pc ≈ 2  10 dyn / cm

Now the : 3/2 k T ≈ G mp M / R - 16 - 24 ( k ≈ 1.4  10 erg/K, mp ≈ 1.7  10 g ) 7 Thus: Test ≈ 1.6  10 K 7 True value: Tc ≈ 1.57  10 K -- not bad! Standard Solar Model 8.2 Energy Generation in Stars Energy Production in Stars: Thermonuclear Reactions Mass of nuclei with several and / or does not exactly equal mass of the constituents - slightly smaller because of the binding energy of the nucleus The main process is hydrogen fusion into : 4×1H→4 He + photons and 4 protons, total mass = helium nucleus, mass = 4 × 1.0081 = 4.0324 amu 4.0039 amu Mass difference: 0.0285 amu = 4.7 x 10-26 g € ΔE = ΔMc 2 = 4.3×10−5 erg = 27 MeV Or about 0.7% of the total rest mass

€

The P-P Cycle Binding Energy Per Nucleon vs. Atomic Number

Energetically favorable Energetically unfavorable

Most bound nucleus is iron 56: up to A < 56 fusion releases energy, but A > 56 fusion requires energy input Overcoming the Electrostatic Barrier

Energy Coulomb barrier At r = r0, height of the potential Coulomb barrier is: E(r) 2 Z1Z2e + E = ~ Z1Z2 MeV r0

r r0 € + + + Most particles have too low Nuclear strong Electrostatic thermal energies to overcome interactions repulsion the Coulomb barrier!

What makes the fusion possible is the quantum tunneling effect Probability of tunneling increases steeply with particle energy The Gamow Peak The most energetic nuclei are the most likely to fuse, but very few of them in a thermal distribution of particle speeds:

Number of Gamow P for fusion particles peak

Particle energy

Narrow range of energies around the Gamow peak where significant numbers of particles in the are able to fuse. Energy is >> typical thermal energy, so fusion is slow Thermonuclear Reactions (TNR) • Burning of H into He is the only energy generation process on the , where stars spend most of their lives; all others happen in post-MS evolutionary stages – ~ 4.3 million tons of H into He per second • In addition to the p-p cycle, there is the CNO Cycle, in which the C, N, O, nuclei catalyze the burning of H into He • The rates of TNR are usually very steep functions of temperature, due to high potential barriers

• Generally, more massive stars achieve higher Tc , and can synthesize elements up to Fe; beyond Fe, it happens in SN explosions Self-Regulation in Stars Suppose the fusion rate increases slightly. Then, (1) Temperature increases (2) Pressure increases (3) Core expands (4) Density and temperature decrease (5) Fusion rate decreases So there's a feedback mechanism which prevents the fusion rate from skyrocketing upward This is the inverse of the core collapse mechanism discussed for the 8.3 Energy Transport in Stars Energy Transport Mechanisms in Stars How does the energy get out? 1. Radiatively ( diffusion) 2. Convectively 3. Conduction (generally not important in stars)

… and the reality is fairly complex Radiative Energy Transfer • As the diffuses from the core outwards, the photons are scattered by the dense plasma inside the star • For the Sun, it takes ~ 250,000 years for the energy to reach the surface

• The opacity of the plasma depends on the temperature, density, and chemical composition • If the plasma is too opaque, becomes a more efficient mechanism for the energy transfer When Does the Convection Happen?

The Schwarzschild criterion: Imagine displacing a small mass element vertically upward by a distance dr. Assume that no heat is exchanged with the surrounding, i.e. the process is adiabatic: • Element expands to stay in pressure balance • New density will not generally equal the new ambient density If this mechanical energy transport is more efficient than the radiative case, the medium will be convectively unstable 8.4 The Sun, Our Star Why Study the Sun? • The nearest star - can study it in a greater detail than any others. This can help us understand better the overall stellar and phenomenology – Radiation transfer, convection – Photospheric and chromospheric activity • Kind of important for the life on … – Solar activity has terrestrial consequences • A gateway to the (and physics) – Thermometry of stellar cores and the – Neutrino oscillations A Theoretical Model of the Energy Transfer in the Sun • Hydrogen fusion takes place in a core extending out to about 0.25 • The core is surrounded by a radiative zone extending to about 0.71 solar radius. Energy transfer through radiative diffusion • The radiative zone is surrounded by a rather opaque convective zone of gas at relatively low temperature and pressure. Energy transfer through convection We can probe the solar interior using the Sun’s own vibrations

• is the study of how the Sun vibrates • These vibrations have been used to infer pressures, , chemical composition, and rotation rates within the Sun • Major contributions from Caltech (Libbrecht et al.) Rotation of the Solar Interior Solar / Surface Layers

• The Sun’s atmosphere has three main layers: 1. the 2. the 3. the corona • The visible surface of the Sun, the photosphere, is the lowest layer in the solar atmosphere Emergent spectrum of the photosphere is close to a blackbody T ~ 5800 K Convection in the photosphere produces granules are low-temperature regions in the photosphere Sunspots come in groups, and follow the Sun’s . They start from the higher latitudes and migrate towards the Solar equator.

They correlate well with other manifestations of the Solar activity. • The Sun’s surface features (including numbers) vary in an 11-year cycle; it is really a 22-year cycle in which the surface increases, decreases, and then increases again with the opposite polarity • There are probably also longer period cycles

Note the scale: Earth’s diameter is < 1% of the Solar diameter Activity in the corona includes coronal mass ejections and coronal holes The Sun’s magnetic field also produces other forms of solar activity

• A is a brief eruption of hot, ionized from a sunspot group • A is a much larger eruption that involves immense amounts of gas from the corona The Sun is Now Being Monitored by a Number of Satellites and Ground-Based Observatories

SOHO far-UV composite Yohkoh soft X-ray  Solar Neutrinos and the Birth of Neutrino Astronomy

• Detection of Solar neutrinos offers a unique probe of deep stellar interiors - a fundamental test of our understanding of stars and their energy production

• For many years, there was a factor of 3 discrepancy between the theoretical predictions and the experiment (the “ problem”). The resolution of it provided a fundamental physical insight (neutrino oscillations) Solar Neutrino Flux Neutrinos from the main p-p chain are of very low energy. Less important reactions (energetically) yield a smaller flux of higher energy neutrinos: − 2 p + e + p→ H + ν e `pep’ - 1.4 MeV neutrino 3 4 + He + p→ He + e + ν e 0 - 18.8 MeV − 7 7 e + Be→ Li + ν e 0.383, 0.861 MeV 8 8 + 0 - 15 MeV B→ Be + e + ν e Since we get 2 neutrinos for each 28 MeV of energy, we can use Solar luminosity to calculate neutrino flux at Earth:

2Lsun 1 10 2 € Neutrino flux = × 2 ~ 6 × 10 neutrinos/s/cm 28 MeV 4πd

€ Homestake Mine Detector First attempt to detect Solar neutrinos began in the 1960s:

Detector is a large tank containing 600 tons of

C2Cl4, situated at 1500m depth in a mine in South Dakota. Neutrinos interact with the to produce a radioactive isotope of : 37 - 37 ν e + Cl → e + Ar

+ an which is not observed. € Super Kamiokande look for Cherenkov radiation from high energy electrons in Measure: − − ν e + e → ν e + e water

€