Coupled-Channel Quark Cluster Model

• The Coupled-Channel Quark Cluster Model is a hadron-structure framework that treats observed hadrons as mixtures of compact quark cores and explicit hadronic channels.
• It employs coupling mechanisms such as ³P₀ pair creation, quark exchange, and dispersion relations to induce mass shifts, continuum admixtures, and threshold effects.
• The framework is applied to heavy-quarkonium, multiquark, and hexaquark systems, providing insights into resonance formation and decay patterns near thresholds.

The coupled-channel quark cluster model is a class of hadron-structure frameworks in which an observed hadron is treated not as a single isolated quark-model eigenstate but as a superposition of compact quark configurations and explicit hadronic channels. In the formulations surveyed here, the coupled basis may contain a bare (q\bar q) or (c\bar c) core plus open-flavor meson-meson continua, meson-meson and diquark-antidiquark cluster rearrangements, or six-quark subamplitudes associated with different pair clusterings. The common organizing principle is that channel coupling—implemented through (^{3}P_0) pair creation, string breaking, quark exchange, or dispersion-relation kernels—modifies masses, induces continuum admixtures, generates widths above threshold, and can produce threshold-dominated molecular or hybridized states rather than pure valence configurations [2109.02539][1412.7406][2208.14075][2306.03526][1003.0257].

1. Conceptual scope and channel content

Across the cited literature, the term denotes a family of related constructions rather than a single canonical Hamiltonian. Some realizations are “unquenched” quark models for heavy quarkonia, where a confined (c\bar c) state is dressed by open-charm loops; others are multichannel cluster calculations in which meson-meson, hidden-color, and diquark-antidiquark partitions are solved simultaneously; still others are relativistic six-quark coupled-subamplitude systems formulated through dispersion relations [2402.02765][2306.03526][1003.0257].

Realization	Channel basis	Characteristic coupling
Unquenched quarkonium model	bare (c\bar c) + open-charm continuum	(^{3}P_0) pair creation
Multichannel Schrödinger model	confinement (q\bar q) + meson-meson channels	delta-shell string breaking
Multiquark cluster calculation	(1!\otimes!1), (8!\otimes!8), (\bar 3!\otimes!3), (6!\otimes!\bar 6)	channel mixing in full Hamiltonian
Relativistic hexaquark formalism	pair-cluster subamplitudes (A_{ij})	coupled dispersion kernels

In the quarkonium versions, the physical state is explicitly decomposed into a compact core and hadronic continua. In the multiquark versions, the “cluster” language refers to alternative quark rearrangements with the same total quantum numbers. In the six-quark case, the cluster concept is promoted to a relativistic coupled-subamplitude hierarchy, formally analogous to a relativistic Faddeev–Yakubovsky construction [2109.02539][1412.7406][1003.0257].

A recurring physical theme is threshold sensitivity. When a bare state lies near an (S)-wave or (P)-wave–(S)-wave threshold, the continuum component can become numerically important, alter spectroscopy by tens of MeV, and control decay or production patterns. This is the central mechanism behind the heavy-quarkonium and multiquark applications discussed below [2109.02539][2402.02765][2604.23078].

2. Hilbert-space construction and coupling mechanisms

A standard unquenched formulation expands the physical hadron as
[
|A\rangle = c_0|\psi_0\rangle + \sum_{BC}\int d^3p\, c_{BC}(p)\, |BC;p\rangle ,
]
with (|\psi_0\rangle) a bare charmonium state from a quenched potential model and (|BC\rangle) charmed-meson continuum states. The Hamiltonian is written as
[
H = H_0 + H_{BC} + H_I,
]
where (H_0) is a Cornell-model charmonium Hamiltonian, (H_{BC}) is the free two-meson continuum, and (H_I) is a pair-creation interaction. In these implementations the interaction is commonly taken in the standard (^{3}P_0) form
[
H_I = 2m_q\gamma \int d^3x\, \bar\psi_q \psi_q ,
]
so that light (q\bar q) creation couples the compact core to open-flavor channels [2109.02539].

A closely related multichannel Schrödinger realization writes
[
|\psi\rangle = \sum_c |\psi_{q\bar q}^{\,c}\rangle + \sum_j |\psi_{MM}^{\,j}\rangle ,
]
and solves a coupled radial equation in which the confined (q\bar q) sector is governed by a harmonic-oscillator potential,
[
V_{Q\bar Q}(r)=\frac12 \mu_c \omega^2 r^2,
]
while the meson-meson sector is treated as free relative motion and the two are linked by a delta-shell transition potential,
[
V_{cj}(r)=\frac{\lambda\, g_{cj}}{2\mu_c}\,\delta(r-a).
]
In that model the same string-breaking radius (a) is used for (J/\psi), (\psi(2S)), and (X(3872)), with (m_c=1.562~\text{GeV}), (\omega=0.190~\text{GeV}), (\lambda_\psi=2.53), and (a=1.95~\text{GeV}^{-1}) after fitting the (J/\psi) and (\psi(2S)) masses [1412.7406].

In multiquark cluster calculations, the basis is enlarged further. For (T_{cc}), for example, the coupled-channel basis contains eight channels: six meson-meson channels and two diquark-antidiquark channels,
[
[DD^*]_{1\otimes 1},\ [D^*D]_{1\otimes 1},\ [D^D^]_{1\otimes 1},\ [DD^*]_{8\otimes 8},\ [D^*D]_{8\otimes 8},\ [D^D^]_{8\otimes 8},\ [(cc)(\bar q\bar q)^*]_{6\otimes \bar 6},\ [(cc)^*(\bar q\bar q)]{\bar 3\otimes 3}.
]
The coupled problem is then reduced to a generalized eigenvalue equation of the form
[
\sum_j \left( H{ij} - E\,N_{ij} \right) c_j = 0,
]
implemented with the Gaussian expansion method [2306.03526].

The six-quark dispersion-relation formalism adopts a different but related decomposition:
[
A = \sum_{i<j} A_{ij}, \qquad i,j=1,\dots,6,
]
so that the full six-body amplitude is represented by 15 pair subamplitudes. For the illustrative (uuuuuu) system, a representative pair amplitude is further decomposed into three subamplitudes (A_1), (A_2), and (A_3) depending on progressively larger cluster invariants. This cluster hierarchy is the dynamical backbone of the coupled-channel construction in the hexaquark sector [1003.0257].

3. Self-energies, effective kernels, and analytic structure

In unquenched quarkonium models, the coupled-channel effect enters spectroscopy through self-energy integrals. A typical mass shift is
[
\Delta M = \sum_{BC}\int d^3p\, \frac{\left|\langle BC;p|H_I|\psi_0\rangle\right|^2}{M-E_{BC}-i\epsilon},
]
and the continuum diagnostic
[
P_{BC} = \int d^3p\, \frac{\left|\langle BC;p|H_I|\psi_0\rangle\right|^2}{(M-E_{BC})^2}
]
is used to quantify how strongly a state couples to a given intermediate channel. In one implementation, because (\gamma) is not fixed from a full spectrum fit, these channel “probabilities” are explicitly described as unnormalized coupling strengths rather than absolute probabilities [2109.02539].

A systematic charmonium model uses the same (^{3}P_0) amplitudes to generate masses, widths, and core fractions simultaneously. For states below threshold the pole condition is
[
M-M_A-\Pi_{BC}(M)=0,
]
with
[
\Pi_{BC}(M) = \sum_{BC}\int_0^\infty \frac{1}{M-E_B-E_C} \sum_{JL}|\mathcal M_{JL}|^2\,P^2\,dP.
]
For states above threshold, (\Pi_{BC}) acquires an imaginary part; the real mass is extracted from the principal value of the integral, while the total open-charm width is
[
\Gamma_\text{total}=-2\,\operatorname{Im}\Pi_{BC}(M).
]
In that framework the approximate (c\bar c) fraction is
[
P_{c\bar c}= \left[ 1+\sum_{BC}\int_0^\infty \frac{1}{(M-E_B-E_C)^2} \sum_{JL}|\mathcal M_{JL}|^2\,P^2\,dP \right]^{-1},
]
though the paper states that for broad resonances it is not a strict probability [2402.02765].

The Feshbach viewpoint makes the effective interaction structure explicit. For a two-channel (\bar c c)-(\bar D D) system,
[
H = \begin{pmatrix} T^{\bar c c} & 0 \ 0 & T^{\bar D D} + \Delta \end{pmatrix} + \begin{pmatrix} V^{\rm conf} & V^t \ V^t & V^{\bar D D} \end{pmatrix},
]
eliminating one channel yields an effective Hamiltonian in the other. Eliminating the hadron channel gives
[
H^{\bar c c}{\rm eff}(E) = T^{\bar c c} +V^{\rm conf} +V^t \left(E-T^{\bar D D}-\Delta-V^{\bar D D}+i0^+\right)^{-1} V^t.
]
The resulting effective potential is therefore generically non-local and energy dependent, irrespective of the properties of the transition potential. When the eliminated channel contains scattering states, the effective interaction acquires an imaginary part above threshold, encoding decay into the open hadronic channel. A local derivative expansion,
[
\langle \mathbf r'{\bar c c}| V^{\bar c c}{\rm eff}(E) |\mathbf r{\bar c c}\rangle = V^{\rm conf}(\mathbf r)\delta(\mathbf r'-\mathbf r) + \frac{[V^t(\mathbf r)]^2}{E-\Delta}\delta(\mathbf r'-\mathbf r) +O(\nabla^2),
]
is possible formally, but the cited analysis stresses that such a truncation loses the imaginary part above threshold and can therefore distort the analytic structure [2208.14075].

The multichannel Schrödinger approach with a delta-shell transition potential exhibits the same logic in a different representation. There the channel wave functions are continuous at the string-breaking radius (r=a), while their derivatives have discontinuities fixed by the delta interaction; the matching conditions determine the secular equation for the eigenenergy and the normalized channel composition [1412.7406].

4. Heavy-quarkonium realizations

The heavy-quarkonium sector is the most developed arena for coupled-channel quark-cluster methods. In the charmoniumlike vector region, one coupled-channel quark model treats (Y(4260)) as predominantly a (D_1\bar D) molecule mixed with a short-distance charmonium core, and then uses heavy-quark spin symmetry (HQSS) to identify partner states. The relevant thresholds are (D_1\bar D) at (4292) MeV, (D_2^*\bar D) at (4331) MeV, (D_1\bar D^*) at (4432) MeV, and (D_2^*\bar D^*) at (4471) MeV. Within that analysis, the dominant coupled channel near the (Y(4360)) region is (D_1\bar D^*), while for (\psi(4415)) the strongest coupling is (D_2^*\bar D^*). The resulting HQSS multiplet picture is (Y(4260)) as mostly (D_1\bar D), (Y(4360)) as mostly (D_1\bar D^*), and (\psi(4415)) as mostly (D_2^*\bar D^*) [2109.02539].

A broader charmonium coupled-channel study uses the same (^{3}P_0) mechanism to analyze masses, open-charm widths, and (S)-(D) mixing. There the vector assignments are (\psi(3770)) as (1^3D_1)-dominated, (\psi(4040)) as (3^3S_1)-dominated, (\psi(4160)) as (2^3D_1)-dominated, (\psi(4360)) as (4^3S_1)-dominated, and (\psi(4415)) as (3^3D_1)-dominated. The extracted mixing angles are (-10.5^\circ) for (\psi(3770)), (-23.4^\circ) for (\psi(4040)), (-25.9^\circ) for (\psi(4160)), (-27.1^\circ) for (\psi(4360)), and (-24.6^\circ) for (\psi(4415)). This illustrates a central feature of the formalism: loop-induced mixing can be as important as the mass shifts themselves [2402.02765].

The (X(3872)) is a benchmark case for channel hybridization. In the multichannel Schrödinger model with HO confinement and a delta-shell transition potential, the fitted wave function contains (26.8\%) (c\bar c), (65.0\%) (D^0D^{*0}), (7.0\%) (D^\pm D^{*\mp}), and (1.2\%) (D^D^). In the same model the radiative widths are (\Gamma(X\to J/\psi\gamma)\approx 24.6~\text{keV}) and (\Gamma(X\to\psi(2S)\gamma)\approx 28.9~\text{keV}), giving (\mathcal R_\psi \approx 1.17) [1412.7406]. A different charmonium coupled-channel calculation instead treats (\chi_{c1}(3872)) as a (\chi_{c1}(2P)) state with significant continuum contribution ((\sim 57\%)), a physical mass (M_{\rm cou}=3870.1~\text{MeV}), and a dominant (D\bar D^*) self-energy contribution (\delta m_{D\bar D^*}=-57.5~\text{MeV}) [2402.02765]. This suggests that the coupled-channel framework accommodates strongly mixed core-plus-continuum interpretations while leaving the precise partition between “quarkonium” and “molecule” implementation dependent.

Coupled-channel effects are equally important in bottomonium transitions. For hindered M1 processes between (P)-wave bottomonia, virtual (B^{(*)}\bar B^{(*)}) loops can dominate the quenched spin-flip amplitude. The triangle-loop power counting gives
[
A_{\text{triangle}} \propto \frac{E_\gamma}{m_b\,v},
]
so the amplitude is enhanced near open-bottom thresholds. The same study emphasizes that two-loop pion-exchange contributions can be order one relative to the triangle contribution, which limits precision but reinforces the conclusion that the coupled-channel mechanism can overwhelm the direct quenched amplitude [1604.00770].

5. Multiquark and hadron-level cluster realizations

The multiquark sector extends the coupled-channel quark cluster concept beyond quarkonium dressing. In the hidden-local-symmetry extension of the chiral quark model, the (T_{cc}) calculation combines confinement, one-gluon exchange, (\sigma), (\pi), (\omega), and (\rho) exchange and solves an eight-channel coupled problem. In a single-channel calculation (T_{cc}) is not bound, but after channel coupling the binding becomes
[
E_B = M - E_{\rm th} = -4.9\ \mathrm{MeV},
]
relative to the threshold (E_{\rm th}(DD^*)=3903.1\ \mathrm{MeV}). The same work reports RMS distances
[
r_{cc}=1.56\ \mathrm{fm},\quad r_{\bar q c}=1.24\ \mathrm{fm},\quad r_{\bar q\bar q}=1.70\ \mathrm{fm},
]
consistent with a loose molecular state. For (T_{bb}), by contrast, full channel coupling gives (-88.2\ \mathrm{MeV}) relative to (E_{\rm th}(BB^*)=10549.4\ \mathrm{MeV}), with much smaller RMS distances, indicating a more compact configuration [2306.03526].

The (Z_c(3900)) channel provides a hadron-plus-quark realization of the same logic. The coupled basis contains
[
J/\psi\,\pi,\quad \eta_c\,\rho,\quad \psi(2S)\,\pi,\quad D\bar D^*,\quad D^*\bar D^*,
]
for (J^{PC}=1^{+-}), (I=1). In this model the diagonal meson-exchange potentials are weak, the off-diagonal one-meson-exchange transitions are also weak, but short-distance quark exchange generates strong off-diagonal transitions, particularly (D\bar D^*\leftrightarrow J/\psi\pi). The quark-exchange interaction is effectively only off-diagonal at leading order, and after a global reduction factor (\beta\sim 0.3) the coupled-channel amplitudes reproduce the qualitative HAL QCD pattern: weak diagonal (D\bar D^*), strong open-charm to hidden-charm mixing, cusp-like threshold structures when more channels are added, and no bound state or resonance pole in the main coupled-channel analysis [2604.23078].

The relativistic hexaquark program pushes the cluster idea to six-body spectroscopy. There the total amplitude is built from 15 pair subamplitudes, the quark-quark interaction is modeled by an effective four-fermion interaction with gluon quantum numbers, and two-body diquark amplitudes are written in (N/D)-type form
[
a_n(s_{ik}) = \frac{G_n^2(s_{ik})}{1 - B_n(s_{ik})}.
]
By extracting leading singularities and evaluating reduced amplitudes at the central Dalitz point
[
s_0 = \frac{s + 4\sum_i m_i^2}{6},
]
the integral system is reduced to algebraic equations with a common denominator (D(s)), and the hexaquark masses are obtained from the pole condition
[
D(s)=0.
]
Using (\Lambda=11), (m_{u,d}=410\ \text{MeV}), (m_s=557\ \text{MeV}), and two fitted gluon-coupling sets, that study reports 39 dibaryons, including a deuteron-like state near (1865\ \text{MeV}) and an (I=3), (J^P=0^+,2^+) (\Delta\Delta)-type state near (2379\ \text{MeV}) [1003.0257].

6. Observables, diagnostics, and methodological limits

A notable strength of the coupled-channel quark cluster framework is that it connects spectroscopy to channel-specific observables. In the HQSS analysis of the (Y) and (\psi) states, decay modes are proposed as discriminants of long- and short-distance structure. If (Y(4260)) is a (D_1\bar D) molecule, the dominant (D_1\to D^*\pi) decay implies a natural (D^*\bar D^*\pi) final state but not (D\bar D^*\pi). If (\psi(4415)) is dominantly (D_2^*\bar D^*), then both (D\bar D^*\pi) and (D^*\bar D^*\pi) should be populated, and a broad enhancement around (4.40) GeV in (e^+e^-\to D\bar D^*\pi) is noted as consistent with that assignment [2109.02539].

In the charmonium applications, the same transition amplitudes that produce continuum dressing also determine open-charm widths and leptonic observables. For vector states the loop-induced (S)-(D) admixture is essential because pure (D)-wave vectors would otherwise have too-small (\Gamma_{ee}). For (\chi_{c1}(3872)), a quasi-three-body treatment using the finite width of (D^{*0}) gives (\Gamma(\chi_{c1}(2P))\simeq 0.94~\text{MeV}), consistent with the observed width (1.19\pm0.21) MeV within that model’s assumptions [2402.02765]. For the multichannel (X(3872)) model, unquenching strongly suppresses the quenched prediction (\Gamma_e(X\to\psi(2S)\gamma)=158~\text{keV}) and yields the more balanced radiative pattern quoted above [1412.7406].

The formalism, however, has persistent methodological limits. One coupled-channel quark model states explicitly that its (P_{BC}) values are unnormalized coupling strengths because the (^{3}P_0) strength (\gamma) is not fitted globally [2109.02539]. Another notes that (P_{c\bar c}) is only an approximate measure for broad resonances [2402.02765]. The (X(3872)) multichannel Schrödinger calculation assumes no direct meson-meson interaction in the (MM) sector [1412.7406]. The Feshbach analysis shows that eliminating channels produces non-local and energy-dependent kernels and that a naive local approximation can remove the imaginary part above threshold, thereby misrepresenting decay physics [2208.14075]. In bottomonium hindered M1 transitions, absolute widths depend on unknown couplings (g_1), (g_1'), and (g_1''), and some two-loop effects are estimated to be order one, so ratios of widths are identified as more robust discriminants than absolute partial widths [1604.00770].

These limitations do not weaken the central conclusion shared by the cited work. Rather, they define the technical frontier of the field: coupled-channel quark cluster models are most informative when the basis captures the relevant nearby thresholds, the transition operators preserve the correct analytic structure, and spectroscopy is constrained simultaneously by masses, widths, and channel-specific decay or scattering observables [2208.14075][2402.02765][2604.23078].