In astrophysics, stellar nucleosynthesis is the creation of chemical elements by nuclear fusion reactions within stars. Stellar nucleosynthesis has occurred since the original creation of hydrogen, helium and lithium during the Big Bang. As a predictive theory, it yields accurate estimates of the observed abundances of the elements. It explains why the observed abundances of elements change over time and why some elements and their isotopes are much more abundant than others. The theory was initially proposed by Fred Hoyle in 1946, [1] who later refined it in 1954. [2] Further advances were made, especially to nucleosynthesis by neutron capture of the elements heavier than iron, by Margaret and Geoffrey Burbidge, William Alfred Fowler and Fred Hoyle in their famous 1957 B2FH paper, [3] which became one of the most heavily cited papers in astrophysics history.
Stars evolve because of changes in their composition (the abundance of their constituent elements) over their lifespans, first by burning hydrogen (main sequence star), then helium (horizontal branch star), and progressively burning higher elements. However, this does not by itself significantly alter the abundances of elements in the universe as the elements are contained within the star. Later in its life, a low-mass star will slowly eject its atmosphere via stellar wind, forming a planetary nebula, while a higher–mass star will eject mass via a sudden catastrophic event called a supernova. The term supernova nucleosynthesis is used to describe the creation of elements during the explosion of a massive star or white dwarf.
The advanced sequence of burning fuels is driven by gravitational collapse and its associated heating, resulting in the subsequent burning of carbon, oxygen and silicon. However, most of the nucleosynthesis in the mass range A = 28–56 (from silicon to nickel) is actually caused by the upper layers of the star collapsing onto the core, creating a compressional shock wave rebounding outward. The shock front briefly raises temperatures by roughly 50%, thereby causing furious burning for about a second. This final burning in massive stars, called explosive nucleosynthesis or supernova nucleosynthesis, is the final epoch of stellar nucleosynthesis.
A stimulus to the development of the theory of nucleosynthesis was the discovery of variations in the abundances of elements found in the universe. The need for a physical description was already inspired by the relative abundances of the chemical elements in the solar system. Those abundances, when plotted on a graph as a function of the atomic number of the element, have a jagged sawtooth shape that varies by factors of tens of millions (see history of nucleosynthesis theory). [4] This suggested a natural process that is not random. A second stimulus to understanding the processes of stellar nucleosynthesis occurred during the 20th century, when it was realized that the energy released from nuclear fusion reactions accounted for the longevity of the Sun as a source of heat and light. [5]
In 1920, Arthur Eddington, on the basis of the precise measurements of atomic masses by F.W. Aston and a preliminary suggestion by Jean Perrin, proposed that stars obtained their energy from nuclear fusion of hydrogen to form helium and raised the possibility that the heavier elements are produced in stars. [6] [7] [8] This was a preliminary step toward the idea of stellar nucleosynthesis. In 1928 George Gamow derived what is now called the Gamow factor, a quantum-mechanical formula yielding the probability for two contiguous nuclei to overcome the electrostatic Coulomb barrier between them and approach each other closely enough to undergo nuclear reaction due to the strong nuclear force which is effective only at very short distances. [9] : 410 In the following decade the Gamow factor was used by Atkinson and Houtermans and later by Edward Teller and Gamow himself to derive the rate at which nuclear reactions would occur at the high temperatures believed to exist in stellar interiors.
In 1939, in a Nobel lecture entitled "Energy Production in Stars", Hans Bethe analyzed the different possibilities for reactions by which hydrogen is fused into helium. [10] He defined two processes that he believed to be the sources of energy in stars. The first one, the proton–proton chain reaction, is the dominant energy source in stars with masses up to about the mass of the Sun. The second process, the carbon–nitrogen–oxygen cycle, which was also considered by Carl Friedrich von Weizsäcker in 1938, is more important in more massive main-sequence stars. [11] : 167 These works concerned the energy generation capable of keeping stars hot. A clear physical description of the proton–proton chain and of the CNO cycle appears in a 1968 textbook. [12] : 365 Bethe's two papers did not address the creation of heavier nuclei, however. That theory was begun by Fred Hoyle in 1946 with his argument that a collection of very hot nuclei would assemble thermodynamically into iron. [1] Hoyle followed that in 1954 with a paper describing how advanced fusion stages within massive stars would synthesize the elements from carbon to iron in mass. [2] [13]
Hoyle's theory was extended to other processes, beginning with the publication of the 1957 review paper "Synthesis of the Elements in Stars" by Burbidge, Burbidge, Fowler and Hoyle, more commonly referred to as the B2FH paper. [3] This review paper collected and refined earlier research into a heavily cited picture that gave promise of accounting for the observed relative abundances of the elements; but it did not itself enlarge Hoyle's 1954 picture for the origin of primary nuclei as much as many assumed, except in the understanding of nucleosynthesis of those elements heavier than iron by neutron capture. Significant improvements were made by Alastair G. W. Cameron and by Donald D. Clayton. In 1957 Cameron presented his own independent approach to nucleosynthesis, [14] informed by Hoyle's example, and introduced computers into time-dependent calculations of evolution of nuclear systems. Clayton calculated the first time-dependent models of the s-process in 1961 [15] and of the r-process in 1965, [16] as well as of the burning of silicon into the abundant alpha-particle nuclei and iron-group elements in 1968, [17] [18] and discovered radiogenic chronologies [19] for determining the age of the elements.
The most important reactions in stellar nucleosynthesis:
Hydrogen fusion (nuclear fusion of four protons to form a helium-4 nucleus [20] ) is the dominant process that generates energy in the cores of main-sequence stars. It is also called "hydrogen burning", which should not be confused with the chemical combustion of hydrogen in an oxidizing atmosphere. There are two predominant processes by which stellar hydrogen fusion occurs: proton–proton chain and the carbon–nitrogen–oxygen (CNO) cycle. Ninety percent of all stars, with the exception of white dwarfs, are fusing hydrogen by these two processes. [21] : 245
In the cores of lower-mass main-sequence stars such as the Sun, the dominant energy production process is the proton–proton chain reaction. This creates a helium-4 nucleus through a sequence of reactions that begin with the fusion of two protons to form a deuterium nucleus (one proton plus one neutron) along with an ejected positron and neutrino. [22] In each complete fusion cycle, the proton–proton chain reaction releases about 26.2 MeV. [22] Proton-proton chain with a dependence of approximately T^4, reaction cycle is relatively insensitive to temperature; a 10% rise of temperature would increase energy production by this method by 46%, hence, this hydrogen fusion process can occur in up to a third of the star's radius and occupy half the star's mass. For stars above 35% of the Sun's mass, [23] the energy flux toward the surface is sufficiently low and energy transfer from the core region remains by radiative heat transfer, rather than by convective heat transfer. [24] As a result, there is little mixing of fresh hydrogen into the core or fusion products outward.
In higher-mass stars, the dominant energy production process is the CNO cycle, which is a catalytic cycle that uses nuclei of carbon, nitrogen and oxygen as intermediaries and in the end produces a helium nucleus as with the proton–proton chain. [22] During a complete CNO cycle, 25.0 MeV of energy is released. The difference in energy production of this cycle, compared to the proton–proton chain reaction, is accounted for by the energy lost through neutrino emission. [22] CNO cycle is highly sensitive to temperature, with rates proportional to T^{16-20}, a 10% rise of temperature would produce a 350% rise in energy production. About 90% of the CNO cycle energy generation occurs within the inner 15% of the star's mass, hence it is strongly concentrated at the core. [25] This results in such an intense outward energy flux that convective energy transfer becomes more important than does radiative transfer. As a result, the core region becomes a convection zone, which stirs the hydrogen fusion region and keeps it well mixed with the surrounding proton-rich region. [26] This core convection occurs in stars where the CNO cycle contributes more than 20% of the total energy. As the star ages and the core temperature increases, the region occupied by the convection zone slowly shrinks from 20% of the mass down to the inner 8% of the mass. [25] The Sun produces on the order of 1% of its energy from the CNO cycle. [27] [a] [28] : 357 [29] [b]
The type of hydrogen fusion process that dominates in a star is determined by the temperature dependency differences between the two reactions. The proton–proton chain reaction starts at temperatures about 4×106 K , [30] making it the dominant fusion mechanism in smaller stars. A self-maintaining CNO chain requires a higher temperature of approximately 1.6×107 K, but thereafter it increases more rapidly in efficiency as the temperature rises, than does the proton–proton reaction. [31] Above approximately 1.7×107 K, the CNO cycle becomes the dominant source of energy. This temperature is achieved in the cores of main-sequence stars with at least 1.3 times the mass of the Sun. [32] The Sun itself has a core temperature of about 1.57×107 K. [33] : 5 As a main-sequence star ages, the core temperature will rise, resulting in a steadily increasing contribution from its CNO cycle. [25]
Main sequence stars accumulate helium in their cores as a result of hydrogen fusion, but the core does not become hot enough to initiate helium fusion. Helium fusion first begins when a star leaves the red giant branch after accumulating sufficient helium in its core to ignite it. In stars around the mass of the Sun, this begins at the tip of the red giant branch with a helium flash from a degenerate helium core, and the star moves to the horizontal branch where it burns helium in its core. More massive stars ignite helium in their core without a flash and execute a blue loop before reaching the asymptotic giant branch. Such a star initially moves away from the AGB toward bluer colours, then loops back again to what is called the Hayashi track. An important consequence of blue loops is that they give rise to classical Cepheid variables, of central importance in determining distances in the Milky Way and to nearby galaxies. [34] : 250 Despite the name, stars on a blue loop from the red giant branch are typically not blue in colour but are rather yellow giants, possibly Cepheid variables. They fuse helium until the core is largely carbon and oxygen. The most massive stars become supergiants when they leave the main sequence and quickly start helium fusion as they become red supergiants. After the helium is exhausted in the core of a star, helium fusion will continue in a shell around the carbon–oxygen core. [20] [24]
In all cases, helium is fused to carbon via the triple-alpha process, i.e., three helium nuclei are transformed into carbon via 8Be. [35] : 30 This can then form oxygen, neon, and heavier elements via the alpha process. In this way, the alpha process preferentially produces elements with even numbers of protons by the capture of helium nuclei. Elements with odd numbers of protons are formed by other fusion pathways. [36] : 398
The reaction rate density between species A and B, having number densities nA,B, is given by: where k is the reaction rate constant of each single elementary binary reaction composing the nuclear fusion process: here, σ(v) is the cross-section at relative velocity v, and averaging is performed over all velocities.
Semi-classically, the cross section is proportional to , where is the de Broglie wavelength. Thus semi-classically the cross section is proportional to .
However, since the reaction involves quantum tunneling, there is an exponential damping at low energies that depends on Gamow factor EG, giving an Arrhenius equation: where S(E) depends on the details of the nuclear interaction, and has the dimension of an energy multiplied for a cross section.
One then integrates over all energies to get the total reaction rate, using the Maxwell–Boltzmann distribution and the relation: where is the reduced mass.
Since this integration has an exponential damping at high energies of the form and at low energies from the Gamow factor, the integral almost vanished everywhere except around the peak, called Gamow peak, [37] : 185 at E0, where:
Thus:
The exponent can then be approximated around E0 as:
And the reaction rate is approximated as: [38]
Values of S(E0) are typically 10−3 – 103 keV·b , but are damped by a huge factor when involving a beta decay, due to the relation between the intermediate bound state (e.g. diproton) half-life and the beta decay half-life, as in the proton–proton chain reaction. Note that typical core temperatures in main-sequence stars give kT of the order of keV. [39] : ch. 3
Thus, the limiting reaction in the CNO cycle, proton capture by 14
7N
, has S(E0) ~ S(0) = 3.5 keV·b, while the limiting reaction in the proton–proton chain reaction, the creation of deuterium from two protons, has a much lower S(E0) ~ S(0) = 4×10−22 keV·b. [40] [41] Incidentally, since the former reaction has a much higher Gamow factor, and due to the relative abundance of elements in typical stars, the two reaction rates are equal at a temperature value that is within the core temperature ranges of main-sequence stars. [42]
The CNO cycle is one of the two known sets of fusion reactions by which stars convert hydrogen to helium, the other being the proton–proton chain reaction, which is more efficient at the Sun's core temperature. The CNO cycle is hypothesized to be dominant in stars that are more than 1.3 times as massive as the Sun.
In astronomy, the main sequence is a classification of stars which appear on plots of stellar color versus brightness as a continuous and distinctive band. Stars on this band are known as main-sequence stars or dwarf stars, and positions of stars on and off the band are believed to indicate their physical properties, as well as their progress through several types of star life-cycles. These are the most numerous true stars in the universe and include the Sun. Color-magnitude plots are known as Hertzsprung–Russell diagrams after Ejnar Hertzsprung and Henry Norris Russell.
Nuclear fusion is a reaction in which two or more atomic nuclei, combine to form one or more atomic nuclei and neutrons. The difference in mass between the reactants and products is manifested as either the release or absorption of energy. This difference in mass arises as a result of the difference in nuclear binding energy between the atomic nuclei before and after the fusion reaction. Nuclear fusion is the process that powers active or main-sequence stars and other high-magnitude stars, where large amounts of energy are released.
The proton–proton chain, also commonly referred to as the p–p chain, is one of two known sets of nuclear fusion reactions by which stars convert hydrogen to helium. It dominates in stars with masses less than or equal to that of the Sun, whereas the CNO cycle, the other known reaction, is suggested by theoretical models to dominate in stars with masses greater than about 1.3 solar masses.
Stellar evolution is the process by which a star changes over the course of its lifetime and how it can lead to the creation of a new star. Depending on the mass of the star, its lifetime can range from a few million years for the most massive to trillions of years for the least massive, which is considerably longer than the current age of the universe. The table shows the lifetimes of stars as a function of their masses. All stars are formed from collapsing clouds of gas and dust, often called nebulae or molecular clouds. Over the course of millions of years, these protostars settle down into a state of equilibrium, becoming what is known as a main-sequence star.
In physical cosmology, Big Bang nucleosynthesis is the production of nuclei other than those of the lightest isotope of hydrogen during the early phases of the universe. This type of nucleosynthesis is thought by most cosmologists to have occurred from 10 seconds to 20 minutes after the Big Bang. It is thought to be responsible for the formation of most of the universe's helium, along with small fractions of the hydrogen isotope deuterium, the helium isotope helium-3 (3He), and a very small fraction of the lithium isotope lithium-7 (7Li). In addition to these stable nuclei, two unstable or radioactive isotopes were produced: the heavy hydrogen isotope tritium and the beryllium isotope beryllium-7 (7Be). These unstable isotopes later decayed into 3He and 7Li, respectively, as above.
Nucleosynthesis is the process that creates new atomic nuclei from pre-existing nucleons and nuclei. According to current theories, the first nuclei were formed a few minutes after the Big Bang, through nuclear reactions in a process called Big Bang nucleosynthesis. After about 20 minutes, the universe had expanded and cooled to a point at which these high-energy collisions among nucleons ended, so only the fastest and simplest reactions occurred, leaving our universe containing hydrogen and helium. The rest is traces of other elements such as lithium and the hydrogen isotope deuterium. Nucleosynthesis in stars and their explosions later produced the variety of elements and isotopes that we have today, in a process called cosmic chemical evolution. The amounts of total mass in elements heavier than hydrogen and helium remains small, so that the universe still has approximately the same composition.
The triple-alpha process is a set of nuclear fusion reactions by which three helium-4 nuclei are transformed into carbon.
The alpha process, also known as alpha capture or the alpha ladder, is one of two classes of nuclear fusion reactions by which stars convert helium into heavier elements. The other class is a cycle of reactions called the triple-alpha process, which consumes only helium, and produces carbon. The alpha process most commonly occurs in massive stars and during supernovae.
The carbon-burning process or carbon fusion is a set of nuclear fusion reactions that take place in the cores of massive stars (at least 8 at birth) that combines carbon into other elements. It requires high temperatures (> 5×108 K or 50 keV) and densities (> 3×109 kg/m3).
In astrophysics, silicon burning is a very brief sequence of nuclear fusion reactions that occur in massive stars with a minimum of about 8–11 solar masses. Silicon burning is the final stage of fusion for massive stars that have run out of the fuels that power them for their long lives in the main sequence on the Hertzsprung–Russell diagram. It follows the previous stages of hydrogen, helium, carbon, neon and oxygen burning processes.
Stellar structure models describe the internal structure of a star in detail and make predictions about the luminosity, the color and the future evolution of the star. Different classes and ages of stars have different internal structures, reflecting their elemental makeup and energy transport mechanisms.
The Oddo–Harkins rule holds that an element with an even atomic number is more abundant than the elements with immediately adjacent atomic numbers. For example, carbon, with atomic number 6, is more abundant than boron (5) and nitrogen (7). Generally, the relative abundance of an even atomic numbered element is roughly two orders of magnitude greater than the relative abundances of the immediately adjacent odd atomic numbered elements to either side. This pattern was first reported by Giuseppe Oddo in 1914 and William Draper Harkins in 1917. The Oddo–Harkins rule is true for all elements beginning with carbon produced by stellar nucleosynthesis but not true for the lightest elements below carbon produced by big bang nucleosynthesis and cosmic ray spallation.
Supernova nucleosynthesis is the nucleosynthesis of chemical elements in supernova explosions.
Nuclear astrophysics studies the origin of the chemical elements and isotopes, and the role of nuclear energy generation, in cosmic sources such as stars, supernovae, novae, and violent binary-star interactions. It is an interdisciplinary part of both nuclear physics and astrophysics, involving close collaboration among researchers in various subfields of each of these fields. This includes, notably, nuclear reactions and their rates as they occur in cosmic environments, and modeling of astrophysical objects where these nuclear reactions may occur, but also considerations of cosmic evolution of isotopic and elemental composition (often called chemical evolution). Constraints from observations involve multiple messengers, all across the electromagnetic spectrum (nuclear gamma-rays, X-rays, optical, and radio/sub-mm astronomy), as well as isotopic measurements of solar-system materials such as meteorites and their stardust inclusions, cosmic rays, material deposits on Earth and Moon). Nuclear physics experiments address stability (i.e., lifetimes and masses) for atomic nuclei well beyond the regime of stable nuclides into the realm of radioactive/unstable nuclei, almost to the limits of bound nuclei (the drip lines), and under high density (up to neutron star matter) and high temperature (plasma temperatures up to 109 K). Theories and simulations are essential parts herein, as cosmic nuclear reaction environments cannot be realized, but at best partially approximated by experiments.
The standard solar model (SSM) is a mathematical model of the Sun as a spherical ball of gas. This stellar model, technically the spherically symmetric quasi-static model of a star, has stellar structure described by several differential equations derived from basic physical principles. The model is constrained by boundary conditions, namely the luminosity, radius, age and composition of the Sun, which are well determined. The age of the Sun cannot be measured directly; one way to estimate it is from the age of the oldest meteorites, and models of the evolution of the Solar System. The composition in the photosphere of the modern-day Sun, by mass, is 74.9% hydrogen and 23.8% helium. All heavier elements, called metals in astronomy, account for less than 2 percent of the mass. The SSM is used to test the validity of stellar evolution theory. In fact, the only way to determine the two free parameters of the stellar evolution model, the helium abundance and the mixing length parameter, are to adjust the SSM to "fit" the observed Sun.
The B2FH paper was a landmark scientific paper on the origin of the chemical elements. The paper's title is Synthesis of the Elements in Stars, but it became known as B2FH from the initials of its authors: Margaret Burbidge, Geoffrey Burbidge, William A. Fowler, and Fred Hoyle. It was written from 1955 to 1956 at the University of Cambridge and Caltech, then published in Reviews of Modern Physics in 1957.
Beryllium-8 is a radionuclide with 4 neutrons and 4 protons. It is an unbound resonance and nominally an isotope of beryllium. It decays into two alpha particles with a half-life on the order of 8.19×10−17 seconds. This has important ramifications in stellar nucleosynthesis as it creates a bottleneck in the creation of heavier chemical elements. The properties of 8Be have also led to speculation on the fine tuning of the Universe, and theoretical investigations on cosmological evolution had 8Be been stable.
A stellar core is the extremely hot, dense region at the center of a star. For an ordinary main sequence star, the core region is the volume where the temperature and pressure conditions allow for energy production through thermonuclear fusion of hydrogen into helium. This energy in turn counterbalances the mass of the star pressing inward; a process that self-maintains the conditions in thermal and hydrostatic equilibrium. The minimum temperature required for stellar hydrogen fusion exceeds 107 K (10 MK), while the density at the core of the Sun is over 100 g/cm3. The core is surrounded by the stellar envelope, which transports energy from the core to the stellar atmosphere where it is radiated away into space.
Georgeanne (Jan) Caughlan was an American astrophysicist known for her work on stellar energy generation. Her compilation of experimental data of the rates of nuclear reactions was instrumental in establishing the theory of nucleosynthesis that led to a Nobel Prize for William A. Fowler.