Modified Godfrey–Isgur Model

• Modified Godfrey–Isgur Model is a relativized quark model that introduces color screening to soften long-distance confinement, addressing overpredictions in high excitation states.
• It retains key features of the original GI approach such as relativistic kinematics, smeared one-gluon exchange, and spin-dependent forces while modifying the confining potential.
• The model effectively simulates vacuum polarization, quark-pair creation, and channel-coupling effects, leading to improved spectral predictions in meson, baryon, and tetraquark studies.

Searching arXiv for recent and foundational papers on the modified Godfrey–Isgur model to support the article. The Modified Godfrey–Isgur model is a relativized quark model derived from the Godfrey–Isgur framework by modifying the long-distance confining interaction, most commonly through color screening, while retaining the GI treatment of relativistic kinematics, smeared one-gluon exchange, and spin-dependent forces. In the standard spectroscopy usage, the central replacement is

$br \;\to\; \frac{b(1-e^{-\mu r})}{\mu},$

or the corresponding form with an added constant $c$, so that the interaction remains approximately linear at short distance and saturates at large distance. Across the literature, this model is used to compute hadron masses and wave functions for mesons, baryons, and diquark-based tetraquarks, with the screening interpreted as an effective representation of vacuum polarization, quark-pair creation, string breaking, or coupled-channel effects that become important for higher excitations (Wang et al., 2022, Weng et al., 2024).

1. Definition and terminological scope

In the light- and heavy-hadron spectroscopy literature, the Modified Godfrey–Isgur model usually denotes the relativized quark model of Godfrey and Isgur with color-screening effects introduced into the confining interaction. The motivation stated repeatedly is that the original GI model, with a strictly linear long-distance potential, works well for many low-lying hadrons but tends to overpredict higher radial and orbital excitations, especially for light mesons and high-spin states (Feng et al., 2022, Wang et al., 2018).

The standard modification is therefore not a wholesale replacement of the GI framework. Rather, it preserves the relativized structure of the model and alters the confinement sector to soften the force at large interquark separation. In meson applications this screening is explicitly tied to vacuum polarization from dynamical light-quark pairs, string breaking, or channel coupling; in baryon applications it is presented as an effective way to encode $q\bar q$ or $qq$ pair-creation effects without solving a coupled-channel problem explicitly (Weng et al., 2024).

The phrase is, however, not perfectly uniform across the broader potential-model literature. In heavy-light Isgur–Wise studies, closely related Cornell-type models with an added constant $c$, or with different perturbative organizations of the linear and Coulomb terms, are described as “modified Godfrey–Isgur-type” in spirit rather than as the standard screened-GI construction. Those works are better regarded as phenomenological relatives than as the canonical MGI spectroscopy model (Hazarika et al., 2011, Hazarika et al., 2011).

2. Hamiltonian structure and relativization

The mesonic MGI Hamiltonian keeps the relativized GI kinetic term,

$$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$

and decomposes the effective interaction into the same operator classes used in GI: a short-range Coulombic one-gluon-exchange term, confinement, contact hyperfine, tensor, and spin-orbit contributions. Representative formulations write

$$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$

or, in the nonrelativistic limit,

$$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$

The hyperfine sector contains both the contact spin-spin and tensor interactions, while the spin-orbit part is split into vector color-magnetic and scalar Thomas-precession pieces (Feng et al., 2022, Chen et al., 19 Oct 2025).

A defining feature inherited from GI is relativization by nonlocal smearing and momentum-dependent operator dressing. The smearing is implemented through a Gaussian kernel,

$$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$

with a mass-dependent width $\sigma_{ij}$. The Coulomb piece is further modified by factors such as

$c$0

and the spin-dependent operators by analogous powers of $c$1 involving fitted exponents $c$2. These ingredients are essential to the model’s relativized character, especially for light quarks (Wang et al., 2018, Pang et al., 2018).

In several implementations the running coupling is represented by a frozen sum of exponentials. One explicit form is

$c$3

with corresponding coordinate-space kernels expressed through error functions after smearing. This is GI-like rather than unique to MGI, but it is part of the standard operator content retained by screened versions of the model (Weng et al., 2024).

3. Screened confinement as the distinguishing modification

The central modification is the replacement of linear confinement by a screened form,

$c$4

or equivalently $c$5 when the constant term is written separately. This replacement has two explicit limits emphasized throughout the literature: $c$6 Thus the model agrees with GI at short distance but flattens at large distance (Wang et al., 2022, Pang et al., 29 Aug 2025).

The physical interpretation given in the spectroscopy papers is consistent across sectors. For higher excited states, large $c$7 or $c$8 separations make vacuum polarization, light-quark pair creation, string breaking, and channel-coupling effects more important. The screened confinement is therefore intended to mimic unquenched dynamics phenomenologically rather than derive them from an explicit continuum calculation (Weng et al., 2024, Feng et al., 2022).

An important point of model identity is that the rest of the GI machinery is usually kept intact. In the $c$9 study, the authors state explicitly that only the confinement term is replaced by the screened form and the rest of the GI framework is retained. In baryon work, the scalar spin-orbit term is updated consistently with the screened central confinement, but the relativized kinetic energy, smeared OGE interaction, and GI-style momentum factors remain in place (Pang et al., 2018, Weng et al., 2024).

This also clarifies a common misconception. The MGI model is not itself a coupled-channel calculation. It is an effective screened-potential approximation to some of the same physics. That distinction becomes important near thresholds, where explicit continuum self-energies can induce state-dependent mass shifts and altered splittings beyond what a static screened potential alone can capture (Ferretti et al., 2013, Wang et al., 2018).

4. Parameterization and numerical realization

MGI studies do not use a single universal parameter set. Instead, the literature summarized here shows two common practices: refitting the screened model to a sector of interest, or importing a parameter set established in earlier MGI work for nearby systems. This suggests a pragmatic rather than strictly universal use of the framework.

For higher excited $q\bar q$0 mesons, one frequently used light-meson parameter set is

$q\bar q$1

$q\bar q$2

with $q\bar q$3, $q\bar q$4, $q\bar q$5, and $q\bar q$6 (Feng et al., 2022). Closely related light-meson studies use

$q\bar q$7

$q\bar q$8

with $q\bar q$9, $qq$0, $qq$1, and $qq$2 (Pang et al., 29 Aug 2025).

In bottomonium, a refit to 18 established states gives

$qq$3

and yields $qq$4, compared with $qq$5 for the original GI model on the same selected data (Wang et al., 2018). In the kaon sector, a screened refit to 11 established states gives $qq$6 versus $qq$7 for the original GI parameters (Pang et al., 2017). In the heavy-baryon study, the MGI and GI descriptions are numerically close, with $qq$8 for MGI and $qq$9 for GI (Weng et al., 2024).

The bound-state problem is solved with standard basis-expansion techniques. Several light-meson papers use a simple-harmonic-oscillator basis with

$c$0

taking $c$1 and determining $c$2 variationally from

$c$3

(Pang et al., 29 Aug 2025, Chen et al., 19 Oct 2025). Other works use the Gaussian expansion method for solving the relativized eigenvalue equation (Lü et al., 2016, Feng et al., 2022).

5. Spectroscopic applications

The model’s most visible impact appears in higher excitations. In higher bottomonium, the screened model is reported to be nearly identical to GI for low-lying states but substantially different for higher $c$4, $c$5, $c$6, and $c$7 states. The original GI model overshoots $c$8 by about $c$9 MeV, while the screened model substantially reduces that mismatch; the same study extends the spectrum through $$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$0, $$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$1, $$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$2, $$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$3, and $$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$4 and uses the resulting wave functions for radiative, annihilation, hadronic, and open-bottom decay calculations (Wang et al., 2018).

For excited light vectors, screening lowers masses dramatically relative to GI. In the higher-$$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$5 analysis, the MGI predictions

$$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$6

$$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$7

are lower than the original GI values by roughly $$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$8–$$\tilde H=\sqrt{m_1^2+\mathbf p^2}+\sqrt{m_2^2+\mathbf p^2}+\tilde V_{\mathrm{eff}}(\mathbf p,\mathbf r),$$9 MeV, and the authors interpret this as strong evidence that screening is crucial for higher excitations (Feng et al., 2022). In a related $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$0-assignment study, the model supports

$$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$1

after combining MGI masses with $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$2 decay widths (Li et al., 2021).

The same mass-plus-decay strategy is used across many light-meson families. In the $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$3 family, the model gives $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$4 MeV, and the corresponding strong-decay analysis leads to the conclusion that $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$5, $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$6, and $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$7 may be the same resonance and are most likely the $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$8 state (Wang et al., 2022). In the $$\tilde{V}^{\mathrm{eff}} = \tilde G_{12} +\tilde V^{\mathrm{cont}} +\tilde V^{\mathrm{tens}} +\tilde V^{\mathrm{so(v)}} +\tilde S_{12}(r) +\tilde V^{\mathrm{so(s)}},$$9 family, the MGI masses

$$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$0

$$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$1

are combined with QPC widths to argue that $$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$2 is more likely $$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$3 than $$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$4 (Pang et al., 29 Aug 2025). For the missing $$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$5 family, the screened model predicts $$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$6 and $$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$7 at $$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$8 GeV and $$V_{\mathrm{eff}}(r)=H^{\mathrm{conf}}+H^{\mathrm{hyp}}+H^{\mathrm{so}}.$$9 at $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$0 GeV, with the screening systematically lowering masses relative to GI; that study reports $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$1 for the screened fit versus $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$2 for the original GI model on 41 established meson states (Pang et al., 2018). In the $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$3 sector, MGI predicts

$$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$4

supporting $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$5 as the ground state and $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$6 as the first radial excitation, while treating $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$7 and $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$8 as likely the same ground $$\rho_{ij}(\mathbf r-\mathbf r') = \frac{\sigma_{ij}^3}{\pi^{3/2}} e^{-\sigma_{ij}^2(\mathbf r-\mathbf r')^2}, \qquad \tilde f(r)=\int d^3r'\,\rho_{ij}(\mathbf r-\mathbf r')f(r'),$$9 state (Chen et al., 19 Oct 2025).

Kaon spectroscopy is another major application. A systematic kaon study treats the $\sigma_{ij}$0-$\sigma_{ij}$1 system in MGI with screening and reports a substantially improved global fit over GI, then uses QPC decays to organize the observed kaon family and to predict missing states such as $\sigma_{ij}$2, $\sigma_{ij}$3, and $\sigma_{ij}$4 (Pang et al., 2017). In the highly excited scalar-kaon sector, the screened spectrum gives

$\sigma_{ij}$5

which is used, together with a predicted width of $\sigma_{ij}$6 MeV, to identify $\sigma_{ij}$7 as $\sigma_{ij}$8 (Pang et al., 16 Mar 2025). A related high-spin-kaon study applies the same framework to $\sigma_{ij}$9 states and finds that the screening effect has a bigger influence on high-spin kaons, motivating the use of MGI over the original GI model in that sector (Wang et al., 16 Jul 2025).

6. Extensions, comparative performance, and limitations

The MGI idea has also been extended beyond ordinary mesons. In the heavy-baryon relativized quark model with chromodynamics, the linear string potential is replaced by a screened one, and the resulting MGI and GI spectra are reported to be similar overall. All heavy baryons observed so far are described as three-quark states in both models, with MGI slightly better in the global fit but not dramatically so. The authors therefore present MGI as a physically motivated refinement rather than a radical restructuring of heavy-baryon spectroscopy (Weng et al., 2024).

Diquark-based tetraquark work uses two related strategies. One is a screened relativized diquark model, where GI-type dynamics with screened confinement are applied first to the diquark and antidiquark and then to the diquark–antidiquark bound state. In the $c$00 study, the MGI model gives lower masses than GI for the same diquark-antidiquark configuration, and for the $c$01 isoscalar $c$02 state one quoted MGI value is $c$03 MeV for $c$04, especially close to the observed $c$05 mass (Dong et al., 2024). The other strategy is a GI-based diquark model without the screened modification in the main calculation; in the open-charm/open-bottom tetraquark study, the effective change from mesons to diquarks is

$c$06

and the resulting $c$07 mass near $c$08 MeV disfavors identifying $c$09 as a compact tetraquark in that framework (Lü et al., 2016).

The model’s strengths are therefore specific rather than universal. It generally improves the description of higher radial and orbital excitations, provides numerical wave functions usable in $c$10, radiative, annihilation, and hadronic-transition calculations, and often lowers masses into better agreement with observed high-lying states (Wang et al., 2018). At the same time, the literature is explicit about several limitations. Screening is phenomenological, not a first-principles solution of continuum dynamics; low-lying states can be nearly unchanged relative to GI; and near-threshold systems may still require explicit coupled-channel or self-energy treatments. The comparison between screened-potential bottomonium and the unquenched bottomonium calculation based on GI bare masses makes this point especially clear: threshold effects can reorganize splittings in ways that are not reducible to a universal static flattening of the confining potential alone (Ferretti et al., 2013, Wang et al., 2018).

In that sense, the Modified Godfrey–Isgur model occupies a specific methodological niche. It is a GI-like relativized quark model in which screened confinement encodes some unquenched-QCD effects economically, making it particularly useful for the spectroscopy and decay phenomenology of highly excited hadrons, while remaining distinct from explicit coupled-channel, continuum, or lattice-QCD treatments.