Coupled Channels Formalism

• Coupled Channels Formalism is a quantum scattering framework that coherently mixes multiple reaction channels to ensure exact two-body unitarity.
• It utilizes integral equations like the Lippmann–Schwinger equation and matrix representations (T-matrix and S-matrix) to extract resonance properties and phase shifts.
• Employing dispersion relations, N/D parameterization, and analytic continuation, the formalism effectively models threshold effects and multi-Riemann sheet structures.

The coupled channels formalism provides a unified, rigorous framework for describing quantum scattering and reaction dynamics wherein multiple reaction or decay channels are coherently mixed by the underlying interactions. It underpins modern analyses of hadronic and nuclear spectra, resonance formation and decay, threshold effects, and multi-channel final-state interactions. Central to this approach is the requirement of exact two-body unitarity, encoded via a matrix S-matrix or T-matrix formalism, often supplemented by dispersion-theoretic or integral-equation representations to guarantee analyticity. Across hadronic, atomic, and nuclear physics, the formalism enables the extraction of resonance properties, cross sections, phase shifts, and observable amplitudes even in strongly coupled or nonperturbative regimes (Oller, 17 Jan 2025).

1. S-matrix, T-matrix, and Unitarity in Coupled Channels

In a system with $n$ open channels (e.g., two-body states of varying internal structure), the partial-wave T-matrix $T(s)$ is represented as an $n \times n$ matrix, with elements $T_{ij}(s)$ encoding the on-shell transition amplitude from channel $j$ to channel $i$ at invariant energy $s$. The phase-space matrix $\rho(s)$ is diagonal, with

$$\rho_i(s) = \frac{\sqrt{\lambda(s, m_{1,i}^2, m_{2,i}^2)}}{16\pi s}\,,$$

where $\lambda(x, y, z)$ is the Källén function. The coupled-channel S-matrix is then

$$S(s) = I + 2i\, \rho^{1/2}(s)\, T(s)\, \rho^{1/2}(s)\,,$$

ensuring $S S^\dagger = I$. The unitarity condition for the physical S-matrix above all open thresholds requires

$\Im T(s) = T(s)\, \rho(s)\, T(s)^*\,,$

which restricts the analytic structure and imbues the amplitudes with the correct cut discontinuities (Oller, 17 Jan 2025).

2. The Lippmann–Schwinger Equation and Integral Representations

The coupled-channel Lippmann–Schwinger (LS) equation resums all multiple scatterings, starting from a potential $V_{ij}(E)$ coupling initial and final states. In operator or matrix notation: $$T(E) = V(E) + V(E)\, \frac{1}{E - H_0 + i0^+}\, T(E),$$ where $H_0$ is the free Hamiltonian, or in explicit (partial wave/projected momentum) basis: $$T_{ij}(E) = V_{ij}(E) + \sum_{k=1}^n V_{ik}(E)\, G_k(E)\, T_{kj}(E)\,,$$ with the loop function $G_k(E)$ incorporating the appropriate free propagation and phase space in channel $k$. This framework realizes the analytic continuation to complex energies required for rigorous resonance analysis (Oller, 17 Jan 2025).

3. Analytic Structure: $N/D$ Parameterization and CDD Poles

The partial-wave amplitude can be decomposed as

$T(s) = D(s)^{-1} N(s)\,,$

with $D(s)$ encoding only right-hand cuts (unitarity) and $N(s)$ containing only left-hand (crossed channel) singularities. Along the unitarity cut ($s > s_{\rm th}$), one has

$\Im D(s) = -N(s) \rho(s)\,,$

while $\Im N(s) = 0$ in this domain. CDD poles are implemented in $D(s)$ as additional resonance phenomena not reducible to rescattering alone: $$D(s)=\ldots + \sum_{i=1}^{N_c} \frac{\gamma_i}{s-s_i}\,,$$ and introduce necessary analytic flexibility to match experimental data in the presence of, e.g., “elementary” resonances (Oller, 17 Jan 2025).

The general matrix parameterization

$T(s) = (N(s)^{-1} + G(s))^{-1}\,,$

with $N(s)$ polynomials and CDD poles, and $G(s)$ a matrix of right-hand-cut functions with $\Im G_{ii}(s) = \rho_i(s)$ for $s > s_{{\rm th},i}$, guarantees unitary and analytic coupled-channel amplitudes (Oller, 17 Jan 2025).

4. Resonances and Riemann Sheets

Each channel threshold creates a branch point, so for $n$ channels there are $2^n$ Riemann sheets, distinguished by the sign conventions for $\Im k_i$ across cuts. Resonances appear as poles on unphysical sheets closest to the physical region; their manifestation in observable line-shapes depends decisively on both sheet and proximity to relevant thresholds (Oller, 17 Jan 2025). The analytic continuation of the loop function enables exact placement of the resonance pole: $g_i^{(\text{II})}(s) = g_i(s) + 2i\,\rho_i(s)\,,$ where sheet II corresponds to $\Im k_1 < 0$ for the lightest channel.

5. Scattering Near Thresholds: Model Examples

The two-potential formalism separates a nonresonant background $V$ from an explicit $s$-channel “seed” pole: $$V_T(p,p';E) = V(p,p') + \frac{f(p) f(p')}{E - E_0}\,,$$ with the full solution

$$T = T_V + \frac{\Theta(E)\, \Theta(E)}{E - E_0 + \Sigma(E)}\,,$$

where the self-energy $\Sigma(E)$ and dressing $\Theta(E)$ account for loop-induced mixing and mass shifts. This decomposition allows for quantitative assessment of compositeness and structure of near-threshold states (Oller, 17 Jan 2025).

The Khuri–Treiman approach for three-body decays (e.g., $\eta \to 3\pi$) employs subenergy dispersion relations with unitarity-imposed discontinuities, realizing exact resummation of two-body rescattering in multi-channel environments (Oller, 17 Jan 2025).

6. Dispersion-Theoretic Final State Interaction Solutions

Omnès Solution

In single-channel final-state interaction problems with known phase $\varphi(s)$, the form factor is constructed as

$$F(s) = \frac{P(s)}{Q(s)} \exp\left\{ \frac{s}{\pi} \int_{s_{\rm th}}^\infty \! \frac{\varphi(s')}{s'(s'-s)}\,ds' \right\}\,,$$

where $P,Q$ capture subtractions/poles/zeros. This guarantees $F(s)$ has the correct unitarity phase along the right-hand cut (Watson’s theorem) (Oller, 17 Jan 2025).

Muskhelishvili–Omnès Matrix Solution

For general coupled channels, an $n\times n$ matrix $D(s)$ is sought such that

$S(s) = D(s)^{-1} D(s)^*\,.$

The determinant has the Omnès-type integral representation

$$\det D(s) = \frac{P(s)}{Q(s)} \exp\left\{ \frac{s - s_{{\rm th},1}}{2\pi} \int_{s_{{\rm th},1}}^\infty \! \frac{\Phi(s')}{(s' - s_{{\rm th},1})(s' - s)}\,ds' \right\},$$

with $\Phi(s)$ the sum of eigenphases—each column of $D^{-1}(s)$ yields one basis solution for any form factor saturating the correct unitarity (Oller, 17 Jan 2025).

7. Broader Applicability and Conceptual Implications

The coupled channels formalism—encompassing the $N/D$ method, explicit CDD pole handling, Lippmann–Schwinger integral equations, and rigorous dispersion techniques—forms the analytic bedrock for interpreting multi-channel quantum scattering. This ensures that theoretical amplitudes are not only unitary and analytic, but also include the full multi-sheet resonance structure observed in experimental spectroscopy. The formalism applies equally to multi-hadron and multi-nucleon systems, atomic collisions, and resonant event distributions in complex quantum systems, demonstrating its foundational character across quantum dynamics (Oller, 17 Jan 2025).