--- title: Modified Godfrey–Isgur Model url: https://www.emergentmind.com/topics/modified-godfrey-isgur-model type: topic --- # Modified Godfrey–Isgur Model 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 [2211.09023][2405.19039]. ## 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 [2206.10132][1802.04938]. 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 [2405.19039]. 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 [1102.4970][1101.1584]. ## 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 [2206.10132][2510.16935]. 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 \[ \tilde G(r)\to \left(1+\frac{p^2}{E_1E_2}\right)^{1/2} \tilde G(r) \left(1+\frac{p^2}{E_1E_2}\right)^{1/2}, \] and the spin-dependent operators by analogous powers of \(m_i/E_i\) involving fitted exponents \(\epsilon_i\). These ingredients are essential to the model’s relativized character, especially for light quarks [1802.04938][1810.02694]. In several implementations the running coupling is represented by a frozen sum of exponentials. One explicit form is \[ \alpha_s(Q^2)=0.25 e^{-Q^2}+0.15 e^{-Q^2/10}+0.20 e^{-Q^2/1000}, \] 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 [2405.19039]. ## 3. Screened confinement as the distinguishing modification The central modification is the replacement of linear confinement by a screened form, \[ S(r)=br+c \quad\longrightarrow\quad S(r)=\frac{b(1-e^{-\mu r})}{\mu}+c, \] or equivalently \(br\to b(1-e^{-\mu r})/\mu\) when the constant term is written separately. This replacement has two explicit limits emphasized throughout the literature: \[ V^{\mathrm{scr}}(r)\approx br \quad (r\to 0), \qquad V^{\mathrm{scr}}(r)\to \frac{b}{\mu} \quad (r\to\infty). \] Thus the model agrees with GI at short distance but flattens at large distance [2211.09023][2508.21442]. The physical interpretation given in the spectroscopy papers is consistent across sectors. For higher excited states, large \(q\bar q\) or \(qqq\) 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 [2405.19039][2206.10132]. An important point of model identity is that the rest of the GI machinery is usually kept intact. In the \(5^{++}\) 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 [1810.02694][2405.19039]. 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 [1306.2874][1802.04938]. ## 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 \(\rho\) mesons, one frequently used light-meson parameter set is \[ m_u=m_d=0.163~\mathrm{GeV},\quad m_s=0.387~\mathrm{GeV},\quad b=0.221~\mathrm{GeV}^2,\quad c=-0.240~\mathrm{GeV}, \] \[ \sigma_0=1.799~\mathrm{GeV},\quad s=1.497,\quad \mu=0.0635~\mathrm{GeV}, \] with \(\epsilon_{\rm c}=-0.138\), \(\epsilon_{\rm sov}=0.157\), \(\epsilon_{\rm sos}=0.9726\), and \(\epsilon_{\rm t}=0.893\) [2206.10132]. Closely related light-meson studies use \[ m_{u(d)}=0.162~\mathrm{GeV},\quad m_s=0.377~\mathrm{GeV},\quad b=0.222~\mathrm{GeV}^2,\quad c=-0.228~\mathrm{GeV}, \] \[ \sigma_0=1.791~\mathrm{GeV},\quad s=0.711,\quad \mu=0.0779~\mathrm{GeV}, \] with \(\epsilon_c=-0.137\), \(\epsilon_t=0.493\), \(\epsilon_{\rm so(v)}=0.0550\), and \(\epsilon_{\rm so(s)}=0.366\) [2508.21442]. In bottomonium, a refit to 18 established states gives \[ m_b=5.027~\text{GeV},\quad b=0.21355,\quad c=-0.36804,\quad \mu=0.07426~\text{GeV}, \] and yields \(\chi^2/n=11.3\), compared with \(31.4\) for the original GI model on the same selected data [1802.04938]. In the kaon sector, a screened refit to 11 established states gives \(\chi^2/n=12.6\) versus \(90.2\) for the original GI parameters [1705.03144]. In the heavy-baryon study, the MGI and GI descriptions are numerically close, with \(\chi^2=1558/k^2\) for MGI and \(\chi^2=1776/k^2\) for GI [2405.19039]. The bound-state problem is solved with standard basis-expansion techniques. Several light-meson papers use a simple-harmonic-oscillator basis with \[ R_{nL}(r)=\sum_{n=1}^{n_{\max}} C_n\,R_{nL}^{\rm SHO}(r,\beta), \] taking \(n_{\max}=21\) and determining \(\beta\) variationally from \[ \frac{\partial E_{nL}}{\partial \beta_i}=0, \qquad \frac{\partial^2 E_{nL}}{\partial \beta_i^2}>0 \] [2508.21442][2510.16935]. Other works use the Gaussian expansion method for solving the relativized eigenvalue equation [1603.06417][2206.10132]. ## 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 \(D\), \(F\), \(G\), and \(n\ge 4\) states. The original GI model overshoots \(\Upsilon(6S)\) by about \(100\) MeV, while the screened model substantially reduces that mismatch; the same study extends the spectrum through \(8S\), \(6P\), \(5D\), \(4F\), and \(3G\) and uses the resulting wave functions for radiative, annihilation, hadronic, and open-bottom decay calculations [1802.04938]. For excited light vectors, screening lowers masses dramatically relative to GI. In the higher-\(\rho\) analysis, the MGI predictions \[ \rho(5^3S_1): 2542~\mathrm{MeV},\quad \rho(6^3S_1): 2774~\mathrm{MeV},\quad \rho(7^3S_1): 2967~\mathrm{MeV}, \] \[ \rho(4^3D_1): 2624~\mathrm{MeV},\quad \rho(5^3D_1): 2840~\mathrm{MeV},\quad \rho(6^3D_1): 3020~\mathrm{MeV} \] are lower than the original GI values by roughly \(300\)–\(500\) MeV, and the authors interpret this as strong evidence that screening is crucial for higher excitations [2206.10132]. In a related \(\rho\)-assignment study, the model supports \[ Y(2040)\to \rho(2^3D_1),\qquad \rho(1900)\to \rho(3^3S_1),\qquad \rho(2150)\to \rho(4^3S_1), \] after combining MGI masses with \(^{3}P_0\) decay widths [2102.05356]. The same mass-plus-decay strategy is used across many light-meson families. In the \(\omega\) family, the model gives \(\omega(3D)=2284.1\) MeV, and the corresponding strong-decay analysis leads to the conclusion that \(X(2232)\), \(X(2200)\), and \(X(2222)\) may be the same resonance and are most likely the \(\omega(3D)\) state [2211.09023]. In the \(a_4\) family, the MGI masses \[ a_4(1^3F_4)=1928~\mathrm{MeV},\quad a_4(2^3F_4)=2243~\mathrm{MeV},\quad a_4(3^3F_4)=2466~\mathrm{MeV},\quad a_4(4^3F_4)=2640~\mathrm{MeV}, \] \[ a_4(1^3H_4)=2405~\mathrm{MeV},\quad a_4(2^3H_4)=2589~\mathrm{MeV} \] are combined with QPC widths to argue that \(a_4(2610)\) is more likely \(2^3H_4\) than \(4^3F_4\) [2508.21442]. For the missing \(5^{++}\) family, the screened model predicts \(a_5(1H)\) and \(f_5(1H)\) at \(2.492\) GeV and \(f_5'(1H)\) at \(2.679\) GeV, with the screening systematically lowering masses relative to GI; that study reports \(\chi^2/n=82\) for the screened fit versus \(2638\) for the original GI model on 41 established meson states [1810.02694]. In the \(3^{++}\) sector, MGI predicts \[ a_3(1F)=1955~\mathrm{MeV},\quad a_3(2F)=2234~\mathrm{MeV},\quad f_3(1F)=1955~\mathrm{MeV},\quad f_3(2F)=2234~\mathrm{MeV}, \] supporting \(f_3(2050)\) as the ground state and \(f_3(2300)\) as the first radial excitation, while treating \(a_3(1875)\) and \(a_3(2030)\) as likely the same ground \(a_3\) state [2510.16935]. Kaon spectroscopy is another major application. A systematic kaon study treats the \(u/d\)-\(s\) 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 \(K_3(2075)\sim1F\), \(K_4(2310)\sim1G\), and \(K_2(1990)\sim2D\) [1705.03144]. In the highly excited scalar-kaon sector, the screened spectrum gives \[ K_0^*(4P)=2451.2~\text{MeV},\quad K_0^*(5P)=2658.9~\text{MeV},\quad K_0^*(6P)=2831.8~\text{MeV}, \] which is used, together with a predicted width of \(432\) MeV, to identify \(\kappa(2600)\) as \(K_0^*(5P)\) [2503.12393]. A related high-spin-kaon study applies the same framework to \(J^P=3^\pm,4^\pm,5^\pm\) 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 [2507.12072]. ## 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 [2405.19039]. 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 \(T_{cc}\) study, the MGI model gives lower masses than GI for the same diquark-antidiquark configuration, and for the \(1S\) isoscalar \(I(J^P)=0(1^+)\) state one quoted MGI value is \(3877\) MeV for \(\mu=70\), especially close to the observed \(T_{cc}(3875)^+\) mass [2407.19383]. 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 \[ \tilde V_{qq}(\mathbf p,\mathbf r)=\frac12 \tilde V_{q\bar q}(\mathbf p,\mathbf r), \] and the resulting \(sq\bar b\bar q\) mass near \(6150\) MeV disfavors identifying \(X(5568)\) as a compact tetraquark in that framework [1603.06417]. 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 \(^{3}P_0\), radiative, annihilation, and hadronic-transition calculations, and often lowers masses into better agreement with observed high-lying states [1802.04938]. 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 [1306.2874][1802.04938]. 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. Source: https://www.emergentmind.com/topics/modified-godfrey-isgur-model