Loading [MathJax]/jax/output/HTML-CSS/jax.js

Strangeonium-like hybrids on the lattice

  • The strangeonium-like sˉsg hybrids are investigated from lattice QCD in the quenched approximation. In the Coulomb gauge, spatially extended operators are constructed for 1 and (0,1,2)+ states with the color octet sˉs component being separated from the chromomagnetic field strength by the spatial distance r, whose matrix elements between the vacuum and the corresponding states are interpreted as Bethe-Salpeter (BS) wave functions. In each of the (1,2)+ channels, the masses and the BS wave functions are reliably derived. The 1+ ground state mass is approximately 2.1-2.2 GeV, and that of 2+ is approximately 2.3-2.4 GeV, whereas the mass of the first excited state is approximately 1.4 GeV higher. This mass splitting is much larger compared to that expected based on the phenomenological flux-tube model or constituent gluon model for hybrids, which is usually a few hundred MeV. The BS wave functions with respect to rexhibit clear radial nodal structures of a non-relativistic two-body system, which imply that r is a meaningful dynamical variable for these hybrids and motivate a color halo picture of hybrids, in which the color octet sˉs is surrounded by gluonic degrees of freedom. In the 1 channel, the properties of the lowest two states are consistent with those of ϕ(1020) and ϕ(1680). We did not obtain convincing information with respect to ϕ(2170). However, we argue that regardless of whether ϕ(2170) is a conventional sˉs meson or a sˉsg hybrid in the color halo scenario, the ratio of partial decay widths Γ(ϕη) and Γ(ϕη) observed by BESIII can be understood based on the mechanism of hadronic transition of a strangeonium-like meson in addition to ηη mixing.
  • [1] N. Brambilla, S. Eidelman, C. Hanhart et al., Phys. Rept. 873, 1-154 (2020), arXiv:1907.07583
    [2] P. Lacock, C. Michael, P. Boyle et al. (UKQCD), Phys. Lett. B 401, 308 (1997), arXiv:heplat/9611011 doi: 10.1016/S0370-2693(97)00384-5
    [3] J. J. Dudek, Phys. Rev. D 84, 074023 (2011), arXiv:1106.5515[hep-ph doi: 10.1103/PhysRevD.84.074023
    [4] L. Liu, G. Moir, M. Peardon et al. (Hadron Spectrum), JHEP 07, 126 (2012), arXiv:1204.5425[hep-ph
    [5] M. Tanabashi et al. (Particle Data Group), Phys. Rev. D 98, 030001 (2018)
    [6] S.-L. Zhu, Phys. Lett. B 625, 212 (2005), arXiv:hepph/0507025[hep-ph
    [7] B. Aubert et al. (BaBar), Phys. Rev. D 74, 091103 (2006), arXiv:hep-ex/0610018 doi: 10.1103/PhysRevD.74.091103
    [8] M. Ablikim et al. (BES), Phys. Rev. Lett. 100, 102003 (2008), arXiv:0712.1143[hep-ex doi: 10.1103/PhysRevLett.100.102003
    [9] C. Shen et al. (Belle), Phys. Rev. D 80, 031101 (2009), arXiv:0808.0006[hep-ex doi: 10.1103/PhysRevD.80.031101
    [10] D. Horn and J. Mandula, Phys. Rev. D 17, 898 (1978) doi: 10.1103/PhysRevD.17.898
    [11] N. Isgur, R. Kokoski, and J. Paton, Phys. Rev. Lett. 54, 869 (1985) doi: 10.1103/PhysRevLett.54.869
    [12] K. Juge, J. Kuti, and C. Morningstar, Nucl. Phys. B Proc. Suppl. 63, 326 (1998), arXiv:hep-lat/9709131 doi: 10.1016/S0920-5632(97)00759-7
    [13] K. Juge, J. Kuti, and C. Morningstar, Phys. Rev. Lett. 82, 4400 (1999), arXiv:hep-ph/9902336 doi: 10.1103/PhysRevLett.82.4400
    [14] E. Braaten, C. Langmack, and D. H. Smith, Phys. Rev. D 90, 014044 (2014), arXiv:1402.0438[hep-ph doi: 10.1103/PhysRevD.90.014044
    [15] P. Guo, A. P. Szczepaniak, G. Galata et al., Phys. Rev. D 77, 056005 (2008), arXiv:0707.3156[hep-ph doi: 10.1103/PhysRevD.77.056005
    [16] P. Guo, A. P. Szczepaniak, G. Galata et al., Phys. Rev. D 78, 056003 (2008), arXiv:0807.2721[hep-ph doi: 10.1103/PhysRevD.78.056003
    [17] P. Guo, T. Ypez-Martnez, and A. P. Szczepaniak, Phys. Rev. D 89, 116005 (2014), arXiv:1402.5863[hep-ph doi: 10.1103/PhysRevD.89.116005
    [18] P. Lacock, C. Michael, P. Boyle et al. (UKQCD), Phys. Rev. D 54, 6997 (1996), arXiv:heplat/9605025[hep-lat
    [19] C. W. Bernard et al. (MILC), Phys. Rev. D 56, 7039 (1997), arXiv:hep-lat/9707008[hep-lat
    [20] X. Liao and T. Manke, (2002), arXiv: hep-lat/0210030 [hep-lat]
    [21] C. Bernard, T. Burch, E. B. Gregory et al., Phys. Rev. D 68, 074505 (2003), arXiv:heplat/0301024[hep-lat
    [22] Z.-H. Mei and X.-Q. Luo, Int. J. Mod. Phys. A 18, 5713 (2003), arXiv:hep-lat/0206012[hep-lat
    [23] J. J. Dudek, R. G. Edwards, M. J. Peardon et al., Phys. Rev. Lett. 103, 262001 (2009), arXiv:0909.0200[hep-ph doi: 10.1103/PhysRevLett.103.262001
    [24] Y.-B. Yang, Y. Chen, G. Li et al., Phys. Rev. D 86, 094511 (2012), arXiv:1202.2205[hep-ph
    [25] J. J. Dudek and E. Rrapaj, Phys. Rev. D 78, 094504 (2008), arXiv:0809.2582[hep-ph
    [26] J. J. Dudek, R. Edwards, and C. E. Thomas, Phys. Rev. D 79, 094504 (2009), arXiv:0902.2241[hep-ph
    [27] Y. Ma, W. Sun, Y. Chen et al., (2019), arXiv: 1910.09819 [hep-lat]
    [28] C. J. Morningstar and M. J. Peardon, Phys. Rev. D 60, 034509 (1999), arXiv:hep-lat/9901004[hep-lat
    [29] Y. Chen et al., Phys. Rev. D 73, 014516 (2006), arXiv:hep-lat/0510074[hep-lat
    [30] S.-q. Su, L.-m. Liu, X. Li et al., Int. J. Mod. Phys. A 21, 1015 (2006), arXiv:hep-lat/0412034[hep-lat
    [31] Y.-B. Yang, Y. Chen, L.-C. Gui et al. (CLQCD), Phys. Rev. D 87, 014501 (2013), arXiv:1206.2086[hep-lat
    [32] C. Davies, E. Follana, I. Kendall et al. (HPQCD), Phys. Rev. D 81, 034506 (2010), arXiv:0910.1229[hep-lat doi: 10.1103/PhysRevD.81.034506
    [33] G. S. Bali, Phys. Rept. 343, 1 (2001), arXiv:hepph/0001312 doi: 10.1016/S0370-1573(00)00079-X
    [34] G.-J. Ding and M.-L. Yan, Phys. Lett. B 650, 390 (2007), arXiv:hep-ph/0611319 doi: 10.1016/j.physletb.2007.05.026
    [35] N. Akbar and S. Noor, (2020), arXiv: 2005.02626 [hepph]
    [36] J. Lees et al. (BaBar), Phys. Rev. D 86, 012008 (2012), arXiv:1103.3001[hep-ex doi: 10.1103/PhysRevD.86.012008
    [37] M. Ablikim et al. (BESⅢ), Phys. Rev. D 100, 032009 (2019), arXiv:1907.06015[hep-ex doi: 10.1103/PhysRevD.100.032009
    [38] M. Ablikim et al. (BESⅢ), Phys. Rev. Lett. 124, 112001 (2020), arXiv:2001.04131[hep-ex doi: 10.1103/PhysRevLett.124.112001
    [39] M. Ablikim et al. (BESⅢ), Phys. Rev. D 102, 012008 (2020), arXiv:2003.13064[hep-ex doi: 10.1103/PhysRevD.102.012008
    [40] M. Ablikim et al. (BESⅢ), Phys. Rev. D 99, 032001 (2019), arXiv:1811.08742[hep-ex doi: 10.1103/PhysRevD.99.032001
  • [1] N. Brambilla, S. Eidelman, C. Hanhart et al., Phys. Rept. 873, 1-154 (2020), arXiv:1907.07583
    [2] P. Lacock, C. Michael, P. Boyle et al. (UKQCD), Phys. Lett. B 401, 308 (1997), arXiv:heplat/9611011 doi: 10.1016/S0370-2693(97)00384-5
    [3] J. J. Dudek, Phys. Rev. D 84, 074023 (2011), arXiv:1106.5515[hep-ph doi: 10.1103/PhysRevD.84.074023
    [4] L. Liu, G. Moir, M. Peardon et al. (Hadron Spectrum), JHEP 07, 126 (2012), arXiv:1204.5425[hep-ph
    [5] M. Tanabashi et al. (Particle Data Group), Phys. Rev. D 98, 030001 (2018)
    [6] S.-L. Zhu, Phys. Lett. B 625, 212 (2005), arXiv:hepph/0507025[hep-ph
    [7] B. Aubert et al. (BaBar), Phys. Rev. D 74, 091103 (2006), arXiv:hep-ex/0610018 doi: 10.1103/PhysRevD.74.091103
    [8] M. Ablikim et al. (BES), Phys. Rev. Lett. 100, 102003 (2008), arXiv:0712.1143[hep-ex doi: 10.1103/PhysRevLett.100.102003
    [9] C. Shen et al. (Belle), Phys. Rev. D 80, 031101 (2009), arXiv:0808.0006[hep-ex doi: 10.1103/PhysRevD.80.031101
    [10] D. Horn and J. Mandula, Phys. Rev. D 17, 898 (1978) doi: 10.1103/PhysRevD.17.898
    [11] N. Isgur, R. Kokoski, and J. Paton, Phys. Rev. Lett. 54, 869 (1985) doi: 10.1103/PhysRevLett.54.869
    [12] K. Juge, J. Kuti, and C. Morningstar, Nucl. Phys. B Proc. Suppl. 63, 326 (1998), arXiv:hep-lat/9709131 doi: 10.1016/S0920-5632(97)00759-7
    [13] K. Juge, J. Kuti, and C. Morningstar, Phys. Rev. Lett. 82, 4400 (1999), arXiv:hep-ph/9902336 doi: 10.1103/PhysRevLett.82.4400
    [14] E. Braaten, C. Langmack, and D. H. Smith, Phys. Rev. D 90, 014044 (2014), arXiv:1402.0438[hep-ph doi: 10.1103/PhysRevD.90.014044
    [15] P. Guo, A. P. Szczepaniak, G. Galata et al., Phys. Rev. D 77, 056005 (2008), arXiv:0707.3156[hep-ph doi: 10.1103/PhysRevD.77.056005
    [16] P. Guo, A. P. Szczepaniak, G. Galata et al., Phys. Rev. D 78, 056003 (2008), arXiv:0807.2721[hep-ph doi: 10.1103/PhysRevD.78.056003
    [17] P. Guo, T. Ypez-Martnez, and A. P. Szczepaniak, Phys. Rev. D 89, 116005 (2014), arXiv:1402.5863[hep-ph doi: 10.1103/PhysRevD.89.116005
    [18] P. Lacock, C. Michael, P. Boyle et al. (UKQCD), Phys. Rev. D 54, 6997 (1996), arXiv:heplat/9605025[hep-lat
    [19] C. W. Bernard et al. (MILC), Phys. Rev. D 56, 7039 (1997), arXiv:hep-lat/9707008[hep-lat
    [20] X. Liao and T. Manke, (2002), arXiv: hep-lat/0210030 [hep-lat]
    [21] C. Bernard, T. Burch, E. B. Gregory et al., Phys. Rev. D 68, 074505 (2003), arXiv:heplat/0301024[hep-lat
    [22] Z.-H. Mei and X.-Q. Luo, Int. J. Mod. Phys. A 18, 5713 (2003), arXiv:hep-lat/0206012[hep-lat
    [23] J. J. Dudek, R. G. Edwards, M. J. Peardon et al., Phys. Rev. Lett. 103, 262001 (2009), arXiv:0909.0200[hep-ph doi: 10.1103/PhysRevLett.103.262001
    [24] Y.-B. Yang, Y. Chen, G. Li et al., Phys. Rev. D 86, 094511 (2012), arXiv:1202.2205[hep-ph
    [25] J. J. Dudek and E. Rrapaj, Phys. Rev. D 78, 094504 (2008), arXiv:0809.2582[hep-ph
    [26] J. J. Dudek, R. Edwards, and C. E. Thomas, Phys. Rev. D 79, 094504 (2009), arXiv:0902.2241[hep-ph
    [27] Y. Ma, W. Sun, Y. Chen et al., (2019), arXiv: 1910.09819 [hep-lat]
    [28] C. J. Morningstar and M. J. Peardon, Phys. Rev. D 60, 034509 (1999), arXiv:hep-lat/9901004[hep-lat
    [29] Y. Chen et al., Phys. Rev. D 73, 014516 (2006), arXiv:hep-lat/0510074[hep-lat
    [30] S.-q. Su, L.-m. Liu, X. Li et al., Int. J. Mod. Phys. A 21, 1015 (2006), arXiv:hep-lat/0412034[hep-lat
    [31] Y.-B. Yang, Y. Chen, L.-C. Gui et al. (CLQCD), Phys. Rev. D 87, 014501 (2013), arXiv:1206.2086[hep-lat
    [32] C. Davies, E. Follana, I. Kendall et al. (HPQCD), Phys. Rev. D 81, 034506 (2010), arXiv:0910.1229[hep-lat doi: 10.1103/PhysRevD.81.034506
    [33] G. S. Bali, Phys. Rept. 343, 1 (2001), arXiv:hepph/0001312 doi: 10.1016/S0370-1573(00)00079-X
    [34] G.-J. Ding and M.-L. Yan, Phys. Lett. B 650, 390 (2007), arXiv:hep-ph/0611319 doi: 10.1016/j.physletb.2007.05.026
    [35] N. Akbar and S. Noor, (2020), arXiv: 2005.02626 [hepph]
    [36] J. Lees et al. (BaBar), Phys. Rev. D 86, 012008 (2012), arXiv:1103.3001[hep-ex doi: 10.1103/PhysRevD.86.012008
    [37] M. Ablikim et al. (BESⅢ), Phys. Rev. D 100, 032009 (2019), arXiv:1907.06015[hep-ex doi: 10.1103/PhysRevD.100.032009
    [38] M. Ablikim et al. (BESⅢ), Phys. Rev. Lett. 124, 112001 (2020), arXiv:2001.04131[hep-ex doi: 10.1103/PhysRevLett.124.112001
    [39] M. Ablikim et al. (BESⅢ), Phys. Rev. D 102, 012008 (2020), arXiv:2003.13064[hep-ex doi: 10.1103/PhysRevD.102.012008
    [40] M. Ablikim et al. (BESⅢ), Phys. Rev. D 99, 032001 (2019), arXiv:1811.08742[hep-ex doi: 10.1103/PhysRevD.99.032001
  • 加载中

Cited by

1. Akbar, N., Zeeshan, A., Arshad, S. et al. Spectrum of hybrid charmonium, bottomonium, and Bc mesons by power series method[J]. European Physical Journal Plus, 2024, 139(10): 954. doi: 10.1140/epjp/s13360-024-05729-4
2. Rinaldi, M., Vento, V. Hybrid spectroscopy within the graviton soft-wall model[J]. Physical Review D, 2024, 109(11): 114030. doi: 10.1103/PhysRevD.109.114030
3. Ablikim, M., Achasov, M.N., Adlarson, P. et al. Measurement of the e + e − → KS0KL0 π 0 cross sections from √s = 2.000 to 3.080 GeV[J]. Journal of High Energy Physics, 2024, 2024(1): 180. doi: 10.1007/JHEP01(2024)180
4. Wang, Q.-N., Lian, D.-K., Chen, W. Predictions of the hybrid mesons with exotic quantum numbers JPC=2+-[J]. Physical Review D, 2023, 108(11): 114010. doi: 10.1103/PhysRevD.108.114010
5. Zhu, W., Wang, X. Determination of the resonant parameters of excited vector strangenia with data[J]. Chinese Physics C, 2023, 47(11): 113003. doi: 10.1088/1674-1137/acf034
6. Malabarba, B.B., Khemchandani, K.P., Torres, A.M. φ (2170) decaying to φη and φη′[J]. Physical Review D, 2023, 108(3): 036010. doi: 10.1103/PhysRevD.108.036010
7. Chen, F., Jiang, X., Chen, Y. et al. 1-+ hybrid meson in J /ψ radiative decays from lattice QCD[J]. Physical Review D, 2023, 107(5): 054511. doi: 10.1103/PhysRevD.107.054511
8. Zhu, W.J., Gao, B.S., Wang, X.L. et al. Study of e+e- →ηφ via initial state radiation at Belle[J]. Physical Review D, 2023, 107(1): 012006. doi: 10.1103/PhysRevD.107.012006
9. Ablikim, M., Achasov, M.N., Adlarson, P. et al. Measurement of e+e − → K+K − π 0 cross section and observation of a resonant structure[J]. Journal of High Energy Physics, 2022, 2022(7): 45. doi: 10.1007/JHEP07(2022)045
10. Huang, L.. Study of ϕ(2170) at BESIII[J]. Proceedings of Science, 2022.
11. Ma, Y., Sun, W., Chen, Y. et al. Color halo scenario of charmonium-like hybrids[J]. Chinese Physics C, 2021, 45(9): 093111. doi: 10.1088/1674-1137/ac0ee2
12. Ablikim, M., Achasov, M.N., Adlarson, P. et al. Study of the process e+e- →φη at center-of-mass energies between 2.00 and 3.08 GeV STUDY of the PROCESS ... ABLIKIM M. et al.[J]. Physical Review D, 2021, 104(3): 032007. doi: 10.1103/PhysRevD.104.032007
13. Liu, D.. Study of f(2170) at BESIII[J]. Proceedings of Science, 2021.

Figures(6) / Tables(5)

Get Citation
Yunheng Ma, Ying Chen, Ming Gong and Zhaofeng Liu. Strangeonium-like hybrids on the lattice[J]. Chinese Physics C. doi: 10.1088/1674-1137/abc241
Yunheng Ma, Ying Chen, Ming Gong and Zhaofeng Liu. Strangeonium-like hybrids on the lattice[J]. Chinese Physics C.  doi: 10.1088/1674-1137/abc241 shu
Milestone
Received: 2020-07-30
Article Metric

Article Views(2557)
PDF Downloads(37)
Cited by(13)
Policy on re-use
To reuse of Open Access content published by CPC, for content published under the terms of the Creative Commons Attribution 3.0 license (“CC CY”), the users don’t need to request permission to copy, distribute and display the final published version of the article and to create derivative works, subject to appropriate attribution.
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Email This Article

Title:
Email:

Strangeonium-like hybrids on the lattice

    Corresponding author: Ying Chen, cheny@ihep.ac.cn
  • 1. Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China
  • 2. School of Physics, University of Chinese Academy of Sciences, Beijing 100049, China

Abstract: The strangeonium-like sˉsg hybrids are investigated from lattice QCD in the quenched approximation. In the Coulomb gauge, spatially extended operators are constructed for 1 and (0,1,2)+ states with the color octet sˉs component being separated from the chromomagnetic field strength by the spatial distance r, whose matrix elements between the vacuum and the corresponding states are interpreted as Bethe-Salpeter (BS) wave functions. In each of the (1,2)+ channels, the masses and the BS wave functions are reliably derived. The 1+ ground state mass is approximately 2.1-2.2 GeV, and that of 2+ is approximately 2.3-2.4 GeV, whereas the mass of the first excited state is approximately 1.4 GeV higher. This mass splitting is much larger compared to that expected based on the phenomenological flux-tube model or constituent gluon model for hybrids, which is usually a few hundred MeV. The BS wave functions with respect to rexhibit clear radial nodal structures of a non-relativistic two-body system, which imply that r is a meaningful dynamical variable for these hybrids and motivate a color halo picture of hybrids, in which the color octet sˉs is surrounded by gluonic degrees of freedom. In the 1 channel, the properties of the lowest two states are consistent with those of ϕ(1020) and ϕ(1680). We did not obtain convincing information with respect to ϕ(2170). However, we argue that regardless of whether ϕ(2170) is a conventional sˉs meson or a sˉsg hybrid in the color halo scenario, the ratio of partial decay widths Γ(ϕη) and Γ(ϕη) observed by BESIII can be understood based on the mechanism of hadronic transition of a strangeonium-like meson in addition to ηη mixing.

    HTML

    I.   INTRODUCTION
    • The naive quark model describes hadrons as qˉq mesons and qqq baryons. Since quarks and gluons are the fundamental degrees of freedom of QCD, if gluons can act as building blocks similar to quarks to build up hadrons, from a phenomenological perspective, there may exist glueballs that are purely made up of gluons and hybrids that are composed of quarks and gluons. Glueballs, hybrids, and multiquark states (tetraquarks, pentaquarks, etc.) are usually called exotic hadrons in contrast to conventional qˉq mesons and qqq baryons. Exotic hadrons are long-standing topics of interest in theoretical and experimental studies of particle physics, especially with the discovery of numerous XYZ particles in the current era based on experimental investigations. These particles exhibit exotic properties in their production and decay processes and are potential candidates for exotic hadrons [1]. Given that qˉqg hybrids are made up of a quark-antiquark pair and a gluon, the JPC=1+ states are most interesting since this quantum number is prohibited in conventional qˉq mesons. There have been many theoretical studies on hybrids from a phenomenological perspective and the lattice QCD approach. It was determined that the lowest 1+ hybrid usually has a mass of approximately 1 GeV higher than the ground state vector meson with the same qˉq component. For example, the mass of the 1+ hybrid with light flavors was estimated to be approximately 1.9 GeV, whereas those of the strangeonium-like and the charmonium-like counterparts are approximately 2.1-2.3 GeV [2, 3] and 4.1-4.3 GeV [4], respectively. Reliable candidates for 1+ hybrids have not been experimentally determined to date. The vector charmonium-like state Y(4260) (or ψ(4230) named by PDG 2018 [5]), due to its very different properties compared to the conventional charmonia and the closeness of its mass to that of the 1+ cˉcg hybrid, has a possible assignment of 1 hybrid [6]. ϕ(2170) [5], also known as Y(2175), was first observed as a result of the BABAR Collaboration in the initial-state-radiation process e+eγISRϕf0(980) in 2006 [7] and was confirmed later by BES and Belle [8, 9]. The similarity of its property to Y(4260) also suggests a sˉsg hybrid interpretation of ϕ(2170).

      Phenomenologically, qˉqg hybrids are usually studied in the constituent gluon model [10], in which the gluon acts as an effective degree of freedom, similar to constituent quarks in the quark model, or the flux tube model in which the gluon is taken as a transverse vibration mode of the flux-tube that binds the qˉq pair [11]. For the hybrids with a heavy quark-anti-quark pair QˉQ, the gluonic excitations along the flux-tube are fast objects, such that in the Oppenheimer approximation [12-14], their distribution obeys cylinder symmetry along the QˉQ-axis, and their motion effects on the QˉQ can be taken as a centrifugal barrier, apart from the binding linear potential. Based on the hybrid potentials simulated from the lattice QCD, one can solve the Schrödinger equation for the QˉQ system to obtain predictions regarding the spectrum of hybrids with properly tuned parameters. There have also been phenomenological studies on heavy quarkonium-like hybrids in which gluonic excitations are treated in a mean field manner [15-17].

      Lattice QCD is an ab initio non-perturbative approach for the study of strong interactions in the low energy scale and is applied extensively to the investigation of hybrids [18-26]. The mass of the hybrids can be derived from the correlation functions of hybrid-like operators ˉqΓqB, where ˉqq is the color octet, B is the chromomagnetic field strength, Γ represents specific combinations of γ matrices, and the symbol represents any possible summation of the spatial indices of Γ and B. A recent lattice calculation [4] revealed that there exists a {1,(0,1,2)+} charmonium-like supermultiplet with nearly degenerate masses of approximately 4.2-4.4 GeV, which overlaps strongly with the hybrid-like operators. This observation implies that these states may have similar internal dynamics, whereas the spin-spin coupling of the ˉqq and B yields the different quantum numbers. In our previous work [27], the internal structure of this supermultiplet was investigated by calculating their Bethe-Salpeter (BS) wave functions based on the lattice QCD in the quenched approximation, in which the spatially extended interpolating field operators ˉqΓq(x,t)B(x+r,t) are introduced in the Coulomb gauge, whose matrix element between the vacuum and a state is defined as the BS wave function. It was determined that the BS wave functions of the states in this multiplet are very similar and exhibit interesting nodal structures, which imply that the distance between the cˉc and the B operator is a meaningful dynamic variable for hybrids.

      In this work, we extend the aforementioned study strategy to strangeonium-like hybrids and focus on the {1,(0,1,2)+} states, to determine if a situation similar to that of the cˉcg hybrids can also occur for the sˉsg states. In contrast, since the quantum numbers 1 and 0+ are permitted by the qˉq mesons, in these channels, we will also use the spatially extended sˉs operators for the quark fields with spatial separations to extract the related BS wave functions, from which we can investigate the internal structure of these states. By comparison of these two kinds of BS wave functions, we may obtain useful information on the possible different formation pattern of hybrids from conventional mesons. For ϕ(2170), since it can be 33S1, 23D1 sˉs, or a candidate for the vector sˉsg hybrid, its properties will be discussed based on the results of this study.

      This work is organized as follows. Section II gives a detailed description of our lattice setup and the numerical strategy, including the construction of the spatially extended operators, the data analysis procedure, and the results of the spectrum and BS wave functions. The discussion and the comparison of our results with those of relevant phenomenological studies are presented in Section III. Section IV is an overall summary.

    II.   NUMERICAL DETAILS
    • The pure gauge configurations are generated via the tadpole-improved gauge action [28, 29] on anisotropic lattices with an aspect ratio of ξ=as/at=5, where as and at are the spatial and temporal lattice spacing, respectively. Two lattices L3×T=163×160(β=2.4) and 243×192(β=2.8) with different lattice spacings are used to check the discretization artifacts. The parameters of the gauge ensembles are listed in Table 1, where the as values are determined from r10=410(20) MeV. For the strange valence quark, we use the tadpole-improved clover action, whose parameters are carefully tuned by requiring the dispersion relations of the vector and pseudoscalar mesons to be reproduced [30]. As will be addressed in the following sections, we will use spatially extended operators to calculate the relevant correlation functions. Therefore, the configurations are first fixed to the Coulomb gauge based on the standard gauge fixing procedure used in lattice QCD studies before the quark propagators are computed. In this investigation, we work in the quenched approximation for which the effects of dynamical quarks are omitted. As such, two-particle scattering states will not appear in the spectrum.

      βξas/fmLas/fmL3×TNconf
      2.450.222(2)3.55163×160500
      2.850.138(1)3.31243×192200

      Table 1.  The input parameters for the calculation. Values of the coupling β, anisotropy ξ, lattice spacing as, lattice size, and number of measurements are listed. as/r0 is determined by the static potential; the first error of as is the statistical error, and the second one is due to the uncertainty of the scale parameter r10=410(20) MeV.

    • A.   Interpolation field operators

    • The major goal of this study is to investigate the inner structure of sˉsg hybrids of the quantum numbers JPC=1,(0,1,2)+. We introduce two types of spatially extended operators as sink operators. The first type includes the following sˉsg operators:

      Ok1(r,t)=x,|r|=rˉsa(x,t)γ5sb(x,t)Babk(x+r,t),O0+(r,t)=x,|r|=rˉsa(x,t)γisb(x,t)Babi(x+r,t),Ok1+(r,t)=x,|r|=rˉsa(x,t)γisb(x,t)Babj(x+r,t)εijk,Ok2+(r,t)=x,|r|=rˉsa(x,t)γisb(x,t)Babj(x+r,t)|εijk|,

      (1)

      where i,j,k=1,2,3 are the spatial indices, a,b=1,2,3 are color indices, and Babi=12εijkFabjk is the chromomagnetic field strength. The summation over r with the same distance r leads to the operators having correct quantum numbers. It should be noted that on a hypercubic spacetime latice, the spin J=2 corresponds to T2E where T2 and E are the irreducible representations of the lattice symmetry group O, so the Ok2+(r,t) are the three components of T2. Generally speaking, the masses of the T2 state and the E state, which correspond to the same J=2 state, are different by O(ans) terms (here, the power n depends on the specific lattice setups). This difference vanishes in the continuum limit as0. In this work, we ignore the mass splitting due to the discretization effect and treat the mass of the T2 state as an approximation of the J=2 state, so we do not give the explicit expression of the E operators in Eq. (1). The possible ground hybrids in these four quantum numbers form a supermultiplet as expected. Obviously, the two constituent quarks are localized at the same space-time point, and the gluon component is placed at another space point. The BS wave function that we attempt to extract reflects the dynamics of these two parts.

      Since the quantum numbers JPC=0+,1 are conventional ones for qˉq mesons, we also introduce the second type of spatially extended sˉs operators by splitting the strange quark field s and its conjugate ˉs by the spatial separation r. This is explicitly expressed as

      P0+(r,t)=x,|r|=rˉs(x,t)γ5s(x+r,t),Pk1(r,t)=x,|r|=rˉs(x,t)γks(x+r,t),

      (2)

      where the summation over r with |r|=r is performed to guarantee the correct JPC.

      In practice, we calculate the wall-source correlation functions of these operators. For example, we use the following wall-source operators for 1+ states:

      OW,k(τ)=y,zˉsa(y,τ)γiBabj(z,τ)sb(z,τ)εijk,

      (3)

      where τ indicates the source time slice. The wall-source operators for other JPC states vary accordingly. The wall source operators for the sˉs operators are PW(τ)=y,zˉs(y,τ)Γs(z,τ) with Γ=γ5,γi for 0+ and 1, respectively. Finally, we calculate the correlation functions (for simplicity, we set τ=0 and omit the subscripts and superscripts that refer to specific symmetry channels and different spatial components) as

      C(r,t)=Ok(r,t)OW,k(0).

      (4)

      After the intermediate state insertion, the correlation function C(r,t) can be parameterized as

      C(r,t)=1Ncn12mnL30|O(r,t)|nn|OW(0)|0=1Ncn12mnL30|O(r,0)|nn|OW(0)|0emntnΦn(r)emnt,

      (5)

      where Nc is the degenerate degree of r=|r|, mn is the mass of the n-th state, and Φn(r) is defined as the corresponding Bethe-Salpeter wave function up to an irrelevant constant factor. It should be noted that mn is independent of r, so that if we simultaneously fit C(r,t) with different r using Eq. (5), we can obtain mn and Φn(r). To be specific, if nr different C(r,t) values with different r values are considered, and N mass terms are involved in the fitting model, the number of parameters to be fitted is Nnr+N. Since we usually have 20-30 statistically meaningful data points for each C(r,t), the number of degrees of freedom is large enough in the fitting procedure.

    • B.   Results of 1+ and 2+ states

    • We start with the 1+ channel, since JPC is a typical exotic quantum number which cannot be assigned to a qˉq meson in the quark model. After the correlation function C(r,t) is calculated for the β=2.4 and β=2.8 lattices, Eq. (5) is used for data analysis wherein we use N=3 mass terms. For the β=2.4 lattice, the r range is from 0 to 0.9 fm (converted using the lattice spacing as in Table 1), and the upper limit of the fit window [tmin,tmax] is uniformly set to be tmax=20 for all the C(r,t), while the tmin varies from 6 to 3. For the β=2.8 lattice, r has a value up to 0.8 fm, and tmax is set to tmax=30, with tmin varying from 9 to 6. For each lattice, we perform a simultaneous correlated fit for all the different C(r,t)values using the jackknife covariance matrix. Table 2 summarizes the fit results for the mass mn with n=1,2,3 for different time windows [tmin,tmax] as well as χ2 per degree of freedom (χ2/dof), which is approximately one and indicates that the fits are reasonable. The mass of the lowest three states is stable to some extent for different tmin and as such is a reliable parameter. For the lowest two states, the mass values for the β=2.4 lattice are slightly larger (approximately 100 MeV larger) than those for β=2.8 lattice. This difference might be attributed to the finite lattice spacing effect, and the strange quark mass parameters on the two lattices are not tuned to be exactly the same in the sense of their physical meaning. Combining the results from the two lattices, we can determine that the mass of the lowest 1+ sˉsg state is approximately 2.2 GeV. However, the mass splitting of the lowest two states is approximately 1.4 GeV, which is almost the same as that of the 1+ charmonium-like hybrids, and thus exhibits quark mass independence to some extent.

      tminχ2/dofm1 /GeVm2 /GeVm3 /GeV
      β=2.4tmax=20(sˉsg)
      61.122.232(22)3.56(21)7.5(2.7)
      51.232.228(22)3.61(26)4.9(7)
      41.362.248(13)3.71(11)5.7(5)
      31.382.255(09)3.65(07)5.4(2)
      β=2.8tmax=30(sˉsg)
      91.532.099(16)3.55(07)7.7(7)
      81.412.168(15)3.78(13)5.1(3)
      71.472.100(15)3.40(07)5.8(2)
      61.262.110(13)3.47(06)5.5(1)

      Table 2.  The fitted masses mn of the 1+ states with n=1,2,3 from different time windows [tmin,tmax] and χ2 per degree of freedom (χ2/dof) for the β=2.4 and β=2.8 lattices. All the masses are converted to values in physical units using the lattice spacing as in Table 1.

      Along with the masses, the BS wave functions of the states of Φn(r) can be extracted from the joint fit to C(r,t) as shown in Fig. 1. The radial separation r is converted to the value in physical units, and the wave functions on the two lattices are compatible with each other. The BS wave functions exhibit a clear nodal structure along the r direction: the BS wave function Φ1(r) of the ground state has no radial node, and that of the first excited state (Φ2(r) ere) has one node, whereas that of the third state has two nodes. These nodal structures are very similar to the non-relativistic two-body Schrödinger wave functions in a central potential. It should be noted that r is the spatial separation between the sˉs component and the color chromomagnetic field strength B. The r behavior of the wave functions of the excitations may imply that within the 1+ sˉsg hybrid, the relative movement between the sˉs and the gluonic degrees of freedom can be viewed qualitatively as a two-body system in which r is a physically meaningful dynamic variable. The same data analysis strategy is applied to the 2+ channel. The fitted masses of the lowest three states are listed in Table 3, and the wave functions are shown in Fig. 2. In comparison with the case of 1+, the masses of the 2+ states are slightly higher (100-200 MeV higher for the ground states) than their 1+ counterparts, but the pattern of the spectrum is similar. For the BS wave functions, the r behavior of the 1+ and 2+ states are similar. These observations support the assertion that the 1+ and 2+ states have almost the same inner structure and dynamics, whereas the small mass difference can be attributed to the different coupling between the spin of the sˉs subsystem and the gluonic degrees of freedom. This meets our expectation that 1+ and 2+ states with the nearly degenerate mass can be in the same supermultiplet. Of course, the possibility exists that these 2+ states can be the conventional sˉs mesons since JPC=2+ is permitted for a qˉq system. However, the masses we obtain are much higher than those of the 1D2 sˉs states in the quark model. However, a previous lattice study [31] on charmonium states found that the 2+ cˉcg operator couples almost exclusively to a state with a mass of 4.4 GeV instead of the expected 11D2 charmonium state ηc2, with a mass of approximately 3.8 GeV. If this is also the case for the sˉs states, the hybrid assignment is favorable for the 2+ states we obtain in this work.

      Figure 1.  (color online) The BS wave functions Φn(r) (normalized as Φn(0)=1) of the lowest two 1+ states. r is the spatial separation between the sˉs component and the chromomagnetic operator Bi and is converted to the value in physical units. Open and filled data points are the results for β=2.4 and β=2.8, respectively.

      Figure 2.  (color online) The BS wave functions Φn(r) (normalized as Φn(0)=1) of the lowest two 2+ states. r is the spatial separation between the sˉs component and the chromomagnetic operator Bi and is converted to the value in physical units. Open and filled data points are the results for β=2.4 and β=2.8, respectively.

      tminχ2/dofm1 /GeVm2 /GeVm3 /GeV
      β=2.4tmax=20(sˉsg)
      61.272.416(41)3.65(23)7.2(2.4)
      51.372.406(41)3.60(25)5.5(8)
      41.472.442(26)3.69(19)4.9(3)
      31.562.426(19)3.60(10)5.0(1)
      β=2.8tmax=30(sˉsg)
      91.472.361(23)3.89(09)12.2(2.4)
      81.242.321(34)3.64(18)5.5(5)
      71.372.341(27)3.71(12)6.0(3)
      61.402.359(19)4.00(08)6.2(2)

      Table 3.  The fitted masses mn of the 2+ states with n=1,2,3 for different time windows [tmin,tmax] and χ2 per degree of freedom (χ2/dof) for the β=2.4 and β=2.8 lattices. All the masses are converted to the values in physical units using the lattice spacing as in Table 1.

    • C.   Results of 0+ and 1 states

    • The 0+ and 1 states have conventional quantum numbers for qˉq mesons, and the mesons with these quantum numbers are usually assigned to the n1S0 and n3S1 states in the quark model. Therefore, we start with the analysis of the wall-source correlation function C(r,t) of the sˉs operators with the s and the ˉs field separated by a spatial distance r. We also use the function form of Eq. (5) with N=3 mass terms. The upper bound of the fit window is fixed to tmax=40 and 30, and the lower bound gradually decreases to tmin=5 and 7 for β=2.4 and β=2.8, respectively. The fitted masses of the 1 states are listed in Table 4 and those of the 0+ states are listed in Table 5. As previously indicated, we use the physical mass of ϕ(1020) to set the mass parameters of the strange quark in the fermion action on the two lattices with smaller gauge ensembles. However, the fitted mass of the ground state deviates from the physical mass slightly, which means that the strange quark masses are not tuned with sufficient precision. Therefore, one should consider this slight deviation when examining the data in the tables.

      tminχ2/dofm1 /GeVm2 /GeVm3 /GeV
      β=2.4tmax=40(sˉs)
      80.841.013(1)1.753(77)2.16(22)
      70.811.013(1)1.787(93)2.08(19)
      60.831.014(1)1.732(47)2.13(12)
      50.921.015(1)1.709(36)2.11(08)
      β=2.4tmax=45(sˉsg)
      161.401.011(1)1.72(12)
      151.421.011(1)1.66(10)
      141.421.009(1)1.77(08)
      131.671.007(1)1.83(06)
      β=2.8tmax=30(sˉs)
      101.931.001(5)1.634(114)2.17(29)
      91.971.003(4)1.633(86)2.02(18)
      82.110.999(3)1.665(81)2.30(20)
      72.390.998(3)1.668(52)2.34(17)
      β=2.8tmax=45(sˉsg)
      191.171.006(4)1.51(6)
      181.391.003(3)1.55(6)
      171.501.005(3)1.51(5)
      161.450.998(2)1.58(4)

      Table 4.  The fitted masses mn of the 1 states for two different types of operators and different time windows [tmin,tmax] and χ2 per degree of freedom (χ2/dof) for the β=2.4 and β=2.8 lattices. All the masses are converted to values in physical units using the lattice spacing as in Table 1.

      tminχ2/dofm1 /GeVm2 /GeVm3 /GeV
      β=2.4tmax=40(sˉs)
      80.700.7012(2)1.690(34)2.21(17)
      70.690.7010(2)1.698(39)2.12(14)
      60.730.7010(2)1.699(30)2.12(10)
      50.890.7014(2)1.669(22)2.10(07)
      β=2.4tmax=45(sˉsg)
      161.300.7008(3)1.711(73)
      151.400.7007(3)1.680(56)
      141.420.7007(3)1.672(46)
      131.560.7009(3)1.659(36)
      β=2.8tmax=30(sˉs)
      102.590.6483(8)1.736(62)2.84(34)
      92.510.6505(8)1.703(76)2.33(25)
      82.370.6512(8)1.679(66)2.23(20)
      72.350.6516(7)1.620(30)2.46(11)
      β=2.8tmax=45(sˉsg)
      191.090.6508(10)1.621(85)
      181.030.6510(10)1.557(67)
      171.040.6490(06)1.711(50)
      161.060.6491(05)1.740(39)

      Table 5.  The fitted masses mn of the 0+ states for two different types of operators and different time windows [tmin,tmax] and χ2 per degree of freedom (χ2/dof) for the β=2.4 and β=2.8 lattices. All the masses are converted to values in physical units using the lattice spacing as in Table 1.

      In the 1 channel, the masses of the ground state and the first excited state are approximately 1 GeV and 1.7 GeV, respectively, which are compatible with those of ϕ(1020) and ϕ(1680). For the β=2.4 lattice, the fitted mass m3 of the third state is also stable with respect to tmin, and the value is approximately 2.1 GeV, which is close to the expected mass of the 33S1 state for the quark model. The m3 on the β=2.8 lattice is also in this mass range but fluctuates more strongly with tmin.

      In the 0+ channel, the ground state mass can be precisely determined with m10.701 GeV at β=2.4 and 0.651 GeV at β=2.8. Since the sˉs pseudoscalar meson (labeled as ηs) is not a physical state, we cannot directly compare our result to the physical value . A previous calculation has been performed for the Nf=2+1 full-QCD lattice formalism, which gives the prediction mηs=0.686(4) GeV [32], lying between our values for the two lattices. The deviation is small and can be attributed to our less precise tuning of the strange quark mass parameter and the other systematic uncertainties. The mass of the first excited state is approximately 1.6-1.7 GeV, which is almost degenerate with that of the first excited 1 state. There is no physical correspondence of this state at present, but it can be compared to the pseudoscalar η(1295)/η(1475). However, it should be noted that this state is a pure sˉs state, which results in a higher mass. As indicated in the previous section, the BS wave functions Φn(r) of the 0+ and 1 states can be derived simultaneously with the masses, as shown in Figs. 3 and 4. It should be noted that r represents the separation between the s and ˉs field. For the ground state and the first excited state in each channel, the Φn(r) values exhibit the expectation of the quark model in that the wave function of the ground state has no radial node whereas that of the first excited state has one node. Therefore, the two states can be assigned to be the 1S and 2S states of a non-relativistic sˉs system. The behavior of the wave function of the third state is strange in that it has no radial nodes even though it has two inflection points. Since we only use three mass terms to fit the correlation functions, the third state may have substantial contamination from higher states, which may result in this phenomenon. As such, the results of the third state are not seriously considered.

      Figure 3.  (color online) The BS wave functions Φn(r) (normalized as Φn(0)=1) of the lowest two 0+ states. r is the spatial separation between the quark fields s and ˉs and is converted to the value in physical units. Open and filled data points are the results for β=2.4 and β=2.8, respectively.

      Figure 4.  (color online) BS wave functions Φn(r) (normalized as Φn(0)=1) of the lowest two 1 states. r is the spatial separation between the quark fields s and ˉs and is converted to the value in physical units. Open and filled data points are the results for β=2.4 and β=2.8, respectively.

      We also use the sˉsg-type operators (in Eq. (1)) to explore the properties of the 0+ and 1 states. We use Eq. (5) with N=2 mass terms to fit the correlation functions in large time ranges (tmax=45 for both lattices). The masses are listed in Table 4 and Table 5, in which it is apparent that they are consistent with those from the sˉs operators. Figures 5 and 6 show the BS wave functions with r as the spatial separation between the sˉs components and the chromomagnetic operator. The results for the two lattices are compatible with each other. It is interesting to see that for the 0+ and 1 channels, such wave functions of the ground state (1S state) and the first excited state (2S state) lie almost upon each other and there are no significant differences. In these two channels, given the masses that are compatible with the masses of the states obtained using the sˉs operators, the ground states and the first excited states can be assigned to the 1S and 2S sˉs mesons. As such the similarity of the wave functions of 1S and 2S states with respect to the distance r between the gluonic component and the sˉs component can be interpreted as follows: the s (or ˉs) field along with the chromomagnetic field can be viewed as a dressed s (or ˉs) field in the fundamental representation of the color SU(3) group, which annihilates the s (anti)quark of the sˉs meson. In this sense, r represents the spatial size of the dressed quark field and the r fall-off does not have a dynamic significance. This is in contrast to the wave functions of the 1+ and 2+ states whose nodal structures imply that r is a dynamic variable for hybrids.

      Figure 5.  (color online) The BS wave functions Φn(r) (normalized as Φn(0)=1) of the lowest two 1 states. r is the spatial separation between the sˉs component and the chromomagnetic operator Bi and is converted to the value in physical units. Open and filled data points are the result for β=2.4 and β=2.8, respectively. Since the two states can be assigned to be the 1S and 2S sˉs mesons, the similar r-behavior of their BS wave functions imply that this r is not a typical dynamical variable for the sˉs states.

      Figure 6.  (color online) Similar to Fig. 5 but for 0+ states.

    III.   DISCUSSION
    • We discuss the obtained results in this section. The BS wave functions of the 1+ and 2+ states show the typical behaviors of non-relativistic two-body Schrödinger wave functions with a central potential, in terms of the correspondence of the spectrum and the nodal structure of the wave functions. We emphasize that the variable r is the spatial distance between the sˉs component and the chromomagnetic field strength of the operators. Since 1+ is an exotic quantum number for qˉq mesons, the states with this quantum number must be a hybrid meson with additional gluonic degrees of freedom. The similarity of the spectrum and the wave functions of the 2+ states to those of the 1+ states indicates that they are also hybrid states. In this sense, the wave functions imply that r can be a meaningful dynamical variable for sˉsg hybrid mesons. In a previous lattice study on cˉcg hybrids [27], the same behaviors of the wave functions and the spectrum pattern were observed for the (0,1,2)+ and 1 supermultiplets, based on which a "color halo" concept has been proposed, which states that a hybrid meson can be viewed as a relatively compact color octet qˉq pair surrounded by color octet gluonic degrees of freedom, such that the wave functions depict the relative motion between the qˉq pair and the gluonic excitation. In a non-relativistic model, the binding mechanism can be the potential between two effective color octet charges. Previous lattice studies [33] revealed that the potential of two static color charges has the feature of a Casimir scaling

      VD(r)=VD,0CDαr+σDr,

      (6)

      where D is the color SU(3) representation of the charge with CD as the eigenvalue of the second order Casimir operator of the color SU(3), VD,0 is the potential constant, α is the coefficient of the Coulomb part, and σD is the string tension, which is related to the conventional string tension σ between a static quark and antiquark pair by σD=34CDσ. For the octet charge in this work, this relation is σD=9/4σ, which means that the interaction between the color octet objects is stronger than that between the color triplet ones. This explains the observation that the mass splitting (approximately 1.2 GeV) between the ground state and the first excited hybrid state is larger than the 1S2S mass splitting of the qˉq states (approximately 0.6 GeV).

      The "color halo" idea is conceptually different from the flux-tube framework of hybrids in the market [12-14, 34, 35], whereby the quark and anti-quark are bound by an effective potential induced by the excitation of gluonic degrees of freedom. In the leading Born-Oppenheimer approximation, the QˉQ of a heavy quarkonium-like hybrid can be viewed as static color sources; the excited gluonic degrees of freedom are distributed along the QˉQ axis and obey the cylinder symmetry, whose effect can be treated as an excited static potential denoted by Λϵη, where Λ=0,1,2, is the projected total angular momentum of the gluons with respect to the QˉQ axis and is labeled as Σ,Π,Δ for Λ=0,1,2, and so on, η represents the combined parity (P) and the charge conjugate (C) of gluon excitations with η=g,u for PC=±, respectively, and ϵ is the P parity of the glue state. Therefore, the quantum number of a QˉQ state with this potential is

      P=ϵ(1)L+Λ+1,C=ηϵ(1)L+S+Λ,

      (7)

      where ˆL=ˆLQˉQ+ˆJgwith ˆLQˉQ being the orbital angular momentum of QˉQ with respect to the midpoint of the QˉQ axis and ˆJg being the total angular momentum of the gluons. The ground Σ+g potential has Λ=0,ϵ=+ and η=+ are the conventional static potential of QˉQ of the Cornell type, and the P and C quantum numbers reproduce the conventional quantum number. The 1 and (0,1,2)+ hybrid supermultiplet is associated with the Π+u(L=1) potential, such that the radial Shrödinger equation is

      d2dr2u(r)+2μ[EVeff(r)]u(r)=0,

      (8)

      where r is the distance between Q and ˉQ, μ is the reduced mass of the QˉQ pair, and u(r) is related to the radial wave function ϕ(r) by u(r)=rϕ(r). The effective potential Veff is

      Veff=VQˉQ(r)+ˆL2QˉQ2μr2

      (9)

      with ˆL2QˉQ=L(L+1)2Λ2+ˆJ2g and ˆJ2g=2. Obviously, the eigenvalues of E are independent of the total spin S of the QˉQ pair. One can use the lattice results to determine VQˉQ and then solve the preceding equation to obtain the masses of the hybrids. We do not wish to delve into such details in this work but mention that ϕn(r) behaves as a P-wave wave function in a central potential, and the mass splitting of the ground state and the first radial excited state is only a few hundred MeV for the 1 and (0,1,2)+ hybrid states [14]. Even though the preceding deduction is based on the heavy quarkonium-like hybrids, this concept has been also applied to the phenomenological studies of strangeonium hybrids [34].

      In contrast to the flux-tube concept, we observe that for (1,2)+ strangeonium hybrids, the mass splitting of the ground and the first excited states is approximately 1.2-1.4 GeV, which is much larger than the prediction of the flux-tube model, and the nodal structure is present with respect to the spatial distance between the sˉs and the chromomagnetic field strength. A similar phenomenon also appears for charmonium-like hybrids without a clear quark mass dependence. It should be emphasized that even though the interpretation of the wave functions is debatable, the pattern of the spectrum should be reliable and model-independent since it is derived directly from the lattice QCD calculation.

      For the 0+ and 1 channels, we obtain consistent results for the masses of the ground states and the first excited state using the sˉs and sˉsg type operators. Since we use the physical mass of the ϕ(1020) meson to set the strange quark mass parameters, it is natural to almost reproduce the physical value of the mass of the vector ground state. The ground state mass of the pseudoscalar is approximately 650700 MeV, which is compatible with the previous lattice result of ηs. The masses of the first excited states in both channels are closely degenerate at 1.7 GeV, and the mass of the first excited sˉs vector meson is in agreement with that of ϕ(1680). However, the BS wave functions in both channels, defined based on the dependence of spatial distance between the s and ˉs quark field, show the expected radially nodal behavior of the non-relativistic sˉs two-body system. Therefore, the ground and the first excited states can be assigned to be the 1S and 2S sˉs mesons, respectively. We also obtain some information on the third state based on the sˉs type operator in each channel, whose mass is approximately 2.1-2.3 GeV. For the vector channel, this mass value is close to the mass of ϕ(2170). However, since we only use three mass terms for the data fitting, the third state may have substantial contaminations from higher states, and the result is not reliable. When we use the sˉsg operator to study these two channels, we can only obtain information on the lowest two states. At present, there is no definitive conclusion to whether there is a 1 and (0,1,2)+ supermultiplet of the strangeonium hybrids.

      Finally, we present some arguments related to ϕ(2170). Its mass is in the range of the 33S1 and 23D1 sˉs predicted by the quark model. If a 1 and (0,1,2)+ sˉsg hybrid multiplet exists with nearly degenerate masses at approximately 2.12.3 GeV, ϕ(2170) can also be a candidate for the 1 member. However, the assignment of its characteristic is still an open question. Till now, ϕ(2170) has been observed in many final states including ϕ(1020), such as ϕf0(980), ϕππ, ϕη, and ϕη. In the K+Kππ and K+KK+K final states [36-38], there are also sizable components including ϕ(1020). If ϕ(2170) is a candidate for the 1 sˉsg hybrid, this decay pattern can be understood based on the color halo concept of the hybrids: the binding between the color octet sˉs and the gluonic degrees of freedom can easily break up such that the sˉs component is neutralized to ϕ(1020) and the gluons are hadronized to light hadrons, which are in the flavor singlet. Furthermore, in contrast to the hadronic transition of conventional excited strangeonium states, these decays are less OZI suppressed due to the existing gluons within the strangeonium hybrids. Recently, the BESIII Collaboration reported the observation of ϕ(2170) in the process e+eηϕ with the resonance parameters MR=2177.5±4.8(stat)±19.5(syst) MeV and ΓR=149.0±15.6(stat)±8.9(syst) MeV, and Br(ϕ(2170)ηϕ)Γe+e is measured to be 7.1±0.7(stat)±0.7(syst) eV [39]. Combining the result of Br(ϕ(2170)ηϕ)Γe+e=1.7±0.7(stat)±1.3(syst) eV, one has

      Br(ϕ(2170)ηϕ)Γe+eBr(ϕ(2170)ηϕ)Γe+e=0.23±0.10(stat)±0.18(syst).

      (10)

      This ratio is much larger than the predictions of phenomenological studies based on the flux tube model or the constituent gluon model of hybrids, with the mechanism whereby the flux tube or the constituent gluon breaks up into a light qˉq pair that reorganizes into two mesons with the original constituent sˉs. However, this ratio can be explained directly from the flavor octet-singlet mixing and the kinetics. If ϕ(2170) is a sˉsg hybrid in the 'color halo' picture, then the decay ϕ(2170)ϕη(η) can take place as follows: a gluon is emitted by a constituent strange quark (or antiquark) and the original gluon(s) couple to the flavor singlet component of the η(η) meson. If ϕ(2170) is a higher excited sˉs meson, then the η(η) is generated by two gluons emitted by the sˉs pair. It should be noted that this process can be enhanced by the QCD axial anomaly. Since the decay dynamics is expected to be the same for the ϕη and ϕη decay modes, the ratio of the partial widths can be attributed to the ηη mixing and the kinetic factors

      Γ(ϕ(2170)ϕη)Γ(ϕ(2170)ϕη)=tan2θ(kηkη)3,

      (11)

      where θ is the flavor octet-singlet mixing angle of the ηη system, and kη() is the magnitude of the decay momentum. If we take the physical masses mη=547 MeV and mη=958 MeV and the mixing angle θas varying between 10 and 20, this ratio is estimated to be between 0.14 and 0.58 and compatible with the experimental value (the mixing angle θ is derived to be approximately |θ|13 using the central value 0.23). As such, for ϕ(2170), this ratio may not be an ideal criterion to distinguish a hybrid assignment from a conventional sˉs meson. Apart from the decay modes involving ϕ(1020), BESIII have reported KˉK decay modes of resonances X observed in the procesess e+eXKˉK with a peak of X at approximately 2.2 GeV [38, 40]. For the characteristics of ϕ(2170) to be revealed, all the observed decay modes should be considered jointly in scrutinized theoretical discussions.

    IV.   SUMMARY
    • The strangeonium-like hybrids were investigated based on lattice QCD in the quenched approximation. Two anisotropic lattices with different lattice spacings were used to examine finite as effects. We constructed spatially extended sˉsg operators with the sˉs component separated from the chromomagnetic field strength operator by a spatial distance r. We investigated the 1 and (0,1,2)+ channels and calculated the corresponding correlation functions based on these operators in the Coulomb gauge. The ground state mass of the 1+ states was determined to be 2.1-2.2 GeV and that of the 2+ states was approximately 200 MeV higher. These results are consistent with previous lattice calculations and phenomenological studies. The masses of the first excited state were approximately 3.6 GeV in these two channels, such that the mass splitting of the first excited state and the ground state was approximately 1.2-1.4 GeV. This is much higher than the predictions obtained based on the flux-tube model, which is only a few hundred MeV. The BS wave functions of these states, defined by the matrix elements of the aforementioned operators between the vacuum and the states, were extracted and exhibit clear nodal structures in the r direction. This indicates that r is a meaningful dynamical variable reflecting the relative motion of the center-of-mass of the sˉs against the gluonic degrees of freedom. Both the spectrum and the wave functions of these sˉsg states have features similar to those of their cˉcg counterparts and are consistent with the "color halo" concept of the hybrids in that the color octet qˉq pair is surrounded by gluons.

      In the 0+ and 1 channels, we used both the spatially extended sˉs and sˉsg operators to perform the calculations. The ground state mass of the vector sˉs meson almost reproduced the mass of ϕ(1020) (it should be noted that we use the mass of the ϕ(1020) to set the strange quark mass parameters), and the ground state mass of the pseudoscalar was 650-700 MeV, which is in agreement with the ηs mass determined by previous lattice calculations. In both channels, the masses of the first excited states were almost degenerate at approximately 1.7 GeV and compatible with the mass of ϕ(1680). The BS wave functions of these states with respect to the distance between s and ˉs were qualitatively similar to the nonrelativistic wave function of a two-body system in that the BS wave function of the first excited state had a radial node. Therefore, the first excited state can be a 2S sˉs meson. In contrast, in each channel, the BS wave functions with respect to the spatial distance of the octet sˉs and the gluonic degrees of freedom had similar profiles for the ground and the first excited state, which meant that this distance was less significant for sˉs mesons.

      We did not obtain reliable results for the 3S sˉs mesons and the possible vector sˉsg hybrids. Therefore, we are unable to present a convincing explanation for ϕ(2170). Since the mass of ϕ(2170) is compatible with the quark model prediction of the 3S sˉs meson and the predicted mass of the lowest sˉsg hybrids, both assignments of ϕ(2170) are possible. We argue that if ϕ(2170) is either the 3S sˉs meson or a vector sˉsg hybrid within the "color halo" picture discussed in the preceding section, the ratio Γ(ϕη)/Γ(ϕη)=0.23±0.10(stat)±0.18(syst) can be understood based on the hadronic transition of a strangeoium-like meson in addition to ηη mixing. Nevertheless, the nature of ϕ(2170) is still an open question to be investigated based on further experimental and theoretical studies.

Reference (40)

目录

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return