Stellar Interiors. Hydrostatic Equilibrium 2. Mass continuity 3. Equation of State. The pressure integral 4. Stellar Energy Sources. Where does it come from? 5. Intro to Nuclear Reactions. Fission 2. Fusion 3. Binding Energies 4. Bringing them together 5. Energy / Tunneling 6. Cross Section 7. p-p chain 8. CNO Cycle 6. Luminosity Gradient 7. Energy Transport. Radiative 2. Convective Transport 8. Stellar Model Building Fig. The basic structure of a star. Figure shows the basic structure of a star. (It is not scale!) The light we see comes from photons in the stellar atmosphere. Hydrostatic Equilibrium Page
P+ P r F grav P Fig. 2 A small volume of gas has three force vectors acting in the vertical directions: gravity, pressure from below and pressure from above. Let's imagine a small cylindrical volume of gas that is located in a larger region. There will be three forces acting on it: F grav is the force of gravity, and two forces from the pressure, one from the pressure above: P + ΔP and one from pressure below: P. The net force from the pressure will be : where A is the area (whose surface normal points in the y direction) of the gas volume. The gravitation force will be given by: where M r is the mass enclosed by a sphere of radius r centered about the earth's center. If this volume of gas is in equilibrium, then the pressure forces and the gravitational forces will be equal: Solving for ΔP: F pres = A[P (P + ΔP)] () GM r ρaδr F grav = (2) r 2 GM r ρaδr A[P (P + ΔP)] = (3) r 2 GM r ρ ΔP = Δr (4) r 2 If we consider this cylinder to be infinitely thin: Δr =, then we have a differential form: dp GM r ρ = (5) dr r 2 Eq. (5) is called the equation of hydrostatic equilibrium (for a spherical symmetric body). To apply this to a star, we should acknowledge the limiting assumptions we made:. The star is spherical 2. The star is not rotating 3. The star is not expanding or contracting 4. Gravity and Pressure gradients are the only forces acting. A better way to write this would be to include the dependencies: Page 2
dp dr GM(r)ρ(r) = (6) r 2 We can see that there are 3 quantities that depend on r: Pressure ( P), Mass ( M) and Density ( ρ). This makes a unique solution not possible. We need more equations! Mass continuity dm dr = 4πr 2 ρ(r) (7) The equation of mass continuity relates density and radial position and total mass. In short it says that the total mass is just the sum of all the infinitesimal spherical shells (i.e. onion skins) Equation of State PV = nrt (8) The equation of state (8) relates density, temperature, and pressure for a gas. For many stars, the ideal gas law will work just fine. Mean Molecular Weight μ m m H where m H =.6735... 27 kg, and m is the average mass of a gas particle in kg. μ therefore gives the the average mass of a free particles in the gas, in units of hydrogen. And, let's say that n = ρ/m. This allows us to write the ideal gas law as: P g = ρkt μm H (9) The pressure integral P = pvdp 3 n p () Derivation of the pressure integral: Assuming a particle in the x direction, normal to a surface: fδt = Δp = 2p x i^ Since it must travel the distance Δx twice in the time Δt: Δt = 2 Δx v x Page 3
Thus the average force: 2p f = x = Δt p x v x Δx We can sub: pv p 3 x v x f(p) = pv 3 Δx The number of particles with a given momentum between p and p + dp is dp, then the total: Contribution to the total force by all particles within that momentum range: Divide by A N = N p dp N p df(p) = f(p) N p dp = pv dp 3 Δx n p dp N p ΔV dp N p P = pv dp 3 n p () For non-relativistic particles, p = mv: P = m dv 3 n v v 2 (2) Recall the Maxwell-Boltzmann distribution: Putting this into Eq. (2) we'll get a simple expression: (which is just the ideal gas law since n = N/V ) Stellar Energy Sources m 3/2 n v dv = n( ) 4πv 2 e mv2 2kT dv 2πkT P g = nkt (3) Where does it come from? The first suggestion to explain where stars got the energy was to look at the gravitational potential. Let's say we have two particles, then the gravitational Page 4
potential energy between them is: U = G Mm r (4) If the distance r between these particles decreases, then the potential energy must get more negative, and we would expect to find energy in other forms, kinetic for example. How much energy is there in gravity? Of course there are more than two particles in a star, so we'll have to expand that equation a bit. If a point mass dm is located outside a spherically symmetric mass M r, then the force on that point mass will be directed to the center and will have a magnitude: The potential energy of the point mass is thus: d F g,i = G M r dm i r 2 Now we'll consider a shell of thickness dr and mass dm: d U g,i = G M r d r m i (5) where ρ is the density and 4πr 2 dr is the volume of the shell. Thus, eq (5) becomes We can integrate over the radius of the star R like this: U g Of course we probably don't know ρ very well at this point, but we can guess based on it's average value: with M being the total mass of the star. Also: dm = 4π r 2 ρdr M r 4πr 2 ρ d U g = G dr r R = 4πG M r ρr dr (6) ρ ρ = M 4 π 3 R3 Putting this into eq. (6), we get: 4 M r π 3 r3 ρ 6π U 2 g Gρ 2R 5 3 GM 2 5 5 R Applying the virial theorum, we have: E 3 GM 2 R (7) (8) Page 5
If we plug some sun values into this: Now if the sun luminosity was constant at today's value: 3.84 W, then this would mean there is enough energy to last for about 7 years (about million years). This is notably too short given other observations from geology that suggest the earth and moon and things have been around for billions of years. There must be another source of energy contributing to the star's luminosity. Intro to Nuclear Reactions 3 GM 2. 4 J R 26 Fission Fig. 3 http://www.nuclear-power.net/nuclear-power/fission/ When an atom breaks apart. Power plants... Fusion Fig. 4 http://www.nuclearpower.net/nuclear-power/nuclearfusion/ When you smash two atoms together. Suns... Binding Energies Page 6
Fig. 5 Carroll & Ostlie, figure -9. Take four Hydrogen atoms ( proton & electron), each have a mass of.7825324 u. Together, they therefore have a combined mass of 4.333 u. One Helium-4 atom has a mass of 4.263 u. This is obtained by adding the mass of the neutrons and the protons together. Thus, let's say we were able to mash 4 Hydrogen atoms together and make a Helium. Now there would be a change in mass between the initial and final constitutions: Δm = 4.333 4.263 =.28697 u. Since E = mc 2, we can see that this amount of mass is associated with an energy of 26.73 MeV. This is essentially the amount of energy in the formation of a Helium nucleus, aka the binding energy. (This is also the amount of energy that would be required to take apart a Helium nucleus.) The trick is to figure out exactly how we can combine 4 Hydrogen atoms to make a Helium. (that's the study of High energy physics earlier last century) Bringing them together Page 7
U(r) F(r) PHY 454 - stellar-structures - J. Hedberg - 27 One impediment is the Coulomb repulsion between two protons, both of which are positively charged. r Fig. 6 The plots of the Coulomb force (top) and the Coulomb potential (bottom). r As the separation distance between two charged (alike) particles decreases, the force of repulsion increases. The Coulomb potential between two protons located fm apart would be: (.6 9 C) 2 U.4 6 ev.4 MeV 4π(8.8 2 C 2 J m )(.6 9 J ev )( 5 m) If two protons wanted to get closer than fm, they would need to have a kinetic of energy greater than this value to overcome the potential barrier. The energy of a proton at the center of the sun is about: 3kT E = c 3(.38 23 J K )(.47 7 K) 2 kev (2) 2 2(.6 9 J ev ) This would be insufficient energy to get over the barrier. This is where quantum tunneling comes in. Tunneling Play/Pause barrier. Fig. 7 A -d simulation of a wave packet incident on a potential barrier. Some of the wave tunnels through the Simulation based on code from https://github.com/jakevdp/pyschrodinger Criteria for tunneling: a) the distance between protons should be comparable to the de Broglie wavelength: λ DB h = = 7 3 m (2) p and b) the kinetic energy ( E) of the proton needs to be comparable to the electrostatic potential energy at that separation. (2 E /2 p U = = E Page 8
U (r) [A.U] PHY 454 - stellar-structures - J. Hedberg - 27 e 2 e 2 (2m p E) /2 U = = E 4πϵ λ DB 4πϵ h (22) Solving for the kinetic energy: e E ( 2 2m p ) αm 3keV 4πϵ h 2 2π 2 p c 2 (23) Since this energy is comparable to the kinetic energy of a proton in the core of the sun, we can expect tunneling through the Coulomb barrier to be possible. The next step is to figure out just how likely such events are. Energy / Tunneling Coulomb repulsion (/ r) 2 4 6 8 r (fm) Strong Nuclear potential well A qualitative graph showing the what happens when two protons get close engouh for the strong force to dominate. Cross Section The cross section gives a probability of interaction. It is the number of reactions per target nucleus per time divided by the flux of incident particles. number of reactions/nucleus/time σ(e) number of incident particles/area/time Cross Section Page 9
ds = v(e) dt σ (E) Fig. 8 Imagine a cylinder of length ds = v(e)dt and area σ(e). All the particles within the volume of the cylinder will reach the wall in the time dt. If a particle hits the wall, then there is a reaction. This is an abstraction of the nuclear collision process. (there is no cylinder, no wall...) n : just a particle n i : an incident particle n ie de : # of incident particles having energies between E and E + de per unit volume dn E : # of particles that "hit the wall" / time or # of reactions / time n E de : # of particles with an energy between E and E + de, (might not be incident) We have a target particle ( x) and an incident particle ( i). If the incident particle strikes the area σ(e), then there will be a nuclear reaction. Take the number of incident particles per unit volume having energy between E and E + de to be: n ie de. (This is like a density) To find the number of such particles inside the cylinder: d N E = n ie de density σ(e)v(e)dt volume The number of particles with the appropriate kinetic energy is some fraction of the whole: n ie n ie de n de = de = de n i E n E de n n E The number of reactions per target nucleus per time dt having energies between E and E + de is reactions per nucleus time interval dn = E n = σ(e)v(e) i de dt n n E Lastly, if there are n x targets per unit volume, then the total number of reactions per unit volume per unit time, given all possible energies is: ix = σ(e)v(e) de x i Page
r ix = n x n i σ(e)v(e) de n n E (24) But, what exactly is σe? σ(e) σ(e) πλ 2 DB π( h ) p E 2 (25) (Since K = E = μ m v 2 /2 = p 2 /2μ m ) Also, for a barrier U and particle with energy E, the probability for tunneling through will go as: In this case the U is the Coulomb potential: = and E is just the kinetic energy: E = μ. 4πϵ r m v 2 /2 We can assume the distance r is basically the de Broglie wavelength: r λ = h/p. In which case we obtain: where b is given by σ(e) e 2π2 U/E (26) U c Z Z 2 e 2 σ(e) e be /2 (27) b πμ/2 m Z Z 2 e 2 2 /2 ϵ h We can combine eq (25) and (27) and define a general function S(E) to get a functional form for σe: σ(e) = S(E) E e be /2 (28) Finally, we can write a reaction rate as: 3/2 2 n r ix = ( ) i n x S(E) de kt ( μ m π) e be /2 e E/kT (29) /2 Page
likelihood [A.U.] PHY 454 - stellar-structures - J. Hedberg - 27. e -E/k T e -b E -/2 ( 6 ).8.6.4 e -E/k T e -b E -/2 ( 3 ).2 2 4 6 8 Energy (kev) Fig. 9 The likelihood that a nuclear reaction will occur. The peak will be at bkt E = ( ) 2 2/3 (3) p-p chain 4 4 H 2He + 2 e + + 2 ν e + 2γ H H H H ν e ν e e + e + 2 H H H 2 H γ γ 3 He 3 He H H 4 He Fig. The full p-p chain One way to convert Hydrogen to Helium is through the proton-proton chain (PP). The overall reaction combines 4 protons to create Helium-4 plus some other stuff. The first step in this reaction take two protons and produces a deuterium, a positron, and an electron neutrino. H + H 2 H + + + e Page 2
H + H 2 H + e + + ν e (3) Then the deuterium combines with another proton to make a Helium-3 atom plus a gamma ray. 2 H + H 3 2 He + γ (32) Lastly, two Helium-3s combine to create Helium-4 and 2 protons 3 He + 3He 4He + 2 H (33) 2 2 2 CNO Cycle 4 He H H 2 C γ 5 N 3 N ν e e + ν e e + 3 C 5 O H γ H 4 N γ Fig. The CNO cycle In stars that more than.3m, a different nuclear fusion process will dominate. This is the CNO cycle. It's more complicated, but the basic idea is that Carbon, Nitrogen, and Oxygen act a catalysts to speed the fusion of hydrogen. The net result is: 4 4p HE + 2 e + + 2 ν e + 3γ (34) Luminosity Gradient dl = ϵ dm (35) Here, ϵ is the total energy released per kilogram per second by all the star stuff: gravity and nuclear reactions together. Page 3
For spherically symmetric star, the mass of a thin shell of thickness dr will be: which leads to: dm = d M r = ρdv = 4πr 2 ρdr (36) dl r dr = 4πr 2 ρϵ (37) Energy Transport Radiative Radiation Pressure Energy density of Blackbody radiation: We can integrate this over all wavelengths/frequencies Thus we obtain: 4π dv = dv = c B v ( a is the radiation constant given by: 7.565767 6 J m 3 K 4 ) u v u λ 8πh ν 3 / c 3 u = dλ = dν The radiation pressure can be shown to be one-third of the energy density: (38) e hν/kt u ν 4π u = (T)dλ = = a c 4σT B 4 λ T 4 (39) c at 4 P rad = 3 (4) Page 4
spherica shell star center T T + dt r r + dr Fig. 2 A thin spherical shell centered on the star center. Let dr r The radiation pressure on the inner surface: at P rad (r) = 4 3 (4) The radiation pressure on the outer surface: a(t + dt) 4 a P rad (r + dr) = = ( + 4 ) 3 3 T 4 dt a ( + ) T 3 T 4 dt T 4 (42) (assuming dt/t ) Thus, the net radiation force will be given by: which using eqs. (4) and (42), we can get: F rad F rad = [ P rad (r) P rad (r + dr)]4πr 2 (43) a 4 4π = a dt 3 T 4 dt r 2 6π r 2 T 3 (44) T 3 Page 5
The optical depth is given by: dτ = ρ(r)κ(r)dr (45) where κ is the opacity. The probability that a photon will be absorbed while passing through the shell is dp dτ. The rate at which photons carry energy through the shell is the luminosity, L(r). Momentum is carried through the shell is just L(r)/c, since p = E/c. The force on the shell, i.e. the rate at which photon momentum is transferred to the shell will be F rad L(r) L(r) = dτ = ρ(r)κ(r)dr c c Combing eq. (46) equation with eq. (44), we get: 6π ar 2 T(r) 3 L(r) dt = ρ(r)κ(r)dr 3 c (46) (47) This can be rearranged to yield the equation of radiative energy transport: dt dr 3ρ(r)κ(r)L(r) = (48) 6πac T(r) 3 r 2 This equation related the luminosity to the temperature gradiant of a star. Convective Transport This is the equation of convective energy transport. dt dr T(r) dp = ( ) γ P(r) dr (49) It shows that the temperature gradient is proportional to the pressure gradient. ( γ is the adiabatic index and depends on the gas in question. ) Page 6
Fig. 3 The different transport modes are separated into regions. Stellar Model Building dp dr dm dr dl dr dt dr = GM r ρ r 2 = 4πr 2 ρ(r) = 4πr 2 ρ(r)ϵ(r) 3ρ(r)κ(r)L(r) = [for radiative] 6πac T(r) 3 r 2 T(r) dp = ( ) [for convective] γ P(r) dr These are the basic equations that show the relationship between the main attributes of a star. They are 4 differential equations that relate: pressure, mass, density, temperature, energy generation, luminosity, and the radius of the star. Boundary Conditions In order to solve these equations, we need boundary conditions: At the center we should expect: and at the surface of the star: M r L r } as r T P ρ as r R Page 7
. Build a star mesa http://mesa-web.asu.edu/index.html Page 8