Approximating the $S$ matrix for solving the Marchenko equation: the case of channels with different thresholds

This paper extends Marchenko inverse scattering theory to multi-channel systems with different thresholds by approximating the $S$-matrix with a rational and sinc series expansion, demonstrating that closed-channel submatrices can be reconstructed from open-channel data and validating the method through both synthetic potential tests and the analysis of $\pi N$ scattering data.

The Gist

Imagine you are a detective trying to figure out what a mysterious machine looks like inside, but you can't open it. All you can do is throw balls at it from the outside and watch how they bounce back. In physics, this is called scattering. The "balls" are particles (like pions and protons), and the "machine" is the force field or potential that holds them together.

The goal of this paper is to solve the Inverse Problem: taking the data from the bouncing balls (the scattering data) and working backward to reconstruct the shape and strength of the invisible machine (the potential) that caused the bounce.

Here is a breakdown of the paper's achievements using simple analogies:

1. The Old Way vs. The New Way

The Old Problem:
Previously, scientists tried to guess the machine's shape by assuming it looked like a specific type of curve (like a bell curve) and just tweaking the numbers until the bouncing balls matched the data.

• The Flaw: This is like trying to fit a square peg in a round hole. If the data is very precise, the math gets messy and creates "ghosts"—fake, impossible features in the machine that don't actually exist. To fix this, scientists had to blur the data, which made the final picture fuzzy.

The New Solution (The "Sinc" Trick):
The author, N. A. Khokhlov, developed a new mathematical "translator." Instead of forcing the data into a rigid shape, he uses a flexible tool made of two parts:

1. A Rational Term: A smooth, basic curve that gets the general shape right.
2. A Truncated Sinc Series: Think of this as a set of tiny, precise "pixels" or "stitches" that fill in the gaps.

By combining these, the author can recreate the exact path of the bouncing balls without creating any "ghosts" or fake features. It's like drawing a perfect portrait using a smooth outline and then adding thousands of tiny, precise dots to capture every detail, rather than guessing the whole face at once.

2. The "Threshold" Puzzle (The Doorway Problem)

The paper deals with a specific complication: Different Thresholds.
Imagine a building with two doors.

• Door A is low; you can walk through it easily (low energy).
• Door B is high; you need to jump to get through (high energy).

In physics, particles can only interact through "Door B" if they have enough energy. Below that energy, Door B is "closed."

• The Challenge: Scientists can easily measure what happens when particles go through Door A. But they can't directly measure what happens at Door B when it's closed.
• The Breakthrough: The author proves that even though Door B is closed, the way particles behave at Door A contains hidden clues about Door B. By using his new mathematical method, he can "extrapolate" (guess the shape of) the closed door based entirely on the data from the open door. He shows that the "closed" part of the machine is mathematically locked to the "open" part, so if you know one, you can reconstruct the other.

3. Relativistic Corrections (Speeding Up)

The paper also updates the math to handle particles moving at near-light speeds (relativistic effects).

• The Analogy: Imagine you are driving a car. At slow speeds, the rules of driving are simple. But as you approach the speed limit of the universe, the rules change (time slows down, mass increases).
• The Fix: The author rewrote the "driving manual" (the Marchenko equation) so it works correctly even when the particles are zooming fast. This ensures the reconstructed machine looks right even at high energies.

4. The Test Drive

To prove this works, the author did two things:

1. The Simulation: He created a fake machine with a known shape, generated fake scattering data, and then used his new method to try and rebuild the machine. Result: He successfully rebuilt the original shape perfectly.
2. The Real World: He applied this to real data from Pion-Nucleon scattering (a common interaction in nuclear physics). He managed to reconstruct the forces between these particles, successfully identifying a "resonance" (a specific vibration or state of the machine) that other methods struggled to see clearly.

Summary

In short, this paper provides a better, more flexible toolkit for physicists to reverse-engineer the invisible forces of nature.

• It removes "ghosts" from the data.
• It allows scientists to see "closed doors" by looking at "open doors."
• It works correctly even when particles are moving incredibly fast.

It's like upgrading from a blurry, guesswork-based map to a high-definition GPS that can navigate complex, multi-level terrain without getting lost.

Technical Summary

1. Problem Statement

The paper addresses the Inverse Scattering Problem (IP) in theoretical nuclear physics, specifically the challenge of reconstructing interparticle interaction potentials from scattering data.

• Context: Existing methods based on Marchenko, Gelfand-Levitan, and Krein theories provide unique local potentials but require complete knowledge of the S-matrix from zero to infinite energy.
• Specific Challenges:
  1. Coupled Channels with Different Thresholds: Most previous implementations of Marchenko theory were limited to cases with equal thresholds or single channels. Real-world scenarios (like $\pi N$ scattering) involve multiple channels opening at different energy thresholds ($W_1 < W_2 < \dots$).
  2. Relativistic Kinematics: Standard non-relativistic approximations fail at higher energies where relativistic effects become significant.
  3. Approximation Artifacts: Traditional methods approximate the S-matrix using rational-fraction expansions. While analytically convenient, these often introduce spurious (unphysical) poles near the physical momentum axis, degrading the accuracy of the reconstructed potential.
  4. Data Limitations: Experimental data is only available for open channels (energies above specific thresholds). The behavior of closed channels (below their thresholds) is inaccessible, yet required for a complete solution.

2. Methodology

The author proposes a generalized framework combining relativistic corrections with a novel numerical approximation technique for the S-matrix.

A. Relativistic Generalization of Marchenko Theory

• The work starts from a relativistic two-particle system where the wavefunction is an eigenfunction of the mass operator.
• By rearranging the mass operator equation, the author derives a formal equation identical to the Schrödinger equation:
  $$(\hat{q}^2 + \hat{V})\chi = Q^2\chi$$
  where $Q$ is a diagonal matrix of channel momenta derived from relativistic kinematics ($q_\alpha^2 = \frac{M^2}{4} - m_{\alpha}^2 + \dots$).
• This allows the application of the standard Marchenko integral equation formalism, provided the input kernel $F(x,y)$ is modified to account for the relativistic momentum-energy relation and the matrix structure of coupled channels.

B. Analytic Structure and Threshold Analysis

• The paper rigorously analyzes the analytic structure of the S-matrix in the inter-threshold region (where some channels are open and others are closed).
• Key Finding: The submatrix corresponding to closed channels ($S_-$) can be reconstructed from the experimentally accessible open-channel submatrix ($S_+$).
• The author derives explicit relations showing that the closed-channel elements depend on the open-channel parameters and specific coupling terms (e.g., $\gamma$), which behave predictably near thresholds (scaling as $q^{L+1/2}$).

C. S-Matrix Approximation Strategy

To avoid spurious poles, the author generalizes a method previously used for single channels [Ref. 21] to the multichannel case:

1. Rational Term ($S^{(0)}$): An initial approximation is constructed using rational functions (ratios of polynomials) to capture the global asymptotic behavior and unitarity above thresholds.
2. Truncated Sinc Series ($\Delta S$): To correct the rational approximation and fit high-precision data without introducing poles, a truncated series of sinc functions is added:
   $$S(q) \approx S^{(0)}(q) - \sum_{k=-N}^{N} s_k \lambda_k(q)$$
   where $\lambda_k(q)$ are sinc functions with finite support.
3. Separable Kernel: This specific approximation transforms the input kernel of the Marchenko equation into a separable (finite-rank) form. This allows the integral equation to be solved analytically via a system of linear equations, yielding potentials explicitly in terms of Riccati-Hankel functions.

3. Key Contributions

• Generalization to Unequal Thresholds: The first successful application of the Marchenko inversion method to coupled channels with distinct thresholds, including the derivation of the necessary relativistic corrections.
• Reconstruction of Closed Channels: A theoretical demonstration that the S-matrix for closed channels (below threshold) can be uniquely determined from open-channel data, resolving a major gap in inverse scattering for multi-channel systems.
• Pole-Free Approximation: The introduction of the rational-plus-sinc series method for multichannel S-matrices, which avoids the unphysical poles common in pure rational expansions.
• Algorithmic Implementation: A complete algorithm for converting experimental phase shifts and mixing parameters into a separable kernel for the Marchenko equation.

4. Results

The method was validated through two distinct tests:

A. Synthetic Data Test (Two-Channel Model)

• Setup: A known Gaussian potential ($V = V_0 e^{-ar^2}$) was used to generate "exact" S-matrix data for a two-channel system with different thresholds.
• Outcome:
  • The approximation method successfully reconstructed the S-matrix, showing significant improvement over the initial rational approximation, particularly near the threshold.
  • The inverse problem solution successfully recovered the original potential with high fidelity.
  • Minor oscillations were observed at large $r$ (an artifact of the finite range approximation), but the core potential structure was preserved.

B. Experimental Data Application ($\pi N$ Scattering)

• Setup: The method was applied to experimental Partial Wave Analysis (PWA) data for the $S_{31}$ $\pi N$ scattering state up to $E_{lab} \approx 2.5$ GeV.
• Model: The $\pi N$ channel was coupled to a $\pi^2 N$ channel (treated as a quasiparticle with mass $3m_\pi$).
• Outcome:
  • The reconstructed potentials were used to solve the direct scattering problem.
  • Cut-off Sensitivity: Potentials cut off at $r=8$ fm provided a good description of the data, including the resonance region. Potentials cut off at $r=4$ fm failed to reproduce the resonance.
  • The method successfully interpolated and extrapolated the S-matrix behavior across the inelastic threshold.

5. Significance

• Bridging Theory and Experiment: This work provides a robust, non-parametric tool to extract interaction potentials directly from scattering data without relying on phenomenological fitting, specifically for complex multi-channel systems.
• Relativistic Consistency: By incorporating relativistic kinematics directly into the Marchenko framework, the method is applicable to high-energy nuclear physics (e.g., nucleon-nucleon and pion-nucleon scattering) where non-relativistic approximations break down.
• Future Directions: The author identifies the extension of this formalism to three or more coupled channels and the inclusion of explicit three-body channels as the next critical step for describing complex nuclear interactions.

In conclusion, Khokhlov presents a mathematically rigorous and numerically stable extension of Marchenko theory that overcomes the limitations of previous methods regarding threshold differences and relativistic effects, offering a viable path for high-precision potential reconstruction in nuclear physics.