PHYSICAL REVIEW D 102, 116014 (2020)

Discovery potential for the LHCb fully charm Xð6900Þ state via pp¯ annihilation reaction

† ‡ Xiao-Yun Wang,1,* Qing-Yong Lin,2, Hao Xu,3 Ya-Ping Xie,4 Yin Huang,5 and Xurong Chen4,6,7, 1Department of physics, Lanzhou University of Technology, Lanzhou 730050, China 2Department of Physics, Jimei University, Xiamen 361021, China 3Department of Applied Physics, School of Science, Northwestern Polytechnical University, Xi’an 710129, China 4Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China 5School of Physical Science and Technology, Southwest Jiaotong University, Chendu 610031, China 6University of Chinese Academy of Sciences, Beijing 100049, China 7Guangdong Provincial Key Laboratory of Nuclear Science, Institute of Quantum , South China Normal University, Guangzhou 510006, China

(Received 2 August 2020; accepted 16 November 2020; published 18 December 2020)

Inspired by the observation of the fully-charm tetraquark Xð6900Þ state at LHCb, the production of Xð6900Þ in pp¯ → J=ψJ=ψ reaction is studied within an effective Lagrangian approach and Breit-Wigner formula. The numerical results show that the cross section of Xð6900Þ at the c.m. energy of 6.9 GeV is much larger than that from the background contribution. Moreover, we estimate dozens of signal events can be detected by the D0 experiment, which indicates that searching for the Xð6900Þ via - scattering may be a very important and promising way. Therefore, related experiments are suggested to be carried out.

DOI: 10.1103/PhysRevD.102.116014

I. INTRODUCTION Γ ¼ 168 33 69 MeV; ð2Þ In recent decades, more and more exotic states have been observed [1]. These states are not respectively. Actually, based on various theoretical models, only conducive to the development of the hadron spectrum a lot of studies have been conducted on fully-heavy but also provide an important opportunity for us to better tetraquark states [5–32], which are important to reveal the understand the multiquark states and strong interactions structure and properties of the fully-heavy tetraquark states. [2,3]. In the past years, the experimental explorations focus After its observation, Xð6900Þ has inspired extensive studies on the multiquark states composed of light and heavy on its underlying properties, where it is interpreted as a P- , while studies for the states composed entirely of wave tetraquark state in a nonrelativistic model [27], heavy quarks are still very limited [1–3]. the first radial excitation of ccc¯ c¯ in an extended relativized Very recently, the LHCb experiment reported a narrow [28]. Besides, in the QCD sum rule framework, structure, named Xð6900Þ, in the J=ψ pair invariant mass the narrow structure Xð6900Þ can be interpreted as a second spectrum with a significance level of more than 5σ [4]. radial excited S-wave tetraquark state [29] with JPC ¼ 0−þ Based on the simple model with interference, the mass and or a P-wave state [30] with 1−þ, respectively. In width of the Xð6900Þ resonance are measured to be Refs. [31,32], the theoretical results indicate that there may exist the resonance states, with masses ranging between M ¼ 6886 11 11 MeV; ð1Þ 6.3 to 7.4 GeV, and the quantum numbers JP ¼ 0þ; 1þ, and 2þ. In Ref. [33], the spin parity of the state formed by the

* two-vector system was discussed, and a possible method for Corresponding author. – [email protected] determining its was given. In Refs. [34 † [email protected] 36], the spectrum or mass of fully-heavy was ‡ [email protected] discussed in model with different approaches. In addition to the study of the mass spectrum and the Published by the American Physical Society under the terms of decay behavior of Xð6900Þ, studying the production of the Creative Commons Attribution 4.0 International license. ð6900Þ Further distribution of this work must maintain attribution to X through different scattering reactions is also very the author(s) and the published article’s title, journal citation, important and will help us to understand its nature as and DOI. Funded by SCOAP3. genuine states. For example, in Ref. [37], the production of

2470-0010=2020=102(11)=116014(5) 116014-1 Published by the American Physical Society WANG, LIN, XU, XIE, HUANG, and CHEN PHYS. REV. D 102, 116014 (2020) the ground ccc¯ c¯ with JPC ¼ 0þþ; 2þþ in pp collisions was The production cross section of Xð6900Þ structure calculated, and the upper limit of the cross section is about relative to that of all J=ψ pair, times the branching fraction 40 fb for the four channel. In addition to pp BrðX → ψψÞ, R, is determined as ½2.6 0.6 ðstatÞ collisions, determining the tetraquark states via the annihi- 0.8 ðsystÞ% by the LHCb experiment [4]. Therefore, in lation of positive and negative is usually an this work, we roughly take the value of the branching important and effective way [38–44]. In this work, the ratio BrðX → ψψÞ ≃ 2%. discovery potential for the tetraquark state Xð6900Þ via pp¯ Having fixed the branch ratio BrðX → ψψÞ, we turn to annihilation reaction will be investigated. One will estimate the value of BrðX → pp¯ Þ, which is indispensable in Eq. (3). the cross section of Xð6900Þ production via pp¯ → J=ψJ=ψ However, it has never been determined in experiments. reaction and analyze the corresponding background con- Meanwhile, we notice in Refs. [39,40] that the branching tribution. The theoretical results obtained will be an impor- fractions of BrðYð4260Þ → pp¯ Þ and Brðψð4040Þ → pp¯ Þ tant theoretical basis for the pp¯ annihilation experiment. were estimated using the branching ratio of J=ψ state, This paper is organized as follows. After the Introd- multiplying by the ratio of the width of the Yð4260Þ or uction, we present the formalism for the production of ψð4040Þ to the width of J=ψ. Thus, in this work, we will Xð6900Þ in Sec. II. The numerical results of the Xð6900Þ adopt the same method as in Refs. [39,40] to naively production follow in Sec. III. Finally, the paper ends with a determine the branching ratio of X decaying into pp¯ .Due brief summary. to the lack of deep understanding of the fully-charm tetraquark Xð6900Þ, its inner structure and quantum num- bers are still unknown. Different JPC assignments were II. FORMALISM assumed to evaluate the mass of a full-charm tetraquark, as A. Xð6900Þ production in pp¯ annihilation discussed in Sec. I. To carry out the estimation of the cross section of the s channel, we assume that the quantum The Feynman diagram of the pp¯ → J=ψJ=ψ reaction via number JPC of Xð6900Þ is 0þþ or 2þþ, the same as the the s channel Xð6900Þ exchange is depicted in Fig. 1(a). assumption in Ref. [37]. In the known charmonium family, For convenience, Xð6900Þ and J=ψ will be abbreviated as þþ þþ the spin parities of χ 0 and χ 2 are 0 and 2 , X and ψ, respectively. c c respectively. Therefore, we plan to simply replace J=ψ For the production of the resonance Xð6900Þ in the s with χ 0 and χ 2 to estimate the branching ratio of X → pp¯ , channel of pp¯ → ψψ reaction, the cross section can be c c namely, calculated by the standard Breit-Wigner formula [1] Γ χc0 BrðX → pp¯ Þ ≃ Brðχc0 → pp¯ Þ × ; 2J þ1 4π Γ2 BrðX → pp¯ ÞBrðX → ψψÞ Γ σ ¼ ; X X ð2 þ1Þð2 þ1Þ 2 4 ð − Þ2 þΓ2 4 PC ¼ 0þþ ð Þ S1 S2 kin W W0 = for X with J 4 ð Þ 3 and

Γχ where W is the c.m. energy, J is the spin of the resonance ð → ¯ Þ ≃ ðχ → ¯ Þ c2 Br X pp Br c2 pp × Γ ; Xð6900Þ, and S1 and S2 are the spins of the initial X PC þþ antiproton and proton, respectively. Moreover, kin is the for X with J ¼ 2 ; ð5Þ c.m. momentum in the initial state, W0 is the c.m. energy Γχ ¼ 10 8 Γχ ¼ 1 97 Γ ¼ corresponding to the mass threshold of Xð6900Þ, and Γ is where c0 . MeV, c2 . MeV, and X 168 χ χ ð6900Þ the full width of Xð6900Þ. MeV are the total width of c0, c2, and X states, respectively. By taking the branching ratio −4 Brðχc0 → pp¯ Þ¼2.24 × 10 and Brðχc2 → pp¯ Þ¼ 7.33 × 10−5, one gets BrðX → pp¯ Þ¼1.44 × 10−5 for X with JPC ¼ 0þþ and BrðX → pp¯ Þ¼8.6 × 10−7 for X with JPC ¼ 2þþ. Here, it must be noted that if we use the (a) branching ratio of J=ψ for estimation BrðX → pp¯ Þ¼ 1.17 × 10−6 is obtained, which is just between the value estimated by the branching ratio of χc0 and χc2.

B. Background analysis ¯ → ψψ (b) (c) Figures 1(b) and 1(c) present the pp reaction with the t and u channels by exchanging a , which FIG. 1. Feynman diagrams for the reaction pp¯ →Xð6900Þ→ can be considered as the main background contributions for ψψ. the production of Xð6900Þ. By employing the effective

116014-2 DISCOVERY POTENTIAL FOR THE LHCB FULLY-CHARM … PHYS. REV. D 102, 116014 (2020)

Lagrangian approach, the cross section of the pp¯ → ψψ III. NUMERICAL RESULTS reaction can be calculated. After the above preparations, we calculated the total The Lagrangian density for the vertex of ψNN is written cross section for the reaction pp¯ → ψψ from threshold to as [41] 12 GeV of the center of mass energy, as depicted in Fig. 2. ¯ μ As can be seen from the Fig. 2, the cross section from the LψNN ¼ −gψNNNγ Nψ μ; ð6Þ Xð6900Þ contribution has a distinct peak near the center of where ψ and N denote the fields of J=ψ and the nucleon, mass energy of 6.9 GeV. Whatever JPC assignment (0þþ or respectively. 2þþ) is chosen, the cross section of X production through The values of coupling constant gψNN can be derived the s channel can reach dozens of picobarns, which is much from the corresponding decay width, larger than the cross section from the background contri- qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi bution. In addition, we also calculated the background 2 2 cross section without the form factor (FF). Although the mψ − 4m Γ ¼ð Þ2 N ð 2 þ 2 2 Þ background cross section without FF is an order of ψ→pp¯ gψNN 2 mψ mN : 12πmψ magnitude higher than that with the FF added, it is still much lower than that of the s channel. The main reason for ≃ 1 6 10−3 Thus, we get gψNN . × , which is calculated by the very small background is that the gψNN coupling the measured branching fractions and total widths of J=ψ constant is very small. (mψ ¼ 3096.916 MeV and Γψ→pp¯ ≃ 0.197 keV) [1]. The above results indicate that the best energy window Based on the above Lagrangians, the scattering ampli- for searching for the fully-charm tetraquark state Xð6900Þ tude for the reactions pp¯ → ψψ can be constructed as via the pp¯ → ψψ process is near the c.m. energy 6.9 GeV, in which the signal can be clearly distinguished from ν μ −iMpp¯ →ψψ ¼ ϵψ ðp2Þu¯ðp1ÞAμνuðk1Þϵψ ðk2Þ; ð7Þ background. The D0 [44] and the forthcoming PANDA [39] experiments are ideal platforms to study new hadronic where u is the Dirac spinor of the nucleon and ϵγ is the states via pp¯ reactions, of which the designed c.m. energy polarization vector of the . is below the double J=ψ threshold for the latter one. By ð6900Þ The reduced amplitude Aμν for the t and u channel taking the cross section of X production calculated background reads above, one finds that the number of events of Xð6900Þ can reach as many as dozens if taking an integrated luminosity =q þ m of 10.4 fb−1 collected with the D0 detector, when a 50% At ¼ð Þγ t N γ F 2ð 2Þ ð Þ μν gψNNgψNN ν 2 − 2 μ t qt ; 8 selection efficiency is adopted. Considering that the back- qt m N ground is very clean, it should be possible for Xð6900Þ to þ u =qu mN 2 2 Aμν ¼ðgψ gψ Þγμ γνF ðq Þ: ð Þ NN NN 2 − 2 u u 9 qu mN

2 For two ψNN vertices, a general form factor F tðqt Þ¼ F ð 2Þ¼ðΛ2 − 2 Þ ðΛ2 − 2 Þ u qu t=u mN = t=u qN is taken into account [41,43]. In Refs. [41,43], it can be found that the cross section of pp¯ → J=ψπ reaction with nucleon exchange was consistent with the E760 and E835 data by taking Λt=u ¼ 1.9 and 3.0 GeV, respectively. In the spirit of estimating the upper limit of background contribution, in this work, we take Λt ¼ Λu ¼ 3 GeV. With the preparation in the previous section, the cross section of the reaction pp¯ → ψψ can be calculated. The differential cross section in the c.m. frame is written as   ⃗c:m: X dσ 1 jk2 j 1 ¼ jMj2 : ð10Þ θ 32π ⃗c:m: 4 d cos s jk1 j λ FIG. 2. The total cross section for the pp¯ → ψψ reaction. The 2 red solid and blue dashed lines are for the cross section of Here, s ¼ðk1 þ p1Þ , and θ denotes the angle of the Xð6900Þ production when the spin parity is 0 þþ and 2 þþ, J=ψ p¯ outgoing relative to beam direction in the respectively. Moreover, the green dot-dashed and violet dotted ⃗c:m: ⃗c:m: c.m. frame. k1 and k2 are the 3-momenta of the initial lines are for the contributions from background with and without photon beam and final J=ψ meson, respectively. the form factor, respectively.

116014-3 WANG, LIN, XU, XIE, HUANG, and CHEN PHYS. REV. D 102, 116014 (2020) be observed by the D0 experiment with a high confidence Xð6900Þ through the pp¯ → ψψ reaction. The results level. It should be noted that the quantum number of indicate that it is accessible to search for Xð6900Þ via Xð6900Þ has not been identified yet. However, to evaluate pp¯ annihilation. the production of Xð6900Þ, the JPC assignments of 0þþ and It should be noted that the pp¯ will be annihilated into þþ 2 are adopted in this work. The results indicate that the first and then those gluons will be converted into the cross sections are at the same order of magnitude for both di-J=ψ final state. The production of the fully-charm cases. It means that it may be an effective way to find tetraquarks at LHC was discussed, and the relative cross Xð6900Þ through the pp¯ annihilation. section for the process gg → J=ψJ=ψ was estimated [45,46]. The cross sections performed are around 10 pb, IV. CONCLUSIONS which have the same order of magnitude as the results In this work, the Xð6900Þ production in the antiproton exhibited here. proton reaction is investigated by employing the effective field theory and the Breit-Wigner formula. We find that ACKNOWLEDGMENTS pp¯ → ψψ can be a suitable process to study the production of Xð6900Þ. One finds that the cross section of X The authors would like to acknowledge Doctor Liupan production can reach dozens of picobarns at the best An for useful discussions. This project is supported by the energy window W ¼ 6.9 GeV, which is much larger than National Natural Science Foundation of China (Grants that of the background. It means that the signal can be No. 12065014 and No. 11705076) and by the Strategic clearly distinguished from the background. Furthermore, Priority Research Program of Chinese Academy of according to our estimation, dozens of Xð6900Þ can be Sciences, Grant No. XDB34030301. This work is partly detected with a data sample of 10.4 fb−1 collected by the supported by HongLiu Support Funds for Excellent Youth D0 detector, which indicates that it is feasible to find Talents of Lanzhou University of Technology.

[1] M. Tanabashi et al. ( Data Group), Review of [13] M. Karliner, S. Nussinov, and J. L. Rosner, QQQ¯ Q¯ states: , Phys. Rev. D 98, 030001 (2018). Masses, production, and decays, Phys. Rev. D 95, 034011 [2] Y. R. Liu, H. X. Chen, W. Chen, X. Liu, and S. L. Zhu, (2017). and tetraquark states, Prog. Part. Nucl. Phys. [14] W. Chen, H. X. Chen, X. Liu, T. G. Steele, and S. L. Zhu, 107, 237 (2019). Hunting for exotic doubly hidden-charm/bottom tetraquark [3] F. K. Guo, X. H. Liu, and S. Sakai, Threshold cusps and states, Phys. Lett. B 773, 247 (2017). triangle singularities in hadronic reactions, Prog. Part. Nucl. [15] Y. Bai, S. Lu, and J. Osborne, Beauty-full tetraquarks, Phys. Phys. 112, 103757 (2020). Lett. B 798, 134930 (2019). [4] R. Aaij et al. (LHCb Collaboration), Observation of [16] Z. G. Wang, Analysis of the QQQ¯ Q¯ tetraquark states with structure in the J=ψ-pair mass spectrum, Sci. Bull. 65, QCD sum rules, Eur. Phys. J. C 77, 432 (2017). 1983 (2020). [17] V.R. Debastiani and F. S. Navarra, A non-relativistic model [5] K. T. Chao, The ðccÞ − ðc¯ c¯Þ (Diquark-Antidiquark) states for the ½cc½c¯ c¯ tetraquark, Chin. Phys. C 43, 013105 (2019). in eþe− annihilation, Z. Phys. C 7, 317 (1981). [18] M. N. Anwar, J. Ferretti, F. K. Guo, E. Santopinto, and B. S. [6] Y. Iwasaki, A possible model for new resonances-exotics Zou, Spectroscopy and decays of the fully-heavy tetra- and hidden charm, Prog. Theor. Phys. 54, 492 (1975). quarks, Eur. Phys. J. C 78, 647 (2018). [7] J. P. Ader, J. M. Richard, and P. Taxil, Do narrow heavy [19] A. Esposito and A. D. Polosa, A bbb¯ b¯ di-bottomonium at multi-quark states exist?, Phys. Rev. D 25, 2370 (1982). the LHC?, Eur. Phys. J. C 78, 782 (2018). [8] L. Heller and J. A. Tjon, On bound states of heavy Q2Q¯ 2 [20] J. Wu, Y. R. Liu, K. Chen, X. Liu, and S. L. Zhu, Heavy- systems, Phys. Rev. D 32, 755 (1985). flavored tetraquark states with the QQQ¯ Q¯ configuration, [9] A. M. Badalian, B. L. Ioffe, and A. V. Smilga, Four quark Phys. Rev. D 97, 094015 (2018). states in the heavy quark system, Nucl. Phys. B281,85 [21] C. Hughes, E. Eichten, and C. T. H. Davies, Searching for (1987). beauty-fully bound tetraquarks using lattice nonrelativistic [10] S. Zouzou, B. Silvestre-Brac, C. Gignoux, and J. M. Richard, QCD, Phys. Rev. D 97, 054505 (2018). Four quark bound states, Z. Phys. C 30, 457 (1986). [22] J. M. Richard, A. Valcarce, and J. Vijande, Few-body quark [11] R. J. Lloyd and J. P. Vary, All charm tetraquarks, Phys. Rev. dynamics for doubly heavy and tetraquarks, Phys. D 70, 014009 (2004). Rev. C 97, 035211 (2018). [12] N. Barnea, J. Vijande, and A. Valcarce, Four-quark spec- [23] G. J. Wang, L. Meng, and S. L. Zhu, Spectrum of the fully- troscopy within the hyperspherical formalism, Phys. Rev. D heavy tetraquark state QQQ¯ 0Q¯ 0, Phys. Rev. D 100, 096013 73, 054004 (2006). (2019).

116014-4 DISCOVERY POTENTIAL FOR THE LHCB FULLY-CHARM … PHYS. REV. D 102, 116014 (2020)

[24] X. Chen, Analysis of hidden-bottom bbb¯ b¯ states, Eur. Phys. [36] R. N. Faustov, V. O. Galkin, and E. M. Savchenko, Masses J. A 55, 106 (2019). of the QQQ¯ Q¯ tetraquarks in the relativistic diquark– [25] M. S. Liu, Q. F. Lü, X. H. Zhong, and Q. Zhao, All-heavy antidiquark picture, arXiv:2009.13237. tetraquarks, Phys. Rev. D 100, 016006 (2019). [37] C. Becchi, A. Giachino, L. Maiani, and E. Santopinto, A [26] C. Deng, H. Chen, and J. Ping, Towards the understanding study of ccc¯ c¯ tetraquark decays in 4 muons and in DðÞD¯ ðÞ of fully-heavy tetraquark states from various models, at LHC, arXiv:2006.14388. arXiv:2003.05154 [Phys. Rev. D (to be published)]. [38] G. Y. Chen and J. P. Ma, Production of Xð3872Þ at PANDA, [27] M. S. Liu, F. X. Liu, X. H. Zhong, and Q. Zhao, Full-heavy Phys. Rev. D 77, 097501 (2008). tetraquark states and their evidences in the LHCb di-J=ψ [39] S. Lange et al. (PANDA Collaboration), New studies of spectrum, arXiv:2006.11952. XYZ states at PANDA, arXiv:1311.7597. [28] Q. F. Lü, D. Y. Chen, and Y. B. Dong, Masses of fully heavy [40] L. Zotti, A. Filippi, S. Marcello, and S. Spataro, The tetraquarks QQQ¯ Q¯ in an extended relativized quark model, ψð4040Þ at the future P¯ ANDA experiment, J. Phys. Conf. Eur. Phys. J. C 80, 871 (2020). Ser. 503, 012007 (2014). [29] Z. G. Wang, Tetraquark candidates in the LHCb’s di-J=ψ [41] Q. Y. Lin, H. S. Xu, and X. Liu, Revisiting the production of mass spectrum, Chin. Phys. C 44, 113106 (2020). charmonium plus a light meson at PANDA, Phys. Rev. D [30] H. X. Chen, W. Chen, X. Liu, and S. L. Zhu, Strong decays 86, 034007 (2012). of fully-charm tetraquarks into di-charmonia, Sci. Bull. 65, [42] X. Y. Wang, J. J. Xie, and X. R. Chen, Production of the 1994 (2020). neutral Z0ð4430Þ in pp¯ → ψ 0π0 reaction, Phys. Rev. D 91, [31] G. Yang, J. Ping, L. He, and Q. Wang, A potential model 014032 (2015). prediction of fully-heavy tetraquarks QQQ¯ Q¯ (Q ¼ c, b), [43] X. Y. Wang and X. R. Chen, Discovery potential for the arXiv:2006.13756. neutral charmonium-like by annihilation, Adv. High Energy [32] X. Jin, Y. Xue, H. Huang, and J. Ping, Full-heavy tetra- Phys. 2015, 918231 (2015). quarks in constituent quark models, arXiv:2006.13745. [44] V. M. Abazov et al. (D0 Collaboration), Inclusive Produc- [33] M. Mikhasenko, L. An, and R. McNulty, The determination tion of the Xð4140Þ State in pp¯ Collisions at D0, Phys. Rev. of the spin and parity of a vector-vector system, arXiv:2007 Lett. 115, 232001 (2015). .05501. [45] A. V. Berezhnoy, A. K. Likhoded, A. V. Luchinsky, [34] J. F. Giron and R. F. Lebed, Simple spectrum of ccc¯ c¯ states and A. A. Novoselov, Production of J=ψ-meson pairs in the dynamical diquark model, Phys. Rev. D 102, 074003 and 4c tetraquark at the LHC, Phys. Rev. D 84, 094023 (2020). (2011). [35] M. A. Bedolla, J. Ferretti, C. D. Roberts, and E. Santopinto, [46] A. V. Berezhnoy, A. V. Luchinsky, and A. A. Novoselov, Spectrum of fully-heavy tetraquarks from a diquark þ Heavy tetraquarks production at the LHC, Phys. Rev. D 86, antidiquark perspective, Eur. Phys. J. C 80, 1004 (2020). 034004 (2012).

116014-5