Introduction To Stellar Evolution
Introduction To Stellar Evolution
This content has been downloaded from IOPscience. Please scroll down to see the full text.
(http://iopscience.iop.org/1742-6596/703/1/012002)
View the table of contents for this issue, or go to the journal homepage for more
Download details:
IP Address: 155.254.207.48
This content was downloaded on 08/08/2016 at 12:17
Abstract. This contribution is meant as a first brief introduction to stellar physics. First I
shortly describe the main physical processes active in stellar structures then I summarize the
most important features during the stellar life-cycle.
1. Introduction
The quantitative study of the physical structures called “stars” developed later with respects
of other physical topics, as electromagnetism or thermodynamics. This because to understand
the operation of stellar structures one must already be acquainted with fundamental physical
mechanisms as photon-matter interactions or energy production through nuclear fusion.
Moreover precise evaluations have been made possible only by the development of electronic
calculators. At the same time telescopes more and more refined have been developed, also
on board of satellites, to observe the whole electromagnetic spectrum without the limitations
induced by the terrestrial atmosphere. Now the characteristics of stars of different mass, chemical
composition and evolutionary phase are known with a quite good precision; the general scenario
is well defined and in agreement with the most of observations of the Sun and of the stars in
our Galaxy.
2. Stellar Observables
The main stellar observables are luminosity and color. The luminosity is the energy emitted
per unit time, but it is not directly measurable for two main reasons. In general, the adopted
detector (telescope) is not sensitive to the whole electromagnetic spectrum, but only to a given
spectral range, which depends on its characteristics. Most important, due to the fact that a star
isotropically emits its energy in all the directions, what is measured is a flux (energy per unit time
and unit surface of the detector) and not a luminosity which, however, is the physical quantity
directly linked to the stellar characteristics. In more detail, if L∗ is the stellar luminosity, a
detector at a distance d from the star measures the flux:
L∗
Φ= (1)
4πd2
The luminosity can be obtained only if a distance estimate is available; it’s thus easy to
understand why the problem of distance determination of astrophysical objects is of paramount
importance. The color of a star, that is the spectral distribution of the emitted energy, can
be directly connected with the temperature of the stellar surface (photosphere). In fact for
Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further distribution
of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Published under licence by IOP Publishing Ltd 1
8th European Summer School on Experimental Nuclear Astrophysics (Santa Tecla School) IOP Publishing
Journal of Physics: Conference Series 703 (2016) 012002 doi:10.1088/1742-6596/703/1/012002
4.0
3.0
RGB
log L/LSun
2.0
1.0
SGB PMS
MS
0.0
PMS
-1.0
3.80 3.75 3.70 3.65 3.60 3.55 3.50
log Teff
Figure 1. Balance between pressure
force (FP ) and gravity force (Fg ) on a Figure 2. Evolution in the Hertzsprung-Russell
given infinitesimal element at distance r diagram of a 1 M star from the Pre-Main
from the center of the star with tickness Sequence to the end of the Red Giant Branch
dr and base area dA. phase.
almost the whole stellar structure photon-matter interactions are so efficient that the star can
be assumed at thermodynamical equilibrium. In this case the gas particles follow a Maxwell-
Boltzmann distribution while the photons a Planck one, with the same temperature for particles
and photons. The temperature clearly varies with the position inside the star; in most cases the
temperature monotonic decreases from the center to the stellar surface. Under these conditions
the energy spectral distribution is the so called “black body distribution”. In reality the emission
of a star with a given surface temperature slightly differs from the one of a black body of the
same temperature, due to emission and absorption phenomena by the atoms/molecules/grains of
the stellar atmosphere. For a precise relation between the observed stellar color and the surface
temperature one needs so solve the photons transport equation in the stellar atmosphere, taking
into account all the possible absorption/emission phenomena. It’s easy to understand that these
calculations are very complex and they become more and more difficult as the photospheric
temperature decreases because, due to the higher presence of molecules and grains, the possible
photon-matter interactions increase more an more. These kind of calculations, which are an
example of how much computing capacity is fundamental for precise stellar physics evaluations,
are performed by few groups in the world and they are still affected by not negligible uncertainties
mainly for stars with relatively low photospheric temperature. In absence of high rotation
velocities and magnetic fields (which is supposed to be the situation of the most of the stars) the
stellar equilibrium configuration is a sphere (as we will better discuss in the following); however,
in general, it’s not possible to directly measure the stellar radius. Stellar surface temperatures
are generally expressed as “effective temperature” , that is the temperature which would have
a star if it would exactly emit as a black body. Thus the relation among luminosity (L), radius
(R) and effective temperature (Tef f ) is:
L = 4πR2 σTef
4
f (2)
where σ is the StefanBoltzmann constant, σ = 5.67 · 10−5 dyne cm−2 K−4 . Usually, observed
stars are shown in a plane which resumes the main stellar observables: effective temperature
(or correspondingly the stellar color) on the abscissa and total luminosity (or correspondingly
the luminosity in a given photometric band) on the ordinate, the so called Hertzsprung-Russell
(Color-Magnitude) diagram.
2
8th European Summer School on Experimental Nuclear Astrophysics (Santa Tecla School) IOP Publishing
Journal of Physics: Conference Series 703 (2016) 012002 doi:10.1088/1742-6596/703/1/012002
where X(r),Y(r) and Z(r) are, respectively, the fractional abundance in mass of hydrogen, helium
and elements heavier than helium, which in astrophysics are called “metals”, at a given radial
coordinate. Clearly X(r)+Y(r)+Z(r)=1.
3
8th European Summer School on Experimental Nuclear Astrophysics (Santa Tecla School) IOP Publishing
Journal of Physics: Conference Series 703 (2016) 012002 doi:10.1088/1742-6596/703/1/012002
where ρ(r) is the corresponding density inside the shell. This is called “continuity equation”.
If at each point of the star one equates the pressure force to the gravity force, as shown in
Fig.1, one obtain the so called “hydrostatic equilibrium equation”:
dP (r) M (r)
= −G 2 ρ(r) (5)
dr r
where P(r) is the pressure and M(r) is the mass inside a sphere of radius r; the above equation
holds at each radial coordinate r.
The hydrostatic equilibrium equation also implies that the most internal regions, which must
“support” the above layers, have an higher pressure. Thus, in general, the pressure is expected
to monotonically decrease from the center to the surface. The two equations written until now,
as a function of the radius, contains three unknowns: P(r), ρ(r) and M(r). P(r) and ρ(r) are
linked through the equation of state, but, except some particular EOS, in this case, also the
temperature (unknown) appears in the EOS. Thus the system is still not resolvable. This is
expected, because the energy production and transport, which play an important role in stellar
structures have still not taken into account.
4
8th European Summer School on Experimental Nuclear Astrophysics (Santa Tecla School) IOP Publishing
Journal of Physics: Conference Series 703 (2016) 012002 doi:10.1088/1742-6596/703/1/012002
regions in which these kind of neutrinos are produces cool down. Summarizing, in each region
of the star energy production through nuclear fusions and/or thermodynamical transformations
and/or neutrino energy losses can happen. Usually we define ε the energy produced/lost per
unit time and unit mass due to the mechanisms described above, that is:
Each process of energy production depends on the physical quantities (density/pressure and
temperature) and chemical composition which vary at each point of the star. In particular
the various nuclear reaction rates can be written as analytical expressions, as a function of
temperature, density and reagent abundances [see e.g. 5]. With these information it’s possible
to write the thermal equilibrium equation in each shell of the star:
dL(r)
= 4πr 2 ρ(r)ε(r); (7)
dr
where L(r) is the luminosity (energy per unit time) which enters in a shell of radius r. Another
equation is added to the system but another unknown, i.e L(r), is also added.
dT (r) 3 k̄(r)ρ(r) 1
=− Φ(r); (8)
dr 4 ac T (r)3
5
8th European Summer School on Experimental Nuclear Astrophysics (Santa Tecla School) IOP Publishing
Journal of Physics: Conference Series 703 (2016) 012002 doi:10.1088/1742-6596/703/1/012002
4. Stellar evolution
In the following I will describe the main stellar features during the different evolutionary
phases, as obtained by solving the stellar structure and evolution equations; some of the main
characteristics can also be roughly obtained by general evaluations discussed above.
6
8th European Summer School on Experimental Nuclear Astrophysics (Santa Tecla School) IOP Publishing
Journal of Physics: Conference Series 703 (2016) 012002 doi:10.1088/1742-6596/703/1/012002
approximation, possible effects due to rotation and magnetic fields, are in equilibrium between
the gravitational force and the pressure force. However, due to stochastic motions, the colder
and more dense regions of the cloud, can start to contract. As soon as the collapse proceeds the
region fragments more and more; until each part in which the original region of the cloud has
been divided will form a star. The way in which the gas of the original cloud falls onto the star
until the achievement of the stellar final mass is still an open problem, because hydrodynamical
protostellar models are still largely debated and not yet settled. The duration of the protostellar
phase too is not known with precision and it varies with the stellar mass, but it is supposed to
range from few 105 to about 106 years, [see e.g. 7, for details].
Anyhow, at the end of the protostellar phase the radiation pressure of the new born star
sweep away the residual gas of the original cloud. The star is now in the “Pre-Main Sequence
phase (PMS)”. Moreover during the first phases of their life stars are supposed to be fully
convective; thus they are completely mixed and possible chemical inhomogeneities of the original
cloud are canceled. During the protostellar phase the star contracted and warmed up until the
achievement of the hydrostatic equilibrium: the gas force exactly counterbalances the gravity
force. The stellar gas pressure is mainly due to its high temperature. However the star is warm
and radiates, loosing energy from the surface; thus, after a while, the star cools down and its
pressure decreases becoming insufficient to completely sustain the star against the gravity. The
star slightly contracts increasing its temperature and reaching again the hydrostatic equilibrium
until the radiation losses from the surface cools down again the star. The PMS phase is
characterized by these consecutive stellar contractions when the star passes from one “quasi-
equilibrium” state to another one, while its internal temperature increases. The timescale of
these processes is the one at which the stars cools down, that is the timescale in which a photon
produced in the stellar interior can reach the surface. This timescale that can be seen also
as the “relaxation time” for departure of a star from the thermal equilibrium (the so called
“Kelvin-Helmholtz” time) varies with the stellar mass but it is of the order of 106 -107 years.
While the star contracts its density increases too and the radius significantly decreases, then
also the luminosity strongly decreases as shown in Fig.2. During the PMS evolution a star reaches
first the temperature for deuterium and light elements (Li, Be, B) burning, then it reaches the
temperature for hydrogen burning into helium in its core, which is the warmest region of the
star. The hydrogen fusion into helium is the first nuclear reaction able to produce an energy
which can replace the radiative losses from the surface. Hydrogen is the first abundant element
to be burned because the coulombian barrier for the fusion of two protons is the lowest one and
H is the most abundant element in the Universe. When the hydrogen fusion is fully active the
star is in equilibrium, with the nuclear reactions which supply to the star the energy lost by
photon emission from the photosphere. Sometimes one hears on TV broadcasts that the stellar
luminosity depends on the efficiency of the nuclear reactions, but actually the exactly opposite
is true. Stars are luminous because are warm and, as all the warm bodies, they irradiate. As
a first approximation, the stellar temperature depends from the fact that the stellar gas must
be warm to have enough pressure to counteract the gravity force; then the surface temperature
also depends on the temperature gradient inside the star (see eq.8). Thus, roughly speaking,
the higher is the mass of the star, the higher is its effective temperature and its luminosity.
When, at the end of the PMS phase, hydrogen fusion ignites, it adjust its efficiency until the
energy losses from the surface, proportional to the effective temperature and the stellar radius,
are completely counterbalanced. This because in stars a feedback mechanism is efficient: if
the produced energy is lower than radiative losses the star still cools down and contract, even
if with a longer timescale. The result is the increase of the temperature of the central regions
together with the burning efficiency, which is very sensitive to the temperature. On the contrary
if the produced energy is higher than the one lost from the surface, the central regions warm
up and expand; this leads to a temperature decrease and thus to a reduction of the burning
7
8th European Summer School on Experimental Nuclear Astrophysics (Santa Tecla School) IOP Publishing
Journal of Physics: Conference Series 703 (2016) 012002 doi:10.1088/1742-6596/703/1/012002
efficiency. At the end of the process the energy produced by nuclear reactions is adjusted to
exactly counterbalance the surface energy losses. Thus, in this phase, the stellar luminosity is
mainly due to the stellar mass (with a smaller dependence on the chemical composition) and the
nuclear reaction efficiency is tuned to exactly produce the energy needed to the star to irradiate.
When nuclear reactions are so efficient to counterbalance the energy losses from the surface,
the star leaves the PMS phase starting its Main Sequence (MS) evolution. The name “Main
Sequence” is due to the fact that stars in this phase occupy in the Hertzsprung-Russell diagram
a mono-parametric sequence (neglecting the secondary effect of the chemical composition) whose
parameter is the stellar mass. It’s worth noticing that stars with masses lower than about 0.08
M do not succeed to reach the central temperatures needed to hydrogen ignition. They remain
in the Pre-Main Sequence phase continuing to contract while the density increases until the star
is sustained by the pressure of the degenerate electrons; these kind of stars are called “brown
dwarfs”. The star is now in equilibrium thanks to the nuclear energy, but the reaction efficiency
decreases according to the reduction of the fuel (hydrogen) abundance. Thus, after a while, the
nuclear energy is no more enough to counterbalance the radiative losses and the star contracts
again, with the timescale of the fuel abundance reduction, that is of the nuclear burning rate.
Thus the temperature of the central regions slightly increases until the energetic equilibrium
is reached again, with a slightly different configuration of the stellar structure; that is the star
evolves with nuclear timescale.
Nuclear timescales strongly vary with the stellar mass; the higher is the mass, the higher
is the luminosity and thus the higher is the rate of the nuclear reactions which counterbalance
the energy losses due to the luminosity. Moreover the higher is the fusion rate, the lower is the
time in which the fuel is exhausted. To be more precise for the most of the stars the luminosity
increases with the mass to a power of 2÷4, while the H abundance linearly increases with the
mass. As an example, the MS time for the Sun is about 10 Gyr while for a 15 M star of the
same chemical composition is about 107 years.
8
8th European Summer School on Experimental Nuclear Astrophysics (Santa Tecla School) IOP Publishing
Journal of Physics: Conference Series 703 (2016) 012002 doi:10.1088/1742-6596/703/1/012002
low density envelope; their luminosity and lifetime are mainly governed by the dense He core and
the surrounding H burning shell, while the temperature gradient in the envelope determines the
stellar effective temperature. In particular the RGB luminosity increases with the He core mass
which grows due to the burning in the H shell which surrounds the He core. The RGB phase
ends when the central regions of the star reach the temperature for helium burning into C and
O (≈ 108 K), thanks to the He core contraction due to its mass growth. Stars with masses lower
than about 2.3 M (the so called “RGB transition” mass), whose interiors did not reach very
high temperatures, remain longtime in the RGB phase because the reaching of the He burning
temperature is more difficult. These stars reaches very high luminosities along the Red Giant
Branch and thus the mass losses by stellar winds, driven by the radiation pressure, increase too.
At the end stars with masses greater than about 0.5 M succeed to ignite central He burning.
The evolutionary theory predict that lower masses will remain in the RGB phase loosing more
and more mass until the only remain will be the dense, degenerate He core, surrounded by a
very thin H shell with all the nuclear burnings switched off. A similar structure, called He
white dwarf (He WD), is maintained in equilibrium by the pressure of degenerate electrons
which counteracts the gravity and it will cools down until it will reach the cosmic background
temperature. However the Main Sequence lifetime of stars with M < 0.5 M is higher than
the present age of the Universe, thus no one of these stars still reached the white dwarf phase.
However He white dwarfs are observed; their progenitors are stars of mass higher that 0.5M ,
which during the RGB phase loose so much mass that the H shell switches off and the stars
leave the RGB phase ending their life as He WDs. Stars which ignite He in the center but did
not explode as supernovae end their life as white dwarfs with a C/O core or with a Ne core if
they succeed to ignite carbon too.
Acknowledgments
I would like to acknowledge the collegues with which I am working in Pisa from many years in
the stellar physics field: M. Dell’Omodarme, P.G. Prada Moroni, E. Tognelli and G. Valle.
References
[1] Cox J P 1968 Principles of stellar structure - Vol.1: Physical principles; Vol.2: Applications
to stars (New York: Gordon and Breach)
[2] Hansen C J, Kawaler S D and Trimble V 2004 Stellar interiors : physical principles, structure,
and evolution (New York: Springer-Verlag)
[3] Rolfs C E and Rodney W S 1988 Cauldrons in the cosmos: Nuclear astrophysics (Chicago,
IL: University of Chicago Press)
[4] Raffelt G G 1996 Stars as laboratories for fundamental physics : the astrophysics of
neutrinos, axions, and other weakly interacting particles (Chicago: University of Chicago
Press)
[5] Fowler W A, Caughlan G R and Zimmerman B A 1967 Annual Review of Astronomy and
Astrophysics 5 525
[6] Salaris M and Cassisi S 2006 Evolution of Stars and Stellar Populations (Wiley-VCH)
[7] Stahler S W and Palla F 2005 The Formation of Stars (Wiley-VCH)