Investigation of exotic state X(3872) in $pp$ collisions at $\sqrt{s}=$7 and 13 TeV

We have used the dynamically constrained phase space coalescence model to study the production of the exotic state $X(3872)$ based on the hadronic final states generated by the parton and hadron cascade model (PACIAE) with $|y|<1$ and $p_T<15.5$ GeV/c in $pp$ collisions at $\sqrt{s}=7$ and 13 TeV, respectively. Here the $X(3872)$ is assumed to consist of bound state $D\bar {D^*}$, which can form three possible structures for the tetraquark state, the nucleus-like state, and the molecular state. The yields of three different structures $X(3872)$ were predicted. The transverse momentum distribution and the rapidity distribution of three different structures $X(3872)$ are also presented. Sizable difference can be found in the transverse momentum and rapidity distributions for the three different $X(3872)$ structures.


I. INTRODUCTION
Hadrons are composite objects formed due to strong interactions formulated by the Quantum Chromodynamics (QCD) theory within the quark model framework. The quark model was formulated in 1964 to classify mesons as bound states made of a quark-antiquark pair, and baryons as bound states made of three quarks. For a long time all known mesons and baryons can be classified within this scheme. However, in principle QCD allows the existence of exotic states: multiquark states, hybrid states and glueball, as was already recognized by Gell-Mann in one of the first publications on the quark model [1], likeqqqq,qqqqqqqqqqq and so on. We will refer to any state which does not appear to fit with the expectations for an ordinaryqq or qqq hadron in the quark model as "exotic". Due to the non-perturbative properties of QCD in the low energy region, it is unavailable for us to use it to study the hadron structures and the hadron-hadron interactions directly. The study of exotic states can provide essential information on low energy QCD, which is absent in the ordinary qqq baryons and qq mesons. Therefore, comprehensive efforts have been made to test their existence and measure their properties.
The first exotic state X(3872) was discovered by the Belle collaboration almost two decades ago in the decay channel B ± → k ± π + π − J/ψ, with a significance of 10σ [2]. The new hadron was then confirmed by CDF [3], D0 [4], BABAR [5], and recently also by LHCb [6]. Forthermore, Belle and other collaborations have reported a lot of "XYZ" particles [7][8][9], which stimulated many researches on hadron spectrum. Some of the new states are unambiguously interpreted as conventional cc states, some are manifestly exotic, while for the others a definite interpretation is still missing. Despite the large amount of experimental data, the nature of the X(3872) state is still unclear. Several interpretations have been * Email:chengang1@cug.edu.cn (Corresponding Author) proposed, such as the conventional χ c1 (2P ) state [10], the molecular state [11][12][13], the tetraquark state [14], the ccg hybrid state [15], the vector glueball [16] or the mixed state [17,18].
In recent years, various theoretical models have been proposed to explain these exotic states. Meanwhile, these efforts have been made with various theoretical approaches like effective field theory [19], QCD sum rule [20], lattice QCD [21], gauge invariant model [22], potential model [23]. The extreme proximity of the X(3872) to the DD * threshold suggests that it could be a D * 0D0 , D 0D * 0 , D − D * + or D + D * − (hadron molecule or binging state).
Exploring the nature of exotic multiquark candidates such as the X(3872) plays a pivotal role in understanding quantum chromodynamics (QCD). In present days the study of exotic hadron has focused on the decay process, but the production process could actually provide more information. Despite significant efforts, there is still a lack of consensus on the process by which it came into being and its internal structure. There are indications that the multi-quark states may directly produced from multiproduction process of high-energy collisions, rather than from hadron decay [24]. Here taking X(3872) as an example, the parton and hadron cascade model (PA-CIAE) [25] plus the dynamically constrained phase space coalescence model (DCPC) [26][27][28] is proposed to simulate the complete evolution process from the initial state of the parton to the final state of the multiparticle in high-energy collisions. Therefore, we predict the yields and properties of the exotic hadron X(3872) and explore the production process and its structure of the exotic hadron X(3872) in high-energy pp collisions. This would provided a new and effective method for us to understand and study exotic hadrons.
In this paper, we will mainly consider three different structures of X(3872), tetraquark state, nucleus-like state, and molecular state. An approach for DCPC model is introduced to connect the multi-particle final state with the PACIAE simulations and to calculate the production of exotic state X(3872) in ultrarelativistic pp collisions. Firstly, the PACIAE models were used to generate the hadronic final states, including D * 0 ,D 0 , D 0 , D * 0 , D − , D * + , D + , and D * − . Then the yields, the transverse momentum distribution, as well as the rapidity distribution for X(3872) as tetraquark state, nucleus-like state, and molecular state were predicted using DCPC model in pp collisions at √ s = 7 and 13 TeV.

II. PACIAE MODEL AND DCPC MODEL
The parton and hadron cascade model PACIAE [25] is based on PYTHIA 6.4 to simulate various collision, such as e + e − , pp, p-A and A-A collisions. In general, PA-CIAE has four main physics stages, consisting of the parton initiation, parton rescattering, hadronization, and hadron rescattering. In the parton initiation, the string fragmentation is switched off temporarily in PACIAE and di(anti)quarks are broken into (anti)quarks. This partonic initial state can be regarded as quark-gluon matter (QGM) formed inside the parton initialization stage in the pp collisions. Then the parton rescattering in QGM is taken into account by the 2 → 2 LO-pQCD partonparton cross sections [29]. Their total and differential cross sections in the parton evolution are computed by the Monte Carlo method. In the hadronization process, the parton can be hadronized by the Lund string fragmentation re gime and/or the phenomenological coalescence model [29]. The final stage is the hadron rescattering process happening between the created hadrons until the hadronic freeze-out.
Theoretically, the production of composite hadronic states like light nuclei are usually calculated in two steps: First, the primordial simple hadrons including mesons and baryons are obtained with the transport model approach. Then, the nuclei or bound states are calculated by the phase-space coalescence model based on the Wigner function [30,31] or by the statistical model [32]. We introduced a dynamically constrained phase-space coalescence (DCPC) model to calculate the yield of bound states after the transport model simulations.
From quantum statistical mechanics [33], one can not simultaneously define both position q ≡ (x, y, z) and momentum p ≡ (p x , p y , p z ) of a particle in sixdimensional phase space because of the uncertainty principle, ∆ q∆ p ∼ h 3 . One can only say this particle lies somewhere within a six-dimensional quantum box or state of volume of ∆ q∆ p volume element in the sixdimensional phase space corresponding to a state of the particle. Therefore, one can estimate the yield of a single particle [33] by where E a , E b , and H denote energy threshold and the energy function of the particle, respectively. The variables q and p are the coordinates and momentum of the particle in the center-of-mass frame of the collision at the moment after hadronization. Furthermore, the yield of a cluster consisting of N particles is defined as following: Therefore, the yield of a X(3872) consisting of DD * cluster in the DCPC model can be calculated by where, The | q 12 | = | q 1 − q 2 | represent the distance between D and D * , m X(3872) denotes the rest mass of X(3872), and ∆m refers to its mass uncertainty. E 1 , E 2 and p 1 , p 2 denote the energies and momenta of the two particles (D and D * ), respectively. Here, the X(3872) is produced by the combination of hadrons D andD * after the mult-particle final states have been produced using the PACIAE model in pp collisions at √ s = 7 and 13 TeV. According to the different distances q 12 between D andD * , the exotic state X(3872) can be separated into three structures, the tetraquark state as D 0 < 0.6 fm, the nuclear-like state as 0.6 < D 0 < 1.36 fm and the molecular state as 1.36 < D 0 < 7 fm [34][35][36][37].

III. RESULTS
In this work, the final states hadrons are produces by the PACIAE model. The model parameters of PA-CIAE were fixed on the default values from the PA-CIAE model, except the parameters of parj(1), parj(2), parj(3), and parj(4). Here, parj(1) is the suppression of diquark-antidiquark pair production compared with the quark-antiquark pair production, parj(2) is the suppression of strange quark pair production compared with up (down) quark pair production, parj (3) is the extra suppression of strange diquark production compared with the normal suppression of a strange quark, and parj(4) is the suppression of spin 1 diquarks compared with spin 0 ones excluding the factor 3 coming from spin counting. These parameters are determined by fitting to the ALICE data [38,39] of D 0 , D + , D * + , π + and K + in midrapidity pp collisions at √ s = 7 TeV. And the comparison of the yields for each final states between the experimental measurement and the calculation from PACIAE model are shown in Table I, which are consistent with each other within uncertainty. The parameters parj(1), parj(2), parj(3 )and parj(4) found to be 0.08, 0.45, 0.40, 0.36, respectively. Notice that the yield of D 0 , D + and D * + for the experimental data in Table I is calculated according to the cross sections in Refs [38]. TeV and comparison with experimental data [38,39], with the 0 < pT < 36 GeV/c for D 0 , 1 < pT < 24 GeV/c for D + and D * + , 0.1 < pT < 3 GeV/c for π + , and 0.2 < pT < 6 GeV/c for K + .
particles Experiment PACIAE  3) and (4). The radius (the distance between the two mesons) and the mass value of X(3872) are given. Here we assume the excited states X(3872) is made up of DD * bound state, which is produced during the hadronization or hadronic evolution period. There are three popular configurations for the production of exotic state X(3872) in multiple production processes: the hadronic molecular state, the nuclear-like states, or the compact tetraquarks. We can determine the structure type of X(3872) from its radius value, and define the radii of the three structures as 8 fm, 1.36 fm and 0.6 fm [34,35], respectively. When the distance D0 between two D mesons is greater than 1.36 fm and less than 8 fm, the molecular state of X(3872) consists of two mesons DD * [35], denided as X m (3872); When D0 ≤ 0.6 fm, four quarks [ccqq] in X(3872) form a compact tetraquarks [34] through strong interaction, denote as X t (3872); When 0.6 fm ≤ D0 ≤ 1.36 fm, X(3872) is considered to consist of two mesons DD * to form a nuclearlike states, denote as X n (3872). The yield distribution of the three exotic states X(3872) with different structures to the uncertain mass parameter ∆m is shown in Fig. 1.
As shown in Fig.1, the yields of the X(3872) increase rapidly with parameter ∆m from 0 to 100 MeV. At the same ∆m and c.m. energy, the yields of molecular state X m (3872) are greater than the nucleus-like state X n (3872), and the yield of the nucleus-like state X n (3872) are greater than that of the tetraquark state X t (3872).
If we take half the decay width of X(3872) as the ∆m parameter, i.e, ∆m = Γ/2, then the yield of the X(3872) would be predicted. When we take the Breit-Wigner width of the X(3872) state that is measured to be (Γ = 3.9 +2.8+0.2 MeV) [40] by the experiment, ie. ∆m = Γ/2 = 1.95 MeV, the yield of exotic states X(3872) from different decay channels and different structures is given by PACIE+DCPC model in pp collision at √ s = 7 and 13 TeV, as shown in Table II.
Results from Table II show that the yields of the X(3872) from D * 0D0 and D 0D * 0 clusters are two orders of magnitude higher than the yields of D − D * + and D + D * − clusters. This indicates that the yield of X(3872) produced in pp collisions mainly comes from the D * 0D0 and D 0D * 0 clusters. While the yields of the X(3872) from D * 0D0 and D 0D * 0 clusters are almost same to each other and the yield of the X(3872) from D − D * + and D + D * − clusters are also close. In the last row of Table II, the total yields of the three different structured X(3872) particles produced in the pp collision are given. Obviously, the yield of X m (3872) (∼ 10 −6 ) is greater than that of X n (3872) (∼ 10 −7 ), and the yield of X n (3872) (∼ 10 −7 ) is greater than that of X t (3872) (∼ 10 −8 ), that is, the yield of X m (3872) is the highest in high-energy pp collisions. In addition, as c.m energy increases from 7 to 13 TeV, the yields of the X(3872) significant increased. The ratio between the different X(3872) structures are shown in Table III.
In fact, the generation and decay of exotic hadron X(3872) had been studied in the non-relativistic wave function method in reference [41]. Employing the factorized formulation with the help of event generators, they investigate the production of exotic hadrons in multi-production processes at high energy colliders. The distributions of rapidity and transverse momentum of X(3872) obtained by our simulation computed are close to their results. In reference [42], Pythia and Herwig models were used to simulate the production of X(3872) as a hadron molecule in the proton-proton/antiproton colliders, and the yields were estimated with the same order of magnitude as those calculated by us.
It should be noted that there are some uncertain factors in the yield of X(3872) given by us. First of all, the determination of the radius parameter for the X(3872) of three structures refers to the results of some theoretical calculations [34,35] rather than the experimental measurement of the radius of X(3872). The selection of different radius values will lead to different yield values. Secondly, we take the half-decay width Γ/2 of the mass spectrum for X(3872) → D 0 * D 0 + c.c. as the mass uncertainty parameter, and the half-decay width (Γ = 3.9 +2.8+0.2 −1.2−1.1 MeV) [40] measured in the experiment has a large systematic error, which also increases the uncertainty of X(3872)yield.
Since this is direct production, and only part of the X(3872) can decay again to D's via strong interaction, this will affect the production ratio of the D's [43]. In fact, the measurements of the X(3872) decays including X(3872) → D 0 * D 0 + c.c., γJ/ψ, π + π − J/ψ, γψ(3686), ωJ/ψ, etc. decay modes are all observed in the Belle, BABAR, BESIII, and LHCB experiments. X(3872) → D * 0D 0 + c.c. decay of branching ratio is about 52.4% [44]. Therefore, the total yield of X(3872) generated in the high-energy pp collision is approximately twice that of our simulation results. Fig.2 shows the transverse momentum pT distributions of the X(3872) as different structures using different bound states in pp collision at √ s = 7 and 13 TeV. Obviously, the shape of the pT distribution of X(3872) produced by different decay final states and different c.m. energies are similar to each other. But the molecular state X m (3872) has a wider pT distribution than tetraquark state X t (3872) and nucleus-like state X n (3872). We also predicted the rapidity distribution of X(3872) are also considered as above three structures in pp collisions at √ s = 7 and 13 TeV, which are symmetrically distributed in the range from −7.5 to 7.5, shown in Fig. 3. From this figure, one can see that the rapidity distributions for the three

IV. CONCLUSIONS
In this paper, D 0 , D + , D * + , π + and K + are generated in pp collision at √ s = 7 and 13 TeV by using the PACIAE model. The productions of each particle are consistent with the experimental data. Then the D andD * are put into the DCPC model to construct the clusters of D * 0D0 , D 0D * 0 ,  D − D * + and D + D * − to produce X(3872). The production and the characteristics of the X(3872) as tetraquark state, nucleus-like state, molecular state are predicted using different bound states. We found that the yields of X(3872) in pp collisions mainly comes from the D * 0D0 and D 0D * 0 clusters. If we take half of its decay width as the mass uncertainty ∆m, the yields of the X(3872) as tetraquark state, nucleus-like state, and molecular state in pp collision are predicted to be 8.22 × 10 −8 , 4.21 × 10 −7 , 1.01 × 10 −6 at √ s = 7 TeV and 9.00 × 10 −8 , 4.72 × 10 −7 , 1.18 × 10 −6 at √ s = 13 TeV, respectively. The transverse momentum pT distribution of X(3872) as tetraquark state, nucleus-like state are similar to each other while the rapidity distribution for the molecular state of X(3872) is slightly wider than that of other two structures.
To further understand the nature of exotic state X(3872), we therefore suggest to measure their yields and properties in pp and heavy-ion collisions by the ALICE and STAR experiments.