Semi Three-Body Unitary Coupled-Channel Model

• The semi three-body unitary coupled-channel model is a framework that reorganizes three-particle systems into coupled spectator–isobar channels to enforce three-body unitarity through nonperturbative rescattering.
• It utilizes approximations such as the isobar treatment, on-shell spectator assumption, finite channel truncation, and simplified contact terms, making it adaptable to diverse hadronic phenomena.
• Analytic continuation techniques within the model reveal resonance pole structures, enabling precise extraction of hadron properties beyond traditional Breit–Wigner fits.

A semi three-body unitary coupled-channel model is a class of hadronic amplitude frameworks in which a physical three-particle system is reorganized into coupled isobar–spectator channels and treated through nonperturbative rescattering equations, while the dynamical space is truncated and some short-distance or unstable-subsystem effects are supplied phenomenologically rather than derived from a full microscopic three-body field theory. In the recent literature, the qualifier “semi” commonly denotes the simultaneous use of the isobar approximation, an on-shell spectator treatment, a finite set of channels and partial waves, and simple parametrizations of short-range production or contact terms, combined with a resummation that preserves the relevant three-body unitarity structure of the reduced channel space (Sakthivasan et al., 23 Jun 2026). Such models have been used to analyze threshold-sensitive and process-dependent structures including the $a_1(1420)$, $\eta(1405/1475)$, $f_1(1285)$, $f_1(1420)$, $T_{cc}^+$, and vector charmonia, especially where Breit–Wigner or conventional isobar descriptions are inadequate (Nakamura et al., 2022).

1. Genealogy and conceptual motivation

The immediate precursor of the semi three-body unitary coupled-channel model is the unitary coupled-channels treatment of heavy-meson decays into three mesons developed to go beyond the conventional isobar model. In that earlier formulation, a heavy meson $M^*$ decays through channels of the type $M^*\to cR\to abc$, where $R$ is a dynamical two-body isobar and $c$ is the spectator. The essential new ingredient is the $Z$-diagram, in which a decay product of $\eta(1405/1475)$0 rescatters with the spectator to form a different $\eta(1405/1475)$1 state; this mechanism generates repeated $\eta(1405/1475)$2 rescattering and enforces the three-body unitarity missing in the conventional isobar approximation (Kamano et al., 2011).

The importance of this extension was quantified in pseudo-data studies of $\eta(1405/1475)$3. There, a unitary coupled-channels model and an isobar model could both reproduce Dalitz-plot distributions reasonably well, yet the extracted masses, widths, and coupling strengths of excited mesons differed significantly, especially for broad states. The same work emphasized that three-body unitarity affects not only line shapes but also the dynamical interpretation of extracted hadron properties (Nakamura, 2012).

A complementary line of development came from dispersion theory. A coupled-channel generalization of the Khuri–Treiman/isobar-type three-body final-state-interaction problem was formulated for spinless particles by writing the decay amplitude as a sum of pairwise partial waves and imposing coupled-channel two-body unitarity in each subchannel. That formalism introduced three equivalent representations, $\eta(1405/1475)$4, $\eta(1405/1475)$5, and $\eta(1405/1475)$6, and treated three-body rescattering through dispersion relations and self-consistent coupled integral equations. This suggests a broader methodological family in which “semi three-body” dynamics is generated from coupled two-body subsystems rather than a fully irreducible three-body Hamiltonian (Guo, 2015).

2. Spectator–isobar construction and dynamical equations

The standard semi three-body construction rewrites the full three-body amplitude in a spectator–isobar basis. In the $\eta(1405/1475)$7 analysis of the $\eta(1405/1475)$8 and $\eta(1405/1475)$9, for example, the full amplitude $f_1(1285)$0 is expressed through dissociation vertices and a reduced spectator–isobar amplitude $f_1(1285)$1, so that the nontrivial dynamics is concentrated in the latter. The reduced amplitude satisfies a Bethe–Salpeter- or Lippmann–Schwinger-type equation,

$f_1(1285)$2

where $f_1(1285)$3 is the one-particle-exchange interaction fixed by three-body unitarity, $f_1(1285)$4 is the dressed isobar propagator encoding the two-body subsystem, and $f_1(1285)$5 is a short-range three-body interaction not fixed by unitarity (Hu et al., 17 Jun 2026).

In the $f_1(1285)$6 analysis the same logic is implemented in an infinite-volume three-body formalism reorganized into isobar + spectator channels. The three-body amplitude is obtained by solving a Lippmann–Schwinger-like integral equation for transitions among these channels. Here the exchange kernel is denoted $f_1(1285)$7, the short-range term $f_1(1285)$8, the isobar propagator $f_1(1285)$9, and the production vertex $f_1(1420)$0. The framework is explicitly described as semi three-body because it uses the isobar approximation, an on-shell spectator treatment, truncation to a finite set of channels and partial waves, and a limited parametrization of short-distance production/contact terms (Sakthivasan et al., 23 Jun 2026).

The same architecture appears in radiative $f_1(1420)$1 decays. The full three-meson amplitude is written as a cyclic sum over spectator assignments and contains a dressed isobar propagator $f_1(1420)$2, dressed decay vertices $f_1(1420)$3, and a dressed propagator matrix for bare resonance states,

$f_1(1420)$4

In this formulation, bare states, quasi-two-body channels $f_1(1420)$5, and rescattering kernels are treated within a single coupled-channel system, rather than as separate ingredients glued together at the amplitude level (Nakamura et al., 2023).

Channel truncation is a defining practical step. The $f_1(1420)$6 study uses a nine-channel $f_1(1420)$7-basis production amplitude,

$f_1(1420)$8

including both nonstrange and strange sectors and repulsive channels such as $f_1(1420)$9 with $T_{cc}^+$0 and $T_{cc}^+$1 with $T_{cc}^+$2 (Sakthivasan et al., 23 Jun 2026). In other applications, the working space is instead built from channels such as $T_{cc}^+$3, $T_{cc}^+$4, $T_{cc}^+$5, $T_{cc}^+$6, $T_{cc}^+$7, $T_{cc}^+$8, and $T_{cc}^+$9 for $M^*$0 (Nakamura et al., 2022), or from $M^*$1 and $M^*$2 partial waves for the $M^*$3 sector (Hu et al., 17 Jun 2026).

3. Unitarity, rescattering, and singularity structure

The central claim of the semi three-body approach is that triangle diagrams and threshold effects should not be inserted as isolated one-loop corrections but embedded into a unitary rescattering series. In the $M^*$4 case, the observed enhancement near $M^*$5 GeV in the $M^*$6 wave is associated with the classic triangle process

$M^*$7

followed by $M^*$8. The singularity appears when the Landau equations admit a real solution, with the Coleman–Norton picture corresponding to the classical alignment in which the kaon from the $M^*$9 decay catches up with the spectator $M^*\to cR\to abc$0. In the unitary formulation, however, this one-loop triangle is only the leading approximation; final-state interactions generated by the coupled-channel integral equation can modify its shape and normalization (Sakthivasan et al., 23 Jun 2026).

The same principle governs the $M^*\to cR\to abc$1 system. There, the large isospin violation in $M^*\to cR\to abc$2 is tied to a triangle singularity in the $M^*\to cR\to abc$3 loop. Because $M^*\to cR\to abc$4, the cancellation between charged and neutral loops becomes incomplete in the narrow interval

$M^*\to cR\to abc$5

producing a sizeable isospin violation and a narrow $M^*\to cR\to abc$6-like structure in the $M^*\to cR\to abc$7 spectrum. In the three-body unitary framework, this effect is generated automatically by the analytic structure of the rescattering amplitude rather than being added by hand (Nakamura et al., 2022).

In the $M^*\to cR\to abc$8–$M^*\to cR\to abc$9 problem, the one-particle-exchange kernel $R$0 is explicitly identified as the mechanism that both enforces the correct three-body unitary structure and produces triangle singularities automatically. The relevant process is

$R$1

so the triangle singularity is an intrinsic part of the coupled $R$2–$R$3 amplitude rather than an optional correction (Hu et al., 17 Jun 2026).

A different manifestation appears for $R$4. There, the three-body content is generated by the instability of the $R$5 and the fact that the exchanged pion can go on shell. The framework is described as self-consistent only if both manifestations of the three-body dynamics, the pion exchange between the $R$6 and $R$7 mesons and the finite $R$8 width, are taken into account simultaneously to ensure that three-body unitarity is preserved. This establishes that semi three-body unitarity is not specific to mesonic isobar problems but also applies to near-threshold hadronic molecules with unstable constituents (Du et al., 2021).

4. Analytic continuation, Riemann sheets, and pole content

A principal advantage of the semi three-body unitary coupled-channel model is that resonances are extracted from poles of analytically continued amplitudes rather than from Breit–Wigner fits. In the $R$9 analyses, poles are found from

$c$0

after analytic continuation of the dressed propagator matrix into the complex plane by contour deformation. The sheet assignment is specified relative to nearby thresholds, notably the $c$1 and $c$2 channels. Two poles are found for $c$3, on different Riemann sheets of the $c$4 channel, and one for $c$5; the two $c$6 poles are explicitly interpreted as a split pole pair on different sheets caused by proximity to the $c$7 threshold, not as two distinct physical states (Nakamura et al., 2022).

The $c$8–$c$9 study emphasizes contour deformation and the handling of moving three-body cuts in the complex spectator-momentum plane. On the relevant unphysical sheets, it finds two robust poles,

$Z$0

$Z$1

together with an additional deeper pole on the same sheet as the $Z$2. Pole trajectories obtained by multiplying the rescattering series by a control parameter $Z$3 show that the $Z$4 evolves from a bare-state pole, whereas the $Z$5 moves toward the physical region as rescattering is turned on, indicating a predominantly dynamically generated origin (Hu et al., 17 Jun 2026).

The $Z$6 analysis illustrates the converse situation: the narrow structure near $Z$7 GeV is reproduced without introducing an additional genuine $Z$8 pole. Instead, the triangle-singularity-driven coupled-channel amplitude suffices, while a pole associated with the nearby ground-state axial vector resonance is extracted. The quoted $Z$9 pole positions are around $\eta(1405/1475)$00 MeV for fixed $\eta(1405/1475)$01, and $\eta(1405/1475)$02 MeV for the integrated $\eta(1405/1475)$03 case, with uncertainties dominated by missing systematic errors (Sakthivasan et al., 23 Jun 2026).

Pole extraction is not limited to light-meson spectroscopy. In the vector-charmonium analysis of $\eta(1405/1475)$04, analytic continuation of a semi three-body unitary coupled-channel amplitude fitted to 20 final states over $\eta(1405/1475)$05–$\eta(1405/1475)$06 GeV yields 14 vector-charmonium poles. The model distinguishes resonances on unphysical sheets, bound states on physical sheets of the nearest-threshold channel, and virtual states on the corresponding unphysical sheets, showing how threshold cusps and near-threshold poles coexist within a single amplitude framework (Nakamura, 16 Sep 2025).

5. Representative physical realizations

$\eta(1405/1475)$07: The nine-channel production amplitude fitted to COMPASS freed-isobar intensities in the $\eta(1405/1475)$08 channel reproduces the narrow enhancement near $\eta(1405/1475)$09 GeV in the $\eta(1405/1475)$10 wave. The model’s key conclusion is that the enhancement is generated by the $\eta(1405/1475)$11 loop dressed by three-body final-state interactions and by the strong $\eta(1405/1475)$12 dynamics of the $\eta(1405/1475)$13, so that an additional genuine $\eta(1405/1475)$14 pole is not required (Sakthivasan et al., 23 Jun 2026).

$\eta(1405/1475)$15: In radiative $\eta(1405/1475)$16 decay, a manifestly three-body unitary coupled-channel model simultaneously fits BESIII Monte-Carlo outputs for the $\eta(1405/1475)$17 component of $\eta(1405/1475)$18 and branching-ratio constraints involving $\eta(1405/1475)$19 and $\eta(1405/1475)$20 final states. Its main dynamical result is a three-pole structure: two poles associated with $\eta(1405/1475)$21 on different sheets of the $\eta(1405/1475)$22 channel and one pole associated with $\eta(1405/1475)$23. The model also makes process-dependent predictions for $\eta(1405/1475)$24, $\eta(1405/1475)$25, $\eta(1405/1475)$26, and $\eta(1405/1475)$27 line shapes, with the $\eta(1405/1475)$28 channel driven by the triangle singularity (Nakamura et al., 2023).

$\eta(1405/1475)$29 and $\eta(1405/1475)$30: In the $\eta(1405/1475)$31 $\eta(1405/1475)$32 system, the coupled $\eta(1405/1475)$33–$\eta(1405/1475)$34 amplitude is constructed in the spectator-isobar representation, and the short-range three-body interaction is constrained by the BESIII $\eta(1405/1475)$35 invariant-mass distribution. The pole analysis supports a dressed-bare-state interpretation for the $\eta(1405/1475)$36 and a predominantly dynamically generated $\eta(1405/1475)$37-wave $\eta(1405/1475)$38 quasi-bound-state interpretation for the $\eta(1405/1475)$39, while the additional deeper pole is traced to the $\eta(1405/1475)$40-wave $\eta(1405/1475)$41 contact interaction and leaves little visible imprint on the physical line shape (Hu et al., 17 Jun 2026).

$\eta(1405/1475)$42: The charged tetraquark is analyzed in a coupled-channel EFT where the $\eta(1405/1475)$43 channels and the explicit $\eta(1405/1475)$44 channels are treated together. Three schemes, contact only, dynamic $\eta(1405/1475)$45 widths without one-pion exchange, and the full contact + OPE + dynamic-width treatment, all fit the $\eta(1405/1475)$46 line shape with $\eta(1405/1475)$47 around $\eta(1405/1475)$48–$\eta(1405/1475)$49, but the pole’s imaginary part is highly scheme dependent. The compositeness parameter is found to be close to unity, implying that the $\eta(1405/1475)$50 is a hadronic molecule generated by the interactions in the $\eta(1405/1475)$51 and $\eta(1405/1475)$52 channels (Du et al., 2021).

Vector charmonia: A semi three-body unitary coupled-channel analysis of BESIII and Belle data for 20 final states identifies several established and near-threshold vector-charmonium poles. The quoted fit quality is $\eta(1405/1475)$53. In the compositeness analysis, $\eta(1405/1475)$54 is assigned a large $\eta(1405/1475)$55 component,

$\eta(1405/1475)$56

while $\eta(1405/1475)$57 and $\eta(1405/1475)$58 are interpreted as substantial mixtures of $\eta(1405/1475)$59, $\eta(1405/1475)$60, $\eta(1405/1475)$61, and $\eta(1405/1475)$62 components (Nakamura, 16 Sep 2025).

Methodological demonstration in $\eta(1405/1475)$63: A relativistic three-body unitary treatment of the $\eta(1405/1475)$64 final-state interaction, restricted to the dominant $\eta(1405/1475)$65 isobar channel and the coupled $\eta(1405/1475)$66 and $\eta(1405/1475)$67 partial waves, was fitted to the ALEPH lineshape. The resulting Dalitz plots showed that rescattering changes the distribution by roughly $\eta(1405/1475)$68, demonstrating the numerical feasibility of solving the manifestly unitary three-body equations with isobars of nonzero spin and coupled partial waves (Sadasivan et al., 2020).

6. Approximations, comparisons, and recurring misconceptions

The word “semi” does not denote a single approximation scheme but a family resemblance. In one usage, the model is semi three-body because it treats three-body rescattering nonperturbatively but within a reduced channel space and with two-body amplitudes supplied as input rather than derived from a full microscopic theory (Sakthivasan et al., 23 Jun 2026). In another, it refers to a quasi-three-body treatment in which the observed final state is genuinely three-body but the dynamics are organized through intermediate $\eta(1405/1475)$69 channels, with some unstable subsystems represented by Breit–Wigner forms; this leads to partial violation of three-body unitarity even though the coupled $\eta(1405/1475)$70 rescattering remains nonperturbative (Nakamura, 16 Sep 2025).

A common misconception is that a narrow enhancement near a threshold must correspond to a new resonance pole. The $\eta(1405/1475)$71 case is a counterexample: within a unitary coupled-channel three-body amplitude, the triangle singularity mechanism can reproduce the observed enhancement without an additional genuine $\eta(1405/1475)$72 pole (Sakthivasan et al., 23 Jun 2026). The converse misconception is that a triangle singularity excludes genuine resonant dynamics. The $\eta(1405/1475)$73 analysis shows that triangle-singularity physics and an actual pole can coexist in the same fully resummed amplitude, with the latter emerging as a predominantly dynamically generated $\eta(1405/1475)$74 state (Hu et al., 17 Jun 2026).

Another recurrent issue concerns multipeak or split-pole patterns. In the $\eta(1405/1475)$75 sector, two poles associated with $\eta(1405/1475)$76 do not imply two distinct physical states; the quoted interpretation is a threshold-induced split pole pair on different Riemann sheets of the $\eta(1405/1475)$77 channel (Nakamura et al., 2022). This illustrates a general lesson of semi three-body coupled-channel spectroscopy: pole multiplicity, process dependence, and threshold cusps are controlled by the analytic structure of the amplitude, not by naive peak counting.

Finally, comparison with the conventional isobar model remains central. Earlier unitary coupled-channels studies showed that the isobar approximation may reproduce Dalitz-plot shapes while still distorting extracted resonance parameters and phases, especially for broad states and in channels relevant to CP-violation analyses. The defining missing ingredient in the isobar model is the rescattering series generated by $\eta(1405/1475)$78-diagrams or one-particle-exchange kernels, which carries the three-body unitarity cut and dresses the production vertices (Kamano et al., 2011). A plausible implication is that semi three-body unitary coupled-channel methods are best understood not as replacements for all phenomenological models, but as the minimal framework required when threshold branch points, unstable constituents, and coupled-channel rescattering materially affect the physical observables.