Stellar nucleosynthesis

Stellar nucleosynthesis is the cosmic alchemy that forges all the chemical elements heavier than hydrogen, helium, and lithium within the fiery hearts of stars. This incredible process not only powers the stars we see but also explains the very composition of the universe around us, from the iron in our blood to the silicon in our planet. Understanding how stars create elements is key to unraveling our own cosmic origins. Stars are the universe's natural foundries, building heavier elements through nuclear fusion from lighter ones. Different types of stars and stellar events, like supernovae, are responsible for synthesizing distinct ranges of elements. The precise conditions within a star, especially temperature, dictate which fusion reactions dominate and how quickly new elements are formed.

Source: Wikipedia

AI Summary

Stars are the universe's natural foundries, building heavier elements through nuclear fusion from lighter ones.
Different types of stars and stellar events, like supernovae, are responsible for synthesizing distinct ranges of elements.
The precise conditions within a star, especially temperature, dictate which fusion reactions dominate and how quickly new elements are formed.

The Universe's Elemental Architects

Long before stars existed, the Big Bang created the lightest elements: hydrogen, helium, and a tiny bit of lithium. But where did everything else come from—the carbon in our bodies, the oxygen we breathe, the iron in Earth's core? The answer lies within the magnificent furnaces of stars, through a process we call stellar nucleosynthesis.

This 'star stuff' creation isn't just a theory; it's a powerful framework that accurately predicts the abundances of elements we observe across the cosmos. It explains why some elements are far more common than others, and why their distribution has changed dramatically over billions of years of cosmic history.

A Story of Cosmic Recycling

Stars aren't static; they constantly evolve, changing their internal composition as they consume fuel. From hydrogen to helium, then carbon, oxygen, and beyond, a star's life is a journey of elemental transformation. But these elements don't stay locked inside forever.

As stars reach the end of their lives, they release their newly forged elements back into space. Less massive stars gently shed their outer layers as stellar winds, forming beautiful planetary nebulae. More massive stars end their lives in spectacular, cataclysmic supernovae, blasting immense quantities of heavy elements across the galaxy.

These violent events—called supernova nucleosynthesis—are responsible for creating many of the heaviest elements, scattering them into interstellar clouds. These enriched clouds then become the raw material for new generations of stars, planets, and even life itself. It's a grand cycle of cosmic creation and rebirth.

Pioneers of Stellar Alchemy

The journey to understanding stellar nucleosynthesis began with a mystery: the highly varied, 'jagged sawtooth' pattern of elemental abundances in our solar system. This non-random distribution hinted at a profound, underlying physical process. Another clue emerged as scientists realized nuclear fusion could power the Sun's incredible longevity.

In 1920, Arthur Eddington first proposed that stars got their energy from fusing hydrogen into helium, and even suggested they might create heavier elements. This was a revolutionary idea! A few years later, in 1928, George Gamow provided a crucial piece of the puzzle: the 'Gamow factor,' explaining how atomic nuclei could overcome their electrostatic repulsion to fuse—even at stellar temperatures.

The Gamow factor describes the probability of quantum tunneling, allowing nuclei to 'tunnel' through the Coulomb barrier rather than needing enough kinetic energy to leap over it. This exponential term, dependent on the Gamow energy (EG) and kinetic energy (E), is why even slight temperature increases drastically boost fusion rates.

P \propto e^{- \sqrt{\frac{E_G}{E}}}

Then, in 1939, Hans Bethe laid out the two main hydrogen fusion pathways in stars: the proton-proton chain and the carbon-nitrogen-oxygen (CNO) cycle. While Bethe focused on energy generation, Fred Hoyle later expanded this, detailing how stars could build elements up to iron. This culminated in the famous 1957 B2FH paper, which synthesized decades of research into the comprehensive theory we know today.

Hydrogen: The Star's First Fuel

The primary energy source for most stars, including our Sun, is the fusion of hydrogen into helium. This process, often called 'hydrogen burning,' is critical for a star's main-sequence life. There are two main ways stars achieve this monumental feat, depending on their mass and internal temperature.

The Proton-Proton Chain

In stars like our Sun, the 'proton-proton chain' dominates. It begins with two protons fusing to form deuterium, releasing a positron and a neutrino. This chain reaction then builds up to a helium-4 nucleus, releasing a total of about 26.2 MeV of energy per complete cycle. This process is highly sensitive to temperature, with a roughly T^4 dependence.

The CNO Cycle

In more massive stars, a different process takes over: the 'CNO cycle.' This is a catalytic process, using carbon, nitrogen, and oxygen nuclei as intermediaries to fuse hydrogen into helium. It's even more temperature-sensitive than the proton-proton chain, with an astounding T^16 to T^20 dependence, meaning a small temperature rise causes a huge jump in energy output.

The Sun's core temperature is about 15.7 million Kelvin, where the proton-proton chain is most efficient. However, in stars about 1.3 times the Sun's mass or larger, where core temperatures exceed roughly 17 million Kelvin, the CNO cycle becomes the dominant energy producer, leading to different stellar structures and evolution paths.

Beyond Hydrogen: Helium Fusion

Once a star exhausts its core hydrogen, it begins to contract and heat up, eventually reaching the extreme temperatures needed to ignite helium fusion. This often occurs when stars become red giants. In Sun-like stars, this ignition can be quite dramatic, starting with a 'helium flash' in a degenerate core.

The Triple-Alpha Process

Helium fusion primarily occurs through the 'triple-alpha process.' Here, three helium nuclei (alpha particles) combine to form one carbon nucleus. This process is tricky, as it requires two helium nuclei to first fuse into unstable beryllium-8, which must then immediately capture a third helium nucleus before it decays. If successful, it's how the universe gets its carbon!

Once carbon is formed, further helium capture can create oxygen, neon, and other heavier elements with even numbers of protons through what's known as the 'alpha process.' This explains why elements like carbon, oxygen, and neon are relatively abundant in the cosmos.

Building the Heaviest Elements

In truly massive stars, fusion doesn't stop at carbon and oxygen. As their cores continue to contract and heat, they sequentially burn heavier and heavier fuels: carbon burning, neon burning, oxygen burning, and finally silicon burning. Each stage produces progressively heavier elements, all the way up to iron and nickel.

But iron is a special element in stellar nucleosynthesis. Fusing elements heavier than iron actually consumes energy rather than releasing it. This energy crisis in the core of a massive star leads to its catastrophic collapse and a supernova explosion, which itself is a final, furious epoch of explosive nucleosynthesis.

Neutron and Proton Capture

Elements heavier than iron are primarily forged not through fusion, but by neutron or proton capture processes. The 's-process' (slow neutron capture) occurs in certain evolved stars, gradually building heavier nuclei. The 'r-process' (rapid neutron capture), however, happens in explosive events like supernovae or neutron star mergers, quickly adding many neutrons to nuclei before they can decay.

Each type of reaction, from slow hydrogen burning to explosive neutron capture, contributes to the rich tapestry of elements that make up our universe. It's a complex, beautiful cosmic ballet that has been ongoing for billions of years, creating everything we see around us.

Key Elemental Pathways: Hydrogen fusion: Proton-proton chain, CNO cycle Helium fusion: Triple-alpha process, Alpha process Heavier element fusion: Carbon, Neon, Oxygen, Silicon burning Beyond Iron: s-process, r-process (neutron capture), rp-process, p-process (proton capture), Photodisintegration

Article

Stellar nucleosynthesis

Logarithmic scales plot of the relative energy output (ε) of the following fusion processes at different temperatures (T): Proton–proton chain (PP) CNO cycle Triple-α process Combined energy generation of PP and CNO within a star The Sun's core temperature (about 1.57×107 K, with ⁠${\displaystyle log{10}T=7.20}$⁠), at which PP is more efficient.

In astrophysics, stellar nucleosynthesis is the creation of chemical elements by nuclear 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, who later refined it in 1954. 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, 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 early in the 20th century of regular variations in the abundances of elements found in the universe according to the Oddo–Harkins rule. 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). This suggested a natural process that is not random. A second stimulus to understanding the processes of stellar nucleosynthesis occurred later in the same century, when it was realized by Arthur Eddington and expanded on by George Gamow and then Hans Bethe that the energy released from nuclear fusion reactions accounted for the longevity of the Sun as a source of heat and light.

History

Stellar nucleosynthesis

In 1920, Arthur Eddington proposed that stars obtained their energy from nuclear fusion of hydrogen to form helium and also raised the possibility that the heavier elements are produced in stars.

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. 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. In the following decade the Gamow factor was used by Robert d'Escourt Atkinson and Fritz 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 a 1939 paper entitled "Energy Production in Stars", Hans Bethe analyzed the different possibilities for reactions by which hydrogen is fused into helium. 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. These works concerned the energy generation capable of keeping stars hot. 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. 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.

Hoyle's theory was extended to other processes, beginning with the publication of the 1957 review paper "Synthesis of the Elements in Stars" by Margaret Burbidge, Geoffrey Burbidge, William Alfred Fowler and Fred Hoyle, more commonly referred to as the B2FH paper. This review paper synthesized information from across fields of nuclear physics, stellar evolution and abundance of the elements to create the foundation of a theory of stellar nucleosynthesis. In 1957 Alastair G. W. Cameron, working independently showed the Oddo–Harkins even-odd abundance rule would follow from processes outlined by the B2FH paper. Clayton calculated the first time-dependent models of the s-process in 1961 and of the r-process in 1965, as well as of the burning of silicon into the abundant alpha-particle nuclei and iron-group elements in 1968, and discovered radiogenic chronologies for determining the age of the elements.

Cross section of a supergiant showing nucleosynthesis and elements formed.

Key reactions

Stellar nucleosynthesis

A version of the periodic table indicating the origins – including stellar nucleosynthesis – of the elements. All elements above 103 are omitted from this table and they are produced only by human synthesis.

The most important reactions in stellar nucleosynthesis:

• Hydrogen fusion: • Deuterium fusion • The proton–proton chain • The carbon–nitrogen–oxygen cycle • Helium fusion: • The triple-alpha process • The alpha process • Fusion of heavier elements: • Lithium burning: a process found most commonly in brown dwarfs • Carbon-burning process • Neon-burning process • Oxygen-burning process • Silicon-burning process • Production of elements heavier than iron: • Neutron capture: • The r-process • The s-process • Proton capture: • The rp-process • The p-process • Photodisintegration

Hydrogen fusion

Hydrogen fusion (nuclear fusion of four protons to form a helium-4 nucleus) 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: the 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.

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. In each complete fusion cycle, the proton–proton chain reaction releases about 26.2 MeV. Proton-proton chain with a dependence of approximately T4, meaning the reaction cycle is highly sensitive 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, 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. 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. 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. CNO cycle is highly sensitive to temperature, with rates proportional to the 16th to 20th power of the temperature; a 10% increase in temperature would result in a 350% increase 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. 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. 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. The Sun produces on the order of 1% of its energy from the CNO cycle.

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, 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. 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. The Sun itself has a core temperature of about 1.57×107 K. As a main-sequence star ages, the core temperature will rise, resulting in a steadily increasing contribution from its CNO cycle.

Helium fusion

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. 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.

In all cases, helium is fused to carbon via the triple-alpha process, i.e., three helium nuclei are transformed into carbon via 8Be. 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.

Reaction rate

Stellar nucleosynthesis

The reaction rate density between species A and B, having number densities nA,B, is given by:${\displaystyle r=n{A}\,n{B}\,k{r}}$where kr is the reaction rate constant of each single elementary binary reaction composing the nuclear fusion process;${\displaystyle k{r}=\langle \sigma (v)\,v\rangle }$where σ(v) is the cross-section at relative velocity v, and averaging is performed over all velocities.

Semi-classically, the cross section is proportional to ${\textstyle \pi \,\lambda ^{2}}$, where ${\textstyle \lambda =h/p}$ is the de Broglie wavelength. Thus semi-classically the cross section is proportional to ${\textstyle {\frac {E}{m}}=c^{2}}$.

However, since the reaction involves quantum tunneling, there is an exponential damping at low energies that depends on Gamow factor EG, given by an Arrhenius-type equation:${\displaystyle \sigma (E)={\frac {S(E)}{E}}e^{-{\sqrt {\frac {E{\text{G}}}{E}}}}.}$Here astrophysical S-factor S(E) depends on the details of the nuclear interaction, and has the dimension of an energy multiplied by a cross section.

One then integrates over all energies to get the total reaction rate, using the Maxwell–Boltzmann distribution and the relation:${\displaystyle {\frac {r}{V}}=n{A}n{B}\int {0}^{\infty }{\Bigl (}{\frac {S(E)}{E}}\,e^{-{\sqrt {\frac {E{\text{G}}}{E}}}}\cdot 2{\sqrt {\frac {E}{\pi (kT)^{3}}}}\,e^{-{\frac {E}{kT}}}\,\cdot {\sqrt {\frac {2E}{m{\text{R}}}}}{\Bigr )}dE}$where k = 86,17 μeV/K, ${\displaystyle m{\text{R}}={\frac {m{A}m{B}}{m{A}+m{B}}}}$ is the reduced mass. The integrand equals

${\displaystyle S(E)\,e^{-{\sqrt {\frac {E{\text{G}}}{E}}}}\cdot 2{\sqrt {2/\pi }}(kT)^{-3/2}\,e^{-{\frac {E}{kT}}}\,/{\sqrt {m{\text{R}}}}.}$

Since this integration of f(E, constant T) has an exponential damping at high energies of the form ${\textstyle \sim e^{-{\frac {E}{kT}}}}$ and at low energies from the Gamow factor, the integral almost vanishes everywhere except around the peak at E0, called Gamow peak. There:${\displaystyle -{\frac {\partial }{\partial E}}\left({\sqrt {\frac {E{\text{G}}}{E}}}+{\frac {E}{kT}}\right)\,=\,0}$

${\displaystyle E{0}=\left({\frac {1}{2}}kT{\sqrt {E{\text{G}}}}\right)^{\frac {2}{3}}}$ and ${\displaystyle {\sqrt {E{\text{G}}}}=E{0}^{\frac {3}{2}}/{\frac {1}{2}}kT}$

The exponent can then be approximated around E0 as:${\displaystyle e^{-({\frac {E}{kT}}+{\sqrt {\frac {E{\text{G}}}{E}}})}\approx e^{-{\frac {3E{0}}{kT}}}e^{{\bigl (}-{\frac {3(E-E{0})^{2}}{4E{0}kT}}{\bigr )}}=e^{-{\frac {3E{0}}{kT}}{\bigl (}1+({\frac {E-E{0}}{2E{0}}})^{2}{\bigr )}}=e^{-{\frac {3E{0}}{kT}}{\bigl (}1+(E/E{0}-1)^{2}/4{\bigr )}}}$

And the reaction rate is approximated as:${\displaystyle {\frac {r}{V}}\approx n{A}\,n{B}\,{\frac {4{\sqrt {(}}2/3)}{\sqrt {m{\text{R}}}}}\,{\sqrt {E{0}}}{\frac {S(E{0})}{kT}}\,e^{-{\frac {3E{0}}{kT}}}}$

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 (the Sun) give kT of the order of 1 keV: ${\textstyle \log {10}kT=-16+\log {10}2.17}$.

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. 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.