Reaction processes of muon-catalyzed fusion in the muonic molecule $ddμ$ studied with the tractable $T$-matrix model

This paper applies a tractable $T$-matrix model to investigate the nuclear reaction processes, fusion rates, and sticking probabilities of the $dd\mu$ muonic molecule, specifically analyzing $p$-wave fusion channels and charge symmetry violations in light of recently reported discrepancies in $p$-wave astrophysical $S(E)$ factors.

The Gist

The Big Picture: The Cosmic "Mousetrap"

Imagine you have a very heavy, stubborn door (the atomic nucleus) that you want to push open to release energy. Usually, you need a massive hydraulic ram (extreme heat and pressure, like in the Sun) to force it open.

Muon-catalyzed fusion is like having a magical, super-lightweight "mousetrap" that shrinks the door down to the size of a postage stamp, making it easy to push open.

In this paper, the authors are studying a specific version of this mousetrap using Deuterium (a heavy form of hydrogen) instead of the usual mix of Deuterium and Tritium. They are trying to figure out exactly how fast this trap snaps shut and what happens to the tiny particles flying out when it does.

The Cast of Characters

1. The Muon ($\mu$): Think of this as a "super-heavy electron." It's about 200 times heavier than a normal electron. Because it's heavy, it orbits the nucleus much closer, acting like a super-strong glue that pulls two nuclei together.
2. The Deuterium Nuclei (d): These are the "doors" we are trying to smash together.
3. The Muonic Molecule ($dd\mu$): This is the mousetrap itself. The muon grabs two deuterium nuclei and holds them so close that they fuse (merge) instantly.
4. The Fusion: When the nuclei merge, they release energy (heat) and a neutron. Sometimes, they turn into Helium-3 and a neutron; other times, they turn into Tritium and a proton.

The Problem: The "Five Different Maps"

The authors faced a tricky problem. To predict how fast this fusion happens, they needed to know the "shape" of the force between the two nuclei.

Imagine trying to drive from New York to London. You ask five different GPS apps for the route.

• App A says it's 3,000 miles.
• App B says it's 3,500 miles.
• App C says it's 2,900 miles.

In the world of physics, these "GPS apps" are five different groups of scientists who have measured the interaction between deuterium nuclei and come up with five different results (called S(E) factors). They disagree significantly on how the nuclei behave at low energies.

The authors asked: "If we use these five different maps, how does it change our prediction for the muon trap?"

The Solution: Three Different Ways to Solve the Puzzle

The authors didn't just pick one map. They used three different mathematical methods (like three different ways of calculating a route) to see if they could get a consistent answer despite the conflicting maps.

1. Method 1 (The Optical Model): Imagine looking at the nuclei through a foggy lens. You can't see the details, but you can see the overall "absorption" of energy. They used a "complex potential" (a mathematical fog) to estimate how fast the fusion happens.
2. Method 2 (The T-Matrix Model - Channel 5 & 8): This is like looking at the specific exit doors. They calculated exactly how the particles fly out into the "continuum" (free space) or get stuck in "bound states" (sticking to the new atom).
3. Method 3 (The T-Matrix Model - Channel 4 & 7): This is the same as Method 2, but looking at the exit from a slightly different angle to check the momentum of the particles.

The Good News: All three methods agreed with each other! Even though the "maps" (the S(E) factors) were different, the math held up.

The Key Findings

1. The Fusion Speed (How fast does the trap snap?)

Because the five "maps" were different, the predicted speed of the fusion varied.

• The Result: The fusion rate could be anywhere between 1.8 and 5.1 (in scientific units).
• The Takeaway: We still don't know the exact speed because the underlying data on how nuclei interact is still debated. We need better measurements of that "low-energy" interaction.

2. The "Sticky" Problem (The Muon's Dilemma)

After the nuclei fuse, they release a lot of energy. The muon is supposed to fly off and start a new fusion reaction (catalyzing the next one).

• The Issue: Sometimes, the muon gets "stuck" to the new Helium-3 atom created in the explosion. It's like a fly getting stuck in honey. If it gets stuck, it can't catalyze more reactions, and the energy chain stops.
• The Result: The authors calculated that the muon gets stuck about 13.3% of the time.
• Why it matters: This is a "showstopper" for using this as a power plant. If the muon gets stuck too often, you can't get enough energy out to pay for the cost of making the muons.

3. The "Charge Symmetry" Mystery

Nature has a rule called "Charge Symmetry," which basically says: "If you swap a proton for a neutron, the physics should look the same."

• In this reaction, the nuclei can split into Helium-3 + Neutron OR Tritium + Proton.
• If Charge Symmetry held perfectly, these two outcomes would happen at exactly the same rate.
• The Result: They don't! One happens about 1.4 times more often than the other. The authors confirmed this "violation" of symmetry using their new calculations, matching what experiments have seen.

4. The "Ultra-Slow" Muon Surprise (The Best Part!)

When the fusion happens, the muon is ejected. The authors calculated the speed and energy of this ejected muon.

• The Analogy: Imagine a cannonball being fired. Usually, you expect it to be moving fast.
• The Surprise: The muon comes out extremely slow.
  • Most of the muons have an energy of about 1 keV (very slow).
  • The average energy is higher (8.2 keV) because a few fly out very fast, but the peak (the most common speed) is very low.
• Why this is exciting: Scientists are trying to build "ultra-slow" muon beams for other experiments. This paper suggests that the dd-muon fusion is a natural factory for these slow muons! It's like finding a machine that naturally produces slow-moving particles without needing a complex brake system.

Summary for the General Audience

This paper is a deep dive into a specific type of nuclear fusion where a heavy electron (muon) acts as a catalyst to smash two hydrogen atoms together.

• The Conflict: Scientists disagree on the exact rules of how these atoms interact at low speeds.
• The Test: The authors used three different math tools to see how this disagreement affects the fusion process.
• The Verdict: The math works consistently, but the final speed of fusion depends on which "rulebook" you use.
• The Catch: The muon gets "stuck" to the new atom about 13% of the time, which is too high for a practical power plant.
• The Bonus: The process naturally spits out very slow muons, which could be a goldmine for future scientific experiments that need gentle, slow-moving particles.

In short: It's a rigorous check of the "rules of the road" for muon fusion, confirming that while it's a fascinating quantum dance, it's not quite ready to power our cities yet—but it might be perfect for building new kinds of scientific tools.

Technical Summary

1. Problem Statement

Muon-catalyzed fusion (µCF) in the deuterium-deuterium (ddµ) molecule has recently regained attention due to potential applications in subcritical thorium reactors and the generation of ultra-slow muon beams. However, a comprehensive understanding of the ddµ fusion process is hindered by significant discrepancies in existing data:

• Inconsistent Astrophysical S(E) Factors: Five different experimental and theoretical groups have reported significantly different p-wave astrophysical S(E) factors for the fundamental nuclear reactions $d+d \to {}^3\text{He} + n$ and $d+d \to t + p$ in the energy range of 1 keV to 1 MeV.
• Lack of Unified Models: Previous studies often relied on specific datasets or simplified approximations (like the sudden approximation for sticking probabilities) without systematically testing how these variations in nuclear input affect the full fusion cycle.
• Charge Symmetry Violation: There is ongoing debate regarding the violation of charge symmetry in p-wave d-d reactions, quantified by the ratio $R_S$ (ratio of S-factors) and $R_Y$ (ratio of fusion yields).

The primary goal of this work is to comprehensively study the nuclear reaction processes in the $(dd\mu)_{J=1}$ molecule using a unified theoretical framework that incorporates all five reported sets of p-wave S(E) factors to calculate fusion rates, sticking probabilities, and emitted particle spectra.

2. Methodology

The authors employ a tractable T-matrix model based on the Lippmann-Schwinger equation, previously validated for the dtµ system, and apply it to the ddµ system. The study proceeds through four distinct steps using three complementary calculation methods to ensure consistency:

• Input Preparation: The authors utilize nuclear interactions (optical potentials and coupling potentials) tuned to reproduce the five distinct sets of p-wave S(E) factors reported by Angulo+, Nebia+, Arai+, Tumino+, and Solovyev.
• Method i: Optical-Potential Model:
  • A complex d-d optical potential (real + imaginary parts) is constructed to reproduce the S(E) factors.
  • The three-body Schrödinger equation for the $(dd\mu)$ molecule is solved using the Gaussian Expansion Method (GEM).
  • The fusion rate is derived from the imaginary part of the complex eigenenergy of the molecular state.
• Method ii: T-matrix Model (Channels 5 & 8):
  • Explicitly calculates transition rates to outgoing channels where the muon sticks to ${}^3\text{He}$ or $t$ (bound states) or remains free (continuum states).
  • Uses a continuum-discretization method (CDCC) to handle the three-body continuum (${}^3\text{He} + n + \mu$ and $t + p + \mu$).
  • Calculates the sticking probability as the ratio of bound-state transitions to total transitions.
• Method iii: T-matrix Model (Channels 4 & 7):
  • Calculates the fusion rate by treating the muon as an emitted particle in the final state.
  • This method is specifically designed to derive the momentum and energy spectra of the emitted muons.
• Consistency Check: The results from all three methods are compared to validate the model's robustness across different potential sets.

3. Key Contributions

• Unified Framework for Discrepant Data: For the first time, the ddµ fusion process is analyzed using all five major sets of p-wave S(E) factors simultaneously, demonstrating how the choice of nuclear input affects macroscopic fusion parameters.
• Validation of the T-matrix Model: The study confirms that the tractable T-matrix model yields results consistent with the more computationally intensive optical-potential model and the coupled-channel approach, validating its use for complex few-body systems.
• First Calculation of Muon Spectra: The paper provides the first theoretical calculation of the absolute momentum and energy spectra of muons emitted from ddµ fusion, a critical step for future experimental detection of ultra-slow muons.
• Charge Symmetry Analysis: A rigorous investigation into the charge symmetry violation ratio ($R_Y$) across different energy regimes and S(E) datasets.

4. Key Results

A. Fusion Rates ($\lambda_{11}$)

• The calculated fusion rates for the $(dd\mu)_{J=v=1}$ state are consistent across the three calculation methods (Optical, T-matrix c=5/8, T-matrix c=4/7).
• However, the absolute values of the fusion rates vary significantly depending on the input S(E) factor set:
  • Range: $1.8 \times 10^8 \text{ s}^{-1}$ to $5.1 \times 10^8 \text{ s}^{-1}$.
  • Effective Rates: Corresponding effective fusion rates ($\tilde{\lambda}_f$) range from $2.0 \times 10^8 \text{ s}^{-1}$ to $5.3 \times 10^8 \text{ s}^{-1}$.
• This spread encompasses most experimental values but highlights the urgent need for more precise low-energy S(E) factor measurements.

B. Sticking Probability ($\omega_d$)

• The initial sticking probability of the muon to ${}^3\text{He}$ is calculated as $\omega_{11}^d = 0.133 \pm 0.001$.
• Crucially, this value is independent of the specific S(E) factor set used.
• The result agrees well with previous literature values based on the sudden approximation and experimental effective sticking probabilities (after accounting for reactivation).

C. Charge Symmetry Violation ($R_Y$)

• The ratio of fusion yields $R_Y = \lambda({}^3\text{He}+n) / \lambda(t+p)$ was calculated.
• Correlation: The study confirms the theoretical prediction that $R_Y \approx R_S$ (the ratio of S-factors) at low energies.
• Discrepancy:
  • Using S(E) factors from Nebia+, Arai+, and Solovyev, $R_Y \approx 1.3 - 1.5$, indicating significant charge symmetry violation (consistent with the latest experimental value of $1.455 \pm 0.011$).
  • Using S(E) factors from Angulo+ and Tumino+, $R_Y \approx 0.9 - 1.0$, suggesting no violation.
• This suggests that the interpretation of charge symmetry violation is highly sensitive to the underlying nuclear data.

D. Emitted Muon Spectra

• Peak Energy: The emitted muons have a peak energy of 1.0 keV, which is significantly lower than the average energy of 8.2 keV due to a long high-energy tail.
• Independence: The shape of the spectrum (peak position and width) is largely independent of the nuclear interaction parameters, making it a robust prediction for experimentalists.
• Comparison: The actual spectra are shifted to lower energies compared to the adiabatic approximation, indicating the muon is less spatially confined by the d-d system during the rapid fusion process.

5. Significance

• Experimental Guidance: The calculated muon energy spectrum (peaking at 1 keV) provides a critical target for future experiments aiming to generate ultra-slow negative muon beams for various applications.
• Nuclear Physics Insight: The study underscores the critical need for precise measurements of p-wave d-d S(E) factors below 100 keV. The current discrepancies in these factors lead to large uncertainties in predicted fusion rates and charge symmetry interpretations.
• Theoretical Benchmark: By successfully applying the T-matrix model to the ddµ system and reproducing results across three different methodologies, the paper establishes a reliable, computationally efficient tool for studying muon-catalyzed fusion in other molecular systems.
• Energy Applications: The findings support the feasibility of using ddµ fusion for thorium subcritical reactors, provided the sticking probability and fusion rates are optimized, while also highlighting the specific challenges posed by the 2.45 MeV neutron production.