Nuclear fusion in stars and laboratories
- We observe stars in the sky. They are (to order of magnitude)
objects like our sun, with characteristic mass of order
kg, characteristic size of
m (well, with a
range of 10
or so either way), and characteristic luminosity of
W (again with a large range).
- Stars are observed to be constant (nearly, anyway). Since
the dynamical time scale is of the order of the free-fall time,
stars are almost always in dynamical equilibrium.
- Stars are supported against gravity (until they grow old) by gas
pressure. From the equation of hydrostatic equilibrium,
where
is pressure,
is the radial coordinate from the centre of
the spherical star, and
is the local gravity due to
the enclosed mass
at
, we may estimate the order of magnitude
of the central pressure. Replace the derivative
and estimate
. Solve for
which is
of order
Pa. Similarly we estimate the mean
density to be of order
kg m
(about the same as condensed matter).
- The support pressure in visible stars, in spite of their rather
high density, is gas pressure (the atoms are highly ionized, and so
they are considerably smaller than
). We may estimate the
temperature using the ideal gas law and the fact that - from
spectroscopy - we find that stars are overwhelmingly composed of H
(i.e.
): using
where
is the number density of ions and
for H, we find
ions m
. From
, and the fact that each ionized H contributes one
electron to
, we find
K as the
characteristic interior temperature for a star. All of these estimates
are confirmed by detailed modelling.
- With such a high internal temperature, a star leaks heat (for
the sun at the rate of
W). Without any
source of energy other than that provided by gravitational contraction,
it must change on a thermal time scale of
which for the sun is of the order of 10
yr. Geologists tell us
that the climate on the earth has changed little over more than
yr, and so there must be some additional source of
energy. By the late 1930's it had become clear that this must be a
nuclear energy source.
- Nuclear energy is quite capable in the case of the sun of
providing the entire energy output. From the typical
binding energy of a nucleon in a nucleus of 8 MeV, compared to the
value of
MeV, we see that combining H into heavier
nuclei releases about 1% of the mass energy as other forms of
energy. Converting 1% of the sun's mass into energy could release
about
J, enough to supply the current luminosity
for
yr.
- We have thus been led to the following view of stellar
evolution.
- Stars form by collapse of interstellar gas clouds under the
influence of gravity. As they shrink, they release gravitational
energy. Once a protostar has become sufficiently dense to be opaque,
this energy cannot be freely radiated away, and the internal
temperature starts to rise. This provides pressure support, so the
protostar shrinks on a thermal time scale.
- As the star shrinks, the internal temperature rises. Eventually
it reaches a value of the order of
K, high enough to begin
fusing H nuclei into He (we will look at the nuclear details
soon). Energy release from this reaction replaces energy that leaks
out of the star, allowing the star to enter a long period of nuclear
equilibrium known as the main sequence phase.
- Eventually the H fuel near the centre (where conditions are
right for nuclear reactions) is exhausted. The star stars to evolve on
a thermal time scale again, with the under-supported core shrinking
(and a low density envelope expanding to a large size). The
temperature at the centre continues to rise, making possible other
nuclear reactions formerly inhibited by the Coulomb barrier, such as
the conversion of He into C and O. These reactions sustain the star for
a while.
- Finally, all useful sources of nuclear energy are exhausted. The
star shrinks or collapses until it finds a structure in which its mass
can be supported by electron degeneracy (a white dwarf) or nuclear
degeneracy (a neutron star), or it collapses completely as a black
hole. If it becomes a neutron star or black hole, the collapse is
accompanied by an enormous gravitational energy release, and the star
briefly becomes a supernova.
- We now look at the fusion reactions that power the main sequence
life of a star. Since a star is (initially) mainly composed of H, we
look for reactions using this nucleus, which has the lowest Coulomb
barrier. We need a reaction that can go in spite of temperatures
corresponding to particle energies of only
keV.
- The big bottleneck to nuclear reactions involving H - apart
from the Coulomb barrier - is the fact that the product of the
collision of two protons,
He, is completely unstable on a
nuclear time scale of
s. The only way around this problem
is for one of the two protons, while they are briefly in contact, to
undergo beta decay and produce a deuteron!
- The initial reaction which makes stellar fusion possible is thus
often written in obvious notation as
H(p,e
)
H.
This reaction is mostly followed by conversion of deuterium into helium-3,
and when the abundance of
He has increased to a significant
level, a final reaction leads to the formation of
He through
This cycle (the ``PPI chain'') results in the release of 6.55 MeV per
proton consumed, or 26.2 MeV per He nucleus formed, plus about 0.5 MeV
that escapes the star directly in the form of neutrinos.
- This is not the only pathway followed by nuclear reactions in
the centre of the sun; as is usually the case inside stars, a number
of reactions take place more or less simultaneously. One side branch
of the chain above is the ``PPII chain''
and in turn a branch from this cycle is the ``PPIII chain''
- To a small extent in the sun, and to a much larger extent in
hotter stars where temperatures and particle energies are hotter,
another cycle catalyzed by C, N and O (called the ``CNO cycle'') leads
to the same basic result, conversion of H into He with release of energy
This reaction is able to compete with the PP chains - in spite of a
considerably higher Coulomb barrier - because the weak decays occur
in nuclei that do not fall apart while waiting for beta-decay to
happen. Note also that this reaction cycle requires the prior presence
of C (or N or O) nuclei (in the sun, about 1 nucleus in 2500 is C),
and does not lead to permanent conversion of C to heavier nuclei - at
the end of the cycle, the original
C nucleus is replaced.
- Stellar nuclear reactions begin at temperatures that correspond
to particle energies very low compared to the Coulomb barriers; even
under these unfavourable conditions enough energy can be released to
replace that leaked out into space. Thus we are in a situation where
most of the reactions will be due to particles in the high energy tail
of the thermal energy distribution, and even for these particles there
is a major tunneling barrier. We now consider how the reaction rates
needed for theoretical stellar models can be calculated.
- We have seen that the reaction rate per target nucleus
is
, where
is the number density of bombarding
particles
,
is the relative velocity, and
is the total reaction
cross section, for example for the first step of the PPI
chain. Because the probability of interaction depends strongly on
particle velocity, we will have to average this rate over the velocity
distribution. Thus the reaction rate
per unit volume of gas is
where
and the factor
is 1 for reactions between dissimilar particles
(for example protons and C nuclei), and
for reactions between
identical particles (this 1/2 prevent counting the same reaction
twice).
- For the cross section we will use the low-energy limiting form
derived previously,
where
in
general; for proton-proton reactions both
's are 1.
- The probability of two nuclei having a relative velocity
, or
a relative kinetic energy
, is given by the Maxwell-Boltzmann
distribution
or
where
is the reduced mass of the interacting pair.
Using the low energy form of the reaction cross section, and changing
to
as the variable, we have
where
.
- This integral cannot be done analytically (although there is no
problem in doing it numerically). But the integrand is strongly peaked
around the value
where
is a minimum; for somewhat smaller
values the integrand drops
quickly to zero because of the Coulomb barrier, while for larger
it goes to zero because of the decrease of the Boltzmann factor. We
may evaluate the integral approximately by using a parabolic Taylor
expansion of
around
,
where
The term with
is zero because we are expanding around
. In this approximation, the integrand is a Gaussian, and the
integration limits may be taken as
. Using
, we find
with
Practical expressions are given by C & G. Notice that the only
dependence on
is through the factor
.
Typically
is of the order of 20 for p-p reactions at
K,
so the exponential factor is very small; nevertheless the number
of collisions is large enough for this reaction to make possible
nuclear release of 10
W.
- For your interest, the factors
for a few reactions (for
energies well below important resonances) are listed in the table
below (these are rather old values, taken from Fowler et al 1967, Ann
Rev Astr Ap 5, 528, which however provides an excellent summary of the
use of nuclear data in stellar computations).
Reaction |
 |
Energy release , MeV |
 |
 |
0.42 |
 |
 |
5.49 |
 |
5.00 |
12.86 |
 |
 |
1.59 |
 |
 |
17.35 |
 |
 |
0.14 |
 |
 |
1.94 |
 |
 |
7.55 |
 |
 |
7.29 |
 |
 |
4.97 |
You can see that
is of the order of 1 for reactions which
involve only the strong interaction, of order
or less for
reactions in which a
must be emitted, and of order
if a weak decay must occur to allow the reaction to proceed. This is a
measure of the weakness of the weak interaction.
- The sun produces enormous numbers of neutrinos (of order
m
s
) mostly from the PPI reactions but also
from all the alternative chains (since it is always necessary to
convert two p's to n's to make each
He nucleus). The PPI neutrinos
have energies of 0.42 MeV or less, but a tiny fraction of the others have
energies above 10 MeV.
- There is much interest in observing solar neutrinos both to test
our models of the sun and to study neutrino physics. Three basic
methods are used.
- Low energy electron neutrinos may be studied using neutrino
capture by neutrons in the elements as
Ga (which changes
into Ge) in the GALLEX experiment in Italy, or
Cl (which
becomes Ar) in the Homestake Mine in the US. The beta-unstable nuclei
produced are detected using x-rays emitted following K-capture.
- The very high energy neutrinos (above about 5 MeV) are studied
by observing them scattering from electrons in very pure
water. Neutrinos entering water are virtually unable to interact with
the nuclei - even an electron neutrino cannot react with a single
proton in the H, and the threshold for reaction in
O is
almost 15 MeV, too high to capture more than a tiny number of
neutrinos. Thus elastic scattering off electrons is all that is
left. This scattering is produced by all flavours of neutrinos,
and is observed via Cerenkov radiation detected by large
photomultipliers. Some 20 kilotonnes of water is used in the Super
Kamiokande facility in Japan (which recently suffered catastrophic
failure).
- The Sudbury Neutrino Observatory (SNO) is a kind of hybrid of
the two methods above. It has a detector with one kilotonne of heavy
water (D
O) which up to now has been used to observe Cerenkov
radiation from energetic electrons produced by
(conversion of the neutron in deuterium into a
proton). Up to now the threshold for detection has been similar to
that in Japan.
- The cross section for (electron) neutrino capture by a neutron
is a lot higher than for scattering by electrons, so neutrino capture
experiments don't need as much target material as a scattering
experiment (30 tonnes at GALLEX, 600 tonnes in the Homestake Mine,
1 kilotonne in Sudbury, compared to 20 kilotonnes at Super
Kamiokande).
- The detection rate of high energy neutrinos by SNO (which
detects only
's) is lower than that observed at Super
Kamiokande (which detects all flavours). It appears that this provides
strong direct evidence that neutrinos oscillate between their flavours
as they travel from the sun to the earth.
- Fusion offers an important potential source of energy on earth,
since it could provide enormous amounts of energy from an almost
inexhaustible energy source, and produce negligible pollution.
- Because of the tiny cross section, proton-proton reactions are
useless. However, deuterium makes up 10
of all H on earth, and
reactions such as
and
which release 3.3 and 4.0 MeV respectively may be
usable. Alternatively, the deuterium-tritium reaction, which
releases more than four times as much energy as a deuterium-deuterium
reaction, has a larger reaction cross section, and uses easily produced
tritium, may become practical.
- The plasmas in which controlled fusion occurs must have a
temperature of the order of 20 keV (
K) and hence must
be confined in some non-mechanical way (e.g. by magnetic fields) and
heated very rapidly. It has turned out to be extremely difficult to do
this in such a way as to extract more fusion energy than the energy
initially expended in making fuel, confining it, and heating
it. Heating and energy extraction have been particularly
problematical. At present - after decades of experimentation - we
are still apparently far from useful power extraction.
- When a star has extracted all the available energy from the
nuclear fuels found in the deep interior (where the temperature is
high enough), it cannot continue to support itself by thermal
pressure. In this case it inexorably shrinks, or even collapses.
- Two structures are available in nature in which the star can
support itself without thermal pressure. These structures are
available only at considerably higher mean densities than those found
in main sequence stars. One, the white dwarf state, in which the
star is supported by electron degeneracy pressure, can be
reached through shrinkage. The other, the neutron star state, in
which support pressure is provided by degenerate neutrons, seems only
to be reachable following a catastrophic collapse. Such a collapse may
be triggered by reaching such high density that capture of
electrons at the top of the Fermi sea by protons to form neutrons
becomes energetically advantageous.
- In both kinds of stars, pressure support is provided by the high
energies forced on fermions at high density by the exclusion
principle. We have already met this effect in action in atomic nuclei,
where the population of successively higher energy levels is forced by
by this effect. We have seen that the total number of fermion states
per unit volume available up to energy
is
Thus to increase the number of particles in a given volume, it is
necessary to increase the maximum
of states occupied
(for non-relativistic particles) as
where
is the number of particles of
a specific type per unit volume. Since higher maximum energy (and also
higher mean energy) means that the particles exert higher pressure -
even if they are at
= 0 - a sufficiently high density may provide
the pressure to support a star against gravity.
- For a white dwarf, minimizing the total energy (non-relativistic
degenerate energy plus gravitational energy) for a given mass, by
varying the radius (C & G do the algebra), gives a radius of roughly
where we assume
. Thus a white dwarf is about the size of the
earth. The corresponding result for a neutron star replaces the 7000
km with 12.6 km; a 1
neutron star is about the size of a
small asteroid.
- Degenerate pressure support is however only conditionally
available. if the density is high enough to make most of the particles
relativistic, the pressure-density relationship changes and the
pressure does not grow fast enough with increased density to support
added weight. There is a limit - about 1.5 solar masses for white
dwarfs, of order 3
for neutron stars - above which the
degenerate structure cannot support a star. If a star collapses to
neutron star dimensions with too large a mass remaining, it becomes a
black hole.
- Returning to stellar evolution: after H fuel in the centre of a
star is exhausted, the core shrinks (and the envelope expands to
create a giant star). Core shrinkage releases gravitational energy,
increasing
to 10
K.
He begins to fuse
Be is unstable, but a tiny equilibrium population is
established, making possible
so the unstable step in the synthesis chain is jumped over. At the same
time alphas begin to react with the carbon:
- When this fuel is exhausted, in a star of solar mass or so,
electrons in the core become degenerate and the star becomes a white
dwarf. In a more massive star contraction - and reactions -
continue in ever more complex paths, for example
at about 10
K.
- As the temperature continues to rise, photons begin to have
enough energy to break up nuclei, and nuclear statistical
equilibrium is established, much like the equilibrium among various
ionization states in a plasma at 10
K. Gammas destroy nuclei which
reform by collisions, so that a whole range of intermediate nuclei
form. Much iron is produced at this stage, but the destruction of
heavy nuclei into lighter ones costs energy, so the pressure
falls below the value needed to support the star. The density
increases, the temperature rises, but at this point there is
essentially no more nuclear energy to be gotten. The core of the star
collapses to become a neutron star or black hole.
- The collapse releases an enormous amount of gravitational
energy, about 10
MeV per particle. This is about 10% of the total
rest mass of a proton or neutron, ten time more energy than has been
liberated during the entire nuclear evolution up to this point. All
the nuclear synthesis done in the core up to this point is
undone. Somehow (it's still not really clear how) a good fraction of
the energy is also transferred to the envelope of the star (which has
not had time to collapse yet). The envelope is blown off at a speed of
10
km s
or more.
- During the collapse, as the electrons are compressed to ever
higher densities and the Fermi energy goes up and up, it becomes
energetically advantageous for the electron capture reaction (inverse beta decay)
to occur, converting most of the protons into neutrons and setting
the stage for creation of a neutron star (if the core mass is not too
high), or a black hole. This reaction releases floods of neutrinos,
enough so that when a supernova occurred in the Large Magellanic
Cloud, a small galaxy at about 160,000 light-yr from us, about 20 of
these neutrinos were detected by Kamiokande and the IMB detector in
the USA.
- Another very interesting aspect of this kind of late nuclear
evolution is that this is when the elements heavier than iron are
created. Already during helium and carbon burning, such reactions as
C(
,n)
O produce neutrons of which many are
absorbed by pre-existing nuclei of iron peak elements, gradually
producing small numbers of heavier elements all the way up to
by slow neutron addition (the s-process). This neutron
addition is slow in the sense that most beta-decays have time to occur
between neutron additions, so as a nucleus evolves it stays close to
the valley of beta stability.
- During a supernova explosion, violent heating of the envelope as
it is ejected leads to many nuclear reactions, some of which release
neutrons as well. In this setting, rapid neutron addition (the r-process) occurs, in which neutrons are added more rapidly than
decays can occur, leading - temporarily - to some nuclei far from
the line of stability. These decay back to the valley, of course, but
some nuclei can only have been formed in this way.
- When the universe first formed, it contained only H, He and
tiny amounts of Li, Be and B. All the heavier elements have formed
either in evolving stars or in supernova explosions. The fact that
elements up to iron form from reactions that occur in large regions of
a star, while heavier elements form only by neutron addition, is the
main reason that the relative abundances of chemical elements drop
abruptly after the iron peak, from elements like Fe and Cr that are
present in the solar system at a level of about 10
to 10
(by number) of H, while significantly heavier elements have relative
number abundances nearer 10
compared to H.
Nuclear fusion in stars and laboratories
This document was generated using the
LaTeX2HTML translator Version 99.2beta8 (1.46)
Copyright © 1993, 1994, 1995, 1996,
Nikos Drakos,
Computer Based Learning Unit, University of Leeds.
Copyright © 1997, 1998, 1999,
Ross Moore,
Mathematics Department, Macquarie University, Sydney.
The command line arguments were:
latex2html lec10-fusion.tex
The translation was initiated by on 2002-04-01
2002-04-01