Detailed Study of the $^{59}$Cu(p,$α)^{56}$Ni Reaction and Constraints on Its Astrophysical Reaction Rate

This paper presents a direct measurement of the $^{59}$Cu$(p,\alpha)^{56}$Ni reaction cross section using the MUSIC detector at FRIB, which, combined with Bayesian model averaging, yields a revised astrophysical reaction rate that is systematically lower than previous evaluations and establishes the competing $(p,\gamma)$ channel as the dominant uncertainty in the NiCu cycle.

The Gist

Imagine the universe as a giant, chaotic kitchen where stars are the chefs, constantly cooking up new elements. Sometimes, these chefs work in a calm, slow oven (like our Sun), but other times, they work in a frantic, explosive kitchen, like during a Type I X-ray burst (a star exploding) or the aftermath of a massive supernova. In these high-pressure, super-hot environments, the "recipe" for creating heavy elements depends entirely on how fast tiny particles collide and react.

This paper is about a specific, crucial "ingredient swap" in that cosmic recipe: the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ reaction.

Here is the story of what the scientists did, explained simply:

1. The Problem: A Traffic Jam in the Cosmic Kitchen

In these explosive stellar events, there is a specific bottleneck called the NiCu cycle. Think of this cycle as a roundabout in a busy city.

• The Goal: The universe wants to build heavier elements (like gold or zinc) by adding protons to atoms.
• The Obstacle: When the atom $^{59}\text{Cu}$ (Copper-59) gets hit by a proton, it has two choices:
  1. Keep the proton: It becomes heavier ($^{60}\text{Zn}$), allowing the recipe to continue toward heavier elements.
  2. Spit out a particle (an alpha particle): It turns back into a lighter atom ($^{56}\text{Ni}$), getting stuck in a loop.

If the "spit out" reaction happens too often, the cosmic traffic jams at the roundabout, and no heavy elements get made. If it happens rarely, the traffic flows, and heavy elements are created. For a long time, scientists didn't know exactly how often this "spit out" reaction happened, so they couldn't predict how the universe cooks heavy elements.

2. The Experiment: Catching the Reaction in Action

To solve this, the team went to the Facility for Rare Isotope Beams (FRIB) in Michigan. They used a massive, high-tech detector called MUSIC (Multi-Sampling Ionization Chamber).

• The Setup: Imagine shooting a stream of unstable Copper-59 atoms (the "bullets") into a tank of gas (methane).
• The Collision: When a Copper atom hits a gas proton, they react. Sometimes the Copper spits out an alpha particle (a helium nucleus) and turns into Nickel-56.
• The Detection: The MUSIC detector is like a super-sensitive 3D camera. It doesn't just take a picture; it tracks the exact path and energy loss of every particle. It can tell the difference between a Copper atom that just bounced off (scattering) and one that actually reacted and changed its identity.
• The Result: They measured this reaction at lower energies than ever before. This is crucial because the "cooking" in stars happens at very specific, lower energy levels than previous experiments could reach.

3. The Analysis: Tuning the Cosmic Recipe Book

Measuring the reaction is only half the battle. To know what happens in a star, they had to predict how the reaction behaves at even lower temperatures (energies) that they couldn't physically test in the lab.

• The Model: They used a computer program called TALYS, which acts like a cosmic recipe book, predicting how particles should behave based on physics rules.
• The Problem: The standard recipe book had been guessing wrong in the past. It was like using a map that said "turn left" when you actually needed to "turn right."
• The Fix: The team used a statistical method called Bayesian Model Averaging. Imagine asking 96 different expert chefs (models) for their opinion on the recipe. Instead of picking just one, they weighed the opinions of all 96 based on how well their predictions matched the new experimental data.
• The Optimization: They fine-tuned the "geometry" of the interaction (how the particles approach each other) until the computer model perfectly matched their new data.

4. The Discovery: The Traffic Jam is Less Severe

The results changed the understanding of the NiCu cycle:

• The Rate is Lower: The new, experimentally confirmed rate of the "spit out" reaction is lower than what was previously thought (specifically lower than the standard REACLIB database).
• The Consequence: Because the "spit out" reaction happens less often than we thought, the traffic jam at the NiCu roundabout is less severe. The "keep the proton" path is more likely to win.
• The New Bottleneck: Since the "spit out" reaction is now well-understood and less of a problem, the main uncertainty in the recipe is no longer this reaction. Instead, the uncertainty lies in the other reaction: $^{59}\text{Cu}(p, \gamma)^{60}\text{Zn}$ (the one where the atom keeps the proton).

Summary

In simple terms, this paper is like a team of mechanics who finally measured exactly how often a specific car engine part fails. They found it fails less often than the manual said. Because of this, they realized the car isn't stuck in traffic as much as we thought. However, now that they know this part works fine, they realize the real problem causing traffic jams is a different part of the engine that hasn't been measured as well yet.

Key Takeaway: The scientists measured a specific nuclear reaction in a lab, proved it happens less frequently than previously estimated, and concluded that this reaction is no longer the main reason heavy elements struggle to form in exploding stars. The focus must now shift to understanding a different reaction to fully solve the mystery of how the universe creates heavy elements.

Technical Summary

Technical Summary of "Detailed Study of the 59Cu(p,α)56Ni Reaction and Constraints on Its Astrophysical Reaction Rate"

Problem Statement
The $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ reaction is a critical component in explosive astrophysical environments, specifically Type I X-ray bursts (XRBs) and the $\nu p$-process in neutrino-driven winds from core-collapse supernovae. This reaction competes with the $^{59}\text{Cu}(p, \gamma)^{60}\text{Zn}$ reaction to regulate the nucleosynthesis flow through the NiCu cycle. At temperatures below $\approx 3$ GK, the $(p, \gamma)$ channel dominates, allowing the synthesis of heavier nuclei. However, at higher temperatures, the $(p, \alpha)$ channel can become dominant, trapping the reaction flow in a closed NiCu cycle and inhibiting the production of elements beyond the iron group. Despite its importance, experimental data for the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ reaction at astrophysically relevant low energies have been scarce. Previous direct measurements were limited to higher energies or relied on angular integration models that introduced significant uncertainties in the total cross-section determination.

Methodology
The authors performed a direct measurement of the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ excitation function using inverse kinematics at the Facility for Rare Isotope Beams (FRIB).

• Experimental Setup: A $^{59}\text{Cu}$ beam (produced via fragmentation of a $^{64}\text{Zn}$ primary beam and reaccelerated to 8.418 MeV/u) was directed into the Multi-Sampling Ionization Chamber (MUSIC), an active-target detector filled with 440 Torr of methane gas. The beam was stopped completely within the detector.
• Event Identification: The MUSIC detector, with its segmented anode, allowed for particle identification via energy loss ($\Delta E$) measurements. The experiment distinguished $(p, \alpha)$ events (producing $^{56}\text{Ni}$ recoils) from unreacted $^{59}\text{Cu}$ beams and other channels like $(p, p)$ and $(p, p')$ based on distinct energy-loss traces and $\Delta E$–$\Delta E$ spectra.
• Data Analysis: Cross sections were measured over a center-of-mass energy range of 2.43–5.88 MeV. Effective energies were calculated accounting for thick-target yield effects. Statistical uncertainties were derived using Poisson statistics or the Feldman–Cousins approach for low-count strips, while systematic uncertainties were evaluated based on event selection criteria and detector resolution.
• Theoretical Extrapolation: To derive the astrophysical reaction rate, the experimental data were extrapolated to lower energies using Hauser–Feshbach statistical-model calculations performed with the TALYS 2.0 code.
  • Optical Model Optimization: The geometry of the Demetriou–3 (DEM-3) $\alpha$-optical model potential (OMP) was systematically optimized to reproduce the measured cross sections without global renormalization.
  • Bayesian Model Averaging (BMA): To quantify model-selection uncertainty, a BMA analysis was conducted over 96 combinations of TALYS parameters, including 8 $\alpha$-OMP parametrizations, 6 level-density models, and variations in the number of discrete states ($msl$) in $^{59}\text{Cu}$.
  • Stellar Rate Calculation: The laboratory rate was converted to a stellar rate by calculating the stellar enhancement factor (SEF), which accounts for transitions from thermally populated excited states in the stellar environment.

Key Results

• Cross Sections: The study provides the first direct, angle- and energy-integrated cross sections for $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ down to $E_{\text{c.m.}} \approx 2.43$ MeV. The results extend previous measurements to lower energies and eliminate the model dependence associated with angular integration in earlier works.
• Model Constraints: The DEM-3 $\alpha$-OMP, with optimized geometry parameters (radius $r_v$ increased by 6% and diffuseness $a_v$ decreased by 6%), provided the best reproduction of the experimental data. The BMA analysis confirmed DEM-3 as the dominant model, accumulating approximately 40% of the total weight.
• Reaction Rate: The recommended stellar reaction rate is systematically lower than the REACLIB evaluation for temperatures $T_9 \lesssim 3$. The rate carries a temperature-dependent uncertainty factor of 1.26–1.63 over $T_9 = 0.2–10$.
• NiCu Cycle Strength: A comparison with the $^{59}\text{Cu}(p, \gamma)^{60}\text{Zn}$ rate (using REACLIB and O'Shea et al. constraints) shows that the $(p, \alpha)$ rate remains consistently below the $(p, \gamma)$ rate across the relevant temperature range ($T_9 \lesssim 3$). This indicates that the NiCu cycle is weak under these conditions, and the reaction flow is not significantly confined to the cycle.

Significance and Claims
The paper claims that these results substantially weaken the inferred strength of the NiCu cycle in explosive astrophysical scenarios. By providing a robust, experimentally constrained reaction rate with quantified uncertainties, the study removes the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ rate as a primary source of uncertainty in the NiCu cycle branching for XRBs and the $\nu p$-process. Consequently, the authors identify the $^{59}\text{Cu}(p, \gamma)^{60}\text{Zn}$ reaction rate as the dominant remaining uncertainty in determining whether nucleosynthesis can proceed beyond the iron group in these proton-rich environments. The work establishes a new benchmark for the $(p, \alpha)$ rate, demonstrating that it is lower than previously assumed in some evaluations, thereby favoring the continuation of the nucleosynthesis flow toward heavier elements rather than a closed cycle.