The carbon challenge

1. Helium burning

Figure 1: Synthesis of carbon in a burning star. Two alpha particles react to form beryllium-8. Beryllium-8 and an alpha particle react to form carbon-12, proceeding via the Hoyle state. This state in carbon-12 is a “resonance” in the beryllium-8 plus alpha-particle system. Once the resonant state is formed, it tends to decay by breaking up into beryllium-8 and an alpha particle. However, approximately four out of ten thousand decays bring the excited carbon-12 nucleus to its stable ground state.

“…In the first-generation stars the ash resulting from hydrogen burning via the p-p chain is entirely helium-4, the creation of heavier elements having been blocked by the instabilities at A=5 and A=8. These are referred to as the mass gaps.
Since carbon-12, the fourth most abundant nuclear species observed in the universe, could not synthesized in its observed abundances in the early universe, the site for its creation has to be in stars. Thus, a major question in the early studies of nucleosynthesis was how the stability gaps were bridged to create carbon-12 using only helium. While the simultaneous interaction of three α-particles to form carbon-12 is energetically possible, the probability for this direct process is much to small to account for the observed carbon-12 abundances.
The solution of this problem was provided in principle by Salpeter and Opik, who proposed that carbon-12 was created via a two-step process. In the first step, two a-particles combine to form beryllium-8 in its ground state.
The ground state of beryllium-8 is known to be unstable against decay into two α-particles with a lifetime of 10-16 s, which is the reason for the mass-8 stability gap. However, as Salpeter pointed out, this lifetime is long compared with the 10-19 s transit time of two α-particles with kinetic energies corresponding to Q.
As a result, a small concentration of beryllium-8 nuclei builds up in equilibrium with the decay products, two α-particles.
The actual equilibrium concentration of beryllium-8 in a helium environment can be calculated using the Saha equation.
For typical values, the Saha equation leads to one beryllium-8 nucleus for every 109 helium-4 nuclei.
In the second step, Salpeter suggested that, since the first step provides an appreciable concentration of beryllium-8 nuclei, these nuclei capture an additional α-particle, thus completing the carbon-12 creation process.
Since the combined effect of these two steps is to transform three α-particles into a carbon-12 nucleus, this set of reactions is referred to as the triple-α process:

3α —> 12C

Assuming the triple-α process to be the correct mechanism for synthesis of 12C, Hoyle showed that the amount of 12C produced in this way is insufficient to expalin the obseved abundance.
To surmount this difficulty, Hoyle suggested that sufficient 12C nuclei could be synthesized if the 8Be(α,γ)12C reaction took place through an s-wave resonance near the 8Be + α threshold, since the existence of such a resonance would greatly accelerate the rate of triple-α process…..
The result of this concentrated effort was a verification of Hoyle’s prediction with regard to both the energetic location and the Jπ =0+ assignment of the state.
This is an impressive example of a prediction of the properties of a nucleus based on purely astrophysical grounds and subsequent verification of the prediction by experiments….”
“Cauldrons in the Cosmos”, Claus E. Rolfs, William S. Rodney

2.The paper
“Ab initio calculation of the Hoyle state
The Hoyle state plays a crucial role in the hydrogen burning of stars heavier than our sun and in the production of carbon and other elements necessary for life. This excited state of the carbon-12 nucleus was postulated by Hoyle as a necessary ingredient for the fusion of three alpha particles to produce carbon at stellar temperatures. Although the Hoyle state was seen experimentally more than a half century ago nuclear theorists have not yet uncovered the nature of this state from first principles. In this letter we report the first ab initio calculation of the low-lying states of carbon-12 using supercomputer lattice simulations and a theoretical framework known as effective field theory. In addition to the ground state and excited spin-2 state, we find a resonance at −85(3) MeV with all of the properties of the Hoyle state and in agreement with the experimentally observed energy. These lattice simulations provide insight into the structure of this unique state and new clues as to the amount of fine-tuning needed in nature for the production of carbon in stars….. http://arxiv.org/PS_cache/arxiv/pdf/1101/1101.2547v2.pdf

3. Analysis
“The carbon challenge
The series of reactions that leads to the production of carbon in stars was predicted more than sixty years ago, but theorists have struggled to reproduce this process from first principles calculations.
Morten Hjorth-Jensen
Carbon, so crucial for the chemistry of all living organisms and life in general, forms at the interior of burning stars in what is called the triple-alpha process (Fig. 1): two alpha particles (helium nuclei) react to form beryllium-8, which in turn reacts with a third alpha particle to form carbon-12. This last step does not produce the carbon-12 as we know it, however, but rather an excited state called the Hoyle state.
This state is a resonance, meaning it is not localized spatially and it has a finite lifetime governed by its distance (in energy) from the particle-emission threshold. The resonance allows approximately four out of ten thousand decays to produce the ground state of carbon-12. Fred Hoyle predicted the resonance in 1954, arguing that without it, stable carbon would not exist [1]. Three years later, an experimental group at Caltech proved he was right [2].
Once carbon is formed, the production of many heavier elements—like the oxygen in our seas and the atmosphere or the silicon needed in cellphones and computers—mostly falls into place. Given its fundamental role, it is natural to ask, Can we understand carbon’s nuclear structure from first principles? And what do we mean by first principles? These are intellectually challenging tasks, which, until recently, have eluded a satisfactory theoretical understanding. Perhaps this is no longer the case. An article appearing in Physical Review Letters by Evgeny Epelbaum of Ruhr University Bochum, Germany, and colleagues [3] aims at a first principles calculation that determines the energies of the ground and Hoyle states of carbon-12, as well as other excited states.
Epelbaum et al.’s work deals with a key intellectual issue in ongoing nuclear physics research: understanding the limits of stability of matter starting from basic building blocks (protons and neutrons) and the theory that describes their interaction, ideally starting from quantum chromodynamics. The complex nature of the nuclear forces between protons and neutrons produces diverse nuclear phenomena, characterized by energies of a few kiloelectronvolts to several billions of electronvolts. Capturing this range of energy scales with computations that are tractable poses a tough challenge to nuclear many-body theory and so-called first principles descriptions of nuclear systems. By first principles, we normally mean methods that allow us to solve exactly, or within controlled approximations, either the nonrelativistic Schrödinger equation or the relativistic Dirac equation for many interacting particles. The input to these methods is typically a given Hamiltonian and the relevant degrees of freedom, such as the number of neutrons and protons.
Within this context, Epelbaum et al. are tackling an extremely complicated and tough theoretical problem. Already, the ground state of carbon-12 is a difficult many-body problem. The Hoyle state lies 7.6 million electronvolts higher in energy than the ground state, but it has the same total angular momentum (J=0) and parity (positive). The fact that the two states share the same quantum numbers, but very different spatial configurations, complicates the calculations. The average density of nuclear matter is much higher than that found in solids or atoms, and the interaction between nucleons is dictated by the strong force. This force, in particular for the important s-waves, is repulsive at short distances [less than half of a femtometer (fm)] but has an attractive component at intermediate and long ranges (up to 2–3 fm). For these reasons, models that assume that nuclei are independent particles can’t fully explain what is observed experimentally, for example, the energy levels of the nuclei. Conversely, the consequences of correlations in nuclear many-particle systems are very difficult to measure experimentally and to interpret theoretically. From a theoretical standpoint, a proper understanding of such many-body forces and correlations is crucial in order to provide a predictive many-body theory. Furthermore, it is central to understand which components of the nuclear forces and correlations are responsible for a particular spectroscopic observable. Lastly, a proper many-body theory should quantify the uncertainty in its predictions.
Years of experience in condensed-matter and high-energy physics have shown the power of using lattice Monte Carlo simulations (where physical properties of a system are discretized on a lattice mesh) to predict the properties of many-body systems, see, for example, a recent review in [4]. Building on this knowledge, Epelbaum et al. combine lattice simulations with a Lagrangian derived from effective field theory [5] to predict, numerically, the energy levels and spatial structure of helium-4, beryllium-8 and carbon-12. (The Lagrangian is basically a perturbative expansion in powers of the exchanged momentum among nucleons and pions). The group demonstrates how different components of the nuclear forces determine the description of the ground states and excited states of these nuclei. They also provide statistical error estimates and show how many-body forces evolve as more nucleons are added to the nucleus. The parameters in effective field theory for the three-nucleon interaction are fixed to reproduce the binding energy of helium-4, while the results for beryllium-8, and carbon-12 are pure predictions. It is important to point out that these calculations would not have been possible without access to modern high-performance computing facilities.
The group’s results agree well with experiment. They also find a resonance for carbon-12, which lies close in energy to the experimentally known Hoyle state and its assumed spatial structure, with three alpha particles arranged almost linearly. Their analysis of this state in terms of different contributions to the nuclear interactions is of particular interest. The inclusion of three-nucleon interactions is essential to match experimental results. Ideally, future work will explore the role played by higher-order corrections, such as the inclusion of four-body interactions. However, it is already rewarding for nuclear structure theory that these calculations are able to reproduce, without an a priori assumption about the structure of a given state, the relevant physics of a complicated many-body system.
Other first principles computations have successfully described the structure of several light nuclei [6, 7], reactions between light ions [8], and heavier isotopes [9]. Such work holds great promise for a truly quantitative description of nuclei and their properties. Recent and future experimental programs in low-energy nuclear physics aim at exploring the limits of stability of nuclear matter, with wide-ranging consequences for our understanding of nuclear structure and the synthesis of the elements. Correlations beyond an independent particle picture that arise in nuclei towards either the proton or neutron driplines—the predicted limits of the chart of nuclei—should then provide us with a better understanding of shell structures and single-particle properties of nuclei. All of these efforts require models that can extract a simple physical picture from complicated many-body systems. Epelbaum et al. have opened a very interesting path toward the achievements of these ambitious goals.
Acknowledgments
I am much indebted to Thomas Papenbrock for a critical read of the manuscript.
References
F. Hoyle, Astrophys. J. Suppl. Ser. 1, 121 (1954).
C. W. Cook et al., Phys. Rev. 107, 508 (1957).
E. Epelbaum, H. Krebs, D. Lee, and U-G. Meißner, Phys. Rev. Lett. 106, 192501 (2011).
D. Lee, Prog. Part. Nucl. Phys. 63, 117 (2009).
E. Epelbaum, H. W. Hammer, and Ulf-G. Meissner, Rev. Mod. Phys. 81, 1773 (2009).
S. C. Pieper and R. B. Wiringa, Ann. Rev. Nucl. Part. Sci. 51, 53 (2001).
P. Maris, A. M. Shirokov, and J. P. Vary, Phys. Rev. C 81, 021301 (2010).
P. Navrátil, S. Quaglioni, I. Stetcu, and B. R. Barrett, J. Phys. G 36, 083101 (2009).
G. Hagen, T. Papenbrock, D. J. Dean, and M. Hjorth-Jensen, Phys. Rev. Lett. 101, 092502 (2008)
http://physics.aps.org/articles/v4/38#c1

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s