$np \leftrightarrow d\gamma$ reactions calculated up to $E_{\gamma}=20$ MeV

This paper calculates the electromagnetic dipole transition cross sections for $np \leftrightarrow d\gamma$ reactions up to 20 MeV using high-order chiral effective field theory interactions and operators, validating results against existing experiments while providing new predictions via an adapted Efros method for future many-body applications.

The Gist

Imagine the atomic nucleus as a tiny, bustling dance floor where particles like protons and neutrons are constantly moving, colliding, and sometimes sticking together to form new pairs. This paper is a detailed report on a very specific dance: the moment a lone proton and a lone neutron meet, grab hands, and form a "deuteron" (a simple two-particle nucleus), all while flashing a burst of light (a photon) to celebrate. The reverse is also studied: what happens when a flash of light hits a deuteron and breaks the pair apart.

Here is a breakdown of what the researchers did, using everyday analogies:

1. The Goal: Mapping the Dance Floor

The scientists wanted to calculate exactly how likely these reactions are to happen across a huge range of energy levels—from the very slow, gentle movements found in the early universe (primordial nucleosynthesis) to much faster, more energetic collisions.

Think of this like trying to predict the outcome of a dance move. If you know the music (the energy) and the dancers' style (the forces between them), you can predict if they will stick together or spin apart. The researchers wanted to create a perfect "score" for this dance that matches what we see in real experiments.

2. The Tools: A New Way to See the Invisible

To do this, they needed a way to describe the "wave functions" of these particles. In quantum physics, particles aren't just solid balls; they are more like ripples in a pond. To calculate how these ripples behave when they crash into each other or break apart, you need a mathematical map.

• The Old Problem: Previous methods were like trying to map a whole ocean by measuring every single drop of water. It was accurate but computationally impossible for complex systems with more than a few particles. Other methods were like using a low-resolution camera; they could see the big picture but missed the fine details needed to calculate the "light flashes" (electromagnetic transitions).
• The New Tool (The Efros Method): The authors adapted a new technique (the "Efros method") that acts like a smart spotlight. Instead of trying to measure the entire ocean, this spotlight focuses only on the most important ripples (the "Short-Range Functions") that actually matter for the calculation. It allows them to get a clear, high-definition picture of the dance without needing to calculate every single drop of water.

3. The Rules of the Dance (The Interaction)

The dancers (protons and neutrons) follow specific rules of movement determined by the "Chiral Effective Field Theory" (χEFT). Think of this as the choreography manual.

• The researchers used a very advanced version of this manual (up to "N4LO"), which includes very subtle, high-level instructions on how the particles interact.
• They also used a specific manual for how the particles emit light (the "electromagnetic operators").

4. The Results: A Perfect Match

The team ran their calculations and compared their "predicted dance scores" against real-world data from experiments.

• The Good News: In most cases, their predictions matched the experimental data almost perfectly. It's as if they predicted exactly how many people would clap at a concert, and the actual crowd clapped at the exact same volume.
• The New Territory: They also calculated results for energy levels where no one had ever measured or predicted anything before. They filled in the blank spots on the map, providing a complete picture from very low energies up to 20 MeV.
• The Small Glitches: At a few very specific, extremely low-energy points, their numbers were slightly off (by a few percent) compared to some experiments. They explain this by saying their "choreography manual" might need a few more pages of instructions (higher-order corrections) to get those specific moves perfect.

5. Why This Matters (For This Paper)

The paper doesn't claim this will immediately cure diseases or build new engines. Instead, its main achievement is proving the new spotlight works.

By successfully using this "Efros spotlight" on a simple two-particle system (the proton and neutron), they have demonstrated that the method is ready to be used on much more complex nuclear systems in the future. It's like successfully testing a new drone on a small park before flying it over a city. They have shown that this new approach can handle the complex math of nuclear reactions accurately and efficiently, paving the way for understanding heavier, more complex atomic nuclei.

In summary: The authors built a new, efficient mathematical "spotlight" to watch how protons and neutrons stick together or break apart. They tested it, found it works beautifully against real-world data, and filled in missing pieces of the puzzle for energies we couldn't see before. This proves the tool is ready for bigger, more complex jobs in the future.

Technical Summary

1. Problem Statement

The paper addresses the calculation of electromagnetic dipole transition cross sections for the radiative capture reaction $np \to d\gamma$ and its inverse, the photodisintegration reaction $d\gamma \to np$.

• Context: These reactions are fundamental to primordial nucleosynthesis and nuclear astrophysics.
• Gap: Experimental data is sparse at low energies, and existing theoretical descriptions often rely on phenomenological potentials or older effective field theory (EFT) orders. Furthermore, there is a lack of modern theoretical calculations for energies where experimental data is missing (specifically in the range $E_\gamma \approx 15\text{--}19$ MeV).
• Challenge: Extending ab initio methods (specifically the No-Core Shell Model, NCSM) to describe continuum scattering states and reactions for systems with $A > 4$ is computationally prohibitive using traditional methods like the Resonating Group Method (RGM) or the standard HORSE method, which require calculating all eigenstates or are limited by convergence issues.

2. Methodology

The authors employ a modern ab initio framework combining high-order interactions with a novel adaptation of the Efros method for continuum calculations.

A. Theoretical Framework

• Interaction: The study uses the LENPIC nucleon-nucleon (NN) interaction derived from Chiral Effective Field Theory ($\chi$EFT) up to Next-to-Next-to-Next-to-Next-to-Leading Order (N4LO). This interaction is regularized using a semi-local coordinate space regulator ($R = 1$ fm).
• Operators:
  • Magnetic Dipole (M1): Calculated using the $\chi$EFT current operator up to N2LO.
  • Electric Dipole (E1): Calculated using the "naive" E1 operator. Higher-order corrections (NLO, N2LO) were deemed negligible based on previous literature.
• Basis: The calculations utilize the Oscillator Basis (Harmonic Oscillator) with a frequency $\hbar\Omega = 28$ MeV.

B. The "Efros Method" Adaptation

The core methodological innovation is the application of an adapted Efros method (based on the Hulthén–Kohn method) to the NCSM framework.

• Wave Function Parameterization: The scattering wave function is parameterized using a set of Short-Range Functions (SRFs) and an infinite series of oscillator functions to describe the long-range behavior.
• K-Matrix Representation: To ensure unitarity for coupled partial waves, the authors use a K-matrix representation ($S = \frac{1+iK}{1-iK}$) rather than the S-matrix directly.
• Truncation Strategy:
  • Unlike the standard HORSE method which requires the complete set of eigenfunctions (computationally expensive for many-body systems), this method uses a limited set of low-lying eigenfunctions (SRFs) of the truncated Hamiltonian.
  • The potential energy matrix is truncated at a maximum oscillator quanta ($N_{max}$), while the kinetic energy matrix is treated as infinite to allow for scattering.
  • Convergence: The authors found that a limited number of SRFs ($v$) and moderate $N_{max}$ values (e.g., $N_{max}=50$ for M1, $N_{max}=20$ for E1) are sufficient to achieve convergence, making the approach feasible for future many-body applications.

C. Bound State Calculation

• The deuteron bound state is obtained by locating the S-matrix pole in the coupled $^3SD_1$ partial waves. This approach ensures the correct asymptotic tail of the wave function, which is critical for low-energy capture cross-sections.

3. Key Contributions

1. Extension of Energy Range: The study provides theoretical cross-sections for $np \to d\gamma$ from astrophysical energies up to $E = 17.78$ MeV (corresponding to $E_\gamma = 20$ MeV).
2. Novel Data Points: The paper reports results for energy ranges where no experimental or modern theoretical data previously existed, particularly for photodisintegration cross-sections between 15 and 19 MeV.
3. Methodological Validation: It serves as a rigorous benchmark for the adapted Efros method. By comparing results against exact solutions (obtained via direct integration or the full HORSE method) for the two-body $np$ system, the authors demonstrate that their truncated SRF approach yields accurate results with significantly reduced computational cost.
4. High-Order Consistency: The results validate the use of N4LO $\chi$EFT interactions and N2LO operators for these reactions, confirming agreement with experimental data across a broad spectrum.

4. Results

• Agreement with Experiment: The calculated cross-sections agree well with existing experimental data for both $np \to d\gamma$ and $d\gamma \to np$ across most energy ranges.
  • Discrepancies: Minor discrepancies (few percent) were observed at very low energies ($E \approx 1.26 \times 10^{-8}$ MeV) and specific higher energies ($E_\gamma = 14.7, 16, 18, 19$ MeV). The authors attribute the low-energy discrepancy to the need for higher-order chiral currents (N3LO) for percent-level precision.
• Partial Wave Contributions:
  • M1 Transition: Dominated by the $^1S_0$ partial wave at energies $E \le 8$ MeV.
  • E1 Transition: Dominated by $^3P_0$, $^3P_1$, and $^3PF_2$ partial waves.
• Convergence:
  • Bound State: Using $N_{max}=120$ and $v=118$ SRFs yielded a binding energy of $E_B = -2.2232$ MeV (close to the experimental $-2.2246$ MeV).
  • Scattering: Reasonable convergence was achieved with $N_{max}=50$ for M1 and $N_{max}=20$ for E1, with uncertainties generally within 1-2%.
• Comparison with Theory: The results align with other modern $\chi$EFT calculations (e.g., Refs. [38]) and older potential models (Paris, Bonn), while resolving discrepancies found in pionless EFT calculations at higher energies.

5. Significance

• Astrophysical Impact: The provided cross-sections fill critical gaps in data required for modeling primordial nucleosynthesis and nuclear astrophysics, particularly in energy regions previously unexplored by theory.
• Pathway to Many-Body Systems: The primary significance lies in the validation of the Efros method adaptation. By demonstrating that accurate continuum wave functions and S-matrix poles can be obtained using a limited set of SRFs and moderate truncations, the authors pave the way for applying this method to many-body nuclear reactions ($A > 4$).
• Computational Efficiency: The study proves that the computationally expensive requirement of calculating all eigenstates (as in HORSE or NCGSM) can be bypassed without sacrificing accuracy, making ab initio reaction theory feasible for complex light nuclei.

In summary, this work successfully bridges the gap between high-precision chiral interactions and continuum reaction theory, providing a robust, computationally efficient framework for future ab initio studies of nuclear astrophysics and scattering phenomena.