A microscopic analysis of sub-barrier photo-induced fission in $^{236}$U$(\gamma,f)$ based on the non-equilibrium Green function method

This paper employs the non-equilibrium Green function method to investigate sub-barrier photo-induced fission in $^{236}$U, successfully reproducing experimental cross-section data and demonstrating that the first eigenchannel dominates the fission probability, thereby providing microscopic support for the Bohr-Wheeler transition-state picture.

The Gist

The Big Picture: Breaking a Heavy Ball with a Flash of Light

Imagine a giant, wobbly ball of clay (the Uranium-236 nucleus). Usually, this ball is very stable and holds its shape together tightly. However, if you hit it just right, it can split in half, releasing a massive amount of energy. This is nuclear fission.

In this specific study, scientists are looking at what happens when you try to split this ball not by hitting it with a hammer (like a neutron), but by shining a very specific, low-energy flash of light (a gamma-ray photon) on it.

The tricky part? The energy of the light is too low to simply smash the ball apart. It's like trying to break a rock with a gentle tap. In classical physics, nothing should happen. But in the quantum world, there's a phenomenon called tunneling. It's as if the rock has a ghostly ability to "phase" through the wall and break apart anyway, even though it didn't have enough energy to climb over the wall.

The author, K. Uzawa, used a sophisticated mathematical tool called the Non-Equilibrium Green Function (NEGF) method to simulate this process. Think of NEGF as a super-advanced "traffic simulator" for subatomic particles.

---

The Analogy: The Quantum Highway

To understand how the paper works, imagine the fission process as a highway system:

1. The Entrance (The Light): The gamma-ray hits the nucleus. This is the "on-ramp" where traffic (energy) enters the system.
2. The Tunnel (The Barrier): The nucleus has a "fission barrier"—a hill it must get over to split. Since the light energy is low, the traffic can't drive over the hill. They have to find a way to tunnel through it.
3. The Exit (The Split): Once through the tunnel, the nucleus splits into two smaller pieces (fission fragments). This is the "off-ramp."

The Problem with Old Simulations:
Previous methods (like Time-Dependent Density Functional Theory) are great at simulating cars driving over the hill, but they are terrible at simulating cars that tunnel through it. They can't easily calculate the odds of a car phasing through a mountain.

The New Solution (NEGF):
The author uses the NEGF method, which was originally designed to simulate electrons moving through computer chips (transistors).

• The Chip Analogy: In a transistor, electrons enter from a source, travel through a channel, and exit at a drain.
• The Nuclear Analogy: Here, the "source" is the light hitting the nucleus, the "channel" is the wobbly, deforming nucleus, and the "drain" is the split fragments.

By treating the nucleus like a quantum circuit, the author can calculate exactly how much "traffic" (probability) flows from the light, through the tunnel, and out the other side as a split nucleus.

---

How the Simulation Was Built

The author didn't just guess; they built a massive digital model:

1. The Map (The Path): They mapped out every possible shape the nucleus could take as it stretches and wobbles before splitting. They used a "fission path" which is like a trail of breadcrumbs showing the easiest route for the nucleus to deform.
2. The Traffic Jams (Excitations): As the nucleus stretches, the particles inside (protons and neutrons) get excited and jump around. The model accounts for these "jumps" (particle-hole excitations), which act like friction or energy dissipation, slowing the process down.
3. The Randomness: To make the simulation realistic, they added a bit of "random noise" to the interactions between particles. This mimics the chaotic nature of the quantum world, ensuring the model doesn't get stuck in a perfect, unrealistic loop.

---

The Results: Did the Light Work?

The author ran the simulation for gamma-ray energies between 5 and 6 MeV. This is the "sub-barrier" zone—the zone where the energy is too low to break the nucleus by force.

• The Prediction: The simulation predicted that the nucleus would split, but the probability would drop sharply as the light energy got lower.
• The Reality Check: They compared their results to real-world experimental data.
  • Above the barrier (High energy): The simulation matched the real data very well.
  • Below the barrier (Low energy): This is the big win. The simulation successfully reproduced the "suppression" (the drop in splitting probability) seen in experiments. It showed that even with low energy, the nucleus can split via tunneling, but it's much harder.

The "One Lane" Discovery

One of the most fascinating findings in the paper comes from an "eigenchannel analysis."

Imagine the highway has many lanes. You might expect that traffic could flow through any of them. However, the analysis showed that almost all the traffic flows through just one single lane (the first eigenchannel).

• What this means: Even though the nucleus is a complex mess of billions of particles, the path to splitting is surprisingly simple. It's like a crowded city where, despite having thousands of streets, 99% of the traffic to a specific destination funnels through a single, narrow tunnel.
• Why it matters: This supports an old theory called the Bohr-Wheeler transition-state picture, which suggests that fission happens through a specific "gateway" state. The author's microscopic simulation proves this theory is correct, even when looking at the quantum details.

Why Should We Care?

1. Understanding the Basics: This helps us understand the fundamental rules of how matter behaves at the smallest scales.
2. Nuclear Energy & Safety: Better models mean we can predict how nuclear reactors behave more accurately, especially in rare or extreme conditions.
3. The Universe: Understanding how heavy elements split helps astrophysicists understand how elements are created in neutron star collisions (the r-process), which is how the universe gets heavy elements like gold and uranium.

Summary

In simple terms, this paper is a success story of using a "traffic simulator" (NEGF) to predict how a heavy nucleus splits when hit by weak light.

The author showed that:

1. The method works for low-energy "tunneling" scenarios where older methods fail.
2. The results match real-world experiments.
3. Despite the complexity of the nucleus, the splitting process is dominated by a single, primary pathway, confirming long-held theories about how nuclear fission works.

It's a bridge between the messy, chaotic quantum world and the clean, predictable laws of physics that govern our universe.

Technical Summary

1. Problem Statement

The paper addresses the challenge of describing nuclear fission on a microscopic level, particularly for photo-induced fission in the sub-barrier energy region (incident gamma-ray energies below the fission barrier).

• Limitations of Existing Methods:
  • Time-Dependent Density Functional Theory (TDDFT): While effective for large-amplitude deformation and dissipation, it cannot describe quantum tunneling through the fission barrier, making it unsuitable for calculating fission probabilities below the barrier.
  • Time-Dependent Generator Coordinate Method (TDGCM): Capable of describing tunneling but often lacks a consistent treatment of the initial state preparation for induced reactions (i.e., defining the wave packet from the reaction entrance channel).
• Specific Challenge: Calculating fission cross-sections for $^{236}\text{U}(\gamma, f)$ requires a framework that naturally handles the transition from the photo-absorption entrance channel to the fission exit channel, including quantum tunneling and particle-hole excitations, without arbitrary initial state assumptions.

2. Methodology

The author employs the Non-Equilibrium Green Function (NEGF) method, a framework originally developed for electron transport in mesoscopic systems, adapted here for nuclear fission.

A. Theoretical Framework

• Model Space Construction: The many-body basis is constructed by superposing Skyrme-Hartree-Fock (SHF) wave functions along a fission path. These basis states, denoted as $|Q, E_\mu\rangle$, include:
  • $Q$: Deformation parameter (along the fission path in the $Q_{20}$-$Q_{30}$ space).
  • $E_\mu$: Particle-hole excitation energy relative to the reference state at deformation $Q$.
• Hamiltonian and Interactions:
  • The Hamiltonian includes the deformation potential energy and particle-hole excitation energies.
  • Residual Interactions: Three types are included to model dissipation and mixing:
    1. Monopole pairing interaction: To reproduce low-lying $0^+$ states.
    2. Random residual interaction: To induce chaotic mixing and dissipative shape evolution.
    3. Diabatic interaction: Derived from the Gaussian Overlap Approximation (GOA) to connect diabatically connected configurations.
• Boundary Conditions (Decay Widths):
  • Entrance Channel ($\gamma_{in}$): Photo-absorption width is determined empirically using the E1 gamma strength function and level density.
  • Exit Channel ($\gamma_{out}$): Photo-emission width calculated via the Brink-Axel hypothesis.
  • Fission Channel ($fis$): The fission width matrix is defined at the scission point ($Q_{max}$). Since empirical fission transmission coefficients are unavailable for states beyond the barrier, a fixed value ($\langle\Gamma_{fis}\rangle = 100$ keV) is adopted, as the fission probability is known to be insensitive to its exact magnitude in this regime.

B. Numerical Implementation

• Green Function Calculation: The retarded Green function $G(E) = (NE - H + i\Gamma/2)^{-1}$ is computed.
• Eigenchannel Decomposition: To analyze the transmission mechanism, the transmission coefficient $T_{ab}$ is decomposed using Singular Value Decomposition (SVD) into eigenchannels ($T_n = |t_n|^2$). This identifies the dominant pathways for wave propagation from the entrance to the exit.
• Computational Strategy: Due to the large dimension of the basis ($N_{dim} \approx 1.45 \times 10^5$), the author uses the shift-invert Arnoldi method to solve the non-Hermitian eigenvalue problem, avoiding the computational prohibitions of solving large linear equations on fine energy meshes repeatedly.

3. Key Contributions

1. Extension of NEGF to Photo-induced Fission: The paper successfully adapts the NEGF formalism (previously applied to neutron-induced fission in $^{235}\text{U}$) to the photon-induced case ($^{236}\text{U}$), handling the specific nature of the photo-absorption entrance channel.
2. Microscopic Sub-Barrier Analysis: It provides a fully microscopic calculation of fission cross-sections in the sub-barrier region ($5.0 \le E_\gamma \le 6.0$ MeV), explicitly including quantum tunneling and particle-hole excitations.
3. Eigenchannel Decomposition in Fission: The paper introduces the application of eigenchannel analysis to nuclear fission, demonstrating that the fission process is dominated by a single transmission channel.
4. Validation of Bohr-Wheeler Picture: The results provide microscopic support for the Bohr-Wheeler transition-state theory, showing that fission in the sub-barrier region proceeds predominantly through the lowest transition state (first eigenchannel).

4. Results

• Fission Cross Sections:
  • The calculated cross-sections for $^{236}\text{U}(\gamma, f)$ in the range $5.0 \le E_\gamma \le 6.0$ MeV reproduce the overall experimental trends.
  • Deep Sub-Barrier ($5.0 - 5.2$ MeV): Excellent agreement with experimental data (within a factor of 2).
  • Near Barrier ($5.3 - 5.7$ MeV): The calculation slightly overestimates the cross-section but captures the decreasing trend.
  • Above Barrier ($> 5.7$ MeV): Agreement is within a factor of 5.
• Photo-Absorption: The calculated photo-absorption cross-section matches empirical values within a factor of 1.2, validating the treatment of the entrance channel.
• Eigenchannel Analysis:
  • The first eigenchannel dominates the fission probability ($|t_1|^2 \gg |t_n|^2$ for $n \ge 2$) across the entire energy range.
  • The singular vectors of the first eigenchannel show a strong correlation with the collective wave functions of the resonant states near the incident energy, confirming that fission occurs via the lowest available transition state.
  • This dominance supports the Bohr-Wheeler transition-state picture, indicating that the number of effective degrees of freedom in the fission process is small (essentially one).

5. Significance and Future Outlook

• Theoretical Advancement: This work demonstrates that the NEGF method is a robust tool for calculating fission cross-sections, particularly in regimes where quantum tunneling is critical and initial state preparation is complex.
• Microscopic Insight: It bridges the gap between macroscopic transition-state models and microscopic many-body dynamics, showing that the "transition state" corresponds to the dominant eigenchannel of the Green function.
• Future Directions:
  • Improved Interactions: Moving beyond simplified random interactions to more realistic residual interactions is crucial, especially for neutron-rich nuclei relevant to the r-process.
  • Dimensionality: Extending the model space beyond the 1D fission path to include multi-dimensional deformation parameters (e.g., triaxiality) to better describe fission fragment mass distributions.
  • Application: The framework is expected to be valuable for predicting fission properties of exotic, neutron-rich nuclei where experimental data is scarce.

In summary, the paper presents a successful microscopic application of the NEGF method to sub-barrier photo-induced fission, validating the dominance of a single transition state and providing a reliable theoretical tool for future nuclear data evaluations.