Download PC3692: Physics of Stellar Structure (and Evolution)

PC3692: Physics of Stellar Structure (and Evolution)
The goal of this course is to apply fundamental physics laws to understand the internal
structure and evolution of stars. In this course, you will employ virtually all the physics you have
learned so far, in particular thermal physics, quantum mechanics and nuclear physics. You will
understand how energy is generated inside stars and transferred outwards. In short, you will
understand how stars shine. The lectures will be divided into six chapters.
1. Observed properties of stars
2. Equations of stellar structure
3. Equation of state
4. Energy transport in stars
5. Nuclear fusion in stars
6. Applications of stellar equations
1.
Observed properties of stars
1.1.
Why do we study stars?
• They are beautiful and seemingly eternal
Twinkle, twinkle, little stars,
how I wonder what you are,
Up above the world so high,
like a diamond in the sky
– Jane Taylor (1806)
A sense of wonder is universal in all civilisations. We wonder what the stars are, as in the
poem. Although we have observed stars for millions of years, they are really physically
understood only in the 20th century, because, as we will see in the course, it requires
knowledges of nuclear physics. Such information has only become available in the 20th
century.
However, the long history in astronomy means that some of the terminologies in astronomy
are peculiar, such as apparent magnitudes and metallicities.
• Stars, in particularly the Sun, plays a crucial role in our lives
Stars → nuclear reaction → Energy+weather (seasons) → life
Stars → synthesise elements (C, O, N) → found in air and our human bodies
• Stars are also interesting astrophysically for at least three reasons:
1. Stars → star clusters → galaxies → universe. Stars to the cosmology is like as atoms in
physics. To understand the whole universe, it is crucial to understand stars.
–2–
2. Stars provide extreme conditions for studying physics. The high density, high temperature
and high pressure are usually not found on Earth. White dwarfs, neutron stars are extremely
dense and are a good testing ground for theories such as general relativity.
3. Stars are fascinating for yet another reason: while the general principle of stars is
well-understood, there are still many puzzles in stars, such as the solar neutrino problem,
star formation and convection; the latter two are long-standing unsolved problems in
astrophysics.
1.2.
Stellar parameters
How do we describe stars? In fact, the internal properties of stars can be primarily described
by just a few parameters, M, R, T and chemical composition (the fractions of different elements,
such as H, He, inside a star), and a related concept luminosity, L.
The stars also have several important external parameters, such as distance and their motions
in space. The distance measurement is essential for determining the internal properties of the
stars, so we discuss it first.
1.2.1.
Distances
Th distance of stars can be measured with a number of different methods. The most
fundamental one, parallax, is also the simplest; it is a purely geometric method (see Fig. 1.2.1).
The viewing angle of the stars from Earth changes due to the Earth motion around the Sun.
The parallax angle is usually quite small. The nearest star has a parallax of about 1 arcsecond.
The distance is then simply
d = 1 AU/π, AU = 1.5 × 1011 m
(1)
When the parallax is one arcsecond, the corresponding distance is called parsec (pc). Since 1
radian=180/π × 60 × 60 ≈ 2 × 105 arcsecond. This means
1 pc = 2 × 105 AU ≈ 3 × 1016 m
(2)
Fig. 1.— Schematic diagram of parallax. The parallax angle is usually denoted as π (labelled as ρ
in the figure)
–3–
1 pc ≈ 3.26 ly
(3)
Example 1.1 A star at 500pc will have a parallax of 1/500=0.002 arcsecond.
The recent satellite Hipparcos can measure parallax as accurate as 0.002 arcsecond, so we can
roughly the star measure distances out to few hundred pc. Future satellite missions, such as GAIA
will be able to measure parallax to even higher accuracies, hence they can measure distances to
all bright stars in our Galaxy.
Astronomers also frequently use kpc=103 pc, Mpc=106 pc, Gpc=109 pc.
1.2.2.
Luminosity of stars and related concepts
Stars radiate and we see the photons from them. The luminosity of a star is defined as the
energy radiated per second,
energy
Luminosity =
(4)
second
Suppose we have a solar panel with perfect efficiency receiving light from stars like the Sun, then
we can measure how much energy is received per square meters per second, which is usually called
flux:
energy
flux =
(5)
second m2
The luminosity is related to the flux by
L = 4πd2 f,
(6)
where d is the distance to the star. And so if we know the distance to a star, then from the
observed flux, its luminosity can be determined.
For our sun, L = 3.8 × 1026 W , including photons in all wavelengths. The solar energy
radiated in the “visual” part of the electromagnetic radiation, L,V , (i.e., at a wavelength of about
5500Å) is about 4.64 × 1025 W , while that for the blue (‘B’) part of the electromagnetic radiation
(i.e., λ ∼ 4000Å) is about 4.67 × 1025 W .
In astronomy, for historical reasons, luminosity is frequently measured in so-called absolute
magnitudes. For example, the B-band and V-band magnitudes of the Sun are ‘defined’ as
M,B = 5.48 and M,V = 4.83, respectively. For other stars, we have
MB = M,B − 2.5 log(LB /L,B ) = 5.48 − 2.5 log(LB /L,B )
(7)
MV = M,V − 2.5 log(LV /L,V ) = 4.83 − 2.5 log(LV /L,V )
(8)
Example 1.2 If LB = 100L,B , then we have
MB = 5.48 − 2.5 log(100) = 0.48
Note that the brighter a star is, the smaller its absolute magnitude!
(9)
–4–
The apparent brightness of a star obviously depends on its distance to us, the closer a star
is to us, the brighter it appears. Astronomers therefore use the so-called apparent magnitude to
describe the brightness of a star.
The apparent magnitude of a star is related to the absolute magnitude by
mB = MB + 5 log(d/10 pc),
(10)
mV = MV + 5 log(d/10 pc),
(11)
mB and mV are also commonly written as B and V . One can show that the apparent magnitudes
are directly related to flux, and therefore it can be measured without knowing the source distance.
The absolute magnitude (luminosity) of a star, in most cases, can only be determined if we know
the distance to the source.
Example 1.3 The distance to the Sun is d = 1 AU, so
B = 5.48 + 5 log(1 AU/10 pc) = −26.08
(12)
Again, the brighter a star is, the smaller its B magnitude is.
1.2.3.
Mass and Radius of stars
Masses and radii of stars are difficult to measure. Most measurements are from studies of
binary stars; most stars (>
∼ 80%) are in binaries. In this subsection, we discuss the procedure for
a special class of binaries, the so-called eclipsing binaries.
In eclipsing binaries, the brightness declines periodically as one stars passes in front of the
other. Of course, for an eclipse to occur, the orbit of the binary must be nearly edge-on. Fig. 3
shows one sample geometrical configuration. The principle of mass determination can be seen the
clearest if the orbits are circular, which we will assume in the following discussion.
From the observations, we can determine the orbital velocities from Doppler shifts, v1 and v2 .
From these two velocities, we can measure the mass ratio.
Fig. 2.— Schematic orbits of a binary in circular orbits. Various quantities are defined. In general
the orbits will be elliptic.
–5–
Fig. 3.— Schematic diagram showing the orbits and corresponding light curves of eclipsing binaries.
M1
v2
M1
=
→
M2
v1
M2
(13)
Exercise 1.1 Treat the solar system as a binary with only the Sun and Jupiter. Estimate the
velocity of the Sun around the centre of mass. Assuming the measurement accuracy in velocity is
5 m s−1 , can the solar motion be detected? From this, discuss what kind of extrasolar planets are
the easiest to be detected in terms of their distances to the parent star and their masses. Compare
your answer with the observed ones at http://exoplanets.org/
Now from the observed period, P , and the two velocities, we can infer the two orbital radii,
a1 and a2 , and the separation, using
v1 =
2πa1
,
P
v2 =
2πa2
, a = a1 + a2 ,
P
(14)
so we can infer the separation, a.
From Kepler’s third law, we have
P2 =
4π 2 a3
→ M1 + M2
G(M1 + M2 )
(15)
So eqs. (13) and (15) can be combined to determine the two masses.
Clearly the time duration of the eclipses is related to the stellar radii. Only a detailed analysis
of the eclipse shape and duration can give the stellar radii. In fact, the exercise is often quite
involved. Notice the depth of the eclipse is related to the relative brightness of the two stars.
Let us suppose a low-mass, tiny star rotates around a very massive giant star, as depicted in
Fig. 3. The separation between the stars is assumed to be much larger than both stellar radii.
The massive star is essentially stationary. Let us denote the relative rotation velocity of the less
–6–
massive star around the giant companion as v; in our case, v = v1 + v2 . For this special case, the
total duration (t) of the eclipse will be related to the period by
t
2R1 + 2R2
=
,
(16)
P
2πa
while the duration of the descent from the maximum light to the minimum light is related to the
diameter of the small star by
2R1
t1 =
(17)
v
We can determine t1 and t2 from the observed light curve, and combined with the measured
velocity v, we can find the stellar radii for both stars from the previous two equations.
1.2.4.
Effective temperature and Chemical Composition
The energy distributions of a star as a function of wavelength (or frequency) is called its
spectrum. The spectrum can be regarded as the finger prints of the star. The spectra of stars can
be crudely approximated as a blackbody. In fact, in astronomy it is common to define an effective
temperature of a star. This is the temperature that a star would have if it were a perfect black
body and had identical luminosity and radius as the real star.
Now the energy per unit time per unit area radiated by a blackbody is given by
4
f = σ Teff
(18)
4
L = 4πR2 × σ Teff
(19)
Therefore, we have
So the radius of a star, the effective temperature and luminosity are related, and if we can measure
two of them, then we can determine the remaining third parameter.
Fig. 4 illustrates the spectrum of the Sun. The top panel shows a low-resolution spectrum
while the right panel shows a high-resolution spectrum between 4660Åand 4990Å. As one can see,
the solar spectrum is approximately a blackbody of temperature 5770 K. The absorption lines in
the spectra are caused by photons being absorbed by atoms in the stellar atmosphere when they
emerge from the solar interior. The number and depths of the absorption lines depend on the
element abundances, i.e., the chemical composition. Detailed analysis of stellar spectra can yield
the effective temperature, the chemical composition of stars and ‘surface gravity’ (= GM/R2 ).
The effective temperature of a star can also be crudely measured using the so-called colour
index. A colour index is simply the magnitude difference in two colours (or filters), e.g., B − V .
Note that the colour index is independent of the source distance. Very hot stars are blue, which
may have B − V = −0.3, while cooler stars are red and may have B − V = 2.
Example 1.4 For the Sun, the B − V colour index is B − V = MB − MV = 5.48 − 4.83 = 0.65
Stars can be divided into various classes according to their spectra. The sequence is historically
termed as O, B, A, F, G, K, M, R, N1 . Each type corresponds to a certain range of effective
1
To be remembered as the initials of Oh, Be A Fine Girl (Guy), Kiss Me, Right Now
–7–
Fig. 4.— The solar spectrum in low-resolution (top) and in high-resolution (bottom). In the top
panel, the smooth dashed curve is a blackbody with a temperature of about 5770 K. There are
many dips in the spectra because photons are absorbed by atoms in the stellar atmosphere. The
continuum in the high-resolution spectrum has been renormalised to unity. It resolves many more
absorption lines. A careful analysis can yield many important stellar parameters, such as the Sun’s
chemical composition, the effective temperature and the ‘surface gravity’ (= GM/R2 ).
–8–
temperature and colour index. From O stars to M stars, they become cooler (temperature from
40,000K drops to 2,800K) and redder (with B − V changes from −0.35 to 1.6). In comparison, the
hottest white dwarfs have Teff ≈ 3 × 105 K, while the hottest neutron stars have Teff ≈ 3 × 107 K.
1.3.
Range of stellar parameters
Now we know how to measure the stellar parameters. What are the typical values of the
stellar parameter values?
Parameter
Radius
Mass
Teff
Luminosity
Chemical composition
Sun
R = 7 × 108 m
M = 2 × 1030 kg
Teff, = 5770K
L = 3.8 × 1026 W
Z = 0.02
Stars
10−2 − 102 R
10−1 − 102 M
103 − 105 K
10−5 − 106 L
∼ 10−3 − 5 Z
Comments:
• Sun is a typical star, with intermediate parameter values.
• Mass has a lower cutoff and upper cutoff. This is because lower mass objects can not initiate
nuclear burning while higher mass stars are unstable.
• The luminosity range varies by about 11 orders of magnitude, much larger than any other
parameters. Some stars are extremely luminous, while others are very faint. We will learn
why this is the case.
1.4.
Properties of Stellar Populations
So far, we have discussed how to measure the parameters for individual stars. In the following,
we will discuss the statistical properties of stellar populations. We can often gain insight into the
astrophysics by examining a large sample at the same time, rather than individually. We shall also
discuss what makes stars to have different properties, such as their temperatures and luminosities.
Spatial distribution of stars
Stars in the universe are not uniformly distributed in space. They are highly clustered. It is
still not completely understood how this highly clustered pattern arises.
Milky Way is a typical disk galaxy (see Fig. 5). It consists a bulge and a thin disk. The
disk is thin because it’s rotationally supported, like an audio CD. You also see black lanes in the
galaxy, these are caused by dust obscuration.
The diameter of our luminous Milky Way is about 30 kpc. The Sun is about 8 kpc off center.
In the picture you can also see some star clusters, called globular clusters. They are spherical
–9–
distributed around the Galactic center. They consist off 105 − 106 stars, and they play an
important role in our understanding of stellar structure and evolution.
1.4.1.
Mass-Luminosity relation
The first relation we examine is that between the mass and luminosity (obtained for nearby
stars). This is shown in Fig. 6. In this figure, both the luminosity and mass are plotted on
logarithmic scales.
Fig. 5.— An image of the Milkyway in the optical (top panel). The bottom panel shows the
approximate linear scales (in light years) of the Milky Way. Notice that the globular clusters are
roughly spherically distributed around the Galactic centre.
Fig. 6.— The Mass-luminosity relation for nearby stars. The scales on both axes are logarithmic.
– 10 –
The data points lie approximately on a straight line. This line is approximated by
2
L ∝ M 3.5
(20)
Since the mass of stars vary from 0.1 to 100 solar masses, we immediately see that the luminosity
of main sequence stars will vary by about 10 or 11 orders of magnitude.
This also has profound consequences for the life times of stars. The total energy available for
a star is, from Einstein’s famous formula, proportional to E = M T 2 . We therefore have
t∝
M c2
M
∝ 3.5 ∝ M −2.5
L
M
(21)
A more massive star has a larger energy reservoir, however, it also consumes its energy at a much
faster rate. As a result, the life time of a massive star is shorter.
More precise modelling gives
M
t ≈ 10 yr
M
10
−2.5
(22)
For the Sun,
t ≈ 1010 yr
(23)
For a star with mass below the Sun, the life is then comparable to the age of the universe,
tuniverse ≈ 1.3 × 1010 years. In other words, very low mass stars evolve very slowly, and they do
not evolve beyong hydrogen burning.
However, for a star with M = 10M , we have
t ≈ 1010 × 10−2.5 = 107.5 = 3 × 107 yr tuniverse
(24)
Since the life time is much shorter than the age of the universe, so these stars could not have
formed at the beginning of the universe. Instead, they must have formed very recently, about tens
of millions years ago. So stars are not eternal! they are forming and dying today. In fact, they are
just like human beings, they have births, they grow, evolve, and eventually die.
Unfortunately, how a star is borne is not yet well understood. Star formation is currently one
of the most active research areas in astrophysics.
1.4.2.
Hertzsprung-Russell diagrams
Hertzsprung-Russell diagram for the solar center
In Fig. 7, we show the Hertzsprung-Russell diagram for 41453 nearby stars with accurate
distances measured by the Hipparcos satellite. The Hertzsprung-Russell diagram is also often
called HR digram or colour-magnitude diagram.
In the figure, the horizontal axis plots V − I, a colour index similar to B − V that we discussed
in §1.2.4 The stars become hotter and bluer as you move from the right to the left. The vertical
axis plots the absolute magnitude in the Hipparcos filter. The stars become more luminous as you
move from the bottom to the top.
2
Rigorously, this relation only applies to massive stars; for low-mass stars, this relation has to be modified.
– 11 –
Fig. 7.— HR diagram for 41453 nearby stars with accurate distance measured by the Hipparcos
satellite. The horizontal axis is the V − I colour index, while the vertical axis is the absolute
magnitude in the Hipparcos passband. The I-band is a filter centred around 8000Å.
One striking feature is there is a sequence of stars running from the top left to the bottom
right. This sequence is called main sequence.
You also see a clump of to the right of the main sequence, these stars are called red clump
stars, and the stars further to the right, red giants.
You can also vaguely see some stars in the bottom left; these are white dwarf stars, they are
hot and very faint.
90% of stars belong to the sequence, and it appears that one single parameter, M , determines
the position of a star on the HR diagram. There are about 10% of stars that do not belong to
the main sequence. One important goal of this course is to understand why most stars are on the
main sequence, while other stars do not, and what are the physical differences between different
types of stars.
HR diagram for globular clusters
We already discussed that stars are not uniformly distributed in space, they are highly
clustered. Globular clusters are very dense star clusters with 105 − 106 stars with a size of a few
pc. In Fig. 8 we show one example, M5. The stars in such clusters move with typical speeds of
∼ 5 km s−1 . It is yet not well understood how globular clusters form.
The right panel shows the HR diagram for M5. The x-axis now shows the colour index B − V
while the vertical axis shows the apparent V magnitude. Again the temperature increases from
the right to the left while the luminosity increases from the bottom to the top.
– 12 –
Fig. 8.— Image (left) and HR diagram (right) for globular cluster M5.
There are a few key features, these are main sequence (labelled as A), red giant branch
(labelled as B), and then the horizontal branch stars (labelled as D). The horizontal branch stars
are so named because they are nearly horizontal in the HR digram. Notice in a globular cluster,
we do not have very hot and luminous stars.
The stars in (most) globular clusters are thought to have formed at the same time and a long
time ago, in fact, about 10 billion years ago, about the age of the universe. They are also thought
to have the same chemical compositions.
The fact that they are old offers a crucial clue to understand their HR diagram. As we have
seen, the high mass stars are luminous and blue, so they have short life times. These stars would
have disappeared (more precisely, evolved away) from the main sequence. In fact, as we will
discuss, these stars evolve to form the red giant and horizontal branch stars in the HR diagram.
A key goal in stellar astronomy is to understand the fine features in the HR diagrams of
globular clusters. Different globular clusters have different HR diagrams, and these differences are
thought to arise mainly from the difference in ages. So a detailed study of the HR diagram in
globular clusters can also be used to determine the age of these clusters.
HR diagram for open clusters
There are also another kind of cluster of stars in galaxies, these are called open clusters. They
have only hundreds of stars. They are a few pc across, similar to the size of globular clusters. Fig.
9 shows an example, the Pleiades open cluster. As you can see from this figure, there are bluish
stars in the figure.
The right panel shows the HR diagram for this cluster. This is very different from the globular
cluster M5, nearly all stars are on the main sequence, including some hot and blue giant stars at
the top left.
Fig. 9.— Image (left) and HR diagram (right) for the open cluster Pleiades.
Open clusters are, like globular clusters, thought to have formed essentially at the same time
and with the same chemical composition.
Why are stars in open clusters lie essentially on the main sequence? The key of understanding
is very simple: they are thought to be young. So they still do not have time to evolve off the main
sequence. That’s why all the stars essentially still lie on the main sequence.
1.4.3.
Summary
We have studied the mass-luminosity relation and HR diagrams for the solar neighbourhood,
open and globular clusters. HR diagrams play an essential role in the study of stellar structure
and evolution, and in fact also for cosmologies. So a clear understanding of the basic features in
HR diagrams is very important.
Problem set 1
1. The binding energy per nucleon for 56 Fe is 8.8 MeV per nucleon. Estimate the total energy
released per kilogram of matter by the sequence of reactions which fuse hydrogen to iron.
2. The main sequence of the Pleiades cluster of stars consists of stars with mass less than 6M ;
the more massive stars have already evolved off the main sequence. Estimate the age of the
Pleiades cluster.
3. Given that the luminosity of the sun is 3.8 × 1026 W and that the absolute bolometric
magnitude of the sun is M ≈ 4.7, estimate the distance at which the sun could just be seen
by the naked eye (The naked eye can detect a star of apparent magnitude 6). The apparent
magnitude, m, is related to the absolute magnitude, M , by
m = M + 5 log
d
.
10pc
(25)
Estimate the number of photons incident on the eye per second in this situation.
4. Sketch the HR diagram for an old globular cluster, label the main sequence, the red giant
and horizontal branch. (2001 final exam, 3 marks)
Write down the approximate observed range of stellar masses. What imposes the lower and
upper limits on the mass of stars? [2001 final exam, 5 marks]
Using the observed mass-luminosity relation for nearby stars, argue stars must be forming
at present day. [2002 final exam, 5 marks]