Direct Measurement of the $^{59}$Cu$(p,α)^{56}$Ni Excitation Function to Constrain the Ni--Cu Cycle Strength and Its Impact on Explosive Nucleosynthesis

A new direct measurement of the $^{59}$Cu$(p,\alpha)^{56}$Ni reaction using the MUSIC detector at FRIB reveals a systematically lower stellar rate than previously estimated, which significantly suppresses Ni-Cu cycle recycling in X-ray bursts while enhancing the efficiency of the $\nu$p process in supernova nucleosynthesis.

The Gist

Imagine the universe as a giant, chaotic kitchen where stars are the chefs. Sometimes, these chefs get so hot and energetic that they cook up new ingredients (elements) in a flash. Two of the most dramatic cooking scenarios are Type I X-ray bursts (explosions on the surface of dead stars called neutron stars) and the neutrino-driven winds (hot, fast outflows of gas after a massive star explodes).

In these super-hot kitchens, the chefs try to build heavier elements by smashing protons (hydrogen nuclei) onto existing atoms. But there's a tricky traffic jam waiting to happen at a specific "intersection" involving a rare atom called Copper-59.

The Traffic Jam: The NiCu Cycle

Think of Copper-59 as a busy intersection. When a proton hits it, the atom has two choices:

1. The Exit Ramp (p, γ): It grabs the proton and becomes heavier (Zinc-60), allowing the cooking process to continue building even heavier elements.
2. The U-Turn (p, α): It spits out a chunk (an alpha particle) and turns back into Nickel-56. This is like a car doing a U-turn and going back to the start of the line.

This U-turn is called the NiCu Cycle. If the U-turn happens too often, the heavy elements never get built. If the Exit Ramp is open, the cooking continues. Scientists needed to know exactly how often the U-turn happens to understand how much heavy stuff the universe can make.

The Experiment: Catching the U-Turn

For a long time, scientists had to guess how often this U-turn happened because it's incredibly hard to measure. Previous guesses were like trying to guess the speed of a car by looking at its tire tracks from far away—they had to assume a lot about the road conditions.

In this new study, researchers at the Facility for Rare Isotope Beams (FRIB) decided to measure it directly.

• The Setup: They created a beam of Copper-59 atoms (which are unstable and hard to make) and fired them into a tank of methane gas.
• The Detector: They used a special "active target" detector called MUSIC. Think of this detector as a giant, high-tech honeycomb. When the Copper atoms hit the gas, they sometimes collide with protons in the gas.
• The Measurement: If a U-turn happens (the Copper spits out an alpha particle), the detector sees the specific energy signature of the resulting Nickel atom. By counting these events at different speeds, they mapped out exactly how likely the U-turn is to happen across a wide range of temperatures.

The Big Discovery: The U-Turn is Rarer Than We Thought

The results were surprising. The new measurements showed that the U-turn (p, α) happens much less often than scientists previously thought.

• Old View: We thought the traffic jam was heavy; the NiCu cycle was recycling a lot of material back to the start, stopping the creation of heavy elements.
• New View: The traffic jam is actually light. The "Exit Ramp" is much more open than we expected.

Why This Matters for the Universe

This discovery changes our understanding of two cosmic cooking events:

1. X-ray Bursts (The Neutron Star Explosions):
   In these bursts, the new data suggests that the NiCu cycle recycles less than 0.74% of the material. This means the explosion is more efficient at building heavier elements than we thought, and the "ash" left behind will have a different chemical makeup.

2. The Neutrino-Driven Wind (The Supernova Outflow):
   This is where the universe tries to make elements heavier than Iron. Because the U-turn is weaker, the "Exit Ramp" stays open for longer.

   • The Result: The process can keep building heavier elements at higher temperatures than previously predicted.
   • The Limit: Instead of stopping at a certain point, the process can now push further, potentially creating elements up to a mass number of 109 (instead of stopping around 107). It also shifts the "crossover point" (where the process decides to stop recycling and start building heavy stuff) to a higher temperature, meaning it happens closer to the center of the explosion where the energy is strongest.

The Bottom Line

By directly measuring this specific nuclear reaction, the scientists have removed a huge guess from the recipe of the universe. They found that the "NiCu Cycle" is a much weaker traffic jam than we thought. This means the universe is likely better at cooking up heavy elements in these explosive events than our old models suggested.

The only thing left to figure out is exactly how often the "Exit Ramp" (the proton capture) happens, as that is now the biggest remaining uncertainty in the recipe. But thanks to this experiment, we have a much clearer picture of how the heavy elements in our universe are made.

Technical Summary

Technical Summary: Direct Measurement of the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ Excitation Function

Problem Statement
Explosive nucleosynthesis in proton-rich astrophysical environments, specifically Type I X-ray bursts (XRBs) and the neutrino-driven winds of core-collapse supernovae (the $\nu$p-process), relies on the competition between proton-capture and charged-particle emission reactions. A critical bottleneck in these scenarios is the branching point at $^{59}\text{Cu}$, where the reaction flow is determined by the competition between the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ reaction (which recycles material back to the Ni region) and the $^{59}\text{Cu}(p, \gamma)^{60}\text{Zn}$ reaction (which allows breakout to heavier nuclei). This competition defines the strength of the "NiCu cycle."

Prior to this work, constraints on the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ rate were limited. Previous direct measurements provided only sparse data points at higher energies ($E_{\text{c.m.}} \approx 6.0$ MeV) or relied on indirect spectroscopic information, leading to significant model dependencies and uncertainties. Statistical model predictions (e.g., NON-SMOKER) had been scaled to match limited data, but the energy dependence and absolute magnitude of the cross section remained poorly constrained, particularly at the lower energies relevant to the $\nu$p-process ($T_9 \sim 1\text{--}3$).

Methodology
The authors performed a new 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).

• Beam and Target: An $^{59}\text{Cu}$ beam (8.418 MeV/u, $\sim 94\%$ purity) was generated via fragmentation of a primary $^{64}\text{Zn}$ beam. The beam was stopped within the Multi-Sampling Ionization Chamber (MUSIC), an active-target detector filled with methane gas at 440 Torr.
• Detection: The MUSIC detector, segmented into 18 strips, allowed for the identification of reaction products. The $^{56}\text{Ni}$ recoil nuclei were distinguished from unreacted $^{59}\text{Cu}$ beam particles via their lower energy loss ($\Delta E$). By summing energy deposition in strips downstream of the reaction point, the authors achieved a clean separation of reaction channels.
• Data Extraction: Cross sections were extracted for eight effective center-of-mass energy bins ranging from $E_{\text{c.m.}} = 2.43$ to $5.88$ MeV. Systematic uncertainties were evaluated by varying event-selection criteria.
• Rate Determination: To derive the stellar reaction rate, the experimental cross sections were combined with Hauser-Feshbach statistical model calculations (using the TALYS code). A Bayesian Model Averaging (BMA) analysis was performed over 96 combinations of $\alpha$-optical model potentials ($\alpha$-OMP), level densities, and discrete-level cutoffs to identify the optimal configuration. The geometry of the DEM-3 $\alpha$-OMP was further optimized (adjusting radius and diffuseness) to reproduce the data without post-hoc scaling.

Key Results

1. Excitation Function: The study provides the first direct experimental input for the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ reaction at energies relevant to the $\nu$p-process ($E_{\text{c.m.}} \ge 2.43$ MeV). The measured cross sections are systematically lower than previous statistical model predictions (NON-SMOKER) and differ in energy dependence from the scaled NON-SMOKER curve used in recent literature.
2. Stellar Reaction Rate: The newly derived stellar rate is systematically lower than previous estimates, including those from Ref. [27] and the standard REACLIB library. The uncertainty on the recommended rate is reduced to a factor of 1.26–1.63 over the temperature range $T_9 = 0.2\text{--}10$, a significant improvement over the factors of 10–100 previously used in sensitivity studies.
3. NiCu Cycle Strength: The branching ratio $B_{p\alpha/p\gamma}$ (the ratio of $(p, \alpha)$ to $(p, \gamma)$ rates) is found to be lower than previously estimated.
   • XRBs: At $T_9 = 1$, the recycling of material via the NiCu cycle is constrained to $0.03\text{--}0.74\%$, indicating very weak recycling under XRB conditions.
   • $\nu$p-process: The crossover temperature where the $(p, \gamma)$ channel dominates (breaking out of the cycle) shifts to $T_9 = 3.94^{+0.99}_{-0.85}$, compared to $T_9 \approx 3.2$ predicted by REACLIB.

Astrophysical Impact and Significance
The paper claims that these new constraints significantly alter the understanding of explosive nucleosynthesis in proton-rich environments:

• Enhanced $\nu$p-process Efficiency: The lower $(p, \alpha)$ rate and the higher crossover temperature imply that breakout from the NiCu cycle occurs closer to the proto-neutron star surface. In this region, stronger antineutrino fluxes enhance $(n, p)$ reactions, thereby improving the overall efficiency of the $\nu$p-process.
• Constrained Endpoint: Nucleosynthesis calculations using the new rate show a narrowing of the predicted endpoint of the $\nu$p-process. The final mass number distribution is constrained to $A \sim 106\text{--}109$, compared to the broader range of $A \sim 107\text{--}119$ derived using previous unconstrained rates.
• Reduced Uncertainty: The work demonstrates that the NiCu cycle is substantially weaker than previously inferred. While the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ rate is now well-constrained, the paper notes that the $^{59}\text{Cu}(p, \gamma)^{60}\text{Zn}$ rate remains the dominant source of uncertainty in the NiCu cycle branching ratio.

In summary, this direct measurement provides the strongest constraints to date on the $^{59}\text{Cu}(p, \alpha)^{56}\text{Ni}$ reaction, reducing uncertainties in the NiCu cycle strength and refining predictions for the production of heavy elements in XRBs and core-collapse supernovae.