Unquenched Quark Model Overview

• The Unquenched Quark Model is a framework that extends the constituent quark model by incorporating quark–antiquark pair creation and hadronic continuum mixing, leading to modified state assignments and mass shifts.
• It uses a dressed-state formalism where each hadron is represented as a superposition of a bare valence core and continuum components, with self-energy corrections being crucial near threshold regions.
• The model’s methodologies—ranging from the ³P₀ pair creation mechanism to chiral dynamics—offer systematic insights into heavy quarkonia, baryon observables, and nonvalence effects in spectroscopy.

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 (Santopinto et al., 2015, Tecocoatzi et al., 2015, Santopinto et al., 2015).

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 (Tecocoatzi et al., 2015, Bijker et al., 2017).

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 (Santopinto et al., 2015). 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 (García-Tecocoatzi et al., 2015, Santopinto et al., 2015).

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 (Santopinto et al., 2015, García-Tecocoatzi et al., 2015).

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 $qqq$0 is the relative radial momentum, $qqq$1 the relative orbital angular momentum, and $qqq$2 the coupled total angular momentum of the continuum state (Santopinto et al., 2015, Tecocoatzi et al., 2015).

The denominator

$qqq$3

is the propagator-like factor controlling the virtual admixture of each channel. It makes threshold sensitivity explicit: when $qqq$4 lies close to $qqq$5, continuum components are enhanced. This is the basic reason the UQM is particularly relevant for open-flavor threshold phenomena (Tecocoatzi et al., 2015).

Observables are then evaluated on the dressed state,

$qqq$6

so that every matrix element contains a valence contribution and a continuum contribution absent in the naive quenched model (Santopinto et al., 2015, Santopinto et al., 2015). In the baryon formulation, the same structure is written with an explicit pair-creation strength $qqq$7,

$qqq$8

with

$qqq$9

so the one-loop interpretation is manifest (Bijker et al., 2017).

The literature also contains nonperturbative coupled-basis implementations. In one $q\bar q$0 study, the total state is taken as

$q\bar q$1

and the continuum-induced correction is encoded in an energy-dependent mass-shift matrix $q\bar q$2, rather than in a single perturbative self-energy integral (Ortega et al., 2019). In another implementation for $q\bar q$3, the wavefunction is written as

$q\bar q$4

with the physical mass obtained from a generalized eigenvalue problem,

$q\bar q$5

in a nonorthogonal coupled basis (Tan et al., 2019).

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

The dynamical core of most UQM implementations is the $q\bar q$6 quark–antiquark pair-creation mechanism. In this model, the vacuum creates a $q\bar q$7 pair with $q\bar q$8, i.e. spin triplet, $q\bar q$9-wave, total $q\bar q \to (q\bar q)(q\bar q)$0 (Santopinto et al., 2015, Bijker et al., 2017). In meson spectroscopy, this pair-creation vertex appears both in the dressed wavefunction and in the self-energy correction,

$q\bar q \to (q\bar q)(q\bar q)$1

with

$q\bar q \to (q\bar q)(q\bar q)$2

Here $q\bar q \to (q\bar q)(q\bar q)$3 is the bare mass, often taken from the Godfrey–Isgur relativized quark model, and $q\bar q \to (q\bar q)(q\bar q)$4 is the transition amplitude induced by $q\bar q \to (q\bar q)(q\bar q)$5 (Santopinto et al., 2015, García-Tecocoatzi et al., 2015).

A frequently used refinement is to replace the constant pair-creation strength by an effective one to suppress unphysical heavy-quark pair creation (Santopinto et al., 2015, Santopinto et al., 2015). A stronger modification appears in later coupled-channel studies that found the naive $q\bar q \to (q\bar q)(q\bar q)$6 operator produced unrealistically large negative shifts. These works introduced damping factors such as

$q\bar q \to (q\bar q)(q\bar q)$7

to suppress high-momentum pair creation and pair creation far from the source hadron (Chen et al., 2017, Chen et al., 2023, Chen et al., 2024).

This technical issue is one of the main internal controversies of UQM phenomenology. With the original $q\bar q \to (q\bar q)(q\bar q)$8 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 (Chen et al., 2017, Chen et al., 2023). 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 $q\bar q \to (q\bar q)(q\bar q)$9 model. In the limit of spin-degenerate continuum multiplets, the self-energy operator becomes diagonal in $qqq \to qqq\,q\bar q$0 and $qqq \to qqq\,q\bar q$1, and loop-induced mixing vanishes between valence states of different spin or orbital angular momentum. Expanding around the spin-averaged energy yields the relation

$qqq \to qqq\,q\bar q$2

with

$qqq \to qqq\,q\bar q$3

so physical spin splittings are suppressed relative to bare ones by the valence probability $qqq \to qqq\,q\bar q$4 (Burns, 2014). This explains why several successful quenched-quark-model relations survive unquenching, including the vanishing $qqq \to qqq\,q\bar q$5-wave hyperfine splitting and the approximate stability of the $qqq \to qqq\,q\bar q$6-wave hyperfine–$qqq \to qqq\,q\bar q$7 width relation (Burns, 2014).

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 $qqq \to qqq\,q\bar q$8 in charmonium and about $qqq \to qqq\,q\bar q$9 in bottomonium, with relative mass shifts of a few tens of MeV (Santopinto et al., 2015, García-Tecocoatzi et al., 2015). 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 $A$0 is the canonical example. In one UQM interpretation, it is neither a pure molecular state nor a pure $A$1 state, but a $A$2 $A$3 core strongly mixed with meson–meson continuum components, with approximately $A$4 core and $A$5 continuum (Santopinto et al., 2015, Tecocoatzi et al., 2015). Other implementations produce different decompositions: one nonperturbative coupled-channel treatment finds a predominantly $A$6 molecular state with only $A$7 total $A$8, although that $A$9 piece is overwhelmingly $$\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] ,$$0 (Ortega et al., 2019); another study using a modified $$\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] ,$$1 operator finds about $$\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] ,$$2 $$\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] ,$$3 and $$\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] ,$$4 meson–meson for $$\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] ,$$5 (Tan et al., 2019). 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 $$\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] ,$$6 and $$\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] ,$$7 in a mixed two-quark/four-quark framework found large downward continuum-induced shifts and substantial meson–meson content: $$\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] ,$$8 with $$\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] ,$$9 ${\cal N}$0 and ${\cal N}$1 meson–meson, and ${\cal N}$2 with ${\cal N}$3 ${\cal N}$4 and ${\cal N}$5 meson–meson. The dominant channels are ${\cal N}$6-like for the ${\cal N}$7 state and ${\cal N}$8-like for the ${\cal N}$9 state (Tan et al., 2021). 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 $|A\rangle$0 found that the bare mass

$|A\rangle$1

is shifted downward by

$|A\rangle$2

to

$|A\rangle$3

closer to the PDG value for $|A\rangle$4, and concluded that $|A\rangle$5 should be viewed as a strongly mixed $|A\rangle$6 state (Chen et al., 18 Jul 2025). A more systematic bottomonium analysis finds that high states often contain large non-$|A\rangle$7 components, that hadronic loops induce substantial $|A\rangle$8-$|A\rangle$9 mixing for some vector states, and that continuum fractions can be tens of percent or larger for highly excited levels (Ni et al., 25 Jan 2025). 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 (Sultan et al., 13 Mar 2025). 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 $|BC\,\ell J;\vec K\,k\rangle$0 constituent quark model by explicit $|BC\,\ell J;\vec K\,k\rangle$1 pair creation, producing baryon–meson components in the dressed baryon wavefunction (Bijker et al., 2017, Bijker et al., 2016). 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,

$|BC\,\ell J;\vec K\,k\rangle$2

so the observed $|BC\,\ell J;\vec K\,k\rangle$3 implies $|BC\,\ell J;\vec K\,k\rangle$4. UQM calculations attribute this naturally to baryon–meson continuum components, particularly pion-cloud fluctuations such as $|BC\,\ell J;\vec K\,k\rangle$5, which enhance $|BC\,\ell J;\vec K\,k\rangle$6 over $|BC\,\ell J;\vec K\,k\rangle$7 (Tecocoatzi et al., 2015, García-Tecocoatzi et al., 2015). In a pion-cloud truncation of the baryon UQM, the flavor asymmetry equals the orbital angular momentum contribution,

$|BC\,\ell J;\vec K\,k\rangle$8

and the model yields $|BC\,\ell J;\vec K\,k\rangle$9 or $qqq$00, depending on the experimental input used to fix the asymmetry (Bijker et al., 2016).

The same framework improves weak and electromagnetic transition observables. For the radiative decay $qqq$01, the conventional constituent quark model gives $qqq$02 keV, while experiment gives $qqq$03 keV. In the baryon UQM, pion loops raise the prediction to $qqq$04 keV, adding kaons gives $qqq$05 keV, and including $qqq$06 yields $qqq$07 keV, showing that the pion cloud provides the dominant correction (Bijker et al., 2017). For semileptonic beta decays, the naive quark model predicts

$qqq$08

whereas experiment gives

$qqq$09

The UQM reduces the neutron axial coupling to $qqq$10 while leaving $qqq$11 comparatively stable, again with pionic intermediate states supplying the dominant sea-quark correction (Bijker et al., 2017).

Strangeness observables in the proton are predicted to be very small. The strange magnetic moment operator

$qqq$12

gives

$qqq$13

and the strange radius operator

$qqq$14

gives

$qqq$15

These results are described as negligible but compatible with experiment and recent lattice calculations (Tecocoatzi et al., 2015, García-Tecocoatzi et al., 2015). 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,

$qqq$16

where $qqq$17 denotes a hybrid component with an explicit constituent gluon (Santopinto et al., 2015). The hybrid–quarkonium coupling is written as

$qqq$18

the analogue of the $qqq$19 transition matrix element for hybrid loops (Santopinto et al., 2015). In the Coulomb-gauge QCD treatment reviewed there, the effective adiabatic potential is

$qqq$20

with $qqq$21 for ordinary mesons and $qqq$22 for excited gluonic surfaces, yielding a lowest charmonium hybrid multiplet near $qqq$23 and a lowest bottomonium hybrid multiplet near $qqq$24 (Santopinto et al., 2015). This suggests that at higher excitation energies, especially above about $qqq$25 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 $qqq$26 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,

$qqq$27

and a relativistic correction term in the strong transition amplitude is required to reproduce the observed large widths of radial excitations (Ni et al., 2023). Later work extending this framework to higher $qqq$28 and $qqq$29 excitations finds that most higher states are shifted significantly downward by coupled-channel effects and uses these shifts in spectroscopy assignments such as $qqq$30 and $qqq$31 (Ni et al., 18 May 2026).

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 ($qqq$32, chiral emission, or variants), on channel truncation, and often on phenomenological regulators (Santopinto et al., 2015, Chen et al., 2017). Different implementations can yield very different continuum fractions for the same state, as the case of $qqq$33 illustrates (Santopinto et al., 2015, Ortega et al., 2019, Tan et al., 2019). 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 (Burns, 2014, Sultan et al., 13 Mar 2025).

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.