Resolving anomalous collectivity in the $4_1^+$ to $2_1^+$ transition of $^{58}$Fe

This study resolves a long-standing discrepancy in the $4_1^+$ to $2_1^+$ transition strength of $^{58}$Fe by using Coulomb-excitation measurements to demonstrate that previous lifetime data were overestimated due to outdated electronic stopping power models, thereby confirming the validity of shell-model predictions.

The Gist

The Mystery of the "Overachieving" Iron Nucleus

Imagine the atomic nucleus as a tiny, bustling city. Inside this city, protons and neutrons (the citizens) live in specific neighborhoods called "shells." Physicists have a very successful map of this city called the Shell Model. For decades, this map has predicted exactly how these citizens behave, how they move, and how much energy they use when they jump between neighborhoods.

Most of the time, the map is perfect. But in the city of Iron-58 (a specific type of iron atom), there was a confusing traffic jam.

The Problem: A Speeding Ticket That Didn't Make Sense

In the Iron-58 city, there is a specific event where a citizen jumps from a high-energy neighborhood (the $4^+$ state) down to a lower one (the $2^+$ state).

According to the Shell Model map, this jump should be a casual stroll. However, previous measurements taken in the 1970s and 80s suggested this jump was a sprint. The data said the nucleus was "collective"—meaning all the citizens were moving together in a big, coordinated dance, making the jump much stronger than the map predicted.

This was a problem. If the map is wrong, our whole understanding of how atomic nuclei work might be shaky. Scientists were puzzled: Is Iron-58 actually doing a special dance, or is our measuring tape broken?

The New Investigation: A Fresh Look with Better Tools

The authors of this paper decided to re-measure this jump, but they wanted to avoid the old methods that might have been flawed.

The Old Method (The "Stopwatch" Problem):
Previous scientists used a technique called the Doppler Shift Attenuation Method (DSAM).

• The Analogy: Imagine a race car speeding down a track while honking its horn. As it slows down, the pitch of the horn changes. By measuring how fast the pitch changes, you can calculate how long the car was slowing down (its lifetime).
• The Flaw: To do this math, you need to know exactly how much friction the road (the target material) exerts on the car. The old scientists used an old friction chart from the 1960s (called LSS stopping powers). It turns out, that old chart said the road was much "stickier" than it actually is. Because they thought the road was stickier, they calculated the car was slowing down faster than it really was, leading to a wrong conclusion about how long the race lasted.

The New Method (The "High-Speed Camera"):
Instead of timing how long the nucleus lives, these researchers used Coulomb Excitation.

• The Analogy: Imagine two magnets flying past each other. They don't crash, but their magnetic fields tug on each other, making one spin faster or jump to a new position. By measuring exactly how hard they tugged and how much energy was transferred, you can figure out the properties of the jump without needing to time how long it lasts.
• They fired beams of Iron-56 and Iron-58 at a gold target and watched the "tug-of-war" using a giant array of detectors (like a high-speed camera capturing the event from all angles).

The Results: The Map Was Right All Along

When the team analyzed the new data:

1. Iron-56: The jump matched the Shell Model perfectly. The map was correct here.
2. Iron-58: The new measurement showed the jump was much weaker than the old "sprint" data suggested. It was actually a casual stroll, just as the Shell Model predicted.

The "collective dance" was an illusion caused by a bad friction chart. Iron-58 isn't breaking the rules; it was just measured with a broken ruler.

The "Re-Do": Fixing the Old Data

The team didn't just ignore the old data; they went back and fixed it. They took the original measurements from the 1970s and recalculated them using a modern, accurate friction chart (called SRIM).

• The Result: When they used the correct friction numbers, the old "sprint" turned into a "stroll." The lifetime of the state became longer, and the strength of the jump became weaker. Suddenly, the old data agreed with the new data and the Shell Model.

Why This Matters

This paper is a victory for the Shell Model. It confirms that even in complex atomic cities, the rules of the map hold true. It also serves as a warning to all scientists: Check your friction charts!

If you use outdated data to calculate how fast things slow down (stopping powers), you can create "ghost" phenomena that don't actually exist. In this case, a 40-year-old measurement error made it look like Iron-58 was doing something special, when it was actually just behaving normally.

In short: The Iron-58 nucleus wasn't doing a special dance; we just had the wrong music playing in the background. Once we fixed the music, the dance looked normal again.

Technical Summary

1. Problem Statement

The nuclear structure of iron isotopes ($^{56}\text{Fe}$, $^{58}\text{Fe}$, $^{60}\text{Fe}$) near the $N=Z=28$ shell closure is generally well-described by the nuclear shell model. However, a significant discrepancy existed regarding the reduced electric quadrupole transition strength, $B(E2; 4^+_1 \to 2^+_1)$, in $^{58}\text{Fe}$.

• The Anomaly: Previous experimental values for this transition strength, derived from lifetime measurements using the Doppler Shift Attenuation Method (DSAM) and Doppler Broadened Line Shape (DBLS), yielded a value of approximately 47(7) W.u. (Weisskopf units).
• The Conflict: This experimental value was nearly double the prediction of large-basis shell-model calculations (which predicted $\approx 25$ W.u.) and deviated significantly from the systematics of neighboring isotopes. This suggested an unexplained emergence of collective behavior in $^{58}\text{Fe}$ that contradicted established theoretical models.
• The Goal: To resolve this discrepancy by performing a direct measurement of the $B(E2)$ ratio independent of lifetime-dependent Doppler-shift methods and by re-evaluating the historical lifetime data.

2. Methodology

A. New Coulomb-Excitation Measurement

The authors performed a new experiment at the Australian National University's Heavy Ion Accelerator Facility to determine the $B(E2; 4^+_1 \to 2^+_1)$ value directly via Coulomb excitation, avoiding the uncertainties of lifetime-based DSAM/DBLS analyses.

• Beam and Target: Beams of 220 MeV $^{56}\text{Fe}$ and $^{58}\text{Fe}$ were incident on a $^{197}\text{Au}$ target.
• Detection: The experiment utilized the CAESAR array, consisting of seven Compton-suppressed High-Purity Germanium (HPGe) detectors for $\gamma$-rays and an array of eight silicon photodiodes for particle detection.
• Technique: Particle-$\gamma$ coincidence spectroscopy was used to identify transitions. The relative yields of the $4^+_1 \to 2^+_1$ and $2^+_1 \to 0^+_1$ transitions were measured.
• Analysis: The Gosia code (semi-classical Coulomb-excitation code) was used to extract reduced matrix elements. The analysis involved:
  • Fixing the reference transition $B(E2; 2^+_1 \to 0^+_1)$ based on accepted literature values.
  • Varying matrix elements and testing various sign combinations (phases) of the E2 matrix elements to find the best fit to the experimental yields.
  • Using SRIM (Stopping and Range of Ions in Matter) code to calculate stopping powers for the simulation.

B. Re-evaluation of Historical DSAM Data

The authors re-examined the seminal 1978 DSAM measurement by Bolotin et al., which provided the basis for the anomalous $B(E2)$ value.

• Root Cause Investigation: The original analysis relied on Lindhard-Scharff-Schiøtt (LSS) electronic stopping powers.
• Modern Correction: The authors replaced the LSS stopping powers with modern values calculated using SRIM.
• Re-analysis: They reconstructed the original analysis code (based on logbooks and printed listings) and re-ran the lifetime extraction using the updated stopping powers, accounting for the fact that LSS values significantly overestimate electronic stopping compared to modern data.

3. Key Contributions and Results

A. Coulomb-Excitation Results

• $^{56}\text{Fe}$ (Validation): The measurement for $^{56}\text{Fe}$ yielded $B(E2; 4^+_1 \to 2^+_1) = 23(4)$ W.u., which agreed perfectly with the adopted literature value and shell-model predictions. This validated the experimental setup and analysis procedures.
• $^{58}\text{Fe}$ (Resolution): For $^{58}\text{Fe}$, the extracted $B(E2; 4^+_1 \to 2^+_1)$ values were significantly lower than the previously adopted value.
  • Depending on the sign combinations of the matrix elements (which could not be uniquely determined by the data), the values ranged from 19 to 27 W.u.
  • The maximum value (27 W.u.) is still significantly lower than the historical 47 W.u. but is in excellent agreement with shell-model calculations ($\approx 25$ W.u.).
  • This result eliminates the evidence for "anomalous collectivity" in $^{58}\text{Fe}$.

B. Re-evaluation of DSAM Lifetimes

• Stopping Power Discrepancy: The re-analysis revealed that the LSS electronic stopping powers used in 1978 were approximately twice as large as the values provided by modern SRIM calculations for the energy regime of the experiment.
• Impact on Lifetime: Because the stopping power determines how quickly the recoiling nucleus slows down (and thus how much Doppler shift is observed), overestimating the stopping power leads to an underestimation of the lifetime.
• Corrected Lifetime: When using SRIM stopping powers, the lifetime of the $4^+_1$ state in $^{58}\text{Fe}$ was revised from $\tau \approx 0.34$ ps (historical) to $\tau = 0.66(8)$ ps.
• Consistency: This revised lifetime yields a $B(E2)$ value consistent with the new Coulomb-excitation measurement and shell-model predictions.

C. Shell-Model Comparison

• Using the GXPF1A interaction with effective charges ($e_\pi=1.5, e_\nu=0.5$), the shell-model calculations for $^{58}\text{Fe}$ successfully reproduced the experimental $B(E2)$ values for the yrast states (including the corrected $4^+_1 \to 2^+_1$ transition) without needing to invoke exotic collective structures.
• The study confirmed that the "anomaly" was purely an artifact of outdated stopping power data used in lifetime measurements.

4. Significance

1. Resolution of a Long-standing Discrepancy: The paper definitively resolves a decades-old puzzle regarding the structure of $^{58}\text{Fe}$. It demonstrates that the nucleus behaves as predicted by the shell model, rather than exhibiting the unexpected collective rotation suggested by previous data.
2. Critique of Historical Data: It highlights a critical vulnerability in nuclear physics: lifetime measurements derived from DSAM/DBLS are highly sensitive to the accuracy of stopping power models. The authors warn that historical data relying on LSS stopping powers (pre-1990s) should be treated with extreme caution and likely require re-evaluation.
3. Methodological Validation: The successful application of the new CAESAR array and the multi-step Coulomb excitation technique at ANU demonstrates the capability to measure transition strengths directly, bypassing the systematic errors inherent in lifetime-based methods.
4. Theoretical Confirmation: The results reinforce the validity of large-basis shell-model calculations (specifically the GXPF1A interaction) for describing nuclear structure in the $N=Z=28$ region, confirming that the $fp$-shell model remains robust even for neutron-rich iron isotopes.

In conclusion, the "anomalous collectivity" in $^{58}\text{Fe}$ was a phantom effect caused by outdated stopping power corrections in historical lifetime measurements. The true nuclear structure is well-described by the standard shell model.