Part1 PDF
Part1 PDF
Part1 PDF
Stellar Structure
and
Evolution
Astr 323/Phys 323
Plan:
1. Structure equations
2. Microphysics
∂ρ 2
−2 ∂(ρr v)
= −r
∂t ∂r
∂ρ
→ = −∇(ρv)
∂t
∂r 1
= (1)
∂m 4πr 2 ρ
Gravitational field
∇2Φ = 4πGρ
(G = 6.673 · 10−8 dyn cm2 g−2).
In spherical symmetry:
µ ¶
1 ∂ ∂Φ
r2 = 4πGρ
r 2 ∂r ∂r
∂Φ Gm
g= ∂r
→ g= r2
is solution of Poisson’s equation.
Potential
R∞ Φ vanishes for r → ∞.
− 0 Φdr is the energy required to disperse the com-
plete star to infinity.
Hydrostatic equilibrium
∂P Gm
=− 2 ρ
∂r r
in Lagrangian coordinates:
∂P Gm
=− (2)
∂m 4πr 4
2GM 2
Pc ≈
πR4
(M/2 and R/2 were used for mean mass and radius)
∂P Gm 1 ∂ 2r
+ =−
∂m 4πr 4 4πr 2 ∂t2
5. Examples: τhydro =
• 27 minutes for the Sun
• 18 days for a Red Giant (R = 100R⊙)
• 4.5 seconds for a White Dwarf (R = R⊙ /50)
R
P = ρT
µ
R 2
= cP − c v = cv
µ 3
2. Gravitational energy
Z M
GMr GM 2
Eg = − dMr ≈ −
0 r R
Generally,
−Eg ≈ Et
Why?
The Kelvin-Helmholtz timescale
¯ ¯
|Eg |
L ≈ ¯ dE → τ := ≈ ELt .
¯ g¯
dt ¯ KH L
GM 2 GM 2
|Eg | ≈ 2R
⇒ τKH ≈ 2RL
.
Z M Z M
∂P Gm
4πr 3 dm = − 4πr 3
dm
0 ∂m 0 4πr 4
Z M Z M µ ¶
3 ∂P ∂r
3
¤M
12πr 2
£
4πr dm = 4πr P 0 − P dm
0 ∂m 0 ∂m
So finally:
Z M
P
3 dm = −Eg
0 ρ
Meaning of left hand term:
With this P
ρ
= 32 cv T = 23 u
R
u: specific internal energy ⇒Ei := udm and
Z M
P
3 dm = 2Ei
0 ρ
In summary, we obtain the Virial Theorem
⇒ Eg = −2Ei (3)
ζEi + Eg = 0.
L = − dW
dt
⇒ L = (ζ − 1) dE
dt
i
.
Ėg
ζ=2 ⇒ L=− = Ėi , (4)
2
Interpretation:
star loses energy → decrease in gravitational energy
(contraction) → but same amount goes into increase
in internal energy (heating)
~ : nu-
• in a stationary mass shell: ǫ = ǫn (ρ, T, X)
clear energy generation;
µ ¶
∂Lr
ǫn − dt = dq = du + P dv
∂m
∂Lr ∂u P ∂ρ
= ǫn − + 2
∂m ∂t ρ ∂t
∂T δ ∂P
= ǫn − cP +
∂t ρ ∂t
ǫg : gravothermal energy
µ ¶
∂T δ ∂P 1 ∂T ∇ad ∂P
ǫg = −cP + = −cP T − (5)
∂t ρ ∂t T ∂t P ∂t
∂Lr
= ǫn + ǫg − ǫν . (6)
∂m
Global energy conservation
L = ∂m dm, Lν = ǫν dm, ǫn dm = − dE
R ∂Lr R R
dt
n
.
M
1 ∂P dEG
Z
−3 dMr = 4
0 ρ ∂t dt
Therefore
Z M
ǫg dM = −Ėt − E˙G
0
τn := En/L
⇒ τn = 7 · 1010 yrs
τn ≫ τKH ≫ τhydr
∂T T Gm
=− ∇
∂m P 4πr 4
Transport by radiation
~ (cm2 /g).
Opacity κ(T, ρ, X)
~ = − 1 vlp ∇n
~j = −D∇n ~
3
(D is called the diffusion constant; v the diffusion ve-
locity; lp is the particle free path length and n the
particle density).
∂U ∂T
= 4aT 3
∂r ∂r
and for the radiation flux F (replacing ~j)
4ac T 3 ∂T
F =− ,
3 κρ ∂r
or F = −Krad ∇T .
4ac T 3
Krad = 3 κρ
is the radiative conductivity.
Its value for the solar center is ≈ 3×1015 erg K−1 s−1 cm−1 ,
compared to values of 107 for typical metals or 104
for concrete.
∂T 3 κρLr
=−
∂r 16πac r 2T 3
or, in Lagrangian coordinates
∂T 3 κLr
=− (7)
∂m 64acπ 2 r 4T 3
Stefan-Boltzmann law:
L = 4πσR2 Teff
4
(17)
This is the second relation for the outer b.c., and the
system is closed and can be solved.
2. Adiabatic stratification:
T = const P ∇ad (in case of convection)
For ideal gas, it follows that
1−∇ad
P = const2 ρ1/(1−∇ad ) and → n = ∇ad
.
For ∇ad = 0.4 → n = 3/2.
The left hand side is the Euler version of eq. (9) and
we obtain Poisson’s equation.
Using the polytropic equation, the derivative of P is
dP 1 dρ
= K(1 + )ρ1/n
dr n dn
and using this in Poisson’s equation yields
µ ¶ µ ¶
n+1 1 d 1
2 n −1 dρ
K r ρ = −4πGρ
n r 2 dr dr
d2 y
µ ¶
1 d dy 2 dy
x2 = + = −y n (19)
x2 dx dx dx 2 x dx
Discussion of solutions of Emden’s equation –
Lane-Emden-functions
x0
dP dfP dx Pcx d ln fP 1 Pc P M d ln fP P d ln fP
= PC = = =
dm dx dm fP d ln x M Pc mM d ln x m d ln x
dP P d ln fP Gm P m
= =− → ∼ 4 (23)
dm m d ln x 4πr 4 m r
dr r d ln fr 1 r 1
= = → ∼ 2 (24)
dm m d ln x 4πr 2ρ m r ρ
dT T d ln fT 3κ Lr T Lr
= =− → ∼ 4 3(25)
dm m d ln x 64πac r 4T 3 m r T
dLr Lr d ln fL Lr
= =ǫ → ∼ǫ (26)
dm m d ln x m
From (23) and (24) expressions for P and ρ as func-
tions of r and m can obtained; taking the ratio of these
and using the same ratio from the ideal gas EOS it
follows that
P m T
∼ ∼
ρ r µ
(or rT = µm) and with equation (25) for Lr
⇒ Lr ∼ µ4 m3
in particular, for x = 1 or m = M , we have
L ∼ µ4 M 3 ,