Scattering Faddeev calculations in the double continuum

This paper presents a configuration-space Faddeev formalism that unifies all scattering processes between single and double continua into a single matrix, which is demonstrated through its application to neutron-deuteron scattering.

The Gist

The Big Picture: The Three-Particle Dance

Imagine you are watching a dance floor with three dancers. In the world of quantum physics, these dancers are subatomic particles (like neutrons and protons).

Usually, physicists are good at predicting what happens when:

1. Two dancers hold hands (forming a bound pair, like a deuteron) and the third one bumps into them.
2. All three dancers are free and bump into each other, flying apart in different directions.

The problem arises when you try to calculate the math for the second scenario (all three flying apart) while they are still interacting. It's like trying to describe a chaotic mosh pit where everyone is moving, but you need to know exactly who is holding hands with whom at any given second.

This paper introduces a new, clever way to solve this "three-body problem" using a mathematical framework called Faddeev formalism.

---

The Core Problem: Two Different Languages

The author, Romain Guérout, points out a major headache in these calculations. The mathematical "wave" that describes the particles speaks two different languages depending on what the particles are doing:

1. Language A (The Couple): When two particles stick together and the third one flies away, the math works best using Cartesian coordinates (like a standard X and Y graph). It's like describing a car driving down a straight road.
2. Language B (The Trio): When all three particles fly apart, the math works best using Polar coordinates (like a radar screen with a center point and angles). It's like describing a firework exploding outward from a single point.

The Analogy: Imagine trying to write a story about a group of friends.

• When two friends are walking together and one runs off, you describe it as "Friend A and B are walking, Friend C is running North." (Linear/Cartesian).
• When all three run off in different directions, you describe it as "They are all running away from the center at different angles." (Radial/Polar).

The difficulty is that the computer simulation has to speak both languages simultaneously. If you try to force the "Polar" description into a "Cartesian" grid, the numbers get messy and inaccurate, especially when the particles are far apart.

The Solution: The "Resampling" Trick

The author's breakthrough is a simple but powerful trick: Resampling.

Think of the calculation as a high-resolution photo.

1. The computer first takes the photo using the Polar Grid (the radar screen). This is great for seeing the explosion of three particles.
2. However, to figure out if two particles are still holding hands (the "bound" state), the computer needs to look at that photo through a Cartesian Lens (the X/Y graph).

Instead of trying to force the math to work in one grid, the author's method takes the data from the Polar grid and translates it onto a Cartesian grid. It's like taking a photo of a spinning globe and flattening it onto a map so you can measure the distance between two specific cities accurately.

By doing this "translation," the computer can cleanly separate the "two-particle holding hands" part from the "three-particle flying apart" part without them getting confused.

The Benchmark: The Neutron and the Deuteron

To prove this method works, the author tested it on a famous "benchmark" system: Neutron-Deuteron scattering.

• The Setup: A neutron hits a deuteron (which is a proton and neutron already holding hands).
• The Outcomes:
  • Elastic: The neutron bounces off, and the deuteron stays together.
  • Breakup: The neutron hits hard enough to break the deuteron apart, sending all three particles flying.
  • Recombination: Three free particles crash together and form a deuteron again (very rare).

The author's new method produced results that matched the "gold standard" data from other super-computer calculations perfectly. This proves that the "Resampling" trick is accurate and reliable.

Why This Matters

1. One Matrix to Rule Them All: Before this, physicists often had to run separate calculations for "sticking together" and "flying apart." This method puts all possible outcomes (elastic, breakup, recombination) into one single matrix (a giant table of numbers). It's like having one master switch that controls the entire dance floor, rather than separate switches for every pair of dancers.
2. Checking the Math: The paper also shows how to check if the math is "honest." In physics, probability must be conserved (you can't lose or create energy out of thin air). The author created a new way to check this "conservation" even when mixing the two different coordinate languages, ensuring the results are physically real.

Summary in a Nutshell

Imagine you are trying to predict the outcome of a chaotic three-way collision.

• Old Way: You tried to describe the whole event using only one type of map, which made the details blurry and hard to read.
• New Way (This Paper): You use the best map for the explosion (Polar), then instantly translate the important parts onto a street map (Cartesian) to see exactly who is holding hands.
• Result: You get a crystal-clear picture of every possible outcome, from gentle bounces to total breakups, all calculated with high precision.

This work makes it easier and more accurate for physicists to understand how the fundamental building blocks of our universe interact when they are all moving freely.

Technical Summary

1. Problem Statement

The paper addresses the challenge of performing quantum mechanical scattering calculations for three particles in the double continuum regime, where all three particles are free ($1+1+1$ state).

• Context: While Faddeev formalism is well-established for bound states and single continuum scattering (one particle free, two bound, $1+2$ state), the double continuum presents unique difficulties.
• The Core Difficulty: In the double continuum, the calculated wavefunction contains a superposition of two distinct physical behaviors:
  1. Two-body channels: Where two particles form a bound cluster and the third is free.
  2. Three-body breakup channels: Where all three particles are free.
• These channels are not orthogonal in the standard sense, making it difficult to disentangle the scattering amplitudes for elastic scattering, breakup, and recombination from a single numerical solution. Furthermore, the asymptotic behavior of the wavefunction differs significantly between these regimes (bound plane waves vs. cylindrical waves), complicating the extraction of physical observables.

2. Methodology

The author employs the configuration-space Faddeev formalism using mass-scaled Jacobi coordinates. The methodology is divided into several key components:

A. Mathematical Framework

• Coordinates: The system uses 6D Jacobi vectors $(\mathbf{x}_i, \mathbf{y}_i)$. The total wavefunction $\Psi$ is decomposed into three Faddeev components $\phi_i$.
• Partial Wave Expansion: The angular part is expanded using bipolar spherical harmonics, reducing the problem to coupled integro-differential equations for radial wavefunctions $f_{i,a}(x_i, y_i)$.
• Jacobi Transform: A central mathematical tool is the integral transform $J_{ij}^{ab}$, which couples the different Faddeev components by transforming between different Jacobi coordinate arrangements.

B. Asymptotic Analysis

The paper rigorously defines the asymptotic behavior for the two regimes:

1. Single Continuum ($1+2$): As $y \to \infty$ (with $x$ bounded), the solution behaves as a "bound plane wave" (product of a bound state wavefunction and an outgoing Riccati-Hankel function).
2. Double Continuum ($1+1+1$): As both $x, y \to \infty$, the solution is best described in polar coordinates $(\rho, \alpha)$. The asymptotic form consists of "cylindrical waves" (Hankel functions of order $\nu$) modulated by Delves functions (eigenfunctions of the angular operator).

C. The Novel Extraction Method (Resampling)

The core innovation of this work is a method to extract both two-body and three-body contributions from the numerical solution:

• Grid Strategy: Calculations are performed on a polar grid $(\rho, \alpha)$, which is optimal for describing the cylindrical wave behavior of the breakup channel.
• Resampling: To extract the bound-state (two-body) contribution, the calculated polar wavefunction is resampled onto a Cartesian grid $(x, y)$.
• Fitting Procedure:
  • The wavefunction is projected onto the deuteron bound state to isolate the radial dependence.
  • A non-linear fit is performed to separate the coefficients of the regular/outgoing bound plane waves (elastic channel) and the regular/outgoing cylindrical waves (breakup channel).
  • This allows the simultaneous extraction of the elastic scattering matrix element ($S_{11}$) and the breakup amplitude ($S_{12}, S_{13}$) from a single calculation.

D. Numerical Implementation

• Basis: Cubic Hermite splines are used for radial expansion.
• Solver: The resulting large linear system (dimension $\approx 256,000$) is solved using the iterative GMRES method with a preconditioner based on the kinetic energy operator $(T-E)^{-1}$.
• Boundary Conditions:
  • Below breakup: Outgoing bound plane wave conditions.
  • Above breakup: Outgoing cylindrical wave conditions.

3. Key Contributions

1. Unified Scattering Matrix: The author constructs a single scattering matrix $S$ that encompasses all processes: elastic $1+2$, breakup $1+2 \to 1+1+1$, recombination $1+1+1 \to 1+2$, and elastic $1+1+1$.
2. Coordinate Transformation Technique: The development of a robust resampling technique (Polar $\to$ Cartesian) to disentangle non-orthogonal asymptotic channels.
3. Generalized Unitarity and Reciprocity: The paper defines a custom binary operation $\langle \cdot, \cdot \rangle$ and a matrix product $\ast$ that incorporates the Jacobi transform. This allows the verification of unitarity ($S \ast S^\dagger = I$) and reciprocity ($S = S^T$) for a matrix containing both scalar elements and continuous functions of the polar angle $\alpha$.
4. Benchmark Validation: The method is applied to the neutron-deuteron (n-d) scattering system, a standard benchmark in nuclear physics.

4. Results

The method was tested on the neutron-deuteron ($n-d$) system in the doublet state ($J^\Pi = 1/2^+$) using Yukawa potentials.

• Elastic Scattering ($S_{11}$): The calculated elastic scattering matrix elements show excellent agreement with benchmark data (both below and above breakup) across a wide energy range ($E_{lab} = 4$ to $42$ MeV).
• Breakup Amplitudes ($S_{12}, S_{13}$): The extracted breakup amplitudes for singlet and triplet channels match benchmark data from Friar et al. (1995) very satisfactorily.
  • Note: A slight deviation was observed near $\alpha \approx \pi/2$ due to the finite size of the computational grid, a known limitation in double continuum calculations.
• Three-Body Recombination & Elastic $nnp$: The paper presents the first calculated results for three-body recombination ($n+n+p \to n+d$) and elastic $n+n+p$ scattering, for which no experimental data exists.
• Accuracy: The defects from unitarity and reciprocity are on the order of $10^{-3}$ to $10^{-2}$, which is comparable to previous high-precision calculations, validating the numerical stability of the approach.

5. Significance

• Theoretical Advancement: This work resolves the long-standing difficulty of extracting multiple scattering amplitudes from a single wavefunction in the double continuum without requiring orthogonalization of channels (as done in previous methods).
• Computational Efficiency: By using polar grids for the breakup region and Cartesian resampling for the bound region, the method optimizes the description of both asymptotic behaviors efficiently.
• Benchmarking: The successful reproduction of established n-d benchmark data confirms the reliability of the configuration-space Faddeev approach for complex three-body breakup problems.
• Predictive Power: The ability to calculate recombination and pure three-body elastic scattering provides new theoretical data for nuclear reactions where experimental data is scarce or non-existent.

In conclusion, Guérout presents a robust, unified framework for three-body scattering that seamlessly handles the transition between bound-cluster and fully-free asymptotic states, offering a powerful tool for few-body nuclear and atomic physics.