---
title: Unquenched Quark Model Overview
url: https://www.emergentmind.com/topics/unquenched-quark-model
type: topic
---

# Unquenched Quark Model Overview

Searching arXiv for the topic and core review papers.
I’m going to retrieve relevant arXiv entries on the Unquenched Quark Model and related coupled-channel treatments.
The unquenched quark model (UQM) is an extension of the constituent quark model in which a hadron is treated not as a pure valence configuration—\(q\bar q\) for mesons or \(qqq\) for baryons—but as a dressed state containing explicit higher-Fock components generated by quark–antiquark pair creation and coupling to hadronic continuum channels. In the formulation emphasized in the review literature, the UQM retains the valence quark-model core as the organizing principle for spectroscopy while incorporating meson–meson or baryon–meson continuum mixing, self-energy corrections, threshold effects, and sea-quark contributions to observables [1503.02968][1506.05015][1503.03064].

## 1. Definition and conceptual scope

In the quenched constituent quark model, hadrons are approximated by fixed valence configurations. The UQM “unquenches” this picture by allowing the valence core to fluctuate into intermediate hadronic states through vacuum \(q\bar q\) creation. For mesons this means \(q\bar q \to (q\bar q)(q\bar q)\); for baryons it means \(qqq \to qqq\,q\bar q\), usually organized as baryon–meson components. Above threshold, the same coupling produces strong decays; below threshold, it generates virtual continuum admixtures and mass renormalization [1506.05015][1712.03919].

This framework was presented by Ferretti, Galatà, and Santopinto as a controlled extension of the conventional quark model for mesons, with particular emphasis on charmonium and bottomonium self-energies and on the possibility of extending the formalism to hybrid mesons [1503.02968]. Closely related review treatments describe the UQM as a unified method for hadron structure and spectroscopy, applicable to proton flavor asymmetry and strangeness observables in the baryon sector as well as to continuum-induced mass shifts in heavy quarkonia [1511.05316][1503.03064].

A recurrent theme in the literature is that unquenching does not imply abandoning the constituent-quark description. Rather, the physical hadron is a normalized superposition of a valence seed and continuum components. This suggests a hierarchy: for deeply bound states far from thresholds, continuum effects may remain perturbative; for near-threshold states, the same mechanism can qualitatively alter state assignments and wave-function composition [1503.02968][1511.05316].

## 2. Dressed-state formalism

The standard UQM ansatz for a hadron \(A\) is
\[
\mid \psi_A \rangle = {\cal N} \left[ \mid A \rangle + \sum_{BC \ell J} \int d \vec{K} \, k^2 dk \, \mid BC \ell J;\vec{K} k \rangle  \frac{ \langle BC \ell J;\vec{K} k \mid T^{\dagger} \mid A \rangle } {E_a - E_b - E_c} \right] ,
\]
where \({\cal N}\) is a normalization factor, \(|A\rangle\) is the bare quark-model state, and \(|BC\,\ell J;\vec K\,k\rangle\) is an intermediate two-hadron channel. Here \(k\) is the relative radial momentum, \(\ell\) the relative orbital angular momentum, and \(J\) the coupled total angular momentum of the continuum state [1503.02968][1506.05015].

The denominator
\[
\frac{1}{E_a-E_b-E_c}
\]
is the propagator-like factor controlling the virtual admixture of each channel. It makes threshold sensitivity explicit: when \(E_a\) lies close to \(E_b+E_c\), continuum components are enhanced. This is the basic reason the UQM is particularly relevant for open-flavor threshold phenomena [1506.05015].

Observables are then evaluated on the dressed state,
\[
O=\langle \psi_A|\hat O|\psi_A\rangle ,
\]
so that every matrix element contains a valence contribution and a continuum contribution absent in the naive quenched model [1503.02968][1503.03064]. In the baryon formulation, the same structure is written with an explicit pair-creation strength \(\gamma\),
\[
\mid \psi_A \rangle = {\cal N}_A \left\{ \mid A \rangle + \gamma \sum_{BC\, l\, J} \int d \vec{K}\, k^2 dk \, \mid BC,l,J; \vec{K},k \rangle \, \frac{ \langle BC,l,J; \vec{K},k \mid T^{\dagger}( ^{3}P_{0}) \mid A \rangle } {\Delta E_{A \rightarrow BC}(k)} \right\},
\]
with
\[
\Delta E_{A \rightarrow BC}(k) = m_A - \sqrt{m_B^2 + k^2} - \sqrt{m_C^2 + k^2},
\]
so the one-loop interpretation is manifest [1712.03919].

The literature also contains nonperturbative coupled-basis implementations. In one \(X(3872)\) study, the total state is taken as
\[
|\Psi\rangle = |\Phi_\alpha\rangle + \sum_\beta \chi_\beta(P)\,|\phi_A\phi_B\,\beta\rangle,
\]
and the continuum-induced correction is encoded in an energy-dependent mass-shift matrix \(\mathcal G(E)\), rather than in a single perturbative self-energy integral [1901.02484]. In another implementation for \(X(3872)\), the wavefunction is written as
\[
\Psi_{JM_J} = c_0 \Psi^{(2q)}_{JM_J} + \sum_{i=1}^{N} c_i \Psi^{(4q)}_{i,JM_J},
\]
with the physical mass obtained from a generalized eigenvalue problem,
\[
[(H)-E_n (N)]\, C_n = 0,
\]
in a nonorthogonal coupled basis [1906.09690].

## 3. Pair creation, self-energy, and general coupled-channel structure

The dynamical core of most UQM implementations is the \(^{3}P_0\) quark–antiquark pair-creation mechanism. In this model, the vacuum creates a \(q\bar q\) pair with \(J^{PC}=0^{++}\), i.e. spin triplet, \(P\)-wave, total \(J=0\) [1503.02968][1712.03919]. In meson spectroscopy, this pair-creation vertex appears both in the dressed wavefunction and in the self-energy correction,
\[
M_a = E_a + \Sigma(E_a),
\]
with
\[
\Sigma(E_a) = \sum_{BC\ell J} \int_0^{\infty} k^2 dk \,  \frac{\left| M_{A \rightarrow BC}(k) \right|^2}{E_a - E_b - E_c} .
\]
Here \(E_a\) is the bare mass, often taken from the Godfrey–Isgur relativized quark model, and \(M_{A\to BC}(k)\) is the transition amplitude induced by \(T^\dagger\) [1503.02968][1511.05316].

A frequently used refinement is to replace the constant pair-creation strength by an effective one to suppress unphysical heavy-quark pair creation [1503.02968][1503.03064]. A stronger modification appears in later coupled-channel studies that found the naive \(^{3}P_0\) operator produced unrealistically large negative shifts. These works introduced damping factors such as
\[
e^{-r^2/4f^2}
\quad\text{and}\quad
e^{-R_{AV}^2/R_0^2},
\]
to suppress high-momentum pair creation and pair creation far from the source hadron [1712.04457][2301.12388][2406.00957].

This technical issue is one of the main internal controversies of UQM phenomenology. With the original \(^{3}P_0\) operator, some light-meson and charmonium calculations found mass shifts so large that they challenged the use of the valence quark model as a sensible zeroth-order approximation. In response, improved operators were proposed to reduce the size of the shifts to a phenomenologically acceptable range [1712.04457][2301.12388]. This suggests that quantitative UQM results are highly sensitive to the assumed transition kernel.

More general coupled-channel results were derived by Burns for a broad class of “non-flip, triplet” models, including the \(^{3}P_0\) model. In the limit of spin-degenerate continuum multiplets, the self-energy operator becomes diagonal in \(S\) and \(L\), and loop-induced mixing vanishes between valence states of different spin or orbital angular momentum. Expanding around the spin-averaged energy yields the relation
\[
\delta E_{nSLJ}=Z_{nL}\,\delta M_{nSLJ},
\]
with
\[
Z_{nL}=\frac{1}{1+\langle \omega(E_{nL})\rangle_{nL}},
\]
so physical spin splittings are suppressed relative to bare ones by the valence probability \(Z\) [1411.2485]. This explains why several successful quenched-quark-model relations survive unquenching, including the vanishing \(P\)-wave hyperfine splitting and the approximate stability of the \(S\)-wave hyperfine–\(e^+e^-\) width relation [1411.2485].

## 4. Meson spectroscopy, threshold states, and heavy quarkonia

The best-known applications of the UQM are in charmonium and bottomonium. In the review literature, self-energy effects in heavy quarkonia are described as modest in absolute size but spectroscopically important: about \(2\%-6\%\) in charmonium and about \(1\%\) in bottomonium, with relative mass shifts of a few tens of MeV [1503.02968][1511.05316]. Because heavy-quark spectroscopy often works at the tens-of-MeV level, these shifts can be decisive for state assignments.

Threshold sensitivity is the central phenomenological message. The \(X(3872)\) is the canonical example. In one UQM interpretation, it is neither a pure molecular state nor a pure \(c\bar c\) state, but a \(\chi_{c1}(2^3P_1)\) \(c\bar c\) core strongly mixed with meson–meson continuum components, with approximately \(45\%\) core and \(55\%\) continuum [1503.02968][1506.05015]. Other implementations produce different decompositions: one nonperturbative coupled-channel treatment finds a predominantly \(D^0\bar D^{*0}\) molecular state with only \(3.62\%\) total \(c\bar c\), although that \(c\bar c\) piece is overwhelmingly \(2^3P_1\) [1901.02484]; another study using a modified \(^{3}P_0\) operator finds about \(69.5\%\) \(c\bar c\) and \(30.5\%\) meson–meson for \(X(3872)\) [1906.09690]. The coexistence of these results is not a contradiction of formalism but a reminder that UQM phenomenology depends strongly on the treatment of the continuum sector and on the transition operator.

The same mechanism has been applied to heavy-light threshold states. A systematic study of \(D_{s0}^*(2317)\) and \(D_{s1}(2460)\) in a mixed two-quark/four-quark framework found large downward continuum-induced shifts and substantial meson–meson content: \(D_{s0}\) with \(42.7\%\) \(c\bar s\) and \(57.3\%\) meson–meson, and \(D_{s1}\) with \(35.9\%\) \(c\bar s\) and \(64.1\%\) meson–meson. The dominant channels are \(DK\)-like for the \(0^+\) state and \(D^*K\)-like for the \(1^+\) state [2111.04677]. This supports a mixed-state interpretation rather than a purely conventional or purely molecular one.

In bottomonium, recent UQM studies emphasize that high-lying states above open-bottom thresholds are strongly influenced by continuum channels. One coupled-channel calculation of \(\Upsilon(5S)\) found that the bare mass
\[
10927.6~\text{MeV}
\]
is shifted downward by
\[
\Delta M=-31.4~\text{MeV}
\]
to
\[
10896.2~\text{MeV},
\]
closer to the PDG value for \(\Upsilon(10860)\), and concluded that \(\Upsilon(10860)\) should be viewed as a strongly mixed \(b\bar b-B^{(*)}_{(s)}B^{(*)}_{(s)}\) state [2507.13882]. A more systematic bottomonium analysis finds that high states often contain large non-\(b\bar b\) components, that hadronic loops induce substantial \(S\)-\(D\) mixing for some vector states, and that continuum fractions can be tens of percent or larger for highly excited levels [2501.15110]. By contrast, another 2025 study comparing quenched and unquenched bottomonium concludes that much of the coupled-channel effect can be absorbed into refitted quenched-model parameters, while still leaving detectable consequences in decay constants and continuum probabilities [2503.10178]. This suggests that in bottomonium the UQM can function both as a dynamical threshold tool and as an explanation of the surprising robustness of quenched constituent models.

## 5. Baryons, sea quarks, and nonvalence observables

In the baryon sector, the UQM extends the \(qqq\) constituent quark model by explicit \(q\bar q\) pair creation, producing baryon–meson components in the dressed baryon wavefunction [1712.03919][1608.07629]. The main motivation is that several observables are intrinsically sensitive to sea quarks: the proton flavor asymmetry, the spin decomposition of the proton, strangeness form-factor observables, and some electromagnetic and weak couplings.

The Gottfried-sum-rule violation is a standard example. In the UQM,
\[
S_G = \int_0^1 dx \frac{F_2^p(x)-F_2^n(x)}{x}
= \frac{1}{3} - \frac{2}{3} \int_0^1 dx \left[ \bar{d}(x) - \bar{u}(x) \right],
\]
so the observed \(S_G=0.2281\pm0.0065\) implies \(\bar d \neq \bar u\). UQM calculations attribute this naturally to baryon–meson continuum components, particularly pion-cloud fluctuations such as \(p\to n\pi^+\), which enhance \(\bar d\) over \(\bar u\) [1506.05015][1511.05316]. In a pion-cloud truncation of the baryon UQM, the flavor asymmetry equals the orbital angular momentum contribution,
\[
{\cal A}(p)=\Delta L,
\]
and the model yields \(\Delta L=0.118\) or \(0.158\), depending on the experimental input used to fix the asymmetry [1608.07629].

The same framework improves weak and electromagnetic transition observables. For the radiative decay \(\Delta\to N\gamma\), the conventional constituent quark model gives \(399\) keV, while experiment gives \(703\pm61\) keV. In the baryon UQM, pion loops raise the prediction to \(554\) keV, adding kaons gives \(582\) keV, and including \(\pi,K,\eta,\eta'\) yields \(608\) keV, showing that the pion cloud provides the dominant correction [1712.03919]. For semileptonic beta decays, the naive quark model predicts
\[
g_A(n\to p)=\frac{5}{3},
\qquad
g_A(\Sigma^-\to n)=-\frac{1}{3},
\]
whereas experiment gives
\[
g_A(n\to p)=1.2701\pm0.0025,
\qquad
g_A(\Sigma^-\to n)=-0.340\pm0.017.
\]
The UQM reduces the neutron axial coupling to \(1.34\) while leaving \(\Sigma^-\to n\) comparatively stable, again with pionic intermediate states supplying the dominant sea-quark correction [1712.03919].

Strangeness observables in the proton are predicted to be very small. The strange magnetic moment operator
\[
\vec{\mu}_{s}=\sum_{i} \mu_{i,s}\left[2\vec{s}(q_{i})+\vec{l}(q_{i})- 2\vec{s}(\bar{q}_{i})-\vec{l}(\bar{q}_{i}) \right]
\]
gives
\[
\mu_s=0.0006\,\mu_N,
\]
and the strange radius operator
\[
R^{2}_{s}= \sum^{5}_{i=1}e_{i,s}\left(\vec{r}_{i}-\vec{R}_{\rm cm}\right)^2
\]
gives
\[
\langle R_s^2\rangle = -0.004\ \mathrm{fm}^2.
\]
These results are described as negligible but compatible with experiment and recent lattice calculations [1506.05015][1511.05316]. A plausible implication is that the baryon UQM is most consequential for observables that directly expose nonvalence structure, while leaving some traditional successes of the constituent quark model largely intact.

## 6. Hybrid extensions, later developments, and limitations

The original meson UQM reviews already proposed extending the formalism beyond ordinary meson loops to include hybrid mesons. In that broader Fock-space picture,
\[
|\Psi\rangle = |q\bar q\rangle + |q\bar q q\bar q\rangle + |q\bar q g\rangle + \cdots,
\]
where \(|q\bar q g\rangle\) denotes a hybrid component with an explicit constituent gluon [1503.02968]. The hybrid–quarkonium coupling is written as
\[
\langle \mathcal H | V | Q \rangle = \langle q\bar q g | V | q\bar q \rangle,
\]
the analogue of the \(^{3}P_0\) transition matrix element for hybrid loops [1503.02968]. In the Coulomb-gauge QCD treatment reviewed there, the effective adiabatic potential is
\[
E^\nu(r) = - \frac{\tau}{r} + \beta r + \frac{\nu \pi}{r} + C ,
\]
with \(\nu=0\) for ordinary mesons and \(\nu=1,2,\dots\) for excited gluonic surfaces, yielding a lowest charmonium hybrid multiplet near \(4.476~\text{GeV}\) and a lowest bottomonium hybrid multiplet near \(11.055~\text{GeV}\) [1503.02968]. This suggests that at higher excitation energies, especially above about \(4\) GeV in charmonium, hybrid loops may become as important as ordinary meson loops.

A distinct line of development is the use of chiral dynamics rather than \(^{3}P_0\) pair creation as the effective continuum-coupling mechanism in heavy-light systems. In the “unified unquenched quark model for heavy-light mesons with chiral dynamics,” coupled-channel effects are induced by Goldstone-boson emission, the dressed masses satisfy a once-subtracted self-energy equation,
\[
M= M_A-\operatorname{Re} \sum_{BC} \int_{0}^{\infty} \frac{\left(M_0-M\right)\overline{\left|\langle BC,\boldsymbol{q}|\mathcal{H}_{I} |A \rangle \right|^{2}}{\left(M-E_{BC}\right)\left(M_0-E_{BC}\right)} q^{2} d q ,
\]
and a relativistic correction term in the strong transition amplitude is required to reproduce the observed large widths of radial excitations [2312.04765]. Later work extending this framework to higher \(D\) and \(D_s\) excitations finds that most higher states are shifted significantly downward by coupled-channel effects and uses these shifts in spectroscopy assignments such as \(D_{s1}(2933)^+\) and \(D_{sJ}(3040)^+\) [2605.17832].

Several limitations recur across the literature. The UQM is not a first-principles QCD calculation; it depends on an underlying constituent-quark Hamiltonian, on the chosen transition mechanism (\(^{3}P_0\), chiral emission, or variants), on channel truncation, and often on phenomenological regulators [1503.03064][1712.04457]. Different implementations can yield very different continuum fractions for the same state, as the case of \(X(3872)\) illustrates [1503.02968][1901.02484][1906.09690]. At the same time, broad coupled-channel theorems indicate that some quenched-quark-model regularities survive unquenching and may even explain why simple constituent descriptions remain empirically effective over large parts of the spectrum [1411.2485][2503.10178].

Taken together, these results define the UQM not as a single rigid model but as a family of coupled-channel constituent-quark frameworks. Their common content is the explicit dressing of valence hadrons by continuum degrees of freedom; their main utility is the treatment of thresholds, self-energies, and sea-quark effects; and their main open issue is the quantitative control of the transition kernel and channel space.

Source: https://www.emergentmind.com/topics/unquenched-quark-model