Radiative strength functions from the energy-localized Brink-Axel hypothesis

This paper introduces an energy-localized Brink-Axel hypothesis variant of the shell model Lanczos strength-function method to efficiently compute radiative strength functions for use in Hauser-Feshbach reaction codes, demonstrating its validity on $^{24}$Mg and revealing that while M1 and E1 transitions significantly contribute below the photo-absorption threshold in $^{56}$Fe, current model spaces cannot fully reproduce the low-energy strength observed in Oslo-type experiments.

The Gist

Imagine an atomic nucleus not as a tiny, solid marble, but as a bustling, chaotic city. When this city gets "excited" (heated up or hit by a particle), it tries to cool down by shooting out light particles called photons. Physicists need to predict exactly how much light comes out and at what colors (energies) to understand how stars are born and how nuclear reactors work.

The tool they use to make these predictions is called a Radiative Strength Function (RSF). Think of the RSF as a "traffic report" for the nucleus: it tells you how easy or hard it is for the nucleus to emit light at different energy levels.

For decades, scientists had a rule of thumb called the Brink-Axel hypothesis. It was like saying, "The traffic report for the city center (the ground state) is the same as the traffic report for the suburbs, no matter how hot the day is." This made calculations easy, but the authors of this paper argue it's not quite right.

Here is what this paper actually found and did, explained simply:

1. The Problem with the Old Map

The old way of calculating the RSF was like trying to map a city by looking at a single, frozen snapshot of one specific neighborhood. It worked okay for some things, but it failed to explain what happens when the nucleus is really hot and excited. Also, calculating the full map for every single possible state of a nucleus is like trying to count every single grain of sand on a beach—it takes too much computer power.

2. The New "Local" Map (The Energy-Localized Brink-Axel Hypothesis)

The authors propose a new idea: The traffic report changes depending on where you are in the city.

• If the nucleus is cool (ground state), it emits light in a specific, predictable pattern.
• If the nucleus is hot (highly excited), the pattern changes. Specifically, it starts emitting more low-energy light than the old rules predicted.

They call this the Energy-Localized Brink-Axel (ELBA) hypothesis. Instead of using one master map for the whole city, they suggest using a series of "local maps" that change slightly as the nucleus gets hotter.

3. The Shortcut: The "Lanczos" Flashlight

To prove this, they needed to calculate the light emission for thousands of different excited states. Doing this the old way would take a supercomputer years.

• The Analogy: Imagine trying to see the shape of a dark room. The old way was to turn on a light and take a photo of every single corner individually.
• The New Way: They used a method called the Lanczos Strength-Function (LSF) method. Think of this as a special flashlight that doesn't just show you one corner; it bounces light around the room and uses the echoes to instantly figure out the shape of the whole room without visiting every single spot.
• They combined this flashlight with their "local map" idea. They only needed to shine the light on a few specific excited states (a few "neighborhoods") and could accurately predict the behavior for the whole range of temperatures. This made the calculation 10 times faster and much more efficient.

4. Testing the Theory on Magnesium and Iron

They tested their new method on two elements:

• Magnesium-24: They compared their new "local map" against the old "master map." They found the new method was just as accurate but much simpler to calculate.
• Iron-56: This is the big test. Iron is crucial for understanding how stars explode and how elements are formed.
  • Finding A: They confirmed that as the Iron nucleus gets hotter, the way it emits light changes smoothly. The "low-energy" light (the "Low-Energy Enhancement" or LEE) gets stronger, just as their new hypothesis predicted.
  • Finding B: They found that both magnetic and electric types of light contribute to this glow, not just one type.
  • Finding C (The Limit): Even with their super-fast new method, they hit a wall. When they looked at the very lowest energy light (below 3 MeV) in Iron, their computer model could not fully reproduce what experiments (called Oslo-type experiments) actually see. There is still a "missing piece" of the puzzle that their current model space (the specific set of rules they used for the Iron nucleus) couldn't capture.

Summary

The paper doesn't claim to have solved every mystery of nuclear physics. Instead, it offers a better, faster way to draw the map of how nuclei emit light.

1. They proved that the "traffic report" (RSF) changes as the nucleus gets hotter, not just staying the same.
2. They built a "flashlight" (the Lanczos method) that lets them draw these changing maps quickly without needing to count every single grain of sand.
3. They applied this to Iron and saw the expected changes, but also admitted that for the very lowest energies, their current model still isn't perfect and needs more work.

In short: They made the map more accurate and the drawing process much faster, but they also pointed out exactly where the map is still incomplete.

Technical Summary

1. Problem Statement

Radiative Strength Functions (RSFs) are critical inputs for Hauser-Feshbach (HF) statistical reaction codes, which model nuclear reactions essential for astrophysics (nucleosynthesis) and nuclear technology. However, calculating RSFs from first principles (microscopic models) presents significant challenges:

• Computational Cost: Traditional Large-Scale Shell Model (LSSM) calculations require converging thousands of individual excited-state wavefunctions to capture the statistical average of decay widths, which is computationally prohibitive for heavy nuclei or high excitation energies.
• Theoretical Inconsistency: Standard RSF definitions often rely on the "strong Brink-Axel hypothesis," assuming the RSF is independent of the excitation energy of the initial state. Microscopic evidence suggests this is incorrect; RSFs should evolve with the energy of the de-excited level.
• Low-Energy Enhancement (LEE): Experimental data (e.g., from Oslo method) shows an unexpected enhancement of radiative strength at low energies (below 3 MeV). Existing microscopic models struggle to reproduce the magnitude of this LEE, particularly for electric dipole (E1) transitions in nuclei like $^{56}$Fe.
• Definition Ambiguity: There is a disconnect between the "sum-rule" strength functions used in microscopic theory and the RSFs required by HF reaction codes, leading to confusion regarding level densities and statistical factors.

2. Methodology

The authors propose a unified framework combining a new theoretical definition with an efficient computational algorithm.

A. Theoretical Framework: Level-Density-Free (LDF) Definition

• Reformulation: The authors redefine the RSF to be independent of level density ($\rho$). Instead of averaging over initial states and multiplying by a level density, they define the RSF as an energy-averaged sum of partial decay widths from Compound Nucleus (CN) levels ($c$) to a specific de-excited level ($d$).
• Equation: The new definition is given by:
  $$f^{XL}_{d,j_c}(E_\gamma) = \frac{1}{E_\gamma^{2L+1}} \frac{1}{\Delta E} \sum_{c'} \delta_{j_{c'} j_c} \Gamma^{XL}_{d \leftarrow c'}(E_\gamma)$$
  This formulation relies solely on partial widths ($\Gamma$), which are natural outputs of shell model calculations, avoiding approximations regarding level density distributions.
• Energy-Localized Brink-Axel (ELBA) Hypothesis: The paper moves away from the strong Brink-Axel hypothesis. Instead, it adopts the ELBA hypothesis, which posits that the total strength of a transition operator evolves smoothly with the energy of the de-excited level ($E_d$). This allows the RSF to be calculated for specific energy windows rather than assuming a global average.

B. Computational Technique: Interior-Eigenstate Lanczos Strength-Function (ILSF)

• Lanczos Method: To avoid calculating thousands of wavefunctions, the authors utilize the Lanczos Strength-Function (LSF) method. This algorithm resolves the moments of the strength distribution rather than computing individual transition probabilities.
• Interior Eigenstates: Unlike standard LSF which targets the ground state, this method targets "interior" excited states (pivots) located several MeV above the ground state.
• Averaging: By calculating the strength function for a small number of nearby excited states (e.g., 5 levels) and averaging them, the method approximates the global RSF.
• Efficiency: This approach reduces computational cost by more than an order of magnitude compared to brute-force LSSM, as it avoids storing massive numbers of vectors and converging thousands of states.

3. Key Contributions

1. New RSF Definition: Established a Level-Density-Free (LDF) definition of RSF that is mathematically consistent with Hauser-Feshbach theory and directly compatible with microscopic shell model outputs.
2. ILSF Algorithm: Developed and validated the Interior-Eigenstate Lanczos Strength-Function (ILSF) method, enabling the calculation of RSFs across a wide energy range using a manageable number of pivot states.
3. Validation of ELBA: Demonstrated that the shape of the RSF is not static but evolves smoothly with the excitation energy of the de-excited level, validating the Energy-Localized Brink-Axel hypothesis.
4. Separation of M1 and E1: Successfully separated and calculated both Magnetic Dipole (M1) and Electric Dipole (E1) contributions to the low-energy strength in $^{56}$Fe.

4. Results

The authors applied their method to $^{24}$Mg (for benchmarking) and $^{56}$Fe (for novel results).

Benchmark: $^{24}$Mg

• The LDF-RSF definition showed excellent agreement with the traditional Bartholomew prescription (pixel-averaged method) when sufficient levels were included.
• The method confirmed that the standard deviation of strength functions (Porter-Thomas fluctuations) decreases as the number of included states increases, validating the statistical approach.

Novel Results: $^{56}$Fe

• Evolution of M1 RSF: The M1 strength function shape evolves with excitation energy.
  • Ground State: Shows a standard Giant Dipole Resonance (GDR) shape with no strength below ~4 MeV.
  • Excited States: As the excitation energy of the de-excited level increases, a low-energy enhancement (LEE) appears.
  • High Energy Saturation: At very high excitation energies (above 10-15 MeV), the slope of the LEE flattens, leading to a nearly constant distribution below the GDR. This flattening is attributed to the mixing of valence orbitals at higher energies.
• E1 Contribution: The study calculated the E1 RSF in the $sdpf$ model space.
  • E1 transitions contribute significantly to the radiative strength below the photo-absorption threshold.
  • The E1 strength dominates the total RSF above 3 MeV, while M1 dominates below 3 MeV.
• The "Missing" Strength: Despite including both M1 and E1 contributions in the $sdpf$ model space, the calculated strength below 3 MeV remains a factor of 3 to 10 lower than experimental Oslo-type data.
  • The authors conclude that the current $sdpf$ model space (lacking the $1g_{9/2}$ orbital) is insufficient to reproduce the full magnitude of the LEE.
  • The GDR centroid is predicted at 16 MeV, lower than the experimental value (18.6 MeV), suggesting missing physics in the model space.

5. Significance and Future Outlook

• Microscopic Foundation: This work provides the first coherent microscopic description of RSFs that includes both M1 and E1 contributions over a wide energy range, moving beyond phenomenological models.
• Reaction Code Improvement: The results strongly motivate the use of energy-dependent RSFs in modern Hauser-Feshbach reaction codes. Instead of a single, static RSF, codes should accept tabulated RSFs that depend on both the photon energy ($E_\gamma$) and the excitation energy of the de-excited level ($E_d$).
• Path Forward: The inability to reproduce the sub-3 MeV strength in $^{56}$Fe highlights the need for larger model spaces (including the $g_{9/2}$ orbital) and improved effective interactions. The ILSF method provides the computational feasibility to tackle these larger spaces in the future.
• Efficiency: The ILSF method makes it practical to generate microscopic RSFs for a wide range of nuclei, facilitating more accurate predictions for neutron capture rates in nucleosynthesis and nuclear data evaluations.