-
Hadrons, subatomic particles that are composed of quarks and gluons, occupy a large spectral scope; the lightest hadron, the pion, has a mass of
$ M_\pi \approx 0.14\;\rm{ GeV} $ , while heavy hadrons are heavier than$ 10\; \rm{ GeV} $ [1]. People expect that the underlying theory, quantum chromodynamics (QCD) [2], can illuminate the hadron spectrum and unify the description of light and heavy hadrons. QCD is a non-Abelian local gauge field theory of strong interaction and has been shown to be consistent with experimental observation. Because of the emergent phenomena at the hadronic scale, i.e., confinement and dynamical chiral symmetry breaking (DCSB), non-perturbative QCD is still a relatively unknown part of the standard model.Confinement provides an intrinsic wavelength,
$ \lambda_c\approx 0.5 \;\rm{ fm} $ , for the propagation of quarks and gluons. They behave like the parton at$ r<\lambda_c $ and show different propagation modes at$ r>\lambda_c $ . The propagation of quarks and gluons would certainly be affected by the finite size of hadrons. Surveying hadron physics using QCD requires a non-perturbative method. As a well-established non-perturbative approach, lattice QCD (lQCD) [3-5], a lattice gauge theory formulated on a grid, has led to much progress with respect to hadron physics. While lQCD resorts to brute-force calculation, a functional method like the Dyson-Schwinger equation and Bethe-Salpeter equation (DSBSE) [6-8] approach is complementary to lQCD.In this work, we aim at unifying the description of light, heavy, and heavy-light mesons via the DSBSE approach. In this Poincaré covariant framework, the quark propagator is presented by the Gap equation [6-8]①,
$\begin{split} S_f^{ - 1}(k) =& {Z_2}({\rm i}\gamma \cdot k + {Z_m}{m_f})\\ & + \frac{4}{3}{{\bar g}^2}{Z_1}\int_{dq}^\Lambda {{D_{\mu \nu }}} (l){\gamma _\mu }{S_f}(q)\Gamma _\nu ^f(k,q), \end{split}$
(1) where
$ f = \{u,d,s,c,b,t\} $ represents the quark flavor,$ l = k-q $ ,$ S_{f} $ is the quark propagator,$ m_f $ is the current quark mass,$ \Gamma^f_\nu $ is the quark-gluon vertex,$ D_{\mu\nu} $ is the gluon propagator, and$ \bar{g} $ is the coupling constant.$ Z_1 $ ,$ Z_2 $ , and$ Z_m $ are the renormalization constants of the quark-gluon vertex, the quark field, and the quark mass, respectively.$ \int^\Lambda_{d q} = \int ^{\Lambda} {\rm d}^{4} q/(2\pi)^{4} $ stands for a Poincaré-invariant regularized integration, where$ \Lambda $ is the regularization scale. A meson corresponds to a pole in the quark-antiquark scattering kernel [9], and the Bethe-Salpeter amplitude (BSA),$ \Gamma^{fg}(k;P) $ , where k and P are the relative and the total momentum of the meson, respectively,$ P^2 = -M^2 $ , and M is the meson mass, is solved by the Bethe-Salpeter equation (BSE) [8-10],$[{\Gamma ^{fg}}(k;P)]_\beta ^\alpha = \int_{dq}^\Lambda {[{K^{fg}}(k,q;P)]_{\sigma \beta }^{\alpha \delta }} [{\chi ^{fg}}(q;P)]_\delta ^\sigma ,$
(2) where
$ K^{fg}(k,q;P) $ is the quark-antiquark scattering kernel and$ \alpha $ ,$ \beta $ ,$ \sigma $ , and$ \delta $ are the Dirac indexes.$ \chi^{fg}(q;P) = S_{f}(q_{+}) \Gamma^{fg}(q;P) S_{g}(q_{-}) $ is the wave function,$ q_{+} = q + \iota P/2 $ ,$ q_{-} = q - (1-\iota) P/2 $ , and$ \iota $ is the partitioning parameter describing the momentum partition between the quark and antiquark and does not affect the physical observables.A promising and consistent method to solve the problem of the meson spectrum is building a quark-gluon vertex and constructing a scattering kernel. The forms of the quark-gluon vertex and scattering kernel have been investigated [11], while the most widely used and technically simple one is the rainbow-ladder (RL) approximation,
${\bar g^2}{Z_1}{D_{\mu \nu }}(l)\Gamma _\nu ^f(k,q) \to {[{Z_2}]^2}\tilde D_{\mu \nu }^f(l){\gamma _\nu },$
(3) $[{K^{fg}}(k,q;P)]_{\sigma \beta }^{\alpha \delta } \to - \frac{4}{3}{[{Z_2}]^2}\tilde D_{\mu \nu }^{fg}(l)[\gamma _\mu ^{}]_\sigma ^\alpha [{\gamma _\nu }]_\beta ^\delta ,$
(4) where
$ \tilde{D}^{fg}_{\mu\nu}(l) = \left(\delta_{\mu\nu}-\frac{l_{\mu}l_{\nu}}{l^{2}}\right){\cal G}^{fg}(l^2) $ and$ \tilde{D}^{f}_{\mu\nu}(l) = \left(\delta_{\mu\nu}-\frac{l_{\mu}l_{\nu}}{l^{2}}\right) {\cal G}^f(l^2) $ are the effective quark-antiquark interactions. In the original RL approximation,$ {\cal G}^{fg} = {\cal G}^{f} $ is flavor-symmetrical and modeled by [12]${{\cal G}^f}(s) = {\cal G}_{\rm IR}^f(s) + {{\cal G}_{\rm UV}}(s),$
(5) ${\cal G}_{\rm IR}^f(s) = 8{\pi ^2}\frac{{D_f^2}}{{\omega _f^4}}{{\rm e}^{ - s/\omega _f^2}},$
(6) ${{\cal G}_{\rm UV}}(s) = \frac{{8{\pi ^2}\gamma _m^{}{\cal F}(s)}}{{{\rm{ln}}[\tau + {{(1 + s/\Lambda _{\rm QCD}^2)}^2}]}},$
(7) where
$ s = l^2 $ .$ {\cal G}^f_{\rm IR}(s) $ is the infrared interaction responsible for DCSB, with$ D_f^2\omega_f $ expressing the interaction strength and$ \omega_f $ the interaction width in the momentum space. The form given in Eq. (6) is used because it enables the natural extraction of a monotonic running-coupling and gluon mass [12], whose relationship to QCD can be traced [13].$ {\cal G}_{\rm UV}(s) $ keeps the one-loop perturbative QCD limit in the ultraviolet range.$ {\cal F}(s) = [1 - \exp(-s^2/ [4m_{t}^{4}])]/s $ ,$ \gamma_{m}^{} = 12/(33- 2N_{f}) $ , with$ m_{t} = 1.0\; \rm{ GeV}\, $ ,$ \tau = {\rm e}^{10} - 1 $ ,$ N_f = 5 $ , and$ \Lambda_{\rm{QCD}} = 0.21 \;\rm{ GeV}\, $ . The values of$ m_{t} $ and$ \tau $ are chosen different from Ref. [12], so that$ {\cal G}_{\rm UV}(s) $ is largely suppressed in the infrared region and the dressed function$ {\cal G}_{\rm IR}^{fg}(s) $ is qualitatively right in the limit$ m_f \to \infty $ or$ m_g \to \infty $ .A nontrivial property of
$ \Gamma^f_\nu $ is its dependence on the quark flavor as a result of the dressing effect. By the same token,$ K^{fg} $ depends on the flavors of the scattering quark and antiquark. For a unified description of the system with different quarks, the flavor dependence of$ \Gamma^f_\nu $ and$ K^{fg} $ should be taken into account appropriately, irrespective of the model forms used. The RL approximation is phenomenologically successful for the pseudoscalar and vector mesons [12, 14-16]. The best parameters are$ (D_f^2\omega_f)^{1/3} = 0.8\;\rm{ GeV} $ and$ \omega_f = 0.5\;\rm{ GeV} $ for light mesons [12], and$ (D_f^2\omega_f)^{1/3} \approx [0.6,0.7] \;\rm{ GeV} $ and$ \omega_f = 0.8\;\rm{ GeV} $ for heavy mesons [17, 18]. The strength decreases and$ \omega_f $ increases as the quark mass increases, showing that heavy-flavor quarks probe shorter distances than light-flavor quarks at the corresponding quark-gluon vertexes [19]. The RL approximation fails to describe the heavy-light mesons because of the lack of flavor asymmetry in Eq. (5)-Eq. (7). The spectrum has a larger error than the quarkonia, and the decay constants are entirely false [20, 21]. The heavy-light meson problem has been surveyed for 20 years using this approach [22-27], yet no satisfactory solution has been found. -
To introduce flavor asymmetry, one should involve the axial-vector Ward-Takahashi identity (av-WTI), which guarantees the ground-state pseudoscalar mesons as Goldstone bosons of DCSB [14, 15],
$\begin{split} {P_\mu }\Gamma _{5\mu }^{fg}(k;P) =& S_f^{ - 1}({k_ + }){\rm i}{\gamma _5} + {\rm i}{\gamma _5}S_g^{ - 1}({k_ - })\\ & - {\rm i}[{m_f} + {m_g}]\Gamma _5^{fg}(k;P), \end{split}$
(8) where
$ \Gamma^{fg}_{5\mu} $ and$ \Gamma^{fg}_5 $ are the axial vector and pseudoscalar vertex, respectively. Considering the equations of$ S^{f,g} $ ,$ \Gamma^{fg}_{5\mu} $ , and$ \Gamma^{fg}_5 $ in the RL approximation, Eq. (8) leads to$\begin{split} &\int_{dq}^\Lambda {{{\cal G}^{fg}}} (s){\gamma _\alpha }[{S_f}({q_ + }){\rm i}{\gamma _5} + {\rm i}{\gamma _5}{S_g}({q_ - })]{\gamma _\beta } = \\ &\int_{dq}^\Lambda {{\gamma _\alpha }} [{{\cal G}^f}(s){S_f}({q_ + }){\rm i}{\gamma _5} + {{\cal G}^g}(s){\rm i}{\gamma _5}{S_g}({q_ - })]{\gamma _\beta }. \end{split}$
(9) Eq. (9) tells us that
$ {\cal G}^{fg}(s) $ is some medium value of$ {\cal G}^{f}(s) $ and$ {\cal G}^{g}(s) $ . Considering the scalar part of the propagator,$S^f(q) = -{\rm i}{\not\!{q}}\sigma^f_v(q^2) + \sigma^f_s(q^2)$ , we obtain$ {\cal G}^{fg}(s) = (\sigma^f_s(q^2_+){\cal G}^{f}(s) + \sigma^g_s(q^2_-){\cal G}^{g}(s))/(\sigma^f_s(q^2_+) + \sigma^g_s(q^2_-)) $ . It is well known that the infrared value of$ \sigma^f_s(q^2) $ is proportional to the inverse of the interaction strength, and the width of$ \sigma^f_s(q^2) $ is proportional to$ \omega_f $ . We thus assume$ {\cal G}^{fg}(s) $ to be${{\cal G}^{fg}}(s) = {\cal G}_{\rm IR}^{fg}(s) + {{\cal G}_{\rm UV}}(s),$
(10) ${\cal G}_{\rm IR}^{fg}(s) = 8{\pi ^2}\frac{{{D_f}}}{{\omega _f^2}}\frac{{{D_g}}}{{\omega _g^2}}{{\rm e}^{ - s/({\omega _f}{\omega _g})}}.$
(11) $ {\cal G}_{\rm UV}(s) $ is unchanged from Eq. (7), and as we are dealing with five active quarks,$ {\cal G}_{\rm UV}(s) $ is independent of the quark flavors. The effective interation$ \tilde{D}^{fg}_{\mu\nu} $ represents the total dressing effect of the gluon propagator and the two quark-gluon vertexes. Eq. (11) suggests that the quark and antiquark contribute equally to the interaction strength and width.The preservation of the av-WTI could be checked numerically using the Gell-Mann-Oakes-Renner (GMOR) relation, which is equivalent to the av-WTI [14, 15],
${\tilde f_{{0^ - }}}: = ({m_f} + {m_g}){\rho _{{0^ - }}}/M_{{0^ - }}^2 = {f_{{0^ - }}},$
(12) where
$ M_{0^-} $ is the mass of the pseudoscalar meson and$ f_{0^-} $ is the leptonic decay constant.$ f_{0^-} $ and$ \rho_{0^-} $ are defined by${f_{{0^ - }}}{P_\mu }: = {Z_2}{N_c}\;{\rm{tr}}\int_{dk}^\Lambda \gamma _5^{}\gamma _\mu ^{}{S_f}({k_ + })\Gamma _{{0^ - }}^{fg}(k;P){S_g}({k_ - }),$
(13) ${\rho _{{0^ - }}}: = {Z_4}{N_c}\;{\rm{tr}}\int_{dk}^\Lambda \gamma _5^{}{S_f}({k_ + })\Gamma _{{0^ - }}^{fg}(k;P){S_g}({k_ - }),$
(14) where
$ Z_4 = Z_2 Z_m $ ,$ N_c $ is the color number,$ \rm{tr} $ is the trace of the Dirac index, and$ \Gamma^{fg}_{0^{-}} $ is the BSA of pseudoscalar mesons. The BSA is normalized by [28]$\begin{split} 2{P_\mu } =& {N_c}\frac{\partial }{{\partial {P_\mu }}}\int_{{{dq}}}^\Lambda {} {\rm tr}[\Gamma (q; - K)\\ & \times S({q_ + })\Gamma (q;K)S({q_ - })]{|_{{P^2} = {K^2} = - {M^2}}}, \end{split}$
(15) where
$ N_c = 3 $ is the color number. Before discussing the details and results, we first assure the reader of the preservation of the av-WTI by comparing$ f_{0^-} $ and$ \tilde{f}_{0^-} $ in Fig. 1. They deviate from each other by no more than 3% for all pseudoscalar mesons considered here. We conclude that the av-WTI is perfectly preserved in our approach.Figure 1. (color online) Decay constants of the ground-state pseudoscalar mesons:
$f_{0^-}$ defined by Eq. (13) and$\tilde{f}_{0^-}$ defined by Eq. (12) and Eq. (14) are our results,$f_{\rm{lQCD}}$ are the lattice QCD data given in Table 1. -
In Eq. (11),
$ D_{f,g} $ and$ \omega_{f,g} $ are parameters expressing the flavor-dependent quark-antiquark interaction. However, the flavor dependence of these parameters is a priori unknown. Herein, we treat the$ D_{f} $ and$ \omega_{f} $ of each flavor as free parameters. Working in the isospin symmetry limit, we have 4 independent quarks up to the b quark mass: u (or d), s, c, and b. There are 3 parameters for each flavor:$ D_f $ ,$ \omega_f $ , and$ m_f $ . In total, there are 12 parameters.$ \omega_u $ is treated as an independent variable, while the other 11 parameters are dependent variables, which are fitted by 11 observables: the masses and decay constants of$ \pi $ , K,$ \eta_c $ , and$ \eta_b $ , and the masses of D,$ D_s $ , and B. All the masses and decay constants of the ground-state pseudoscalar mesons (except$ \eta $ and$ \eta^\prime $ ) and all the ground-state vector mesons are predicted.The masses and decay constants of the ground-state pseudoscalar mesons are listed in Table 1. Our outputs are quite stable, with
$ \omega_u $ varying by 10% around$ 0.5 \;\rm{ GeV} $ . With$ \omega_u \in [0.45,0.55] \;\rm{ GeV} $ , the masses are almost unchanged and the decay constants vary within 1.2%. Our output of$ M_{B^{\pm}_s} $ deviates from the experimental value by only$ 0.01 \;\rm{ GeV} $ , which is impossible in the original RL truncated DSBSE. The flavor dependence of the quark gluon interaction even has a significant effect on the$ B_c $ meson.$ M_{B_c} $ produced by the original RL truncated DSBSE is$ 0.11\;\rm{ GeV} $ larger than the experimental value [18]. We reduce the error to less than$ 0.02 \;\rm{ GeV} $ herein. Our output of$ f_D $ ,$ f_{D^{\pm}_s} $ ,$ f_{B} $ ,$ f_{B^{\pm}_s} $ , and$ f_{B_c} $ are all consistent with the lattice QCD results, with deviation of less than 6%. Note that our$ f_{D^{\pm}_s} $ is also in good agreement with a recent experimental measurement [29]. The only absent mesons in Table 1 are$ \eta $ and$ \eta^\prime $ , which are affected by the axial anomaly [30, 31] and beyond our present purpose.herein lQCD expt. herein lQCD $M_{\pi}$ 0.138 $\ast$ 0.138(1) $f_{\pi}$ 0.093 0.093(1) $M_K$ 0.496 $\ast$ 0.496(1) $f_K$ 0.111 0.111(1) $M_D$ 1.867 1.865(3) 1.867(1) $f_D$ 0.151(1) 0.150(1) $M_{D^{\pm}_s}$ 1.968 1.968(3) 1.968(1) $f_{D^{\pm}_s}$ 0.181(1) 0.177(1) $M_{\eta_c}$ 2.984 $\ast$ 2.984(1) $f_{\eta_c}$ 0.278 0.278(2) $M_{B}$ 5.279 5.283(8) 5.279(1) $f_{B}$ 0.141(2) 0.134(1) $M_{B^{\pm}_s}$ 5.377(1) 5.366(8) 5.367(1) $f_{B^{\pm}_s}$ 0.168(2) 0.163(1) $M_{B_c}$ 6.290(3) 6.276(7) 6.275(1) $f_{B_c}$ 0.312(1) 0.307(10) $M_{\eta_b}$ 9.399 $\ast$ 9.399(2) $f_{\eta_b}$ 0.472 0.472(5) Table 1. Masses and decay constants of the ground-state pseudoscalar mesons (in GeV). We use the convention
$f_\pi = 0.093\;\rm{ GeV}$ . The lQCD data are taken from:$M_D$ and$M_{D_s}$ - Ref. [32];$M_B$ and$M_{B_s}$ - Ref. [33];$M_{B_c}$ - Ref. [34];$f_\pi$ and$f_K$ - Ref. [35];$f_D$ ,$f_{D_s}$ ,$f_B$ , and$f_{B_s}$ - Ref. [36];$f_{\eta_c}$ and$f_{\eta_b}$ - Ref. [37];$f_{B_c}$ - Ref. [38].$M_{\pi}$ ,$M_K$ ,$M_{\eta_c}$ , and$M_{\eta_b}$ here and$M_\Upsilon$ in Table 2 are usually used to fit the quark masses in lQCD calculations [39], so there are no lQCD predictions for these quantities. The exprimental data are taken from Ref. [1]. Note that we work in the isospin symmetry limit, so the average value among or between the isospin multiplet is cited for$\pi$ , K, D, and B mesons. All data are given with an accuracy of$0.001 \;\rm{ GeV}$ . In our production, the underlined values are those used to fit the 11 dependent variables, and the others are our output, with the uncertainty corresponding to$\omega_u \in [0.45,0.55]\; \rm{ GeV}$ . The decay constants are fitted to the lQCD data because an accurate and complete experimental estimate of these data is still lacking.A further confirmation of our model is given by the vector mesons. Our predictions of the static vector meson masses and decay constants are listed in Table 2. The decay constant is defined by analogy to Eq. (13)
herein lQCD expt. herein lQCD $M_{\rho}$ 0.724(2) 0.780(16) 0.775(1) $f_{\rho}$ 0.149(1) – $M_{K^*}$ 0.924(2) 0.933(1) 0.896(1) $f_{K^*}$ 0.160(2) – $M_\phi$ 1.070(1) 1.032(16) 1.019(1) $f_\phi$ 0.191(1) 0.170(13) $M_{D^*}$ 2.108(4) 2.013(14) 2.009(1) $f_{D^*}$ 0.174(4) 0.158(6) $M_{D^{*\pm}_s}$ 2.166(7) 2.116(11) 2.112(1) $f_{D^{*\pm}_s}$ 0.206(2) 0.190(5) $M_{J/\psi}$ 3.132(2) 3.098(3) 3.097(1) $f_{J/\psi}$ 0.304(1) 0.286(4) $M_{B^*}$ 5.369(5) 5.321(8) 5.325(1) $f_{B^*}$ 0.132(3) 0.131(5) $M_{B^{*\pm}_s}$ 5.440(1) 5.411(5) 5.415(2) $f_{B^{*\pm}_s}$ 0.152(2) 0.158(4) $M_{B^*_c}$ 6.357(3) 6.331(7) – $f_{B^*_c}$ 0.305(5) 0.298(9) $M_\Upsilon$ 9.454(1) $\ast$ 9.460(1) $f_\Upsilon$ 0.442(3) 0.459(22) Table 2. Masses and decay constants of ground-state vector mesons (in GeV). The lQCD data are taken from:
$M_{\rho}$ - Ref. [40];$M_{K^*}$ - Ref. [41];$M_\phi$ and$f_\phi$ - Ref. [42];$M_{D^*}$ ,$f_{D^*}$ ,$M_{D^{*\pm}_s}$ ,$f_{D^{*\pm}_s}$ ,$M_{B^*}$ ,$f_{B^*}$ ,$M_{B^{*\pm}_s}$ , and$f_{B^{*\pm}_s}$ - Ref. [43];$M_{J/\psi}$ and$f_{J/\psi}$ - Ref. [44];$M_{B^*_c}$ - Ref. [34];$f_{B^*_c}$ - Ref. [38];$f_\Upsilon$ - Ref. [45]. The exprimental data are taken from Ref. [1], and the average value between the isospin multiplet is cited for$M_{D^*}$ . Hitherto, the$B^*_c$ meson has not been discovered experimentally. All data are cited with an accuracy of$0.001 \;\rm{ GeV}$ . The uncertainties of our results correspond to$\omega_u \in [0.45,0.55]\; \rm{ GeV}$ .${f_{{1^ - }}}{M_{{1^ - }}} = {Z_2}{N_c}\;{\rm{tr}}\int_{dk}^\Lambda \gamma _\mu ^{}{S_f}({k_ + })\Gamma _{{1^ - }}^{\mu ,fg}(k;P){S_g}({k_ - }),$
(16) where
$ M_{1^-} $ is the vector meson and$ \Gamma^{\mu,fg}_{1^{-}} $ is the vector meson BSA. The vector mesons also show a weak dependence on$ \omega_u \in [0.45,0.55]\; \rm{ GeV} $ . The deviation from experimental or lQCD values decreases as the mass increases. The mass deviation is approximately 6% for the$ \rho $ meson, decreasing to approximately 1% for the heavy mesons. The decay constant deviation is approximately 12% for the$ \phi $ meson, decreasing to less than 7% for the heavy mesons. This deviation can be attributed to the systematic error of the RL truncation [16]. The success of the pattern of the flavor-dependent interaction, Eq. (10,11,7), is shown by the fact that the deviation is on the same order for both the open-flavor mesons and the quarkonia. We can see again that the flavor dependence has a significant effect on$ B_c $ mesons. While$ M_{B^*_c}\approx 6.54 \;\rm{ GeV} $ and$ f_{B^*_c} \approx 0.43 \;\rm{ GeV} $ in the original RL truncated DSBSE [18], our results, i.e.,$ M_{B^*_c} \approx 6.357\; \rm{ GeV} $ and$ f_{B^*_c} \approx 0.305\; \rm{ GeV} $ , are more consistent with the lQCD predictions.$ B^*_c $ has not been discovered experimentally, so both our and lQCD predictions await experimental verification.Finally, we investigate the flavor dependence of the quark-antiquark interaction. In the heavy quark limit, the dressing of the quark-gluon vertex can be ignored and our adopted interaction is in agreement with QCD [13], so the interaction should saturate
$ {\cal G}^{ff}(l^2) \xrightarrow{m_f\to\infty} 4\pi\alpha_s \frac{Z(l^2)}{l^2} $ , where$ \alpha_s $ is the strong-interaction constant and$ Z(l^2) $ is the dressing function of the gluon propagator, defined by$ \Delta_{\mu\nu}(l) = (\delta_{\mu\nu}-\frac{l_\mu l_\nu}{l^2})\frac{Z(l^2)}{l^2} $ , where$ \Delta_{\mu\nu}(l) $ is the dressed gluon propagator. As we fix$ N_f = 5 $ , both$ \alpha_s $ and$ Z(l^2) $ are independent of the interacting quarks. Phenomenally, the parameters$ D_f $ and$ \omega_f $ should approach a constant as the quark mass increases. In the chiral limit, the interaction is enhanced because of the dressing of the quark-gluon vertex [46-50], which is necessary to trigger chiral symmetry breaking. The potential is properly defined by the Fourier transform of the interaction. For the interesting infrared part of our model, we have$ {\cal V}_{\rm{IR}}^{ff}(\vec{r}) = \int {\rm{d}}^3\vec{l}\, {\cal G}_{\rm{IR}}^{ff}(l^2) {\rm{e}}^{-\vec{l}\cdot\vec{r}/\omega^2_f} \propto {\rm{e}}^{-\vec{r}^2/R_{f}^2} $ , where$ \vec{r} $ is the space coordinate and$ R_{f} = 2/\omega_f $ expresses the radius of the quark-gluon interaction. Additionally, we adopt the following quantity to describe the interaction strength [51]:${\sigma _f} = \frac{1}{{4\pi }}\int_{\Lambda _{{\rm{QCD}}}^{\rm{2}}}^{{{({\rm{10}}{\Lambda _{{\rm{QCD}}}})}^{\rm{2}}}} {} {\rm d}s{\mkern 1mu} {{\cal G}^{ff}}(s)*s.$
(17) The quark mass dependence of
$ \sigma_f $ and$ R_{f} $ is depicted in Fig. 2. The interaction strength and radius decrease as the quark mass increases, which is expected by the fact that the quark-gluon vertex dressing effect should decrease as the quark mass increases [47]. The interaction radius,$ 2/\sqrt{\omega_{f}\omega_{g}} $ , also expresses another fact: the quarks and gluons have a maximum wavelength of the hadron size [52]. -
The fitted parameters corresponding to
$ \omega_u = 0.45, 0.50, 0.55 \;\rm{ GeV} $ are listed in Table A1. The quark mass$ \bar{m}_f^{\zeta} $ is defined by$\bar m_f^\zeta = {\hat m_f}/{\left[ {\frac{1}{2}{\rm{ln}}\frac{{{\zeta ^{\rm{2}}}}}{{\Lambda _{{\rm{QCD}}}^{\rm{2}}}}} \right]^{{\gamma _m}}},\tag{A1}$
(18) ${\hat m_f} = \mathop {\lim }\limits_{{p^2} \to \infty } {\left[ {\frac{1}{2}{\rm{ln}}\frac{{{{\rm{p}}^{\rm{2}}}}}{{\Lambda _{{\rm{QCD}}}^{\rm{2}}}}} \right]^{{\gamma _m}}}{M_f}({p^2}),\tag{A2}$
(19) where
$ \zeta $ is the renormalization scale,$ \hat{m}_f $ is the renormalization-group invariant current-quark mass, and$ M_f(p^2) $ is the quark mass function,$S_f(p) = \dfrac{Z_f(p^2,\zeta^2)}{{\rm i}\gamma\cdot p + M_f(p^2)}$ .
Pattern for flavor-dependent quark-antiquark interaction
- Received Date: 2020-05-09
- Available Online: 2020-11-01
Abstract: A flavor-dependent kernel is constructed based on the rainbow-ladder truncation of the Dyson-Schwinger and Bethe-Salpeter equation approach of quantum chromodynamics. The quark-antiquark interaction is composed of a flavor-dependent infrared part and a flavor-independent ultraviolet part. Our model gives a successful and unified description of the light, heavy, and heavy-light ground pseudoscalar and vector mesons. For the first time, our model shows that the infrared-enhanced quark-antiquark interaction is stronger and wider for lighter quarks.