The nature of $X(3872)$ from high-multiplicity $pp$ collisions

The structure of exotic resonances that do not trivially fit the usual quark model expectations has been a matter of intense scientific debate during the last two decades. A possible way of estimating the size of these states is to study their behavior when immersed in QCD matter. Recently, LHCb has measured the relative abundance of the exotic $X(3872)$ over the ordinary $\psi(2S)$. We use the comover interaction model to study the yield of a compact $X(3872)$. To confirm the reliability of the model in high-multiplicity $pp$ collisions, we describe the suppression of excited over ground $\Upsilon$ states. With this at hand, we show that the size of the compact $X(3872)$ would be slightly larger than that of the $\psi(2S)$. If the $X(3872)$ is instead assumed to be a meson molecule of large size, we argue that its evolution in QCD matter should be described via a coalescence model, as suggested by data on deuteron production. We show that the predictions of this model for the $X(3872)$ are in contrast with data.

The last two decades witnessed a remarkable progress in heavy meson spectroscopy. Several new states, called XYZ, have been observed in the quarkonium sector, close to open flavor thresholds. Their properties are not well described by the conventional quark model/NRQCD, whence they are expected to have an exotic structure. In particular, the X(3872), observed as an unexpected peak in the J/ψ π + π − invariant mass, was the first of the series [1]. Its mass is almost exactly at theD 0 D * 0 threshold, and is remarkably narrow [2,3]. 1 The pion pair is dominated by the ρ meson, thus showing sizable isospin violation, unexpected if the X were an ordinary charmonium. Its structure has been subject of an intense debate [4][5][6][7].
Since the simple cc description cannot account for the observed features of X(3872), more valence quarks are needed. They could be aggregated by color forces in a new kind of hadron, a compact tetraquark of hadronic size ∼ 1 fm (see e.g. [4,[8][9][10][11]). Alternatively, if the coupling to the closest channel is dominant, nuclear forces could bind them in a hadron molecule which, given the extremely small binding energy would have a size of the order of 10 fm [5,[12][13][14][15], or more.
Recently, the LHCb collaboration has presented the production rates of promptly produced X(3872) relative to the ψ(2S), as a function of final state particle multiplicity [16]. This ratio is found to decrease with increasing multiplicity, an effect that has been known for decades to affect the production of ordinary quarkonia in proton- 1 Charge conjugation is understood throughout the paper. nucleus collisions. There is an ample consensus for this to be due to final state breakup interactions of the quarkonia with comoving particles [17].
Moreover, the ALICE collaboration has recently published an analysis for deuteron production in protonproton collisions [18]. The number of deuterons produced increases with multiplicity, hence showing a behavior that is qualitatively different from that of the X(3872). The idea that interactions with comovers could favor the coalescence of a hadron molecule was originally proposed in [19,20] for proton-proton, and in [21][22][23] for nucleusnucleus collisions.
In this work, we show how the combined study of the LHCb and ALICE points to a compact structure for the X(3872).
As a first step, we show that the comover model for compact states [17,28] reproduces the bottomonium yields observed in high-multiplicity pp collisions [29].With this at hand, we apply the same method, which assumes negligible recombination and a breakup cross section of the order of the geometrical one, to a compact tetraquark X(3872). The results are in agreement with the LHCb data. Pushing further the use of this model, we observe that a state of much larger size, would get a severe suppression.
However, the case of a loosely bound, large-size, hadronic molecule is most likely described in a different way. While recombination is known to be irrelevant for compact states [30], data on deuteron production in high multiplicity final states suggest the contrary. Therefore we extend the comover model to i) include the possible recombination of the hadronic pairs, and ii) implement the coalescence mechanism proposed in [19]. We find that,  TABLE I. Fixed values used in our parametrisation of the comover cross sections and the corresponding results. The procedure to compute vσ Q and its uncertainties are described in the text. The "binding energy" is computed with respect to D 0D0 -D + D − for the charmonium, and with respect toD 0 D * 0 [24] for a tetraquark X(3872) (the OZI-favored mode Λ + c Λ − c [25][26][27] being kinematically forbidden at the typical comover's energy). The average for this binding energy is 44 ± 116 keV [2,3], and in our calculations, we use the 1σ error. The radius of the tetraquark is taken from [10,24]. The error on vσ Q depends on the uncertainty of T eff , and on whether considering pionic or gluonic comovers.
while this reproduces well the deuteron data, it fails with the X(3872). Our results are therefore consistent with a compact tetraquark interpretation, and a destruction cross section comparable to that of other compact states.

THE COMOVER INTERACTION MODEL
To include final state interactions for compact states, we follow the comover interaction model (CIM) [17,28]. Within this framework, quarkonia are broken by collisions with comovers-i.e. final state particles with similar rapidities. The density of quarkonium ρ Q , at a given transverse position s and rapidity y, for a collision of impact parameter b, evolves following where vσ Q is velocity times the cross section of quarkonium dissociation, averaged over the momentum distributions of the comoving particles, whose transverse density is ρ c at initial time τ i . The above equation neglects recombination effects which, for a compact object, are irrelevant due to the paucity of heavy quarks produced in the pp environment considered [28,30]. Integrating the equation above from τ i to τ f , we get the quarkonium density for a given position and impact parameter, where the argument of the logarithm comes from τ f /τ i converted into ratio of densities, effectively playing the role of time. Indeed, the interaction stops at τ f , when the densities have diluted down to ρ pp , the value of the minimum-bias pp density at that energy and rapidity, as taken from [31]. The previous equation can already be used for an estimate of the qualitative behavior of the yields as a function of multiplicity. This can be done by neglecting the dependence on the variables and consider average values only. The comover density would simply be related to the number of charged particles by ρ c (b, s, y) 3 2 N ch /σ.
The factor 3/2 accounts for the neutral comovers, while σ is the inelastic pp cross section. The full spatial dependence can be described with an eikonal-Glauber model, as it is standard in heavy ion phenomenology [17,28,32,33]. The basic ingredient is the profile function of the proton, taken as a Fermi function. The comover density ρ c (b, s, y) is proportional to the number of binary parton-parton collisions per unit of transverse area ds and rapidity at a given impact parameter b, which in turns is proportional to the overlap between the protons-see, e.g., [33]. The normalization is fixed to reproduce the minimum-bias pp multiplicity-i.e. the pp multiplicity averaged over all impact parameters. Proceeding this way, the quarkonium yields are obtained weighting Eq. (2) with the pp overlap function.
As seen from Eq. (1), the quarkonium abundance is driven by its interaction cross section with comovers vσ Q . In nucleus-nucleus collisions at lower energies the latter has been fitted from data (independently state by state). However, it has also been effectively related to the geometrical cross section by [17] vσ Q = σ geo Here σ geo Q = πr 2 Q , with r Q the quarkonium radius. Moreover, E thr Q = B Q + m c , where B Q is the distance between the quarkonium mass and the closest open flavor threshold (OZI-favored), and E c = m 2 c + p 2 c the energy of the comovers in the quarkonium center-of-mass frame. For states close to the threshold, the breakup cross section is essentially the geometric one. The average is computed over an isotropic Bose-Einstein distribution for the comover energy, Both n and T eff are phenomenological parameters. Attempts to compute n using the multipole expansion in perturbative QCD at leading order suggest n 4 for pion comovers, making the strong assumption that the scattering is initiated by the gluons inside these pions [34][35][36]. More realistic hadronic models suggest a smaller value [37,38]. The dependence on the comover mass (gluons or pions) is not dramatic. In Table I we quote the average of these two possibilities. A fit on the relative yields of excited-to-ground state Υ data at LHC in pPb collisions gave T eff = 250 ± 50 MeV and n = 1 [17]. We adopt the same values here. 2 We recall two features of the comover approach. First, larger particles are more affected by dissociation, due to larger interaction cross sections. As a consequence, excited states are more suppressed than the ground states. Second, the suppression increases with comover densities, which is proportional to particle multiplicities: it increases with centrality in nucleus-nucleus collisions, and it is stronger in the nucleus direction for proton-nucleus collisions.
To confirm the applicability of the CIM also to compact states in proton-proton collisions, we use it to describe the yields of Υ mesons. These have been measured by CMS at 2.76 TeV, as a function of the number of charged tracks with p T > 400 MeV, reconstructed in the tracker at and |η| < 2.4 [29]. In Figure 1 we show the data together with our results obtained using the breakup cross sections of [17], confirming the validity of the model. The global normalisation corresponds to the experimental value at N ch = 15±2, that we identify with the mean multiplicity. 3 We extend our calculation to charmonia by applying Eq. (3) for the cross sections, using the same T eff and n. Although the non-perturbative value of n could in principle be different from the bottomonium one, we get cross section values-5 mb for ψ(2S)-in the ballpark of those obtained by directly fitting the charmonium data [32]-6 mb for ψ(2S)-, confirming that a unified description is possible.

THE X(3872) IN THE CIM MODEL
The relative production rates of prompt X(3872) over ψ(2S) have been measured by LHCb in pp collisions at 8 TeV, in the forward pseudorapidity region, 2 < η < 5 [16], as a function of the number of charged particle tracks reconstructed in the VELO detector. This ratio is found to decrease with increasing multiplicity.
As mentioned above, the suppression of the state is driven by its interaction cross section with the comovers, as reported in Table I. Following the model described in the previous section, we compute the N X(3872) /N ψ(2S) ratio as a function of N ch -see Figure 2. For a compact tetraquark of typical hadronic size, the same physics described for quarkonia must apply, Relative yields of excited-to-ground state Υ as a function of multiplicity for pp collisions at 2.76 TeV in the central region, as measured by CMS [29]. The bands follow the uncertainties of the six cross sections that contribute via the feed down, and the one of T eff . Our results are normalised to the experimental values corresponding to N ch = 15 ± 2.
FIG. 2. Relative yield of X(3872) vs ψ(2S) as a function of event multiplicity for pp collisions at 8 TeV and forward pseudorapidity, as measured by LHCb [16]. The assumption of a tetraquark of diameter 1.3 fm reproduces well the experimental data. The two error bands are computed as in Figure 1. Applying the geometrical picture to a molecular state predicts a very sharp suppression. The coalescence picture also predicts a behavior that is qualitatively different from data. See text for the definition of error band in this case. with the breakup cross section being dictated by the geometric one. Indeed, assuming a diameter of 1.3 fm [10,24], the CIM gives results which describe well the LHCb data. In particular, it predicts a 20% decrease in the ratio when going from the first to the second multiplicity bin, similar to the Υ states ( Figure 1).
One could apply the same geometrical estimate for the breakup of a molecular X(3872). Clearly, being the size much larger, the corresponding suppression is way too strong to describe the data, as shown in Figure 2. The steep drop of this curve is a qualitative indication that the description of the interaction of a molecule with comover particles must be refined. We explain this in the next section.

COALESCENCE OF HADRON MOLECULES
The implementation of comover interactions discussed above disregards recombination. While negligible for compact states in pp collisions [30], recombination is needed to explain deuteron data, whose yield increases with increasing multiplicity of the final state [18,40]. Were a molecular X(3872) to behave like that, it would be in striking contrast with the data in [16]. Moreover, for a large composite object like a hadron molecule, one expects the interaction with comovers to be dominated by the scattering off the molecule constituents, an idea already put forth in a number of papers [19,[21][22][23]. A popular description of the destruction and recombination of molecules is indeed in terms of coalescence (see, e.g., [41][42][43]). In this picture the constituents are bound/free depending on whether their relative momentum is smaller/larger than some Λ. In [19] it has been proposed that the driving process is given by the scattering πhh πm, where hh are free constituents, m the molecule itself and π the comover.
In what follows we first present a derivation of the corresponding evolution equation. We then obtain the creation and destruction cross sections for the deuteron and the X(3872) in the molecular regime. The effective coupling are extracted from well known hadronic physics, and are compatible with several other results (see, e.g., [44,45]). The particle distributions for the molecules, free constituents and comovers are instead either taken from data, or from Monte Carlo simulations, which have been tuned over the years to reproduce hadronic results (see, e.g., [46]). As a check of the soundness of our result, we note that the destruction cross section for the X(3872) is perfectly compatible with what obtained in [22] using different techniques and a different model. Finally, we employ all this to compute the yields of deuteron and X(3872) as a function of the event multiplicity, which we compare with observations.

Boltzmann equation for hadronic molecules and comovers
In this derivation, we follow [47]. We use the standard noncovariant form of the Boltzmann equation, and the quantities discussed below are defined in the lab frame.
Consider two hadrons, '1' and '2', with positions and momenta (x i , q i ). We adopt the coalescence picture, i.e. that the two hadrons bind if their relative momentum is smaller than some threshold, |q 1 − q 2 | < Λ [41][42][43]. In other words, if the particles are close enough in phase space, their mutual interaction will induce the formation of the molecule. The phase-space density for these two hadrons, where R Λ is the domain where |q 1 − q 2 | < Λ, andR Λ its complement. N m , N hh and N 12 are respectively the number of molecules, free pairs and total pairs. Clearly N 12 = N m (τ ) + N hh (τ ), and it is assumed to be constant in time, hence neglecting creation/annihilation of the constituents themselves. In the lab frame, the variable τ can be taken simply as the time.
We assume that the spatial and momentum distributions factorize as where both f m and f hh are normalized to unity when integrated over their momentum domain. Indeed, in absence of collisions, the particles travel freely, and the spatial and momentum distribution are hence uncorrelated. Our assumption is that, to first order in perturbation theory, this is preserved. Collisions with comovers change the momenta of the constituents, and consequently modify their distribution. Choosing the z-axis to be in the beam direction, when the longitudinal expansion is much faster than the transverse one, the Boltzmann equation reduces to [47]: where the loss and gain terms are, respectively Here W is the nonrelativistic matrix element for the q 2 + q 3 q 2 +q 3 process, and δ P ≡ 2π 4 δ (4) q 2 +q 3 −q 2 −q 3 enforces conservation of energy and momentum. We assumed that the comovers interact with constituents '1' and '2' equally. Without loss of generality, we restrict the interaction to the constituent '2', so that the position of the comover must be x 2 . The interaction with '1' is taken into account later, by a factor of 2 in the cross sections. Moreover, n c is the phase-space distribution of comovers, which again we factorize as n c (x, q, τ ) ρ c (x, τ )f c (q, τ ).
To study the evolution of the density of molecules, we now integrate Eq. (7) over q 1 and q 2 in R Λ , in order to isolate the molecular contribution. We also assume that the momentum distribution of the molecule follows the free-stream (collisionless) distribution [48]. After that, the left hand side of Eq. (7) simply returns ∂ρ m (x i , τ )/∂τ . The loss term instead gives where the average cross section for the destruction of a molecule is defined as The gain term instead requires a bit more care. Imposing momentum conservation, its integral over R Λ gives The relative momentum appearing in the distribution is then q 1 − q 2 + q 3 − q 3 . Now, by construction |q 1 − q 2 | < Λ, while the comovers distribution in the lab frame is dominated by momenta |q 3 | Λ. This means that, barring small integration regions, for most configurations q 1 − q 2 + q 3 − q 3 > ∼ Λ, and we can use Eq. (6) and write with average cross section for the creation of a molecule given by Eqs. (10) and (13) are simply the cross sections for the processes of interest (destruction or creation) summed over all possible final states and averaged over the initial ones, with suitable distributions of the momenta. The Boltzmann equation now reads Finally, if the density of comovers is roughly homogeneous, we can integrate over x 1 and x 2 as well. The final evolution equation then reads In the last line we used the fact that N 12 = N m (τ ) + N hh (τ ), in order to take advantage of the fact that N 12 = constant. In the free-stream approximation, the spatial density of comovers evolves with time as ρ c (τ ) = ρ c /τ , with ρ c the transverse density at formation time τ i . With this approximation, the solution to Eq. (15) for the number of molecules as a function of multiplicity, N ch , is: were N 0 m is the number of molecules generated by hadronization, before any interaction with comovers. The dependence on multiplicity comes from the comovers spatial density. For vσ m = 0 the result above reduces to Eq. (2), upon volume integration.

Effective couplings for comover-constituent interaction
To describe the creation/annihilation of the deuteron and the molecular X(3872) by the interaction with comovers we need to determine their effective couplings. To simplify the computation we describe the comoverconstituent scattering with a constant relativistic matrix element, M πh→πh = g 2 . As our notation suggests, we assume that the comovers are all pions, which are indeed the dominant fraction.
For the case of the deuteron, the effective coupling is obtained by matching the total elastic pion-nucleon cross section from PDG [3], averaged in the [0, 300 MeV] kinetic energy range for the comover (therefore reaching the peak of the intermediate ∆ resonance). One gets g 2 /(4π) 10, which is not far from the threshold value g 2 /(4π) 13.5 [45].
For the X(3872) no data are available for the πD ( * ) → πD ( * ) scattering. To obtain the coupling we then consider the following Lagrangian for the interaction between pions and heavy mesons Here a, b are isospin indices, the traces are taken over Dirac matrices and M is the mass of the heavy mesons. Moreover, H a is the HQET heavy meson multiplet and M ab is the pion matrix. The trilinear coupling is given by y 0.8 [49]. We take the propagator of the D * to be the real part of the pole being at the mass of the D to avoid collinear divergences in the u-channel (note that the D and D * are indeed degenerate at leading order in HQET). Moreover, we choose the projector to be exactly transverse (rather than only on-shell) to make sure that the contribution from the off-shell D * propagation vanishes at threshold. This way, the only contribution to the scattering length is given by the quartic coupling, which is found to be λ 25, by matching with lattice calculation of the πD scattering length [50]. Given the above vertices, the effective coupling g is again obtained matching the total elastic πD cross section averaged over the kinetic energy of the comover in the [0, 300] MeV range. The processes we considered are reported in Figure 3. The result is g 2 /(4π) 5.07.

Creation and destruction average cross sections
The last step we need to perform to implement the evolution equation (16) is to average the cross sections over the particle distributions-see Eqs. (10) and (13). The calculation is conceptually straightforward but rather tedious. We will spare most of it to the reader, and only highlight the main points.
We always work in the approximation of roughly massless comovers, m c 0, and of loosely bound molecule, Λ |q i |. In particular, the latter implies that the momenta of the constituents are approximately equal to each other and to half the momentum of the molecule itself.
To study the momentum distributions we employ Pythia8, and generate 750k pp events at √ s = 7 TeV for deuteron, and 2.5G pp events at 8 TeV for the X(3872), with full-QCD 2 → 2 matrix elements and a cut on the partonic transverse momentum,p ⊥ > 2 GeV. Long-lived particles are prevented from decaying. We select all the events that have at least one constituent pair in the final state and, when more than one pair per event is available, all the combinations are considered. By comovers, we mean all those particles whose momentum lies in a cone ∆R = ∆φ 2 + ∆η 2 < 0.4 from one of the constituents (the other constituent is clearly excluded). We checked that varying this cut within the range ∆R = 0.3 − 0.5 affects the final results by no more than 3%.
In the lab frame, most particles have large momentum component along the beam axis. The polar angles between the constituent momenta, between the comovers and the molecule momentum, and between the comovers and the constituents relative momentum can thus be approximated by cos θ ±1 (see Figure 4). This amounts to replacing the θ-dependence of the corresponding distributions with a sum of δ-functions with support on cos θ = ±1.
After this, the two integrals can be reduced to Where both the sum in the first equation and the Heaviside function in the second arise from requiring for the δ-functions for conservation of energy to have support on the integration region. Here g is the comover-constituent coupling defined in section 3 3.2, P and k the total and relative momenta of the hadron pair and M the mass of the constituents, which we take to be approximately equal. The momentum distributions have been integrated over all variables except those explicitly written. More- , while the minimum value of the momentum q 3 is q min The distributions extracted from Pythia8 are reported in Figure 5 for both the proton-neutron and theD 0 D * 0 pairs. They are well described by the following functional forms: with best fit values given in Table II. The overall constant if fixed by normalization, as described below Eq. (6). The distribution of the total momentum of the molecule can instead be obtained from the experimental distributions in transverse momentum and rapidity. In particular with F m (P ⊥ , y) the experimental distribution, assumed even under y → −y, Y the experimental cut in rapidity, TeV of center-of-mass energy as a function of multiplicity, as reported in [18]. The solid line is our result (16), the uncertainty being determined by varying the coalescence momentum Λ between 50 and 250 MeV [51].
and sinhȳ ≡ (P 2 − P 2 ⊥ )/(4M 2 + P 2 ⊥ ). With all this at hand, we now compute the average cross sections in Eqs. (19). To test the validity of our idea, we first consider the case of the deuteron. In particular, we take M = 938 MeV and a coalescence momentum in the range Λ = 50-250 MeV [51]. The effective coupling g is discussed in section 3 3.2. The distribution in rapidity is approximately uniform, while the one in transverse momentum is well fitted by a Lévy-Tsallis function [18]. We obtain vσ m Λ 150 MeV 3 × 0.51 mb and vσ hh 4.34 mb. The suppression of the creation cross section compared to the destruction one is understood as follows. In a hadronic collision, the vast majority of constituents are produced free, with relative momentum much larger than Λ ( Figure 5). After an interaction with a comover of momentum of order ∼ GeV, it is unlikely for the pair to fall within the small region of phase space with k < ∼ Λ. Technically, this is due to the exponential suppression of the comover distribution at momenta of the order of the free pair ones. Similarly, if a pair is initially bound, it is much more likely for a comover to increase their relative momentum to a value larger than Λ, rather than viceversa.
We estimate the number of initial deuterons with Pythia8, by counting the proton-neutron pairs with relative momentum initially below Λ. We find N 0 m /N 12 = O 10 −4 , which can be neglected. Hence, the dependence of the number of deuterons on multiplicity is fixed up to an overall factor, which we fit to data. We also set ρ pp so that the curve starts at N ch = 1, as in data.
In Figure 6 we compare our results to the ALICE data. For simplicity, we just consider the yield as in Eq. (2), computed on the average comover density, ρ c (b, s, y) 3 2 N ch /σ. The good match with experiment confirms the validity of the coalescence approach, and of the idea proposed in [19], especially the relevance of the comover-constituent interaction to the enhancement of the production of hadron molecules.

THE X(3872) IN THE COALESCENCE MODEL
We now apply the same procedure to the X(3872). The momentum distributions for the comovers and the freeD 0 D * 0 pairs are again given in Eqs. (21) and Table II, while those for the X(3872) are obtained from a NRQCD calculation [52], which reproduces well the prompt production data at high p ⊥ [53]. The coalescence momentum for the X(3872) is unknown, and we take it in the range Λ = 30-360 MeV, as proposed in the literature [19,46,52,53]. This contributes minimally to the error band in Figure 2. With the effective coupling found in section 3 3.2, and taking M (m D +m D * )/2 = 1936 MeV, the creation and annihilation cross section are be vσ m (Λ/30 MeV) 3 × 7.1 × 10 −6 mb and vσ hh 0.50 mb. The creation cross section is much smaller than the deuteron one since the constituents of the X(3872) are produced with harder momenta (Figure 5), and hence the creation process suffers from a stronger suppression from the exponential comover distribution: it is even more unlikely to reduce the relative momentum to a value smaller than Λ. The destruction cross section is in good agreement with what obtained in [22] with different methods and for different processes.
As for the initial number of X(3872) produced by hadronization alone, N 0 m , there is still no consensus. On the one hand, a purely molecular interpretation requires for this number to be very small because of the difficulty in producing the constituents with such a small relative momentum [46,53], as for the deuteron, albeit some controversies on final state interactions [52,53]. On the other hand, it has been suggested that the production of the X(3872) could be dominated by short distance physics, likely associated to a charmonium component of its wave function [54]. Were this to be true, the prompt production cross section could be significantly enhanced. Here we adopt an agnostic viewpoint and let the initial number of molecules vary from N 0 m /N 12 = 0 (the estimate with Pythia8 being 10 −6 ) to N 0 m /N 12 = 1. We consider this to be the main source of uncertainty in our calculation. Either way, the molecular nature of the X(3872) must be manifest when propagating throughout the comovers for distances of O(1 fm).
Our results are reported in Figure 2, and are qalitatively at odds with data, regardless on whether the molecular X is copiously produced by hadronization or not. The number of X(3872) normalized to ψ(2S) always grows, similarly to the deuteron. Indeed, the typical cross sections obtained for the X(3872) are much smaller than the one for ψ(2S) shown in Table I. Thus the increasing behavior of the ratio is dominated by the decreasing ψ(2S) yield, regardless of the details of N 0 m /N 12 .

CONCLUSIONS
The production of X(3872) at high transverse momenta in low multiplicity pp collisions challenges the molecular interpretation [23,46,53], to the extent that it is necessary to assume its hadronization proceeds through a compact cc core. In [19] it was shown that not even the interaction with comoving particles was able to account for the large number of X(3872), if this compact component is not considered. However, the recent high-multiplicity data from the LHCb and ALICE collaborations [16,18], encourages to reconsider the role of comovers.
In this paper we redesigned the molecule-comover interaction model, treating multiple scattering with kinetic theory. This works remarkably well at explaining the deuteron production reported by ALICE [18]. Were a sizeable molecular component to appear in the X(3872) wave function, the same approach should describe its relative yield with respect to ψ(2S) [16]. However, the predicted yield always grows and cannot match the decreasing slope observed by LHCb. The only way to reconcile the results from the coalescence model with experiment is to make the averaged X(3872) molecular destruction cross section, vσ hh , about twenty times larger (effectively that of a compact state), in sharp contradiction with several agreeing determinations of the interaction couplings of pion comovers with D, D * mesons and with the findings in [22], for example. Our analysis concludes that, despite the closeness to theD 0 D * 0 threshold that motivated several studies (see, e.g., [5,[12][13][14][15]), a molecular component cannot dominate the wave function of the X(3872).
The LHCb results are analyzed also with the Comover Interaction Model (CIM) for compact states [17,28], which is well known to describe the quarkonia yields in high multiplicity final states. The yields are determined assuming as effective cross section the geometrical one, i.e. the size of the states. For the first time, we apply the CIM to pp collisions, and match the relative yields of Υ mesons with the ones reported by CMS [29]. These yields decrease with multiplicity, contradicting the statement that such a behavior requires a molecular interpretation, as suggested initially for the X(3872). Conversely, it is perfectly compatible with a compact tetraquark of typical hadronic size, as we showed here.
After the appearing of this work, a modification of the CIM has been effective in reconciling the diffusion of a molecule with the LHCb data [55]. However, we note that such description considers no recombination, which makes it hard to achieve a common understanding of the X(3872) and deuteron data.
We look with interest also to PbPb data [56], which seem to present novelties with respect to pp. It would be useful to have them binned in centrality, to allow a comparison with the deuteron data presented by AL-ICE [57]. Discussions on tetraquarks and molecules in PbPb collisions can be found in [58].
In conclusion, as soon as the loosely-bound molecule description is made concrete, its expected behavior appears different from what observed experimentally. To fill the discrepancy we should tune the interaction cross sections well beyond the range we estimate. In our view, the LHCb data on the X(3872) display the same features characterizing compact states like the Υ mesons.