Study of the doubly charmed tetraquark $${{{{{{\rm{T}}}}}}}

Nature Communications volume 13, Article number: 3351 (2022) Cite this article

3752 Accesses

60 Citations

446 Altmetric

Metrics details

Quantum chromodynamics, the theory of the strong force, describes interactions of coloured quarks and gluons and the formation of hadronic matter. Conventional hadronic matter consists of baryons and mesons made of three quarks and quark-antiquark pairs, respectively. Particles with an alternative quark content are known as exotic states. Here a study is reported of an exotic narrow state in the D0D0π+ mass spectrum just below the D*+D0 mass threshold produced in proton-proton collisions collected with the LHCb detector at the Large Hadron Collider. The state is consistent with the ground isoscalar \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) tetraquark with a quark content of \({{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}\overline{{{{{{\rm{u}}}}}}}\overline{{{{{{\rm{d}}}}}}}\) and spin-parity quantum numbers JP = 1+. Study of the DD mass spectra disfavours interpretation of the resonance as the isovector state. The decay structure via intermediate off-shell D*+ mesons is consistent with the observed D0π+ mass distribution. To analyse the mass of the resonance and its coupling to the D*D system, a dedicated model is developed under the assumption of an isoscalar axial-vector \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state decaying to the D*D channel. Using this model, resonance parameters including the pole position, scattering length, effective range and compositeness are determined to reveal important information about the nature of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state. In addition, an unexpected dependence of the production rate on track multiplicity is observed.

Hadrons with quark content other than that seen in mesons (\({{{{{{\rm{q}}}}}}}_{1}{\bar{{{{{{\rm{q}}}}}}}}_{2}\)) and baryons (q1q2q3) have been actively discussed since the birth of the quark model1,2,3,4,5,6,7,8. Since the discovery of the χc1(3872) state9 many tetraquark and pentaquark candidates, listed in Table 1, have been observed10,11,12,13,14,15,16,17,18,19. For all but the X0(2900) and X1(2900) states the minimal quark content implies the presence of either a \({{{{{\rm{c}}}}}}\overline{{{{{{\rm{c}}}}}}}\) or \({{{{{\rm{b}}}}}}\overline{{{{{{\rm{b}}}}}}}\) quark-antiquark pair. The masses of many tetra- and pentaquark states are close to mass thresholds, e.g., \({{{{{{\rm{D}}}}}}}^{(* )}{\overline{{{{{{\rm{D}}}}}}}}^{(* )}\) or \({{{{{{\rm{B}}}}}}}^{(* )}{\overline{{{{{{\rm{B}}}}}}}}^{(* )}\), where D(*) or B(*) represents a hadron containing a charm or beauty quark, respectively. Therefore, these states are likely to be hadronic molecules16,20,21,22 where colour-singlet hadrons are bound by residual nuclear forces, such as the exchange of a pion or ρ meson23, similar to electromagnetic van der Waals forces attracting electrically neutral atoms and molecules. These states are expected to have a spatial extension significantly larger than a typical compact hadron. Conversely, the only hadron currently observed that contains a pair of c quarks is the \({{{\Xi }}}_{{{{{{\rm{cc}}}}}}}^{++}\) (ccu) baryon, a long-lived, weakly-decaying compact object24,25. The recently observed X(6900) structure in the J/ψJ/ψ mass spectrum26 belongs to both categories simultaneously. Its proximity to the χc0χc1 threshold could indicate a molecular structure27,28. Alternatively, it could be a compact object, where all four quarks are within one confinement volume and each quark interacts directly with the other three quarks via the strong force29,30,31,32.

The existence and properties of \({{{{{{\rm{Q}}}}}}}_{1}{{{{{{\rm{Q}}}}}}}_{2}{\bar{{{{{{\rm{q}}}}}}}}_{1}{\bar{{{{{{\rm{q}}}}}}}}_{2}\) states with two heavy quarks and two light antiquarks have been widely discussed for a long time33,34,35,36,37,38. In the limit of large masses of the heavy quarks the corresponding ground state should be deeply bound. In this limit, the two heavy quarks Q1Q2 form a point-like colour-antitriplet object, analogous to an antiquark, and as a result, the \({{{{{{\rm{Q}}}}}}}_{1}{{{{{{\rm{Q}}}}}}}_{2}{\bar{{{{{{\rm{q}}}}}}}}_{1}{\bar{{{{{{\rm{q}}}}}}}}_{2}\) system has similar degrees of freedom for its light quarks as an antibaryon with a single heavy quark, e.g., the \({\overline{\Lambda }}_{{{{{{\rm{c}}}}}}}^{-}\) or \({\overline{\Lambda }}_{{{{{{\rm{b}}}}}}}^{0}\) antibaryons. The beauty quark is considered heavy enough to sustain the existence of a \({{{{{\rm{b}}}}}}{{{{{\rm{b}}}}}}\overline{{{{{{\rm{u}}}}}}}\overline{{{{{{\rm{d}}}}}}}\) state that is stable with respect to the strong and electromagnetic interactions with a mass of about 200 MeV/c2 below the BB* mass threshold. In the case of the \({{{{{\rm{b}}}}}}{{{{{\rm{c}}}}}}\overline{{{{{{\rm{u}}}}}}}\overline{{{{{{\rm{d}}}}}}}\) and \({{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}\overline{{{{{{\rm{u}}}}}}}\overline{{{{{{\rm{d}}}}}}}\) systems, there is currently no consensus in the literature whether such states exist and if their natural widths are narrow enough to allow for experimental observation. The theoretical predictions for the mass of the \({{\mbox{}}}{{{{{\rm{c}}}}}}{{\mbox{}}}{{{{{\rm{c}}}}}}{{\mbox{}}}\overline{{{{{{\rm{u}}}}}}}{{\mbox{}}}\overline{{{{{{\rm{d}}}}}}}{{\mbox{}}}\) ground state with spin-parity quantum numbers JP = 1+ and isospin I = 0, denoted hereafter as \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\), relative to the D*+D0 mass threshold

lie in the range −300 < δm < 300 MeV/c2 39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66,67,68,69,70, where \({m}_{{{{{{{\rm{D}}}}}}}^{* +}}\) and \({m}_{{{{{{{\rm{D}}}}}}}^{0}}\) denote the known masses of the D*+ and D0 mesons10, with \({{{{{\rm{c}}}}}}\overline{{{{{{\rm{d}}}}}}}\) and \({{{{{\rm{c}}}}}}\overline{{{{{{\rm{u}}}}}}}\) quark content, respectively. The observation of a narrow state in the D0D0π+ mass spectrum near the D*+D0 mass threshold, compatible with being a \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) tetraquark state with \({{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}\overline{{{{{{\rm{u}}}}}}}\overline{{{{{{\rm{d}}}}}}}\) quark content is reported in Ref. 71.

In the work presented here, the properties of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state are studied by constructing a dedicated amplitude model that accounts for the D*+D0 and D*0D+ decay channels. In addition, the mass spectra of other DD(*) and opposite-sign \({{{{{\rm{D}}}}}}{\overline{{{{{{\rm{D}}}}}}}}^{(* )}\) combinations are explored. Furthermore, production-related observables, such as the event multiplicity and transverse momentum (pT) spectra that are sensitive to the internal structure of the state, are discussed. This analysis is based on proton–proton (pp) collision data, corresponding to an integrated luminosity of 9 fb−1, collected with the LHCb detector at centre-of-mass energies of 7, 8 and 13 TeV. The LHCb detector72,73 is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of particles containing b or c quarks and is further described in Methods.

The D0D0π+ final state is reconstructed using the D0 → K−π+ decay channel with two D0 mesons and a pion all produced promptly in the same pp collision. The inclusion of charge-conjugated processes is implied throughout the paper. The selection criteria are similar to those used in Refs. 74,75,76,77 and described in detail in Methods. The background not originating from true D0 mesons is subtracted using an extended unbinned maximum-likelihood fit to the two-dimensional distribution of the masses of the two D0 candidates from selected D0D0π+ combinations, see Methods and Supplementary Fig. 1a. The obtained D0D0π+ mass distribution for selected D0D0π+ combinations is shown in Fig. 1.

Distribution of D0D0π+ mass where the contribution of the non-D0 background has been statistically subtracted. The result of the fit described in the text is overlaid. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

An extended unbinned maximum-likelihood fit to the D0D0π+ mass distribution is performed using a model consisting of signal and background components. The signal component corresponds to the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decay and is described as the convolution of the natural resonance profile with the detector mass resolution function. A relativistic P-wave two-body Breit–Wigner function \({{\mathfrak{F}}}^{{{{{{\rm{BW}}}}}}}\) with a Blatt–Weisskopf form factor78,79 is used in Ref. 71 as the natural resonance profile. That function, while sufficient to reveal the existence of the state, does not account for the resonance being in close vicinity of the D*D threshold. To assess the fundamental properties of resonances that are close to thresholds, advanced parametrisations ought to be used80,81,82,83,84,85,86,87,88,89,90. A unitarised Breit–Wigner profile \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\), described in Methods Eq. (47), is used in this analysis. The function \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) is built under two main assumptions.

The newly observed state has quantum numbers JP = 1+ and isospin I = 0 in accordance with the theoretical expectation for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) ground state.

The \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state is strongly coupled to the D*D channel, which is motivated by the proximity of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) mass to the D*D mass threshold.

The derivation of the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) profile relies on the assumed isospin symmetry for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{* }{{{{{\rm{D}}}}}}\) decays and the coupled-channel interaction of the D*+D0 and D*0D0 system as required by unitarity and causality following Ref. 91. The resulting energy-dependent width of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state accounts explicitly for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\), \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{+}{\pi }^{0}\) and \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{+}\gamma \) decays. The modification of the D* meson lineshape92 due to contributions from triangle diagrams93 to the final-state interactions is neglected. Similarly to the \({{\mathfrak{F}}}^{{{{{{\rm{BW}}}}}}}\) profile, the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) function has two parameters: the peak location mU, defined as the mass value where the real part of the complex amplitude vanishes, and the absolute value of the coupling constant g for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{* }{{{{{\rm{D}}}}}}\) decay.

The detector mass resolution, \({\mathfrak{R}}\), is modelled with the sum of two Gaussian functions with a common mean, and parameters taken from simulation, see Methods. The widths of the Gaussian functions are corrected by a factor of 1.05, which accounts for a small residual difference between simulation and data94,95,96. The root mean square (RMS) of the resolution function is around 400 keV/c2.

A study of the D0π+ mass distribution for selected D0D0π+ combinations in the region above the D*0D+ mass threshold and below 3.9 GeV/c2 shows that approximately 90% of all D0D0π+ combinations contain a true D*+ meson. Therefore, the background component is parameterised with a product of the two-body phase-space function \({{{\Phi }}}_{{{{{{{\rm{D}}}}}}}^{* +}{{{{{{\rm{D}}}}}}}^{0}}\)97 and a positive polynomial function Pn, convolved with the detector resolution function \({\mathfrak{R}}\)

where n denotes the order of the polynomial function, n = 2 is used in the default fit.

The D0D0π+ mass spectrum with non-D0 background subtracted is shown in Fig. 1 with the result of the fit using a model based on the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) signal profile overlaid. The fit gives a signal yield of 186 ± 24 and a mass parameter relative to the D*+D0 mass threshold, δmU of −359 ± 40 keV/c2. The statistical significances of the observed \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) signal and for the δmU < 0 hypothesis are determined using Wilks’ theorem to be 22 and 9 standard deviations, respectively.

The width of the resonance is determined by the coupling constant g for small values of \(\left|g\right|\). With increasing \(\left|g\right|\), the width increases to an asymptotic value determined by the width of the D*+ meson, see Methods and Supplementary Fig. 7. In this regime of large \(\left|g\right|\), the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) signal profile exhibits a scaling property similar to the Flatté function94,98,99. The parameter \(\left|g\right|\) effectively decouples from the fit model, and the model resembles the scattering-length approximation81. The likelihood profile for the parameter \(\left|g\right|\) is shown in Fig. 2, where one can see a plateau at large values. At small values of the \(\left|g\right|\) parameter, \(\left|g\right|\, < \, 1\,{{{{{\rm{GeV}}}}}}\), the likelihood function is independent of \(\left|g\right|\) because the resonance is too narrow for the details of the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) signal profile to be resolved by the detector. The lower limits on the \(\left|g\right|\) parameter of \(\left|g\right| \, > \, 7.7\,(6.2)\,{{{{{\rm{GeV}}}}}}\) at 90% (95%) confidence level (CL) are obtained as the values where the difference in the negative log-likelihood \(-{{\Delta }}\log {{{{{\mathcal{L}}}}}}\) is equal to 1.35 and 1.92, respectively. Smaller values for \(\left|g\right|\) are further used for systematic uncertainty evaluation.

Likelihood profile for the absolute value of the coupling constant g from the fit to the background-subtracted D0D0π+ mass spectrum with a model based on the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) signal profile.

The mode relative to the D*+D0 mass threshold, \(\delta {\mathfrak{m}}\), and the full width at half maximum (FWHM), \({\mathfrak{w}}\), for the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) profile are found to be \(\delta {\mathfrak{m}}=-361\pm 40 \,{{{{{\rm{keV}}}}}}/{c}^{2}\) and \({\mathfrak{w}}=47.8\pm 1.9{{{{{\rm{keV}}}}}}/{c}^{2}\), to be compared with those quantities determined for the \({{\mathfrak{F}}}^{{{{{{\rm{BW}}}}}}}\) signal profile of \(\delta {\mathfrak{m}}=-279\pm 59 \,{{{{{\rm{keV}}}}}}/{c}^{2}\) and \({\mathfrak{w}}=409\pm 163 \,{{{{{\rm{keV}}}}}}/{c}^{2}\). They appear to be rather different. Nonetheless, both functions properly describe the data given the limited sample size, and accounting for the detector resolution, and residual background. To quantify the impact of these experimental effects, two ensembles of pseudoexperiments are performed. Firstly, pseudodata samples are generated with a model based on the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) profile. The parameters used here are obtained from the default fit, and the size of the sample corresponds to the size of data sample. Each pseudodata sample is then analysed with a model based on the \({{\mathfrak{F}}}^{{{{{{\rm{BW}}}}}}}\) function. The obtained mean and RMS values for the parameters δmBW and ΓBW over the ensemble are shown in Table 2. The mass parameter δmBW agrees well with the value determined from data71. The difference for the parameter ΓBW does not exceed one standard deviation. Secondly, an ensemble of pseudodata samples generated with a model based on the \({{\mathfrak{F}}}^{{{{{{\rm{BW}}}}}}}\) profile is analysed with a model based on the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) function. The obtained mean and RMS values for the δmU parameter over an ensemble are also reported in Table 2. These values agree well with the result of the default fit to data. The results of these pseudoexperiments explain the seeming inconsistency between the models and illustrate the importance of an accurate description of the detector resolution and residual background given the limited sample size.

Systematic uncertainties for the δmU parameter are summarised in Table 3 and described in greater detail below. The systematic uncertainty related to the fit model is studied using pseudoexperiments with a set of alternative parameterisations. For each alternative model, an ensemble of pseudoexperiments is performed with parameters obtained from a fit to data. A fit with the baseline model is performed on each pseudoexperiment, and the mean values of the parameters of interest are evaluated over the ensemble. The absolute values of the differences between these mean values and the corresponding parameter values obtained from the fit to data are used to assess the systematic uncertainty due to the choice of the fit model. The maximal value of such differences over the considered set of alternative models is taken as the corresponding systematic uncertainty. The following sources of systematic uncertainty related to the fit model are considered:

Imperfect knowledge of the detector resolution model. To estimate the associated systematic uncertainty a set of alternative resolution functions is tested: a symmetric variant of an Apollonios function100, a modified Gaussian function with symmetric power-law tails on both sides of the distribution101,102, a generalised symmetric Student’s t-distribution103,104, a symmetric Johnson’s SU distribution105,106, and a modified Novosibirsk function107.

A small difference in the detector resolution between data and simulation. A correction factor of 1.05 is applied to account for known discrepancies in modelling the detector resolution in simulation. This factor was studied for different decays94,95,96,108,109,110 and found to lie between 1.0 and 1.1. For decays with relatively low momentum tracks, this factor is close to 1.05, which is the nominal value used in this analysis. This factor is also cross-checked using large samples of D*+ → D0π+ decays, where a value of 1.06 is obtained. To assess the systematic uncertainty related to this factor, detector resolution models with correction factors of 1.0 and 1.1 are studied as alternatives.

Parameterisation of the background component. To assess the associated systematic uncertainty, the order of the positive polynomial function of Eq. (2) is varied. In addition, to estimate a possible effect from a small contribution from three-body D0D0π+ combinations without an intermediate D*+ meson, a more general family of background models is tested

where \({{{\Phi }}}_{{{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}}\) denotes the three-body phase-space function111. The functions B0, B1, B3 and \({B}_{nm}^{\prime}\) with n ≤ 2, m ≤ 1 are used as alternative models for the estimation of the systematic uncertainty.

Values of the coupling constants for the D* → Dπ and D* → Dγ decays affecting the shape of the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) signal profile. These coupling constants are calculated from the known branching fractions of the D* → Dπ and D* → Dγ decays10, the measured natural width of the D*+ meson10,112 and the derived value for the natural width of the D*0 meson66,81,113. To assess the associated systematic uncertainty, a set of alternative models built around the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) profiles, obtained with coupling constants varying within their calculated uncertainties, is studied.

Unknown value of the \(\left|g\right|\) parameter. In the baseline fit the value of the \(\left|g\right|\) parameter is fixed to a large value. To assess the effect of this constraint the fit is repeated using the value of \(\left|g\right|=8.08\,{{{{{\rm{GeV}}}}}}\), that corresponds to \(-2{{\Delta }}\log {{{{{\mathcal{L}}}}}}=1\) for the most conservative likelihood profile for \(\left|g\right|\) that accounts for the systematic uncertainty. The change of 7 keV/c2 of the δmU parameter is assigned as the systematic uncertainty.

The calibration of the momentum scale of the tracking system is based upon large samples of B+ → J/ψK+ and J/ψ → μ+μ− decays114. The accuracy of the procedure has been checked using fully reconstructed B decays together with two-body ϒ(nS) and \({{{{{{\rm{K}}}}}}}_{{{{{{\rm{S}}}}}}}^{0}\) decays and the largest deviation of the bias in the momentum scale of δα = 3 × 10−4 is taken as the uncertainty115. This uncertainty is propagated for the parameters of interest using simulated samples, with momentum scale corrections of \(\left(1\pm \delta \alpha \right)\) applied. Half of the difference between the obtained peak locations is taken as an estimate of the systematic uncertainty.

In the reconstruction the momenta of the charged tracks are corrected for energy loss in the detector material using the Bethe–Bloch formula116,117. The amount of the material traversed in the tracking system by a charged particle is known to have 10% accuracy118. To assess the corresponding uncertainty the magnitude of the calculated corrections is varied by ±10%. Half of the difference between the obtained peak locations is taken as an estimate of the systematic uncertainty due to energy loss corrections.

The mass of D0D0π+ combinations is calculated with the mass of each D0 meson constrained to the known value of the D0 mass10. This procedure produces negligible uncertainties for the δmU parameter due to imprecise knowledge of the D0 mass. However, the small uncertainty of 2 keV/c2 for the known D*+ − D0 mass difference10,112,119 directly affects the values of these parameters and is assigned as a corresponding systematic uncertainty.

For the lower limit on the parameter \(\left|g\right|\), only systematic uncertainties related to the fit model are considered. For each alternative model the likelihood profile curves are built and corresponding 90 and 95% CL lower limits are calculated using the procedure described above. The smallest of the resulting values is taken as the lower limit that accounts for the systematic uncertainty: \(\left|g\right| \, > \, 5.1\,(4.3)\,{{{{{\rm{GeV}}}}}}\) at 90 (95%) CL.

Studying the D0π+ mass distribution for \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decays allows testing the hypothesis that the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decay proceeds through an intermediate off-shell D*+ meson. The background-subtracted D0π+ mass distribution for selected D0D0π+ candidates with the D0D0π+ mass with respect to the D*+D0 mass threshold, \(\delta {m}_{{{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}}\), below zero is shown in Fig. 3. Both D0π+ combinations are included in this plot. The two-dimensional distribution of the mass of one D0π+ combination versus the mass of another D0π+ combination is presented in Supplementary Fig. 10.

Mass distribution for D0π+ pairs from selected D0D0π+ candidates with a mass below the D*+D0 mass threshold with non-D0 background subtracted. The overlaid fit result is described in the text. The background component vanishes in the fit. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

A fit is performed to this distribution with a model containing signal and background components. The signal component is derived from the \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}\) amplitude, see Methods Eq. (48), and is convolved with a detector resolution for the D0π+ mass. This detector resolution function is modelled with a modified Gaussian function with power-law tails on both sides of the distribution101,102 and parameters taken from simulation. Similarly to the correction used for the D0D0π+ mass resolution function \({\mathfrak{R}}\), the width of the Gaussian function is corrected by a factor of 1.06 which is determined by studying large samples of D*+ → D0π+ decays. The RMS of the resolution function is around 220 keV/c2. The shape of the background component is derived from data for \(\delta {m}_{{{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}} > 0.6\,{{{{{\rm{MeV}}}}}}/{c}^{2}\). The fit results are overlaid in Fig. 3. The background component vanishes in the fit, and the D0π+ spectrum is consistent with the hypothesis that the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decay proceeds through an intermediate off-shell D*+ meson. This in turn favours the 1+ assignment for the spin-parity of the state.

Due to the proximity of the observed \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) signal to the D*+D0 mass threshold, and the small energy release in the D*+ → D0π+ decay, the D0D0 mass distribution from the\({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decay forms a narrow peak just above the D0D0 mass threshold. In a similar way, a peaking structure in the D+D0 mass spectrum just above the D+D0 mass threshold is expected from \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{+}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{0}\) and \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{+}{{{{{{\rm{D}}}}}}}^{0}\gamma \) decays, both proceeding via off-shell intermediate D*+D0 and D*0D+ states. The D0D0 and D+D0 final states are reconstructed and selected similarly to the D0D0π+ final state, where the D+ → K−π+π− decay channel is used. The background-subtracted D0D0 and D+D0 mass distributions are shown in Fig. 4 (top), where narrow structures are clearly visible just above the DD thresholds. Fits to these distributions are performed using models consisting of two components: a signal component FDD described in Methods Eqs. (70) and (71) and obtained via integration of the matrix elements for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{\rm{D}}}}}}{{{{{\rm{D}}}}}}\pi /\gamma \) decays with the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) profile, and a background component, parameterised as a product of the two-body phase-space function ΦDD and a positive linear function P1. The fit results are overlaid in Fig. 4 (top). The signal yields in the D0D0 and D+D0 spectra are found to be 263 ± 23 and 171 ± 26, respectively. The statistical significance of the observed \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{{{{{\rm{X}}}}}}\) and \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{+}{{{{{{\rm{D}}}}}}}^{0}{{{{{\rm{X}}}}}}\) signals, where X stands for non-reconstructed pions or photons, is estimated using Wilks’ theorem120 and is found to be in excess of 20 and 10 standard deviations, respectively. The relative yields for the signals observed in the D0D0π+, D0D0 and D0D+ mass spectra agree with the expectations of the model described in Methods where the decay of an isoscalar \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state via the D*D channel with an intermediate off-shell D* meson is assumed.

(Top) DD and DDπ+ mass distributions for selected (left) D0D0 and (right) D+D0 candidates with the non-D background subtracted. The overlaid fit results are described in the text. For visibility the\({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{+}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{0}\) is stacked on top of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{+}{{{{{{\rm{D}}}}}}}^{0}\gamma \) component. (Bottom) Mass distributions for selected (left) D+D+ and (right) D+D0π+ candidates with the non-D background subtracted. The vertical coloured band indicates the expected mass for the hypothetical \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{++}\) state. The overlaid fit results with background-only functions are described in the text. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

The observation of the near-threshold signals in the D0D0 and D+D0 mass spectra, along with the signal shapes and yields, all agree with the isoscalar \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) hypothesis for the narrow signal observed in the D0D0π+ mass spectrum. However, an alternative interpretation could be that this state is the I3 = 0 component of a \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}\) isotriplet (\({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{0}\),\({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\),\({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{++}\)) with \({{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}\overline{{{{{{\rm{u}}}}}}}\overline{{{{{{\rm{u}}}}}}}\), \({{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}\overline{{{{{{\rm{u}}}}}}}\overline{{{{{{\rm{d}}}}}}}\) and \({{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}\overline{{{{{{\rm{d}}}}}}}\overline{{{{{{\rm{d}}}}}}}\) quark content, respectively. Assuming that the observed peak corresponds to the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) component and using the estimates for the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}\) mass splitting from Methods Eqs. (85) and (86), the masses of the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{0}\) and \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{++}\) states are estimated to be slightly below the D0D*0 and slightly above the D+D*+ mass thresholds, respectively:

With these mass assignments, assuming equal production of all three \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}\) components, the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{0}\) state would be an extra narrow state that decays into the D0D0π0 and D0D0γ final states via an off-shell D*0 meson. These decays would contribute to the narrow near-threshold enhancement in the D0D0 spectrum, and increase the signal in the D0D0 mass spectrum by almost a factor of three. The \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{++}\) state would decay via an on-shell D*+ meson \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{++}\,\to {{{{{{\rm{D}}}}}}}^{+}{{{{{{\rm{D}}}}}}}^{* +}\); therefore, it could be a relatively wide state, with a width up to a few Me121. Therefore, it would manifest itself as a peak with a moderate width in the D+D0π+ mass spectrum with a yield comparable to that of the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decays. In addition, it would contribute to the D+D0 mass spectrum, tripling the contribution from the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) decays. However, due to the larger mass of the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{++}\) state and its larger width, this contribution should be wider, making it more difficult to disentangle from the background. Finally, the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{++}\) state would make a contribution to the D+D+ spectrum with a yield similar to the contribution from \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{+}{\pi }^{0}/\gamma \) decays to the D0D+ spectrum, but wider. The mass spectra for D+D+ and D+D0π+ combinations are shown in Fig. 4 (bottom).

Neither distribution exhibits any narrow signal-like structure. Fits to these spectra are performed using the following background-only functions:

The results of these fits are overlaid in Fig. 4 (bottom). The absence of any signals in the D+D+ and D+D0π+ mass spectra is therefore a strong argument in favour of the isoscalar nature of the observed peak in the D0D0π+ mass spectrum.

The interference between two virtual channels for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decay, corresponding to two amplitude terms, see Methods Eq. (35), is studied by setting the term proportional to C in Methods Eq. (39) to be equal to zero. This causes a 43% reduction in the decay rate, pointing to a large interference. The same procedure applied to the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{+}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{0}\) decays gives the contribution of 45% for the interference between the \(\left({{{{{{\rm{D}}}}}}}^{* +}\,\to {{{{{{\rm{D}}}}}}}^{+}{\pi }^{0}\right){{{{{{\rm{D}}}}}}}^{0}\) and \(\left({{{{{{\rm{D}}}}}}}^{* 0}\,\to {{{{{{\rm{D}}}}}}}^{0}{\pi }^{0}\right){{{{{{\rm{D}}}}}}}^{+}\) channels. For \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{+}{{{{{{\rm{D}}}}}}}^{0}\gamma \) decays the role of the interference between the \(\left({{{{{{\rm{D}}}}}}}^{* +}\,\to {{{{{{\rm{D}}}}}}}^{+}\gamma \right){{{{{{\rm{D}}}}}}}^{0}\) and \(\left({{{{{{\rm{D}}}}}}}^{* 0}\,\to {{{{{{\rm{D}}}}}}}^{0}\gamma \right){{{{{{\rm{D}}}}}}}^{+}\) channels is estimated by equating to zero the \({{\mathfrak{F}}}_{+}{{\mathfrak{F}}}_{0}^{* }\) and \({{\mathfrak{F}}}_{+}^{* }{{\mathfrak{F}}}_{0}\) terms in Methods Eqs. (45) and (46). The interference contribution is found to be 33%.

Using the model described earlier and results of the fit to the D0D0π+ mass spectrum, the position of the amplitude pole \(\hat{s}\) in the complex plane, responsible for the appearance of the narrow structure in the D0D0π+ mass spectrum is determined. The pole parameters, mass mpole and width Γpole, are defined through the pole location \(\hat{s}\) as

The pole location \(\hat{s}\) is a solution to the equation

where \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}^{II}(s)\) denotes the amplitude on the second Riemann sheet defined in Methods Eq. (64). For large coupling \(\left|g\right|\) the position of the resonance pole is uniquely determined by the parameter δmU, i.e., the binding energy and the width of the D*+ meson. Figure 5 shows the complex plane of the \(\delta \sqrt{s}\) variable, defined as

All possible positions of the pole for \(\left|g\right|\gg {m}_{{{{{{{\rm{D}}}}}}}^{0}}+{m}_{{{{{{{\rm{D}}}}}}}^{* +}}\) are located on a red dashed curve in Fig. 5. The behaviour of the curve can be understood as follows: with an increase of the binding energy (distance to the D*+D0 mass threshold), the width gets narrower; and when the parameter δmU approaches zero, the pole touches the D0D*+ cut and moves to the other complex sheet, i.e., the state becomes virtual. For smaller values of \(\left|g\right|\), the pole is located between the limiting curve and the ℜ s = 0 line. The pole parameters are found to be

where the first uncertainty is due to the δmU parameter and the second is due to the unknown value of the \(\left|g\right|\) parameter. The peak is well separated from the D*+D0 threshold in the D0D0π+ mass spectrum. Hence, as for an isolated narrow resonance, the parameters of the pole are similar to the visible peak parameters, namely the mode \(\delta {\mathfrak{m}}\) and FWHM \({\mathfrak{w}}\).

Complex plane of the \(\delta \sqrt{s}\) variable. The dashed red line shows the allowed region for large \(\left|g\right|\) values. The filled red circle indicates the best estimate for the pole location and the filled regions show 1σ and 2σ confidence regions. Open blue and green circles show the branch points corresponding to the D*+D0 and D*0D+ channels, respectively, and the corresponding blue and green lines indicate branch cuts. Three other branch points at \(\sqrt{s}\) of \({m}_{{{{{{{\rm{D}}}}}}}^{0}}+{m}_{{{{{{{\rm{D}}}}}}}^{+}}\), \({m}_{{{{{{{\rm{D}}}}}}}^{0}}+{m}_{{{{{{{\rm{D}}}}}}}^{+}}+{m}_{{\pi }^{0}}\) and \(2{m}_{{{{{{{\rm{D}}}}}}}^{0}}+{m}_{{\pi }^{+}}\), corresponding to the openings of the D0D+γ, D0D+π0 and D0D0π+ decay channels, are outside of the displayed region.

The systematic uncertainties quoted here do not account for the possibility that any of the underlying assumptions on which the model is built are not valid. For example, as shown earlier the data are consistent with a wide range of FHWM \({\mathfrak{w}}\) values for the signal profile. Therefore the pole width Γpole is based mainly on the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) amplitude model and the value of the mU parameter determined from the fit to the D0D0π+ mass spectrum.

A study of the behaviour of the \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}(s)\) amplitude in the vicinity of the D*+D0 mass threshold leads to the determination of the low-energy scattering parameters, namely the scattering length, a, and the effective range, r. These parameters are defined via the coefficients of the first two terms of the Taylor expansion of the inverse non-relativistic amplitude122, i.e.,

where k is the wave number. For \(\delta \sqrt{s}\, \lesssim -{{{\Gamma }}}_{{{{{{{\rm{D}}}}}}}^{* +}}\) the inverse amplitude from Eq. (47) matches Eq. (13) up to a scale parameter obtained numerically, see Methods Eq. (74). The value of the scattering length is found to be

Typically, a non-vanishing imaginary part of the scattering length indicates the presence of inelastic channels123; however, in this case the non-zero imaginary part is related to the lower threshold, \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\), and is determined by the width of the D*+ meson. The real part of the scattering length a is negative indicating attraction. This can be interpreted as the characteristic size of the state16,

For the \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}\) amplitude the effective range r is non-positive and proportional to \({\left|g\right|}^{-2}\), see Methods Eq. (76). Its value is consistent with zero for the baseline fit. An upper limit on the −r value is set as

The Weinberg compositeness criterion124,125 makes use of the relation between the scattering length and the effective range to construct the compositeness variable Z,

for which Z = 1 corresponds to a compact state that does not interact with the continuum, while Z = 0 indicates a composite state formed by compound interaction. Using the relation between r and ∣g∣ from Methods Eq. (76), one finds \(Z\propto {\left|g\right|}^{-2}\) for large values of \(\left|g\right|\). The default fit corresponds to large values of \(\left|g\right|\), and thus, Z approaching to zero. A non-zero value of Z would require a smaller value of \(\left|g\right|\), i.e., smaller resonance width, see Supplementary Fig. 7. The following upper limit of the compositeness parameter Z is set:

Another estimate of the characteristic size is obtained from the value of the binding energy ΔE. Within the interpretation of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state as a bound D*+D0 molecular-like state, the binding energy is ΔE = − δmU. The characteristic momentum scale γ16 is estimated to be

where μ is the reduced mass of the D*+D0 system. This value of the momentum scale in turn corresponds to a characteristic size RΔE of the molecular-like state,

which is consistent with the Ra estimate from the scattering length.

For high-energy hadroproduction of a state with such a large size, Ra or RΔE, one expects a strong dependency of the production rate on event multiplicity, similar to that observed for the χc1(3872) state126. The background-subtracted distribution of the number of tracks reconstructed in the vertex detector, Ntracks, is shown in Fig. 6 (left) together with the distributions for low-mass \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs with \({m}_{{{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}}\, < \, 3.87\,{{{{{\rm{GeV}}}}}}/{c}^{2}\) and low-mass D0D0 pairs with mass \(3.75 < {m}_{{{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}}\, < \, 3.87\,{{{{{\rm{GeV}}}}}}/{c}^{2}\). The former is dominated by \({{{{{\rm{p}}}}}}{{{{{\rm{p}}}}}}\,\to {{{{{\rm{c}}}}}}{\overline{{{{{{\rm{c}}}}}}}}{{{{{\rm{X}}}}}}\) production, while the latter is presumably dominated by the double-parton scattering process74,127. The chosen interval for \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs includes the region populated by the \({\chi }_{{{{{{\rm{c}}}}}}1}(3872)\,\to {{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}{\pi }^{0}/\gamma \) decays; however, this contribution is small, see Fig. 7. The χc1(3872) production cross-section is suppressed with respect to the conventional charmonium state ψ(2S) at large track multiplicities126. It is noteworthy that the track multiplicity distribution for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state differs from that of the low-mass \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs, in particular, no suppression at large multiplicity is observed. A p value for the consistency of the track multiplicity distributions for \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) production and low-mass \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs is found to be 0.1%. It is interesting to note that the multiplicity distribution for \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) production and the one for D0D0-pairs with \(3.75\, < \, {m}_{{{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}}\, < \, 3.87\,{{{{{\rm{GeV}}}}}}/{c}^{2}\) are consistent with a corresponding p value of 12%. The similarity between \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) production, which is inherently a single parton scattering process, and the distribution for process dominated by a double-parton scattering is surprising.

(Left) Background-subtracted distributions for the multiplicity of tracks reconstructed in the vertex detector for (red circles) \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) signal, low-mass (blue open squares) \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) and (green filled diamonds) D0D0 pairs. The binning scheme is chosen to have an approximately uniform distribution for \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs. The distributions for the \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) and D0D0 pairs are normalised to the same yields as the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) signal. (right) Background-subtracted transverse momentum spectra for (red circles) \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) signal, (blue open squares) low-mass\({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) and (green filled diamonds) D0D0 pairs. The binning scheme is chosen to have an approximately uniform distribution for \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs. The distributions for the \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) and D0D0 pairs are normalised to the same yields as \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) signal. For better visualisation, the points are slightly displaced from the bin centres. For better visualisation, the points are slightly displaced from the bin centres. Uncertainties on the data points are statistical only and represent one standard deviation.

Background-subtracted D0D0 and \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) mass distributions. The near-threshold enhancement in the D0D0 channel corresponds to partially reconstructed \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decays, while in the\({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) channel the threshold enhancement corresponds to partially reconstructed \({\chi }_{{{{{{\rm{c}}}}}}1}(3872)\to {{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}{\pi }^{0}\) decays. The \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) mass distribution is zero-suppressed for better visualisation. Uncertainties on the data points are statistical only and represent one standard deviation, calculated as a sum in quadrature of the assigned weights from the background-subtraction procedure.

The transverse momentum spectrum for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state is compared with those for the low-mass \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) and D0D0 pairs in Fig. 6 (right). The p values for the consistency of the pT spectra for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state and low-mass \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs are 1.4%, and 0.02% for low-mass D0D0 pairs. More data are needed for further conclusions.

The background-subtracted D0D0 mass distribution in a wider mass range is shown in Fig. 7 together with a similar distribution for \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs. In the D0\({\overline{{{{{{\rm{D}}}}}}}}^{0}\) mass spectrum the near-threshold enhancement is due to \({\chi }_{{{{{{\rm{c}}}}}}1}(3872)\,\to {{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}{\pi }^{0}\) and \({\chi }_{{{{{{\rm{c}}}}}}1}(3872)\,\to {{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\gamma \) decays via intermediate D*0 mesons77. This structure is significantly wider than the structure in the D0D0 mass spectrum from \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decays primarily due to the larger natural width and smaller binding energy for the χc1(3872) state94,95. With more data, and with a better understanding of the dynamics of \({\chi }_{{{{{{\rm{c}}}}}}1}(3872)\,\to {{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}{\pi }^{0}/\gamma \) decays, and therefore of the corresponding shape in the D0\({\overline{{{{{{\rm{D}}}}}}}}^{0}\) mass spectrum, it will be possible to estimate the relative production rates for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) and χc1(3872) states. Background-subtracted D0D0π+ and D0D+ mass distributions together with those for \({\overline{{{{{{\rm{D}}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) and D0D− are shown in Supplementary Figs. 3 and 4.

The exotic narrow tetraquark state \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) observed in Ref. 71 is studied using a dataset corresponding to an integrated luminosity of 9 fb−1, collected by the LHCb experiment in pp collisions at centre-of-mass energies of 7, 8 and 13 TeV. The observed D0π+ mass distribution indicates that the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decay proceeds via an intermediate off-shell D*+ meson. Together with the proximity of the state to the D*D0 mass threshold, this favours the spin-parity quantum numbers JP to be 1+. Narrow near-threshold structures are observed in the D0D0 and D0D+ mass spectra with high significance. These are found to be consistent with originating from off-shell \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{* }{{{{{\rm{D}}}}}}\) decays followed by the D* → Dπ and D* → Dγ decays. No signal is observed in the D+D0π+ mass spectrum, and no structure is observed in the D+D+ mass spectrum. These non-observations provide a strong argument in favour of the isoscalar nature for the observed state, supporting its interpretation as the isoscalar JP = 1+ \({{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}\overline{{{{{{\rm{u}}}}}}}\overline{{{{{{\rm{d}}}}}}}\)-tetraquark ground state. A dedicated unitarised three-body Breit–Wigner amplitude is built on the assumption of strong isoscalar coupling of the axial-vector \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state to the D*D channel. This assumption is supported by the data; however, alternative models are not excluded by the distributions studied in this analysis. Probing alternative models and the validity of the underlying assumptions of this analysis will be a subject for future studies.

Using the developed amplitude model, the mass of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state, relative to the D*+D0 mass threshold, is determined to be

where the first uncertainty is statistic and the second systematic. The lower limit on the absolute value of the coupling constant of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state to the D*D system is

Using the same model, the estimates for the scattering length a, effective range r, and the compositeness, Z are obtained from the low-energy limit of the amplitude to be

The characteristic size calculated from the binding energy is RΔE = 7.49 ± 0.42 fm. This value is consistent with the estimation from the scattering length, Ra = 7.16 ± 0.51 fm. Both RΔE and Ra correspond to a spatial extension significantly exceeding the typical scale for heavy-flavour hadrons. Within this model the resonance pole is found to be located on the second Riemann sheet with respect to the D0D0π+ threshold, at \(\hat{s}={m}_{{{{{{\rm{pole}}}}}}}-\frac{i}{2}{{{\Gamma }}}_{{{{{{\rm{pole}}}}}}}\), where

where the first uncertainty accounts for statistical and systematic uncertainties for the δmU parameters, and the second is due to the unknown value of the \(\left|g\right|\) parameter. The pole position, scattering length, effective range and compositeness form a complete set of observables related to the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) reaction amplitude, which are crucial for inferring the nature of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) tetraquark.

Unlike in the prompt production of the χc1(3872) state, no suppression of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) production at high track multiplicities is observed relative to the low-mass \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\) pairs. The observed similarity with the multiplicity distribution for the low-mass D0D0 production process, that is presumably double-parton-scattering dominated, is unexpected. In the future with a larger dataset and including other decay modes, e.g., D0 → K−π+π+π−, detailed studies of the properties of this new state and its production mechanisms could be possible.

In conclusion, the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) tetraquark observed in D0D0π+ decays is studied in detail, using a unitarised model that accounts for the relevant thresholds by taking into account the D0D0π+ and D0D+π0(γ) decay channels with intermediate D* resonances. This model is found to give an excellent description of the D0π+ mass distribution in the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decay and of the threshold enhancements observed in the D0D0 and D0D+ spectra. Together with the absence of a signal in the D0D+ and D+D0π+ mass distributions this provides a strong argument for interpreting the observed state as the isoscalar \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) tetraquark with spin-parity JP = 1+. The precise \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) mass measurement will rule out or improve on a considerable range of theoretical models on heavy quark systems. The determined pole position and physical quantities derived from low-energy scattering parameters reveal important information about the nature of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) tetraquark. In addition, the counter-intuitive dependence of the production rate on track multiplicity will pose a challenge for theoretical explanations.

The LHCb detector72,73 is a single-arm forward spectrometer covering the pseudorapidity range 2 < η < 5, designed for the study of particles containing b or c quarks. The detector includes a high-precision tracking system consisting of a silicon-strip vertex detector surrounding the pp interaction region, a large-area silicon-strip detector located upstream of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip detectors and straw drift tubes placed downstream of the magnet. The tracking system provides a measurement of the momentum, p, of charged particles with a relative uncertainty that varies from 0.5% at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary pp collision vertex, the impact parameter, is measured with a resolution of (15 + 29/pT) μm, where pT is the component of the momentum transverse to the beam, in GeV/c. Different types of charged hadrons are distinguished using information from two ring-imaging Cherenkov detectors128. Photons, electrons and hadrons are identified by a calorimeter system consisting of scintillating-pad and preshower detectors, an electromagnetic and a hadronic calorimeter. Muons are identified by a system composed of alternating layers of iron and multiwire proportional chambers. The online event selection is performed by a trigger, which consists of a hardware stage, based on information from the calorimeter and muon systems, followed by a software stage, which applies a full event reconstruction.

Simulation is required to model the effects of the detector acceptance, resolution, and the efficiency of the imposed selection requirements. In the simulation, pp collisions are generated using PYTHIA129 with a specific LHCb configuration130. Decays of unstable particles are described by EVTGEN131, in which final-state radiation is generated using PHOTOS132. The interaction of the generated particles with the detector, and its response, are implemented using the GEANT4 toolkit133,134 as described in Ref. 135.

The D0D0, D0D+ and D0D0π+ final states are reconstructed using the D0 → K−π+ and D+ → K−π+π+ decay channels. The selection criteria are similar to those used in Refs. 74,75,76,77. Kaons and pions are selected from well-reconstructed tracks within the acceptance of the spectrometer that are identified using information from the ring-imaging Cherenkov detectors. The kaon and pion candidates that have transverse momenta larger than 250 MeV/c and are inconsistent with being produced at a pp interaction vertex are combined together to form D0 and D+ candidates, referred to as D hereafter. The resulting D candidates are required to have good vertex quality, mass within ±65 and ±50 MeV/c2 of the known D0 and D+ masses10, respectively, transverse momentum larger than 1 GeV/c, decay time larger than 100 μm/c and a momentum direction that is consistent with the vector from the primary to secondary vertex. Selected D0 and D+ candidates consistent with originating from a common primary vertex are combined to form D0D0 and D0D+ candidates. The resulting D0D0 candidates are combined with a pion to form D0D0π+ candidates. At least one of the two D0π+ combinations is required to have good vertex quality and mass not exceeding the known D*+ mass by more than 155 MeV/c2. For each D0D0, D0D+ and D0D0π+ candidate a kinematic fit136 is performed. This fit constrains the mass of the D candidates to their known values and requires both D mesons, and a pion in the case of D0D0π+, to originate from the same primary vertex. A requirement is applied to the quality of this fit to further suppress combinatorial background and reduce background from D candidates produced in two independent pp interactions or in the decays of beauty hadrons74. To suppress background from kaon and pion candidates reconstructed from a common track, all track pairs of the same charge are required to have an opening angle inconsistent with zero and the mass of the combination must be inconsistent with the sum of the masses of the two constituents. For cross-checks additional final states D+D+, D+D0π+, \({{{{{{\rm{D}}}}}}}^{0}{\overline{{{{{{\rm{D}}}}}}}}^{0}\), D0D− and \({\overline{{{{{{\rm{D}}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) are reconstructed, selected and treated in the same way.

Two-dimensional distributions of the mass of one D candidate versus the mass of the other D candidate from selected D0D0π+, D0D0 and D0D+ combinations are shown in Supplementary Fig. 1. These distributions illustrate the relatively small combinatorial background levels due to fake D candidates. This background is subtracted using the sPlot technique137, which is based on an extended unbinned maximum-likelihood fit to these two-dimensional distributions with the function described in Ref. 74.

This function consists of four components:

a component corresponding to genuine D1D2 pairs and described as a product of two signal functions, each parameterised with a modified Novosibirsk function107;

two components corresponding to combinations of one of the D mesons with combinatorial background, described as a product of the signal function and a background function, which is parameterised with a product of an exponential function and a positive first-order polynomial;

a component corresponding to pure background pairs and described by a product of exponential functions and a positive two-dimensional non-factorisable second-order polynomial function.

Based on the results of the fit, each candidate is assigned a positive weight for being signal-like or a negative weight for being background-like, with the masses of the two D0 candidates as discriminating variables. The D0D0π+ mass distributions for each of the subtracted background components are presented in Supplementary Fig. 2. where fit results with background-only functions \({B}_{10}^{\prime}\), defined in Eq. (3) are overlaid.

In the vicinity of the D*+D0 mass threshold the resolution function \({\mathfrak{R}}\) for the D0D0π+ mass is parametrised with the sum of two Gaussian functions with a common mean. The widths of the Gaussian functions are σ1 = 1.05 × 263 keV/c2 and σ2 = 2.413 × σ1 for the narrow and wide components, respectively, and the fraction of the narrow Gaussian is α = 0.778. The parameters α and σ1,2 are taken from simulation, and σ1,2 are corrected with a factor of 1.05 that accounts for a small difference between simulation and data for the mass resolution94,95,96. The RMS of the resolution function is around 400 keV/c2.

Assuming isospin symmetry, the isoscalar vector state \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) that decays into the D*D final state can be expressed as

Therefore, the S-wave amplitudes for the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{* +}{{{{{{\rm{D}}}}}}}^{0}\) and \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{* 0}{{{{{{\rm{D}}}}}}}^{+}\) decays have different signs

where g is a coupling constant, \({\epsilon }_{{{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}}\) is the polarisation vector of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) particle and \({\epsilon }_{{{{{{{\rm{D}}}}}}}^{* }}\) is the polarisation vector of D* meson, and the upper and lower Greek indices imply the summation in the Einstein notation. The S-wave (corresponding to orbital angular momentum equal to zero) approximation is valid for a near-threshold peak. For \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) masses significantly above the D*D threshold, higher-order waves also need to be considered. The amplitudes for the D* → Dπ decays are written as

where f denotes a coupling constant, and pD stands for the momentum of the D meson. The amplitude for the D* → Dγ decays is

where h denotes a coupling constant, μ stands for the magnetic moment for D* → Dγ transitions, \({p}_{{{{{{{\rm{D}}}}}}}^{* }}\) and pγ are the D*-meson and photon momenta, respectively, and ϵγ is the polarisation vector of the photon. The three amplitudes for \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to \pi {{{{{\rm{D}}}}}}{{{{{\rm{D}}}}}}\) and \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to \gamma {{{{{\rm{D}}}}}}{{{{{\rm{D}}}}}}\) decays are

where \({s}_{ij}={p}_{ij}^{2}={({p}_{i}+{p}_{j})}^{2}\) and the \({\mathfrak{F}}\) functions that denote the Breit–Wigner amplitude for the D* mesons are

A small possible distortion of the Breit–Wigner shape of the D* meson due to three-body final-state interactions is neglected in the model. The impact of the energy-dependence of the D* meson self-energy is found to be insignificant. The decays of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state into the D+D+π− final state via off-shell D*0 → D+π− decays are highly suppressed and are not considered here. The last terms in Eqs. (35) and (36) imply the same amplitudes with swapped momenta.

The \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state is assumed to be produced unpolarised; therefore, the squared absolute value of the decay amplitudes with pions in the final state, averaged over the initial spin state are

where

and \(\lambda \left(x,y,z\right)\) stands for the Källén function97. The additional factor of 2! in the denominator of Eq. (39) is due to the presence of two identical particles (D0) in the final state. The squared absolute values of the decay amplitude with a photon in the final state, averaged over the initial spin state is

The coupling constants f and h for the D* → Dπ and D* → Dγ decays are calculated using Eqs. (31)–(34), from the known branching fractions of the D* → Dπ and D* → Dγ decays10, the measured natural width of the D*+ meson10,112 and the derived value for the natural width for the D*0 meson66,81,113. The magnetic moment μ+ is taken to be 1 and the ratio of magnetic moments μ0/μ+ is calculated according to Refs. 138,139,140.

A unitarised three-body Breit–Wigner function is defined as

where \(f\in \left\{{{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+},{{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{+}{\pi }^{0},{{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{+}\gamma \right\}\) denotes the final state. The decay matrix element for each channel integrated over the three-body phase-space is denoted by

where \({\left|{{\mathfrak{M}}}_{f}\right|}^{2}\) is defined by Eqs. (39)–(46) and the unknown coupling constant g is taken out of the expression for \({\left|{{\mathfrak{M}}}_{f}\right|}^{2}\). For large values of s, in excess of s*, such as \(\sqrt{{s}^{* }}-\left({m}_{{{{{{{\rm{D}}}}}}}^{* }}+{m}_{{{{{{\rm{D}}}}}}}\right)\gg {{{\Gamma }}}_{{{{{{{\rm{D}}}}}}}^{* }}\), the functions ϱf(s) are defined as

where \({{{\Phi }}}_{{{{{{{\rm{D}}}}}}}^{* }{{{{{\rm{D}}}}}}}(s)\) denotes the two-body phase-space function, the constants c1,c2 and c3 are chosen to ensure the continuity of the functions ϱf(s), and a value of \(\sqrt{{s}^{* }}=3.9\,{{{{{\rm{GeV}}}}}}/{c}^{2}\) is used. The functions ϱf(s) are shown in Supplementary Fig. 5. The complex-valued width \(\hat{{{\Gamma }}}(s)\) is defined via the self-energy function Σ(s)141

where \({\left|g\right|}^{2}\) is again factored out for convenience. The imaginary part of Σ(s) for real physical values of s is computed through the optical theorem as half of the sum of the decay probability to all available channels142:

The real part of the self-energy function is computed using Kramers–Kronig dispersion relations with a single subtraction143,144,

where the Cauchy principal value (p.v.) integral over ϱtot(s) is understood as

and sf denotes the threshold value for the channel f. The subtraction is needed since the integral ∫ϱtot(s)/s ds diverges. The term \(\xi ({m}_{{{{{{\rm{U}}}}}}}^{2})\) in Eq. (56) corresponds to the choice of subtraction constant such that \(\Re \,{{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}({m}_{{{{{{\rm{U}}}}}}}^{2})=0\). The function ξ(s) is shown in Supplementary Fig. 6.

Alternatively, the isoscalar amplitude \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}\) is constructed using the K-matrix approach145 with two coupled channels, D*+D0 and D*0D+. The relation reads:

where a production vector P and an isoscalar potential K are defined as

The propagation matrix G describes the D*D → D*D rescattering via the virtual loops including the one-particle exchange process91 and expressed in a symbolical way in Supplementary Eq. (1). where suppressed D*0 → D+π− transition is neglected. The D*+D0 ↔ D*0D+ rescattering occurs due to non-diagonal element of the K-matrix (contact interaction) and non-diagonal elements of the G matrix (long-range interaction). The matrix G and the self-energy function Σ(s) from Eqs. (54) and (56), are related as

Similar to the Flatté function98 for large values of the \(\left|g\right|\) parameter, the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) signal profile exhibits a scaling property94,99. For large values of the \(\left|g\right|\) parameter the width approaches asymptotic behaviour, see Supplementary Fig. 7. The unitarised three-body Breit–Wigner function \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}\) for \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{{\rm{D}}}}}}}^{0}{{{{{{\rm{D}}}}}}}^{0}{\pi }^{+}\) decays with parameters mU and \(\left|g\right|\) obtained from the fit to data is shown in Supplementary Fig. 8. The inset illustrates the similarity of the profile with the single-pole profile in the vicinity of the pole

where \({\mathfrak{m}}\) and \({\mathfrak{w}}\) are the mode and FWHM, respectively.

Equation (56) defines Σ(s) and the amplitude \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}(s)\) for real values of s. Analytic continuation to the whole first Riemann sheet is calculated as

where the integral is understood as in Eq. (58). The search for the resonance pole requires knowledge of the amplitude on the second Riemann sheet denoted by \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}^{II}\). According to the optical theorem142, the discontinuity of the inverse amplitude across the unitarity cut is given by \(i{\left|g\right|}^{2}{\varrho }_{{{{{{\rm{tot}}}}}}}(s)\):

For the complex-values s, the analytic continuation of ϱtot(s) is needed: the phase-space integral in Eq. (49) is performed over a two-dimensional complex manifold \({{{{{\mathcal{D}}}}}}\) (see discussion on the continuation in Ref. 146):

where the limits of the second integral represent the Dalitz plot borders111,

The integration is performed along straight lines connecting the end points in the complex plane.

The shapes of the D0D0 and D+D− mass spectra from \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\,\to {{{{{\rm{D}}}}}}{{{{{\rm{D}}}}}}\pi /\gamma \) decays are obtained via integration of the \({|{{\mathfrak{M}}}_{f}|}^{2}\) expressions from Eqs. (39)–(46) over the s and s12 variables with the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) amplitude squared, \({\left|{{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}\left(s\right)\right|}^{2}\), from Eq. (47):

where the lower and upper integration limits for s12 at fixed s and s23 are given in Eq. (66). The function fC(s) is introduced to perform a smooth cutoff of the long tail of the \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) profile. Cutoffs are chosen to suppress the profile for regions \(\left|\sqrt{s}-{\mathfrak{m}}\right|\gg {\mathfrak{w}}\), where \({\mathfrak{m}}\) and \({\mathfrak{w}}\) are the mode and FWHM for the \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}(s)\) distribution. Two cutoff functions fC(s) are studied:

A Gaussian cutoff \({f}_{C}^{G}(s)\) is defined as

A power-law cutoff \({f}_{C}^{P}(s)\) defined as

Fits to the background-subtracted D0D0π+ mass spectrum using a signal profile of the form \({{\mathfrak{F}}}^{{{{{{\rm{U}}}}}}}(s)\times {f}_{C}(s)\) show that the parameter δmU is insensitive to the choice of cutoff function when \({x}_{c}\ge {m}_{{{{{{{\rm{D}}}}}}}^{* 0}}+{m}_{{{{{{{\rm{D}}}}}}}^{+}}\) and σc≥1 MeV/c2. The power-law cutoff function \({f}_{C}^{P}(s)\) with parameters \({x}_{c}={m}_{{{{{{{\rm{D}}}}}}}^{* 0}}+{m}_{{{{{{{\rm{D}}}}}}}^{+}}\) and σc = 1 MeV/c2 is chosen. The shapes for the D0D0 and D+D0 mass distributions are defined as

The unitarized Breit–Wigner amplitude is formally similar to the low-energy expansion given by Eq. (13) once the factor \(\frac{1}{2}{\left|g\right|}^{2}\) is divided out

The function iϱtot(s) matches ik up to a slowly varying energy factor that can be approximated by a constant in the threshold region. The proportionality factor w has the dimension of an inverse mass and is found by matching the decay probability to the two-body phase-space expression:

where c1 is a coefficient computed in Eq. (50). The comparison of \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{NR}}}}}}}^{-1}\) and \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}^{-1}\times 2w/{\left|g\right|}^{2}\) that validates the matching is shown in Supplementary Fig. 9.

The inverse scattering length is defined as the value of the amplitude in Eq. (72) at the D*+D0 threshold:

The imaginary part is fully determined by the available decay channels, while the real part depends on the constant \(\xi ({m}_{{{{{{\rm{U}}}}}}}^{2})\) adjusted in the fit. The quadratic term, k2 in Eq. (72), corresponds to the linear correction ins since k2 = (s − sth)/4 for the non-relativistic case. Hence, the slope of the linear term in the \({{{{{{\mathcal{A}}}}}}}_{{{{{{\rm{U}}}}}}}^{-1}\) amplitude is related to the effective range as follows:

While the degrees of freedom of the light diquark for the isoscalar \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\) state are similar to those for the \({\overline{\Lambda }}_{{{{{{\rm{c}}}}}}}^{-}\) state, for the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}\) isotriplet (\({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{0}\),\({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\),\({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{++}\)) the light diquark degrees of freedom would be similar to those for the \({\overline{\Sigma }}_{{{{{{\rm{c}}}}}}}\) (anti)triplet. Assuming that the difference in the light quark masses, the Coulomb interaction of light quarks in the diquark, and the Coulomb interaction of the light diquark with the c-quark are responsible for the observed mass splitting in the Σc isotriplet, the masses for the Σc states can be written as

where mΣ is a common mass parameter; the second and third terms describe the contribution from the light quark masses, mu and md, into the mass splitting; terms proportional to a describe Coulomb interactions of light quarks in the diquark; terms proportional to b describe the Coulomb interactions of the diquark with the c-quark; and qq denotes the charge of the q-quark. Similar expressions can be written for the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}\) isotriplet:

where \({m}_{{\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}}\) is the common mass parameter, \({q}_{\overline{{{{{{\rm{q}}}}}}}}=-{q}_{{{{{{\rm{q}}}}}}}\) and qcc = 2qc is the charge of a cc diquark. Using the known masses of the light quarks and Σc states10 and taking \({a}^{\prime}=a\) and \({b}^{\prime}=b\), the mass splitting for the \({\hat{{{{{{\rm{T}}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}\) isotriplet is estimated to be

The validity of this approach is tested by comparing the calculated mass splitting between Σbp and Σbm states of −6.7 ± 0.7 MeV/c2 with the measured value of −5.1 ± 0.2 MeV/c210. Based on the small observed difference, an additional uncertainty of 0.8 MeV/c2 is added in quadrature to the results from Eqs. (83) and (84), and finally one gets

These results agree with the assigned uncertainty with results based on a more advanced model from Ref. 147.

LHCb data used in this analysis will be released according to the LHCb external data access policy, which can be downloaded from https://opendata.cern.ch/record/410/files/LHCb-Data-Policy.pdf. The raw data in all of the figures of this manuscript can be downloaded from https://cds.cern.ch/record/2780001, where no access codes are required. In addition, the unbinned background-subtracted data, shown in Figs. 1, 3 and 4 have been added to the HEPData record at https://www.hepdata.net/record/ins1915358.

LHCb software used to process the data analysed in this manuscript is available at GitLab repository. The specific software used in data analysis is available at Zenodo repository https://doi.org/10.5281/zenodo.5595937.

Gell-Mann, M. A schematic model of baryons and mesons. Phys. Lett. 8, 214 (1964).

Article ADS Google Scholar

Zweig, G. An SU3 Model for Strong Interaction Symmetry and its Breaking; Version 1. CERN-TH-401 (CERN, Geneva). http://cds.cern.ch/record/352337 (1964).

Zweig, G. An SU3 Model for Strong Interaction Symmetry and its Breaking; Version 2. CERN-TH-412 (CERN, Geneva). http://cds.cern.ch/record/570209 (1964).

Jaffe, R. L. Multiquark hadrons. I. Phenomenology of \({{{{{{{{\rm{Q}}}}}}}}}^{2}{\bar{{{{{{{{\rm{Q}}}}}}}}}}^{2}\) mesons. Phys. Rev. D15, 267 (1977).

ADS Google Scholar

Jaffe, R. L. Multi-quark hadrons. 2. Methods. Phys. Rev. D15, 281 (1977).

ADS Google Scholar

Rossi, G. C. & Veneziano, G. A possible description of baryon dynamics in dual and gauge theories. Nucl. Phys. B123, 507 (1977).

Article ADS CAS Google Scholar

Jaffe, R. L. \({{{{{{{{\rm{Q}}}}}}}}}^{2}{\bar{{{{{{{{\rm{Q}}}}}}}}}}^{2}\) resonances in the baryon-antibaryon system. Phys. Rev. D17, 1444 (1978).

ADS Google Scholar

Lipkin, H. J. New possibilities for exotic hadrons – anticharmed strange baryons. Phys. Lett. B195, 484 (1987).

Article ADS Google Scholar

Belle Collaboration, Choi, S-K. et al. Observation of a narrow charmoniumlike state in exclusive B± → K±π+π−J/ψ decays. Phys. Rev. Lett. 91262001. https://doi.org/10.1103/PhysRevLett.91.262001 (2003).

Particle Data Group, Zyla, P. A. et al. Review of particle physics. Prog. Theor. Exp. Phys. 2020, 083C01. http://pdglive.lbl.gov/ (2020).

Chen, H.-X., Chen, W., Liu, X. & Zhu, S.-L. The hidden-charm pentaquark and tetraquark states. Phys. Rept. 639, 1 (2016).

Article ADS MathSciNet CAS Google Scholar

Esposito, A., Pilloni, A. & Polosa, A.D. Multiquark resonances. Phys. Rept. 668, 1 (2017).

Article ADS MathSciNet CAS Google Scholar

Ali, A., Lange, J.S. & Stone, S. Exotics: heavy pentaquarks and tetraquarks. Prog. Part. Nucl. Phys. 97, 123 (2017).

Article ADS CAS Google Scholar

Hosaka, A. et al. Exotic hadrons with heavy flavors: X, Y, Z, and related states. PTEP 2016, 062C01 (2016).

Lebed, R.F., Mitchell, E.S. & Swanson, E.S. Heavy-quark QCD exotica. Prog. Part. Nucl. Phys. 93, 143 (2017).

Article ADS CAS Google Scholar

Guo, F.-K. et al. Hadronic molecules. Rev. Mod. Phys. 90, 015004 (2018).

Article ADS CAS Google Scholar

Olsen, S.L., Skwarnicki, T. & Zieminska, D. Nonstandard heavy mesons and baryons: experimental evidence. Rev. Mod. Phys. 90, 015003 (2018).

Article ADS MathSciNet CAS Google Scholar

Brambilla, N. et al. The XYZ states: experimental and theoretical status and perspectives. Phys. Rept. 873, 1 (2020).

Article ADS CAS Google Scholar

Ali, A.L., Maiani, & Polosa, A. D. Multiquark Hadrons. Cambridge University Press. https://doi.org/10.1017/9781316761465 (2019).

Richard, J.M. Exotic hadrons: review and perspectives. Few-Body Systems 57, 1185 (2016).

Article ADS Google Scholar

Oset, E. et al. Tetra and pentaquarks from the molecular perspective. EPJ Web Conf. 199, 01003 (2019).

Article CAS Google Scholar

Martinez Torres, A., Khemchandani, K. P., Roca, L. & Oset, E. Few-body systems consisting of mesons. Few Body Syst. 61, 35 (2020).

Article ADS Google Scholar

Tornqvist, N. A. Possible large deuteronlike meson-meson states bound by pions. Phys. Rev. Lett. 67, 556 (1991).

Article ADS CAS PubMed Google Scholar

LHCb Collaboration, Aaij, R. et al. Observation of the doubly charmed baryon \({{{\Xi }}}_{{{{{{{{\rm{cc}}}}}}}}}^{++}\). Phys. Rev. Lett. 119, 112001. https://doi.org/10.1103/PhysRevLett.119.112001 (2017).

LHCb Collaboration, Aaij, R. et al. Measurement of the lifetime of the doubly charmed baryon \({{{\Xi }}}_{{{{{{{{\rm{cc}}}}}}}}}^{++}\). Phys. Rev. Lett. 121, 052002. https://doi.org/10.1103/PhysRevLett.121.052002 (2018).

LHCb Collaboration, Aaij, R. et al. Observation of structure in the J/ψ-pair mass spectrum. Sci. Bull. 65, 1983 (2020).

Albuquerque, R.M. et al. Doubly-hidden scalar heavy molecules and tetraquarks states from QCD at NLO. Phys. Rev. D102, 094001 (2020).

ADS Google Scholar

Dong, X.-K. et al. Coupled-channel interpretation of the LHCb double-J/ψ spectrum and hints of a new state near the J/ψJ/ψ threshold. Phys. Rev. Lett. 126, 132001 (2021).

Article ADS CAS PubMed Google Scholar

Bedolla, M.A., Ferretti, J., Roberts, C.D. & Santopinto, E. Spectrum of fully-heavy tetraquarks from a diquark+antidiquark perspective. Eur. Phys. J. C80, 1004 (2020).

Article ADS CAS Google Scholar

Karliner, M. & Rosner, J.L. Interpretation of structure in the di-J/ψ spectrum. Phys. Rev. D102, 114039 (2020).

ADS Google Scholar

Lü, Q.-F., Chen, D.Y. & Dong, Y.B. Masses of fully heavy tetraquarks \({{{{{{{\rm{Q}}}}}}}}{{{{{{{\rm{Q}}}}}}}}\bar{{{{{{{{\rm{Q}}}}}}}}}\bar{{{{{{{{\rm{Q}}}}}}}}}\) in an extended relativized quark model. Eur. Phys. J. C80, 871 (2020).

Article ADS CAS Google Scholar

Giron, J.F. & Lebed, R.F. Simple spectrum of \({{{{{{{\rm{c}}}}}}}}\overline{{{{{{{{\rm{c}}}}}}}}}{{{{{{{\rm{c}}}}}}}}\overline{{{{{{{{\rm{c}}}}}}}}}\) states in the dynamical diquark model. Phys. Rev. D102, 074003 (2020).

ADS Google Scholar

Ader, J.P., Richard, J.M. & Taxil, P. Do narrow heavy multiquark states exist? Phys. Rev. D25, 2370 (1982).

ADS Google Scholar

Ballot, J.I. & Richard, J.M. Four quark states in additive potentials. Phys. Lett. B123, 449 (1983).

Article ADS Google Scholar

Zouzou, S., Silvestre-Brac, B., Gignoux, C. & Richard, J.M. Four quark bound states. Z. Phys. C30, 457 (1986).

ADS Google Scholar

Lipkin, H.J. A model-independent approach to multiquark bound states. Phys. Lett. B172, 242 (1986).

Article ADS Google Scholar

Heller, L. & Tjon, J.A. On the existence of stable dimesons. Phys. Rev. D35, 969 (1987).

ADS Google Scholar

Manohar, A.V. & Wise, M.B. Exotic \({{{{{{{\rm{Q}}}}}}}}{{{{{{{\rm{Q}}}}}}}}\bar{{{{{{{{\rm{q}}}}}}}}}\bar{{{{{{{{\rm{q}}}}}}}}}\) states in QCD. Nucl. Phys. B399, 17 (1993).

Article ADS Google Scholar

Carlson, J., Heller, L. & Tjon, J.A. Stability of dimesons. Phys. Rev. D37, 744 (1988).

ADS Google Scholar

Silvestre-Brac, B. & Semay, C. Systematics of L = 0 \({{{{{{{{\rm{q}}}}}}}}}^{2}{\overline{{{{{{{{\rm{q}}}}}}}}}}^{2}\) systems. Z. Phys. C57, 273 (1993).

ADS Google Scholar

Semay, C. & Silvestre-Brac, B. Diquonia and potential models. Z. Phys. C61, 271 (1994).

ADS Google Scholar

Moinester, M.A. How to search for doubly charmed baryons and tetraquarks. Z. Phys. A355, 349 (1996).

Article ADS Google Scholar

Pepin, S., Stancu, F., Genovese, M. & Richard, J.M. Tetraquarks with color blind forces in chiral quark models. Phys. Lett. B393, 119 (1997).

Article ADS Google Scholar

Gelman, B.A. & Nussinov, S. Does a narrow tetraquark \({{{{{{{\rm{c}}}}}}}}{{{{{{{\rm{c}}}}}}}}\overline{{{{{{{{\rm{u}}}}}}}}}\overline{{{{{{{{\rm{d}}}}}}}}}\) state exist?Phys. Lett. B551, 296 (2003).

Article ADS Google Scholar

Vijande, J., Fernandez, F., Valcarce, A. & Silvestre-Brac, B. Tetraquarks in a chiral constituent quark model. Eur. Phys. J. A19, 383 (2004).

Article ADS CAS Google Scholar

Janc, D. & Rosina, M. The Tcc = DD* molecular state. Few Body Syst. 35, 175 (2004).

ADS CAS Google Scholar

Navarra, F. S., Nielsen, M. & Lee, S. H. QCD sum rules study of \({{{{{{{\rm{Q}}}}}}}}{{{{{{{\rm{Q}}}}}}}}-\overline{{{{{{{{\rm{u}}}}}}}}}\overline{{{{{{{{\rm{d}}}}}}}}}\) mesons. Phys. Lett. B649, 166 (2007).

Article ADS CAS Google Scholar

Vijande, J., Weissman, E., Valcarce, A. & Barnea, N. Are there compact heavy four-quark bound states? Phys. Rev. D76, 094027 (2007).

ADS Google Scholar

Ebert, D., Faustov, R. N., Galkin, V. O. & Lucha, W. Masses of tetraquarks with two heavy quarks in the relativistic quark model. Phys. Rev. D76, 114015 (2007).

ADS Google Scholar

Lee, S. H. & Yasui, S. Stable multiquark states with heavy quarks in a diquark model. Eur. Phys. J. C64, 283 (2009).

Article ADS CAS Google Scholar

Yang, Y., Deng, C., Ping, J. & Goldman, T. S-wave \({{{{{{{\rm{QQ}}}}}}}}\overline{{{{{{{{\rm{q}}}}}}}}}\overline{{{{{{{{\rm{q}}}}}}}}}\) state in the constituent quark model. Phys. Rev. D80, 114023 (2009).

ADS Google Scholar

Li, N., Sun, Z.-F., Liu, X. & Zhu, S.-L. Coupled-channel analysis of the possible D(*)D(*), \({\overline{{{{{{{{\rm{B}}}}}}}}}}^{(* )}{\overline{{{{{{{{\rm{B}}}}}}}}}}^{(* )}\) and \({{{{{{{{\rm{D}}}}}}}}}^{(* )}{\overline{{{{{{{{\rm{B}}}}}}}}}}^{(* )}\) molecular states. Phys. Rev. D88, 114008 (2013).

ADS Google Scholar

Feng, G.-Q., Guo, X.-H., & Zou, B.-S. \({{{{{{{{\rm{QQ}}}}}}}}}^{\prime}\overline{{{{{{{{\rm{u}}}}}}}}}\overline{{{{{{{{\rm{d}}}}}}}}}\) bound state in the Bethe-Salpeter equation approach. Preprint at arXiv:1309.7813

Luo, S.-Q. et al. Exotic tetraquark states with the \({{{{{{{\rm{q}}}}}}}}{{{{{{{\rm{q}}}}}}}}\bar{{{{{{{{\rm{Q}}}}}}}}}\bar{{{{{{{{\rm{Q}}}}}}}}}\) configuration. Eur. Phys. J. C77, 709 (2017).

Article ADS CAS Google Scholar

Karliner, M. & Rosner, J.L. Discovery of doubly-charmed \({{{{\Xi }}}_{{{{{{{{\rm{cc}}}}}}}}}}_{{{{{{{{\rm{c}}}}}}}}{{{{{{{\rm{c}}}}}}}}}\) baryon implies a stable \({{{{\Xi }}}_{{{{{{{{\rm{cc}}}}}}}}}}_{{{{{{{{\rm{c}}}}}}}}{{{{{{{\rm{c}}}}}}}}}\) tetraquark. Phys. Rev. Lett. 119, 202001 (2017).

Article PubMed Google Scholar

Eichten, E. J. & Quigg, C. Heavy-quark symmetry implies stable heavy tetraquark mesons \({{{{{{{{\rm{Q}}}}}}}}}_{i}{{{{{{{{\rm{Q}}}}}}}}}_{j}{\bar{{{{{{{{\rm{q}}}}}}}}}}_{k}{\bar{{{{{{{{\rm{q}}}}}}}}}}_{l}\). Phys. Rev. Lett. 119, 202002 (2017).

Article PubMed Google Scholar

Wang, Z.-G. Analysis of the axialvector doubly heavy tetraquark states with QCD sum rules. Acta Phys. Polon. B49, 1781 (2018).

Article ADS Google Scholar

Park, W., Noh, S. & Lee, S. H. Masses of the doubly heavy tetraquarks in a constituent quark model. Acta Phys. Polon. B50, 1151 (2019).

Article ADS MathSciNet Google Scholar

Junnarkar, P., Mathur, N. & Padmanath, M. Study of doubly heavy tetraquarks in lattice QCD. Phys. Rev. D99, 034507 (2019).

ADS Google Scholar

Deng, C., Chen, H. & Ping, J. Systematical investigation on the stability of doubly heavy tetraquark states. Eur. Phys. J. A56, 9 (2020).

Article ADS CAS Google Scholar

Liu, M.-Z. et al. Heavy-quark spin and flavor symmetry partners of the X(3872) revisited: What can we learn from the one boson exchange model? Phys. Rev. D99, 094018 (2019).

ADS Google Scholar

Maiani, L., Polosa, A. D. & Riquer, V. Hydrogen bond of QCD in doubly heavy baryons and tetraquarks. Phys. Rev. D100, 074002 (2019).

ADS Google Scholar

Yang, G., Ping, J. & Segovia, J. Doubly-heavy tetraquarks. Phys. Rev. D101, 014001 (2020).

ADS Google Scholar

Tan, Y., Lu, W. & Ping, J. \({{{{{{{\rm{Q}}}}}}}}{{{{{{{\rm{Q}}}}}}}}\bar{{{{{{{{\rm{q}}}}}}}}}\bar{{{{{{{{\rm{q}}}}}}}}}\) in a chiral constituent quark model. Eur. Phys. J. Plus 135, 716 (2020).

Article Google Scholar

Lü, Q.-F., Chen, D.-Y. & Dong, Y.-B. Masses of doubly heavy tetraquarks \({{{{{{{{\rm{T}}}}}}}}}_{{{{{{{{\rm{Q}}}}}}}}{{{{{{{{\rm{Q}}}}}}}}}^{\prime}}\) in a relativized quark model. Phys. Rev. D102, 034012 (2020).

ADS Google Scholar

Braaten, E., He, L.-P. & Mohapatra, A. Masses of doubly heavy tetraquarks with error bars. Phys. Rev. D103, 016001 (2021).

ADS Google Scholar

Gao, D. et al. Masses of doubly heavy tetraquark states with isospin =\(\frac{1}{2}\) and 1 and spin-parity 1+±. Preprint at arXiv:2007.15213

Cheng, J.-B. et al. Double-heavy tetraquark states with heavy diquark-antiquark symmetry. Chin. Phys. C45, 043102 (2021).

Article ADS CAS Google Scholar

Noh, S., Park, W. & Lee, S. H. The doubly-heavy tetraquarks, \({{{{{{{\rm{q}}}}}}}}{{{{{{{{\rm{q}}}}}}}}}^{\prime}\bar{{{{{{{{\rm{Q}}}}}}}}}{\bar{{{{{{{{\rm{Q}}}}}}}}}}^{\prime}\), in a constituent quark model with a complete set of harmonic oscillator bases. Phys. Rev. D103, 114009 (2021).

ADS Google Scholar

Faustov, R. N., Galkin, V. O. & Savchenko, E. M. Heavy tetraquarks in the relativistic quark model. Universe 7, 94 (2021).

Article ADS CAS Google Scholar

LHCb Collaboration, Aaij, R. et al. Observation of an exotic narrow doubly charmed tetraquark. https://doi.org/10.1038/s41567-022-01614-y (2022).

LHCb Collaboration, Alvers, A. A. Jr. et al. The LHCb detector at the LHC. JINST 3 S08005 (2008).

LHCb Collaboration, Aaij, R. et al. LHCb detector performance. Int. J. Mod. Phys. 30, A301530022 (2015).

LHCb Collaboration, Aaij, R. et al. Observation of double charm production involving open charm in pp collisions at \(\sqrt{s}=\)7TeV. JHEP 06, 141 (2012).

LHCb Collaboration, Aaij, R. et al. Observation of associated production of a Z boson with a D meson in the forward region. JHEP 04, 091 (2014).

LHCb Collaboration, AAij, R. et al. Production of associated ϒres and open charm hadrons in pp collisions at \(\sqrt{s}=\)7 and 8TeV via double parton scattering. JHEP 07, 052 (2016).

LHCb Collaboration, Aaij, R. et al. Near-threshold D\(\overline{{{{{{{{\rm{D}}}}}}}}}\) spectroscopy and observation of a new charmonium state. JHEP 07, 035 (2019).

Blatt, J. M. & Weisskopf, V. F. Theoretical Nuclear Physics. (Springer, New York, 1952).

von Hippel, F. & Quigg, C. Centrifugal-barrier effects in resonance partial decay widths, shapes, and production amplitudes. Phys. Rev. D5, 624 (1972).

ADS Google Scholar

Hanhart, C., Kalashnikova, Y. S., Kudryavtsev, A. E. & Nefediev, A. V. Reconciling the X(3872) with the near-threshold enhancement in the \({{{{{{{{\rm{D}}}}}}}}}^{0}{\overline{{{{{{{{\rm{D}}}}}}}}}}^{* 0}\) final state. Phys. Rev. D76, 034007 (2007).

ADS Google Scholar

Braaten, E. & Lu, M. Line shapes of the X(3872). Phys. Rev. D76, 094028 (2007).

ADS Google Scholar

Braaten, E. & Stapleton, J. Analysis of J/ψ π+π− and \({{{{{{{{\rm{D}}}}}}}}}^{0}{\overline{{{{{{{{\rm{D}}}}}}}}}}^{0}{\pi }^{0}\) decays of the X(3872). Phys. Rev. D81, 014019 (2010).

ADS Google Scholar

Kalashnikova, Y. S. & Nefediev, A. V. Nature of X(3872) from data. Phys. Rev. D80, 074004 (2009).

ADS Google Scholar

Artoisenet, P., Braaten, E. & Kang, D. Using line shapes to discriminate between binding mechanisms for the X(3872). Phys. Rev. D82, 014013 (2010).

ADS Google Scholar

Hanhart, C., Kalashnikova, Y. S. & Nefediev, A. V. Lineshapes for composite particles with unstable constituents. Phys. Rev. D81, 094028 (2010).

ADS Google Scholar

Hanhart, C., Kalashnikova, Y. S. & Nefediev, A. V. Interplay of quark and meson degrees of freedom in a near-threshold resonance: multi-channel case. Eur. Phys. J. A47, 101 (2011).

Article ADS CAS Google Scholar

Hanhart, C. et al. Practical parametrization for line shapes of near-threshold states. Phys. Rev. Lett. 115, 202001 (2015).

Article ADS CAS PubMed Google Scholar

Guo, F.-K. et al. Interplay of quark and meson degrees of freedom in near-threshold states: a practical parametrization for line shapes. Phys. Rev. D93, 074031 (2016).

ADS Google Scholar

Hanhart, C. et al. A practical parametrisation of line shapes of near-threshold resonances. J. Phys. Conf. Ser. 675, 022016 (2016).

Article CAS Google Scholar

Guo, F. K. et al. Phenomenology of near-threshold states: a practical parametrisation for the line shapes. EPJ Web Conf. 137, 06020 (2017).

Article CAS Google Scholar

Mikhasenko, M. et al. Three-body scattering: ladders and resonances. JHEP 08, 080 (2019).

Article ADS MathSciNet Google Scholar

Pasquier, R. & Pasquier, J. Y. Khuri-Treiman-type equations for three-body decay and production processes. Phys. Rev. 170, 1294 (1968).

Article ADS Google Scholar

Aitchison, I. J. R. & Brehm, J. J. Are there important unitarity corrections to the isobar model?. Phys. Lett. B 84, 349 (1979).

Article ADS Google Scholar

LHCb Collaboration, Aaij, R. et al. Study of the line shape of the χc1(3872) state. Phys. Rev. D102, 092005 (2020).

LHCb Collaboration, Aaij, R. et al. Study of the ψ2(3823) and χc1(3872) states in \({{{{{{{{\rm{B}}}}}}}}}^{+}\,\to \left({{{{{{{\rm{J}}}}}}}}/\psi {\pi }^{+}{\pi }^{-}\right){{{{{{{{\rm{K}}}}}}}}}^{+}\) decays. JHEP 08, 123 (2020).

LHCb Collaboration, Aaij, R. et al. Study of \({{{{{{{{\rm{B}}}}}}}}}_{{{{{{{{\rm{s}}}}}}}}}^{0}\,\to {{{{{{{\rm{J}}}}}}}}/\psi {\pi }^{+}{\pi }^{-}{{{{{{{{\rm{K}}}}}}}}}^{+}{{{{{{{{\rm{K}}}}}}}}}^{-}\) decays. JHEP 02, 024 (2021).

Källén, G. Elementary Particle Physics (Addison-Wesley, Reading, Massachusetts, 1964).

Flatté, S. M. Coupled-channel analysis of the πη and \({{{{{{{\rm{K}}}}}}}}\overline{{{{{{{{\rm{K}}}}}}}}}\) systems near \({{{{{{{\rm{K}}}}}}}}\overline{{{{{{{{\rm{K}}}}}}}}}\) threshold. Phys. Lett. B63, 224 (1976).

Baru, V. et al. Flatté-like distributions and the a0(980)/f0(980) mesons. Eur. Phys. J. A23, 523 (2005).

Article ADS CAS Google Scholar

Martínez Santos, D. & Dupertuis, F. Mass distributions marginalized over per-event errors. Nucl. Instrum. Meth. A764, 150 (2014).

Article ADS CAS Google Scholar

Skwarnicki, T. A study of the radiative cascade transitions between the \({\Upsilon }^{\prime}\) and ϒ resonances. PhD thesis. Institute of Nuclear Physics, Krakow. http://inspirehep.net/record/230779/ (1986). DESY-F31-86-02.

LHCb Collaboration, Aaij, R. et al. Observation of J/ψ-pair production in pp collisions at \(\sqrt{s}=\)7TeV. Phys. Lett. B707, 52 (2012).

Gosset, W. S. (Student). The probable error of a mean. Biometrika 6, 1 (1908).

Jackman, S. Bayesian Analysis for the Social Sciences. (John Wiley & Sons, Inc., Hoboken, New Jersey, USA, 2009).

Johnson, N.L. Systems of frequency curves generated by methods of translation. Biometrika 36, 149 (1949).

Article MathSciNet CAS PubMed MATH Google Scholar

Johnson, N.L. Bivariate distributions based on simplie translation systems. Biometrika 36, 297 (1949).

Article MathSciNet CAS PubMed Google Scholar

BaBar Collaboration, Lees, J. P. et al. Branching fraction measurements of the color-suppressed decays \({\overline{{{{{{{{\rm{B}}}}}}}}}}^{0}\) to D(*)0π0, D(*)0η, D(*)0ω, and \({{{{{{{{\rm{D}}}}}}}}}^{(* )0}{\eta }^{\prime}\) and measurement of the polarization in the decay \({\overline{{{{{{{{\rm{B}}}}}}}}}}^{0}\,\to {{{{{{{{\rm{D}}}}}}}}}^{* 0}\omega\). Phys. Rev. D84, 112007 (2011). Erratum https://doi.org/10.1103/PhysRevD.87.039901D87 (2013) 039901(E).

LHCb Collaboration, Aaij, R. et al. χc1 and χc2 resonance parameters with the decays χc1,c2 → J/ψ μ+μ−. Phys. Rev. Lett. 119, 221801 (2017).

LHCb Collaboration, Aaij, R. et al. Observation of a new baryon state in the \({\Lambda }_{{{{{{{{\rm{b}}}}}}}}}^{0}{\pi }^{+}{\pi }^{-}\) mass spectrum. JHEP 06, 136 (2020).

LHCb Collaboration, AAij, R. et al. Updated search for \({{{{{{{{\rm{B}}}}}}}}}_{{{{{{{{\rm{c}}}}}}}}}^{+}\) decays to two open charm mesons. JHEP 12, 117 (2021).

Byckling, E. & Kajantie, K. Particle Kinematics (John Wiley & Sons Inc., New York, 1973).

BaBar Collaboration, Lees, J. P. et al. Measurement of the D*(2010)+ meson width and the D*(2010)+ − D0 mass difference. Phys. Rev. Lett. 111, 111801 (2013).

Guo, F.-K. Novel method for precisely measuring the X(3872) mass. Phys. Rev. Lett. 122, 202002 (2019).

Article ADS CAS PubMed Google Scholar

LHCb Collaboration, Aaij, R. et al. Measurements of the \({\Lambda }_{{{{{{{{\rm{b}}}}}}}}}^{0}\), \({{{\Xi }}}_{{{{{{{{\rm{b}}}}}}}}}^{-}\), and Ωb baryon masses. Phys. Rev. Lett. 110, 182001 (2013).

LHCb Collaboration, Aaij, R. et al. Precision measurement of D meson mass differences. JHEP 06, 065 (2013).

Bethe, H. A. Zur Theorie des Durchgangs schneller Korpuskularstrahlen durch Materie. Annalen der Physik 397, 325 (1930).

Article ADS MATH Google Scholar

Bloch, F. Zur Bremsung rasch bewegter Teilchen beim Durchgang durch Materie. Annalen der Physik 408, 285 (1933).

Article ADS MATH Google Scholar

LHCb Collaboration, Aaij, R. et al. Prompt \({{{{{{{{\rm{K}}}}}}}}}_{{{{{{{{\rm{S}}}}}}}}}^{0}\) production in pp collisions at \(\sqrt{s}=\)0.9TeV. Phys. Lett. B693, 69 (2010).

CLEO Collaboration, Anastassov, A. et al. First measurement of Γ(D*+) and precision measurement of \({m}_{{{{{{{{{\rm{D}}}}}}}}}^{* +}}-{m}_{{{{{{{{{\rm{D}}}}}}}}}^{0}}\). Phys. Rev. D65, 032003 (2002).

Wilks, S. S. The large-sample distribution of the likelihood ratio for testing composite hypotheses. Ann. Math. Stat. 9, 60 (1938).

Article MATH Google Scholar

Del Fabbro, A., Janc, D., Rosina, M. & Treleani, D. Production and detection of doubly charmed tetraquarks. Phys. Rev. D71, 014008 (2005).

ADS Google Scholar

Bethe, H. A. Theory of the effective range in nuclear scattering. Phys. Rev. 76, 38 (1949).

Article ADS CAS MATH Google Scholar

Balakrishnan, N., Kharchenko, V., Forrey, R. C. & Dalgarno, A. Complex scattering lengths in multi-channel atom-molecule collisions. Chemical Physics Letters 280, 5 (1997).

Article ADS CAS Google Scholar

Weinberg, S. Evidence that the deuteron is not an elementary particle. Phys. Rev. 137, B672. https://doi.org/10.1103/PhysRev.137.B672 (1965).

Matuschek, I., Baru, V., Guo, F.-K. & Hanhart, C. On the nature of near-threshold bound and virtual states. Eur. Phys. J. A57, 101 (2021).

Article ADS CAS Google Scholar

LHCb Collaboration, Aaij, R. et al. Modification of χc1(3872) and ψ(2S) production in pp collisions at \(\sqrt{s}=8{{{{{{{\rm{TeV}}}}}}}}\). Phys. Rev. Lett. 126, 092001 (2021).

Belyaev, I. & Savrina, D. Study of double parton scattering processes with heavy quarks, in Multiple Parton Interactions at the LHC (eds Bartalini, P. and Gaunt, J. R.) (World Scientific, Singapore, 2018).

Adinolfi, M. et al. Performance of the LHCb RICH detector at the LHC. Eur. Phys. J. C73, 2431 (2013).

Article ADS CAS Google Scholar

Sjöstrand, T., Mrenna, S. & Skands, P. A brief introduction to PYTHIA 8.1. Comput. Phys. Commun. 178, 852 (2008).

Article ADS MATH CAS Google Scholar

Belyaev, I. et al. Handling of the generation of primary events inGauss, the LHCb simulation framework. J. Phys. Conf. Ser. 331, 032047 (2011).

Article Google Scholar

Lange, D.J. The EvtGen particle decay simulation package. Nucl. Instrum. Meth. A462, 152 (2001).

Article ADS Google Scholar

Davidson, N., Przedzinski, T. & Was, Z.PHOTOS interface in C++: Technical and physics documentation. Comp. Phys. Comm. 199, 86 (2016).

Article ADS MathSciNet CAS Google Scholar

Geant4 Collaboration, Allison, J. et al. Geant4 developments and applications. IEEE Trans. Nucl. Sci. 53, 270 (2006).

Geant4 Collaboration, Allison, J. et al. Geant4: a simulation toolkit. Nucl. Instrum. Meth. A506, 250 (2003).

Clemencic, M. et al. The LHCb simulation application,Gauss: design, evolution and experience. J. Phys. Conf. Ser. 331, 032023 (2011).

Article Google Scholar

Hulsbergen, W.D. Decay chain fitting with a Kalman filter. Nucl. Instrum. Meth. A552, 566 (2005).

Article ADS CAS Google Scholar

Pivk, M. & Le Diberder, F. R. sPlot: a statistical tool to unfold data distributions. Nucl. Instrum. Meth. A555, 356 (2005).

Article ADS CAS Google Scholar

Rosner, J.L. Quark models. NATO Sci. Ser. B 66, 1 (1981).

Google Scholar

Gasiorowicz, S. & Rosner, J. L. Hadron spectra and quarks. Am. J. Phys. 49, 954 (1981).

Article ADS CAS Google Scholar

Rosner, J. L. Hadronic and radiative D* widths. Phys. Rev. D88, 034034 (2013).

ADS Google Scholar

Peskin, M. E. & Schroeder, D. V. An Introduction to Quantum Field Theory (Addison-Wesley, Reading, USA, 1995).

Gribov, V. N. Strong Interactions of Hadrons at High Energies: Gribov Lectures on Theoretical Physics (Cambridge University Press, 2012).

Martin, A. D. & Spearman, T. D. Elementary-Particle Theory (North-Holland, Amsterdam, 1970).

Eden, R. J., Landshoff, P. V., Olive, D. I. & Polkinghorne, J. C. The Analytic S-matrix (Cambridge Univ. Press, Cambridge, 1966).

Aitchison, I. J. R. K-matrix formalism for overlapping resonances. Nucl. Phys. A189, 417 (1972).

Article ADS Google Scholar

JPAC Collaboration, Mikhasenko, M. et al. Pole position of the a1(1260) from τ-decay. Phys. Rev. D98, 096021 (2018).

Karliner, M. & Rosner, J. L. Status of isospin splittings in mesons and baryons. Phys. Rev. D100, 073006 (2019).

ADS Google Scholar

LHCb Collaboration, Aaij, R. et al. Model-independent study of structure in B+ → D+D−K+ decays. Phys. Rev. Lett. 125, 242001 (2020).

LHCb Collaboration, Aaij, R. et al. Amplitude analysis of the B+ → D+D−K+ decay. Phys. Rev. D102, 112003 (2020).

BESIII Collaboration, Ablikim, M. et al. Observation of a charged charmoniumlike structure in e+e− → π+π−J/ψ at \(\sqrt{s}=4.26{{{{{{{\rm{GeV}}}}}}}}\). Phys. Rev. Lett. 110, 252001 (2013).

Belle Collaboration, Liu, Z. Q. et al. Study of e+e− → π+π−J/ψ and observation of a charged charmoniumlike state at Belle. Phys. Rev. Lett. 110, 252002 (2013). Erratum https://doi.org/10.1103/PhysRevLett.111.019901 111, 019901 (2013).

BESIII Collaboration, Ablikim, M. et al. Observation of a charged \({({{{{{{{\rm{D}}}}}}}}{\overline{{{{{{{{\rm{D}}}}}}}}}}^{* })}^{\pm }\) mass peak in \({{{{{{{{\rm{e}}}}}}}}}^{+}{{{{{{{{\rm{e}}}}}}}}}^{-}\,\to \pi {{{{{{{\rm{D}}}}}}}}{\overline{{{{{{{{\rm{D}}}}}}}}}}^{* }\) at \(\sqrt{s}=4.26 {{{{{{{\rm{GeV}}}}}}}}\). Phys. Rev. Lett. 112, 022001 (2014).

BESIII Collaboration, Ablikim, M. et al. Observation of Zc(3900)0 in e+e− → π0π0J/ψ. Phys. Rev. Lett. 115, 112003 (2015).

BESIII Collaboration, Ablikim, M. et al. Observation of a neutral structure near the \({{{{{{{\rm{D}}}}}}}}{\overline{{{{{{{{\rm{D}}}}}}}}}}^{* }\) mass threshold in \({{{{{{{{\rm{e}}}}}}}}}^{+}{{{{{{{{\rm{e}}}}}}}}}^{-}\,\to {({{{{{{{\rm{D}}}}}}}}{\overline{{{{{{{{\rm{D}}}}}}}}}}^{* })}^{0}{\pi }^{0}\) at \(\sqrt{s}=4.226\) and 4.257 GeV. Phys. Rev. Lett. 115, 222002 (2015).

BESIII Collaboration, Ablikim, M. et al. Observation of a charged charmoniumlike structure Zc(4020) and search for the Zc(3900) in e+e− → π+π−hc. Phys. Rev. Lett. 111, 242001 (2013).

BESIII Collaboration, Ablikim, M. et al. Observation of a charged charmoniumlike structure in\({{{{{{{{\rm{e}}}}}}}}}^{+}{{{{{{{{\rm{e}}}}}}}}}^{-}\,\to {\left({{{{{{{{\rm{D}}}}}}}}}^{* }{\overline{{{{{{{{\rm{D}}}}}}}}}}^{* }\right)}^{\pm }{\pi }^{\mp }\) at \(\sqrt{s}=4.26{{{{{{{\rm{GeV}}}}}}}}\). Phys. Rev. Lett. 112, 132001 (2014).

Belle Collaboration, Mizuk, R. et al. Observation of two resonance-like structures in the π+χc1 mass distribution in exclusive \({\overline{{{{{{{{\rm{B}}}}}}}}}}^{0}\,\to {{{{{{{{\rm{K}}}}}}}}}^{-}{\pi }^{+}{\chi }_{{{{{{{{\rm{c}}}}}}}}1}\) decays. Phys. Rev. D78, 072004 (2008).

LHCb Collaboration, Aaij, R. et al. Evidence for a ηc(1S)π− resonance in B0 → ηc(1S)K+π− decays. Eur. Phys. J. C78, 1019 (2018).

Belle Collaboration, Chilikin, K. et al. Observation of a new charged charmoniumlike state in \({\overline{{{{{{{{\rm{B}}}}}}}}}}^{0}\,\to {{{{{{{\rm{J}}}}}}}}/\psi {{{{{{{{\rm{K}}}}}}}}}^{-}{\pi }^{+}\) decays. Phys. Rev. D90, 112009 (2014).

Belle Collaboration, Choi, S. K. et al. Observation of a resonancelike structure in the \({\pi }^{\pm }{\psi }^{\prime}\) mass distribution in exclusive \({{{{{{{\rm{B}}}}}}}}\,\to {{{{{{{\rm{K}}}}}}}}{\pi }^{\pm }{\psi }^{\prime}\) decays. Phys. Rev. Lett. 100, 142001 (2008).

Belle Collaboration, Chilikin, K. et al. Experimental constraints on the spin and parity of the Z(4430)+. Phys. Rev. D88, 074026 (2013).

LHCb Collaboration, Aaij, R. et al. Observation of the resonant character of the Z(4430)− state. Phys. Rev. Lett. 112, 222002 (2014).

LHCb Collaboration, Aaij, R. et al. Model-independent confirmation of the Z(4430)− state. Phys. Rev. D92, 112009 (2015).

BESIII Collaboration, Ablikim, M. et al. Observation of a near-threshold structure in the K+ recoil-mass spectra in \({{{{{{{{\rm{e}}}}}}}}}^{+}{{{{{{{{\rm{e}}}}}}}}}^{-}\,\to {{{{{{{{\rm{K}}}}}}}}}^{+}\left({{{{{{{{\rm{D}}}}}}}}}_{{{{{{{{\rm{s}}}}}}}}}^{-}{{{{{{{{\rm{D}}}}}}}}}^{* 0}+{{{{{{{{\rm{D}}}}}}}}}_{{{{{{{{\rm{s}}}}}}}}}^{* -}{{{{{{{{\rm{D}}}}}}}}}^{0}\right)\). Phys. Rev. Lett. 126, 102001 (2021).

LHCb Collaboration, Aaij, R. et al. Observation of new resonances decaying to J/ψK+ and J/ψ ϕ. Phys. Rev. Lett. 127, 082001 (2020).

CDF Collaboration, Aaltonen, T. et al. Evidence for a narrow near-threshold structure in the J/ψ ϕ mass spectrum in B+ → J/ψ ϕK+ decays. Phys. Rev. Lett. 102, 242002 (2009).

D0 Collaboration, Abazon, V. M. et al. Search for the X(4140) state in B+ → J/ψ ϕK+ decays with the D0 Detector. Phys. Rev. D89, 012004 (2014).

CMS Collaboration, Chartchyan, S. et al. Observation of a peaking structure in the J/ψ ϕ mass spectrum from B+ → J/ψ ϕK+ decays. Phys. Lett. B734, 261 (2014).

LHCb Collaboration, Aaij, R. et al. Observation of exotic J/ψ ϕ structures from amplitude analysis of B+ → J/ψ ϕK+ decays. Phys. Rev. Lett. 118, 022003 (2017).

Belle Collaboration, Bondar, A. et al. Observation of two charged bottomonium-like resonances in ϒ(5S) decays. Phys. Rev. Lett. 108, 122001 (2012).

LHCb Collaboration, AAij, R. et al. Observation of a narrow pentaquark state, Pc(4312)+, and of two-peak structure of the Pc(4450)+. Phys. Rev. Lett. 122, 222001 (2019).

LHCb Collaboration, AAij, R. et al. Observation of J/ψp resonances consistent with pentaquark states in\({\Lambda }_{{{{{{{{\rm{b}}}}}}}}}^{0}\,\to {{{{{{{\rm{J}}}}}}}}/\psi {{{{{{{\rm{p}}}}}}}}{{{{{{{{\rm{K}}}}}}}}}^{-}\) decays. Phys. Rev. Lett. 115, 072001 (2015).

LHCb Collaboration, AAij, R. et al. Evidence for a new structure in the J/ψp and \({{{{{{{\rm{J}}}}}}}}/\psi \overline{{{{{{{{\rm{p}}}}}}}}}\) systems in \({{{{{{{{\rm{B}}}}}}}}}_{{{{{{{{\rm{s}}}}}}}}}^{0}\,\to {{{{{{{\rm{J}}}}}}}}/\psi {{{{{{{\rm{p}}}}}}}}\overline{{{{{{{{\rm{p}}}}}}}}}\) decays. Phys. Rev. Lett. 128, 062001 (2021).

LHCb Collaboration, Aaij, R. et al. Evidence of a J/ψΛ structure and observation of excited Ξ− states in the \({{{\Xi }}}_{{{{{{{{\rm{b}}}}}}}}}^{-}bm\,\to {{{{{{{\rm{J}}}}}}}}/\psi \Lambda {{{{{{{{\rm{K}}}}}}}}}^{-}\) decay. Sci. Bull. 66, 1278 (2021).

Download references

This paper is dedicated to the memory of our dear friend and colleague Simon Eidelman, whose contributions to improving the quality of our papers were greatly appreciated. We express our gratitude to our colleagues in the CERN accelerator departments for the excellent performance of the LHC. We thank the technical and administrative staff at the LHCb institutes. We acknowledge support from CERN and from the national agencies: CAPES, CNPq, FAPERJ and FINEP (Brazil); MOST and NSFC (China); CNRS/IN2P3 (France); BMBF, DFG and MPG (Germany); INFN (Italy); NWO (Netherlands); MNiSW and NCN (Poland); MEN/IFA (Romania); MSHE (Russia); MICINN (Spain); SNSF and SER (Switzerland); NASU (Ukraine); STFC (UK); DOE NP and NSF (USA). We acknowledge the computing resources that are provided by CERN, IN2P3 (France), KIT and DESY (Germany), INFN (Italy), SURF (Netherlands), PIC (Spain), GridPP (UK), RRCKI and Yandex LLC (Russia), CSCS (Switzerland), IFIN-HH (Romania), CBPF (Brazil), PL-GRID (Poland) and NERSC (USA). We are indebted to the communities behind the multiple open-source software packages on which we depend. Individual groups or members have received support from ARC and ARDC (Australia); AvH Foundation (Germany); EPLANET, Marie Skłodowska-Curie Actions and ERC (European Union); A*MIDEX, ANR, IPhU and Labex P2IO, and Région Auvergne-Rhône-Alpes (France); Key Research Program of Frontier Sciences of CAS, CAS PIFI, CAS CCEPP, Fundamental Research Funds for the Central Universities, and Sci. & Tech. Program of Guangzhou (China); RFBR, RSF and Yandex LLC (Russia); GVA, XuntaGal and GENCAT (Spain); the Leverhulme Trust, the Royal Society and UKRI (UK).

Nikhef National Institute for Subatomic Physics, Amsterdam, Netherlands

R. Aaij, J. S. Butter, K. Carvalho Akiba, S. Ferreres Sole, E. Gabriel, R. E. Geertsema, L. M. Greeven, K. Heijhoff, W. Hulsbergen, D. Hynds, E. Jans, S. Klaver, P. Koppenburg, I. Kostiuk, H. S. Kuindersma, M. Lucio Martinez, V. Lukashenko, A. Mauri, M. Merk, A. Pellegrino, C. Sanchez Gras, M. Schubiger, A. Snoch, N. Tuning, A. Usachov, M. van Beuzekom & M. Veronesi

Department of Physics, University of Warwick, Coventry, UK

A. S. W. Abdelmotteleb, F. J. Abudinen Gallego, J. J. Back, A. Beck, T. Blake, A. Brossa Gonzalo, M. F. Cicala, C. M. Costa Sobral, L. M. Garcia Martin, T. Gershon, R. J. Hunter, T. P. Jones, M. Kenzie, M. Kreps, T. Latham, O. Lupton, E. Millard, A. G. Morris, B. Pagare, M. Ramos Pernas, N. Sahoo, M. Vesterinen, J. A. Ward & A. R. Wiederhold

Physik-Institut, Universität Zürich, Zürich, Switzerland

C. Abellán Beteta, M. Atzeni, V. Bellee, R. Bernet, C. Betancourt, Ia Bezshyiko, A. Buonaura, D. De Simone, V. Denysenko, J. Eschle, M. Ferrillo, D. Lancierini, A. Mathad, K. Müller, P. Owen, N. Serra, R. Silva Coutinho, E. Smith, O. Steinkamp & Z. Wang

Oliver Lodge Laboratory, University of Liverpool, Liverpool, UK

T. Ackernley, T. J. V. Bowcock, G. Casse, A. J. Chadwick, C. A. Chavez Barajas, S. Farry, V. Franco Lima, T. Halewood-leagas, J. P. Hammerich, T. Harrison, K. Hennessy, D. Hutchcroft, P. J. Marshall, J. V. Mead, K. Rinnert, E. Rodrigues, S. Scherl, T. Shears, K. A. Thomson, E. Vilella Figueras & L. E. Yeomans

Instituto Galego de Física de Altas Enerxías (IGFAE), Universidade de Santiago de Compostela, Santiago de Compostela, Spain

B. Adeva, P. Baladron Rodriguez, O. Boente Garcia, A. Brea Rodriguez, A. Casais Vidal, V. Chobanova, X. Cid Vidal, I. Corredoira, J. Dalseno, L. Dieste Maronas, A. Fernandez Prieto, A. Gallas Torreira, B. Garcia Plana, A. Gioventù, J. Lomba Castro, S. Lopez Solino, D. Martinez Santos, T. Mombächer, A. Pereiro Castro, M. Plo Casasus, C. Prouve, E. R. R. Rodriguez Rodriguez, M. Romero Lamas, A. Romero Vidal, R. A. Ruiz Fernandez, J. J. Saborido Silva, C. Santamarina Rios, S. Sellam & P. Vazquez Regueiro

H.H. Wills Physics Laboratory, University of Bristol, Bristol, UK

M. Adinolfi, J. L. Amey, S. Bhasin, P. Camargo Magalhaes, M. G. Chapman, R. Lane, S. Maddrell-Mander, L. R. Madhan Mohan, A. M. Marshall, P. Naik, D. P. O’Hanlon, K. Petridis, J. H. Rademacker, A. R. Tanner, J. J. Velthuis, R. Wang, B. D. C. Westhenry & M. Whitehead

Université Clermont Auvergne, CNRS/IN2P3, LPC, Clermont-Ferrand, France

H. Afsharnia, Z. Ajaltouni, E. Cogneras, O. Deschamps, R. Lefèvre, J. Maratas, S. Monteil, P. Perret, H. Sazak, L. Soares Lavra & V. Tisserand

LPNHE, Sorbonne Université, Paris Diderot Sorbonne Paris Cité, CNRS/IN2P3, Paris, France

C. Agapopoulou, E. Ben-Haim, P. Billoir, L. Calefice, M. Charles, L. Del Buono, S. Esen, M. Fontana, V. V. Gligorov, T. Grammatico, F. Polci, A. Scarabotto, D. Y. Tou, P. Vincent & S. G. Weber

University of Michigan, Ann Arbor, USA

C. A. Aidala, D. S. Fitzgerald, S. H. Lee, K. R. Mattioli, C. Nunez, J. D. Roth & D. M. Shangase

INFN Sezione di Milano, Milano, Italy

S. Aiola, P. Gandini, D. Marangotto, A. Merli, N. Neri, M. Petruzzo, P. K. Resmi & E. Spadaro Norella

University of Cincinnati, Cincinnati, OH, USA

S. Akar, T. Boettcher, P. Ilten, B. Meadows, C. Pappenheimer, M. D. Sokoloff, M. Stahl & D. Vieira

Fakultät Physik, Technische Universität Dortmund, Dortmund, Germany

J. Albrecht, A. Battig, M. Becker, M. S. Bieker, J. A. Boelhauve, L. Calefice, E. Dall’Occo, F. De Vellis, H.-P. Dembinski, Q. Fuehring, L. L. Gerken, K. Heinicke, S. E. Hollitt, P. Ibis, V. Jevtic, B. Khanji, J. Langer, V. Lisovskyi, P. Mackowiak, G. Meier, D. S. Mitzel, A. Mödden, M. Saur, A. Seuthe, B. Spaan, H. Stevens & D. Wiedner

European Organization for Nuclear Research (CERN), Geneva, Switzerland

F. Alessio, L. An, C. Bertella, F. Betti, M. P. Blago, N. Bondar, A. Boyer, C. Burr, J. Buytaert, W. Byczynski, M. Cattaneo, G. Cavallero, Ph Charpentier, G. Ciezarek, M. Clemencic, J. Closier, V. Coco, P. Collins, T. Colombo, G. Corti, B. Couturier, C. D’Ambrosio, P. d’Argent, H. Dijkstra, L. Dufour, P. Durante, T. Evans, M. Féo, M. Ferro-Luzzi, A. Fkiaras, R. Forty, M. Frank, C. Frei, S. Gambetta, C. Gaspar, C. Giugliano, E. L. Gkougkousis, E. Govorkova, L. A. Granado Cardoso, T. Gys, C. Haen, J. Haimberger, C. Hasse, M. Hatch, A. M. Hennequin, L. Henry, R. Jacobsson, S. Jakobsen, D. Johnson, B. Jost, N. Jurik, M. Karacson, F. Keizer, R. D. Krawczyk, D. Lacarrere, R. Lindner, B. Malecki, C. Marin Benito, L. Martinazzoli, R. Matev, M. Mazurek, M. Mikhasenko, M. Milovanovic, F. Muheim, M. Mulder, D. Müller, P. Muzzetto, N. Neufeld, P. R. Pais, A. Pearce, M. Pepe Altarelli, M. Petric, F. Pisani, M. Pizzichemi, S. Ponce, D. Popov, L. Promberger, P. Roloff, T. Ruf, M. Salomoni, H. Schindler, B. Schmidt, A. Schopper, R. Schwemmer, B. Sciascia, F. Stagni, S. Stahl, M. Szymanski, F. Teubert, E. Thomas, N. Tuning, A. Valassi, C. Vázquez Sierra, K. Wyllie & O. Zenaiev

School of Physics and Astronomy, University of Glasgow, Glasgow, UK

M. Alexander, D. Bobulska, G. Coombs, L. Douglas, L. Eklund, D. A. Friday, I. Longstaff, N. T. Mchugh, M. Petric, M. Schiller, F. J. P. Soler & P. Spradlin

ICCUB, Universitat de Barcelona, Barcelona, Spain

A. Alfonso Albero, P. Garcia Moreno, L. Garrido, P. Gironella Gironell, S. Gomez Fernandez, E. Graugés, A. Lobo Salvia, J. Mauricio & R. Vazquez Gomez

Department of Physics and Astronomy, University of Manchester, Manchester, UK

Z. Aliouche, R. J. Barlow, A. Bitadze, S. Borghi, A. T. Burke, J. L. Cobbledick, J. E. Davies, A. Davis, O. De Aguiar Francisco, S. De Capua, D. Dutta, C. Fitzpatrick, E. Gersabeck, M. Gersabeck, L. Grillo, E. B. Hansen, M. Hilton, G. Lafferty, J. J. Lane, A. Lupato, T. H. Mcgrath, A. McNab, B. Mitreska, D. Murray, Y. Pan, C. Parkes, F. Reiss, N. Skidmore, P. Svihra, S. Taneja, D. J. White & G. Zunica

Petersburg Nuclear Physics Institute NRC Kurchatov Institute (PNPI NRC KI), Gatchina, Russia

G. Alkhazov, A. Andreianov, N. Bondar, A. Chubykin, V. Chulikov, A. Dzyuba, D. Ilin, A. Inglessi, K. Ivshin, S. Kotriakhova, P. Kravchenko, O. Maev, D. Maisuzenko, N. Sagidova, A. Solovev, I. Solovyev, A. Vorobyev & N. Voropaev

Cavendish Laboratory, University of Cambridge, Cambridge, UK

P. Alvarez Cartelle, F. C. R. Bishop, H. V. Cliff, B. Delaney, D. Foulds-Holt, V. Gibson, C. R. Jones, G. H. Lovell, J. G. Smeaton, I. Williams & S. A. Wotton

Universidade Federal do Rio de Janeiro (UFRJ), Rio de Janeiro, Brazil

S. Amato, L. De Paula, F. Ferreira Rodrigues, M. Gandelman, A. Hicheur, J. H. Lopes, L. Meyer Garcia, I. Nasteva, J. M. Otalora Goicochea, E. Polycarpo, M. S. Rangel, F. L. Souza De Almeida & B. Souza De Paula

Université Paris-Saclay, CNRS/IN2P3, IJCLab, Orsay, France

Y. Amhis, S. Barsuk, J. A. B. Coelho, F. Desse, M. Guittiere, F. Machefert, E. M. Niel, V. Puill, P. Robbe, M. H. Schune, Y. Shimizu, F. C. Volle, G. Wormser & V. Zhovkovska

INFN Sezione di Firenze, Firenze, Italy

L. Anderlini, M. Barbetti, A. Bizzeti, G. Graziani, C. Lucarelli, S. Mariani, G. Passaleva & M. Veltri

INFN Sezione di Ferrara, Ferrara, Italy

M. Andreotti, W. Baldini, C. Bozzi, R. Calabrese, M. Fiorini, E. Franzoso, C. Giugliano, M. Guarise, S. Kotriakhova, E. Luppi, L. Minzoni, I. Neri, S. Okamura, L. L. Pappalardo, B. Passalacqua, B. G. Siddi, I. Slazyk, L. Tomassetti & S. Vecchi

Physikalisches Institut, Ruprecht-Karls-Universität Heidelberg, Heidelberg, Germany

F. Archilli, S. Bachmann, D. Berninghoff, M. Borsato, G. Frau, D. Gerick, J. P. Grabowski, P. A. Günther, X. Han, S. Hansmann-Menzemer, J. Hu, R. Kopecna, B. Leverington, P. Li, H. Malygina, J. Marks, D. J. Unverzagt, U. Uwer, R. Waldi, C. Wang, L. Witola & A. Zhelezov

Institute for High Energy Physics NRC Kurchatov Institute (IHEP NRC KI), Protvino, Russia, Protvino, Russia

A. Artamonov, K. Belous, Y. Guz, S. Kholodenko, V. Obraztsov, S. Poslavskii, V. Romanovskiy, M. Shapkin, O. Stenyakin & O. Yushchenko

Syracuse University, Syracuse, NY, USA

M. Artuso, B. Batsukh, A. Beiter, H. C. Bernstein, S. Blusk, S. Ely, D. M. Keller, Z. Li, X. Liang, R. Mountain, M. E. Olivares, T. T. H. Pham, M. Poliakova, I. Polyakov, M. S. Rudolph, J. D. Shupperd, T. Skwarnicki, S. Stone, A. Venkateswaran, M. Wilkinson, H. Wu, Y. Yao & X. Yuan

Yandex School of Data Analysis, Moscow, Russia

K. Arzymatov, V. Belavin, M. Borisyak, S. Popov, F. Ratnikov & A. Ustyuzhanin

Aix Marseille Univ, CNRS/IN2P3, CPPM, Marseille, France

E. Aslanides, J. Cerasoli, J. Cogan, V. Dedu, D. Gerstel, R. Le Gac, O. Leroy, G. Mancinelli, A. Poluektov & D. Vom Bruch

Laboratoire Leprince-Ringuet, CNRS/IN2P3, Ecole Polytechnique, Institut Polytechnique de Paris, Palaiseau, France

B. Audurier, V. Balagura, F. Fleuret, F. A. Garcia Rosales & E. Maurice

Institute of Physics, Ecole Polytechnique Fédérale de Lausanne (EPFL), Lausanne, Switzerland

M. Bachmayer, A. Bay, F. Blanc, S. A. Bouchiba, S. Celani, S. Cholak, M. De Cian, S. Ek-In, L. Ferreira Lopes, E. Graverini, G. Haefeli, D. Hill, V. S. Kirsebom, V. Macko, M. Marinangeli, T. Nakada, T. Nanut, G. Pietrzyk, R. Quagliani, F. Redi, A. B. Rodrigues, O. Schneider, S. Schulte, L. Shchutska, C. Trippl, A. Tully, M. Van Dijk, M. Vieites Diaz & E. Zaffaroni

Centro Brasileiro de Pesquisas Físicas (CBPF), Rio de Janeiro, Brazil

J. Baptista Leite, I. Bediaga, M. Cruz Torres, J. M. De Miranda, A. C. dos Reis, A. Gomes, A. Massafferri & D. Torres Machado

Università di Firenze, Firenze, Italy

M. Barbetti, C. Lucarelli & S. Mariani

Imperial College London, London, UK

W. Barter, M. Birch, M. J. Bradley, A. Golutvin, M. Hecker, F. Kress, T. Lin, M. McCann, R. D. Moise, R. Newcombe, M. Patel, M. Smith & D. Websdale

INFN Sezione di Genova, Genova, Italy

M. Bartolini, R. Cardinale, F. Fontanelli, A. Petrolini & A. Sergi

Università di Genova, Genova, Italy

M. Bartolini, R. Cardinale, F. Fontanelli, A. Petrolini & A. Sergi

National University of Science and Technology “MISIS”, Moscow, Russia

F. Baryshnikov, S. Dadabaev, S. Didenko, A. Golutvin, S. Gromov, A. Ishteev, A. Kondybayeva, O. Lantwin, S. Legotin, N. Polukhina, D. Saranin, I. Shchemerov, E. Shmanin, O. Steinkamp, D. Strekalina, A. Terentev, E. Trifonova, E. Ursov, E. van Herwijnen & A. Zharkova

I. Physikalisches Institut, RWTH Aachen University, Aachen, Germany

J. M. Basels, L. Carus, S. Escher, A. Guth, J. Heuel, H. Jage, T. Kirn, S. Kretzschmar, C. Langenbruch, M. Materok, S. Nieswand, T. Oeser, G. Sarpis, S. Schael, S. Schmitt & V. Zhukov

AGH - University of Science and Technology, Faculty of Physics and Applied Computer Science, Kraków, Poland

S. Bashir, A. Dendek, M. Firlej, T. Fiutowski, M. Idzik, S. H. Kalavan Kadavath, P. Kopciewicz, W. Krupa, O. Madejczyk, M. W. Majewski, J. Moron, A. Oblakowska-Mucha, B. Rachwal, J. Ryzka, K. Swientek & T. Szumlak

INFN Sezione di Pisa, Pisa, Italy

G. Bassi, F. Bedeschi, N. Kleijne, F. Lazzari, A. Lusiani, M. J. Morello, L. Pica, G. Punzi, M. Rama, R. Ribatti & J. Walsh

INFN Sezione di Cagliari, Monserrato, Italy

S. Belin, D. Brundu, S. Cadeddu, A. Cardini, A. Contu, F. Dettori, F. Dordei, M. Garau, A. Lai, A. Lampis, R. Litvinov, A. Loi, G. Manca, P. Muzzetto, R. Oldeman, B. Saitta & J. Sun

Institute of Nuclear Physics, Moscow State University (SINP MSU), Moscow, Russia

I. Belov, A. Berezhnoy, I. V. Gorelov, M. Korolev, A. Leflat, N. Nikitin, D. Savrina & V. Zhukov

Institute of Theoretical and Experimental Physics NRC Kurchatov Institute (ITEP NRC KI), Moscow, Russia

I. Belyaev, A. Danilina, V. Egorychev, D. Golubkov, P. Gorbounov, A. Konoplyannikov, T. Kvaratskheliya, V. Matiunin, T. Ovsiannikova, D. Pereima, D. Savrina, A. Semennikov & A. Smetkina

INFN Laboratori Nazionali di Frascati, Frascati, Italy

G. Bencivenni, L. Calero Diaz, S. Cali, P. Campana, P. Ciambrone, P. De Simone, A. Di Cicco, P. Di Nezza, M. Giovannetti, G. Lanfranchi, G. Morello, M. Palutan, M. Poli Lener, M. Rotondo, M. Santimaria & B. Sciascia

Universita degli Studi di Padova, Universita e INFN, Padova, Padova, Italy

A. Bertolin, L. Buonincontri, L. Giambastiani, D. Lucchesi, M. Morandin, L. Sestini, G. Simi & D. Zuliani

Henryk Niewodniczanski Institute of Nuclear Physics Polish Academy of Sciences, Kraków, Poland

J. Bhom, J. T. Borsuk, J. Brodzicka, A. Chernov, M. Chrzaszcz, M. W. Dudek, A. Dziurda, M. Goncerz, M. Kucharczyk, T. Lesiak, J. J. Malczewski, K. Prasanth, M. Witek & M. Zdybal

School of Physics and Technology, Wuhan University, Wuhan, China

L. Bian, H. Cai, B. Fang, X. Huang, L. Sun & J. Wang

University of Birmingham, Birmingham, UK

S. Bifani, R. Calladine, G. Chatzikonstantinidis, N. Cooke, J. Plews, M. W. Slater, P. N. Swallow, D. J. D. Thompson & N. K. Watson

Università di Modena e Reggio Emilia, Modena, Italy

A. Bizzeti

Department of Physics, University of Oxford, Oxford, UK

M. Bjørn, D. Cervenkov, K. M. Fischer, A. L. Gilman, F. Goncalves Abrantes, B. R. Gruberg Cazon, T. H. Hancock, N. Harnew, M. John, S. Malde, R. A. Mohammed, C. H. Murphy, T. Pajero, M. Pili, H. Pullen, V. Renaudin, A. Rollings, L. G. Scantlebury Smead, J. C. Smallwood, S. Stanislaus, F. Suljik, M. D. Tat, G. Wilkinson & Y. Zhang

National Research University Higher School of Economics, Moscow, Russia

A. Boldyrev, D. Derkach, M. Hushchyn, M. Karpov, A. Maevskiy, S. Mokhnenko, F. Ratnikov, A. Ryzhikov & A. Ustyuzhanin

Budker Institute of Nuclear Physics (SB RAS), Novosibirsk, Russia

A. Bondar, S. Eidelman, P. Krokovny, V. Kudryavtsev, T. Maltsev, L. Shekhtman & V. Vorobyev

University of Maryland, College Park, MD, USA

S. Braun, A. D. Fernez, M. Franco Sevilla, P. M. Hamilton, A. Jawahery, W. Parker, Y. Sun & Z. Yang

Guangdong Provincial Key Laboratory of Nuclear Science, Guangdong-Hong Kong Joint Laboratory of Quantum Matter, Institute of Quantum Matter, South China Normal University, Guangzhou, China

A. Bursche, J. Hu, H. Li & G. Liu

Institute for Nuclear Research of the Russian Academy of Sciences (INR RAS), Moscow, Russia

A. Butkevich, S. Filippov, E. Gushchin, L. Kravchuk & S. Luchuk

Università di Ferrara, Ferrara, Italy

R. Calabrese, M. Fiorini, C. Giugliano, E. Luppi, L. Minzoni, L. L. Pappalardo, I. Slazyk & L. Tomassetti

INFN Sezione di Milano-Bicocca, Milano, Italy

M. Calvi, S. Capelli, P. Carniti, D. Fazzini, J. García Pardiñas, C. Gotti, M. Martinelli, C. Matteuzzi, S. Meloni, A. Minotti, M. Pizzichemi & E. B. Shields

Università di Milano Bicocca, Milano, Italy

M. Calvi, S. Capelli, P. Carniti, D. Fazzini, J. García Pardiñas, L. Martinazzoli, M. Martinelli, S. Meloni, A. Minotti, M. Pizzichemi & E. B. Shields

DS4DS, La Salle, Universitat Ramon Llull, Barcelona, Spain

M. Calvo Gomez, E. Golobardes, N. Valls Canudas & X. Vilasis-Cardona

University of Chinese Academy of Sciences, Beijing, China

A. F. Campoverde Quezada, Y. Fan, J. Fu, J. He, W. Huang, H. Liu, Q. Liu, X. Lyu, R. Ma, D. Popov, W. Qian, J. Qin, R. I. Rabadan Trejo, G. Tuci, Z. Wang, Z. Xiang, J. Xu, Q. Xu, Z. Xu, S. Yang, Y. Yang, Y. Zheng, X. Zhou, Y. Zhou & Z. Zhu

INFN Sezione di Bologna, Bologna, Italy

L. Capriotti, A. Carbone, A. Falabella, F. Ferrari, D. Galli, S. Maccolini, D. Manuzzi, U. Marconi, C. Patrignani, S. Perazzini, M. Soares, V. Vagnoni, G. Valenti, A. Villa & S. Zucchelli

Università di Bologna, Bologna, Italy

L. Capriotti, A. Carbone, F. Ferrari, D. Galli, S. Maccolini, D. Manuzzi, C. Patrignani & S. Zucchelli

INFN Sezione di Roma Tor Vergata, Roma, Italy

G. Carboni, E. Santovetti & A. Satta

Institute Of High Energy Physics (IHEP), Beijing, China

I. Carli, S. Chen, N. A. Grieser, Y. Li, S. Liu, Y. Lu, L. Ma, M. Tobin, J. Wang & Q. Zou

Univ. Savoie Mont Blanc, CNRS, IN2P3-LAPP, Annecy, France

M. Chefdeville, D. Decamp, A. G. Downes, Y. Hou, J. F. Marchand, M.-N. Minard, C. Normand, B. Pietrzyk, B. Quintana, M. Reboud & S. T’Jampens

Center for High Energy Physics, Tsinghua University, Beijing, China

C. Chen, C. Dong, J. Fan, Y. Gan, G. Gong, C. Gu, X. Hu, F. Jiang, Y. Kang, Y. Luo, H. Mu, H. Qi, N. Qin, Z. Ren, J. Wang, M. Wang, Z. Wang, L. Xu, D. Yang, M. Zeng, L. Zhang & X. Zhu

School of Physics and Astronomy, University of Edinburgh, Edinburgh, UK

P. E. L. Clarke, R. Currie, S. Eisenhardt, S. Gambetta, K. Gizdov, S. E. Mitchell, F. Muheim, M. Needham, F. Oliva, R. H. O’neil, S. Petrucci, G. Robertson, M. R. J. Williams & J. B. Zonneveld

Horia Hulubei National Institute of Physics and Nuclear Engineering, Bucharest-Magurele, Romania

L. Cojocariu, A. Ene, L. Giubega, A. Grecu, F. Maciuc, V. Placinta & M. Straticiuc

INFN Sezione di Bari, Bari, Italy

L. Congedo, M. De Serio, F. Debernardis, R. A. Fini, G. Galati, A. Palano, M. Pappagallo, A. Pastore & S. Simone

Università di Bari, Bari, Italy

L. Congedo, M. De Serio, F. Debernardis, G. Galati, M. Pappagallo & S. Simone

Massachusetts Institute of Technology, Cambridge, MA, USA

D. C. Craik, O. Kitouni, N. S. Nolte, C. Weisser & M. Williams

Los Alamos National Laboratory (LANL), Los Alamos, USA

J. Crkovská, C. L. Da Silva, C. T. Dean, J. M. Durham, E. Epple & G. J. Kunde

Physics and Micro Electronic College, Hunan University, Changsha City, China

L. Dai, J. Yu & S. Zhang

Van Swinderen Institute, University of Groningen, Groningen, Netherlands

K. De Bruyn, C. J. G. Onderwater & M. van Veghel

Universiteit Maastricht, Maastricht, Netherlands

J. A. de Vries, M. Merk, C. J. Pawley & K. Vos

Università di Cagliari, Cagliari, Italy

F. Dettori, G. Manca, R. Oldeman & B. Saitta

Eotvos Lorand University, Budapest, Hungary

B. Dey

Institute for Nuclear Research of the National Academy of Sciences (KINR), Kyiv, Ukraine

V. Dobishuk, S. Koliiev, I. Kostiuk, O. Kot, V. Lukashenko & V. Pugatch

School of Physics, University College Dublin, Dublin, Ireland

A. M. Donohoe, L. Mcconnell, R. McNulty, N. V. Raab & C. B. Van Hulse

NSC Kharkiv Institute of Physics and Technology (NSC KIPT), Kharkiv, Ukraine

A. Dovbnya & S. Kandybei

INFN Sezione di Perugia, Perugia, Italy

V. Duk, G. Martelli & M. Piccini

STFC Rutherford Appleton Laboratory, Didcot, UK

S. Easo, R. Nandakumar, A. Papanestis, S. Ricciardi & F. F. Wilson

School of Physics and Astronomy, Monash University, Melbourne, Australia

U. Egede & T. Hadavizadeh

Novosibirsk State University, Novosibirsk, Russia

S. Eidelman, P. Krokovny, V. Kudryavtsev, T. Maltsev, L. Shekhtman & V. Vorobyev

Department of Physics and Astronomy, Uppsala University, Uppsala, Sweden

L. Eklund

School of Physics State Key Laboratory of Nuclear Physics and Technology, Peking University, Beijing, China

Y. Gao, Y. Shang, Z. Shen, X. Tong, J. Wang, A. Xu, Z. Xu, Z. Yang, S. Zhang, Y. Zhang & T. Zhou

National Center for Nuclear Research (NCBJ), Warsaw, Poland

H. K. Giemza, K. Klimaszewski, M. R. Kmiec, W. Krzemien, D. Melnychuk, A. Szabelski, A. Ukleja & W. Wislicki

Università di Roma Tor Vergata, Roma, Italy

M. Giovannetti & E. Santovetti

Pontifícia Universidade Católica do Rio de Janeiro (PUC-Rio), Rio de Janeiro, Brazil

C. Göbel

Universidade Federal do Triângulo Mineiro (UFTM), Uberaba-MG, Uberaba, Brazil

A. Gomes

Institute of Particle Physics, Central China Normal University, Wuhan, Hubei, China

Q. Han, R. Hou, W. Hu, S. Li, Y. Wang, D. Xiao, Y. Xie, M. Xu, H. Yin, D. Zhang & X. Zhu

Hangzhou Institute for Advanced Study, UCAS, Hangzhou, China

J. He

Instituto de Fisica Corpuscular, Centro Mixto Universidad de Valencia - CSIC, Valencia, Spain

B. K. Jashal, F. Martinez Vidal, A. Oyanguren, C. Remon Alepuz & J. Ruiz Vidal

Nikhef National Institute for Subatomic Physics and VU University Amsterdam, Amsterdam, Netherlands

T. Ketel, H. S. Kuindersma, G. Raven & M. Senghi Soares

National Research Tomsk Polytechnic University, Tomsk, Russia

A. Kharisova, G. Panshin & A. Vagner

Università di Siena, Siena, Italy

F. Lazzari

Università di Padova, Padova, Italy

D. Lucchesi

Scuola Normale Superiore, Pisa, Italy

A. Lusiani, M. J. Morello & L. Pica

National Research Centre Kurchatov Institute, Moscow, Russia

A. Malinin, A. Petrov & V. Shevchenko

Università degli Studi di Milano, Milano, Italy

D. Marangotto & N. Neri

MSU - Iligan Institute of Technology (MSU-IIT), Iligan, Philippines

J. Maratas

INFN Sezione di Roma La Sapienza, Roma, Italy

G. Martellotti, D. Pinci, R. Santacesaria, A. Sarti & C. Satriano

Institut für Physik, Universität Rostock, Rostock, Germany

N. Meinert & H. Viemann

Departamento de Fisica, Universidad Nacional de Colombia, Bogota, Colombia

D. A. Milanes, I. A. Monroy & J. A. Rodriguez Lopez

Universität Bonn - Helmholtz-Institut für Strahlen und Kernphysik, Bonn, Germany

A. B. Morris, S. Neubert, K. O. Padeken & M. Sarpis

P.N. Lebedev Physical Institute, Russian Academy of Science (LPI RAS), Moscow, Russia

N. Polukhina & M. Zavertyaev

Università di Pisa, Pisa, Italy

G. Punzi

Università della Basilicata, Potenza, Italy

C. Satriano

Max-Planck-Institut für Kernphysik (MPIK), Heidelberg, Germany

M. Schmelling & M. Zavertyaev

Università di Urbino, Urbino, Italy

M. Veltri

All contributing authors, as listed at the end of this manuscript, have contributed to the publication, being variously involved in the design and the construction of the detector, in writing software, calibrating sub-systems, operating the detector and acquiring data and finally analysing the processed data.

Correspondence to I. Belyaev.

The authors declare no competing interests.

 Nature Communications thanks Iain Bertram and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. This article has been peer reviewed as part of Springer Nature’s Guided Open Access initiative.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

LHCb collaboration. Study of the doubly charmed tetraquark \({{{{{{\rm{T}}}}}}}_{{{{{{\rm{c}}}}}}{{{{{\rm{c}}}}}}}^{+}\). Nat Commun 13, 3351 (2022). https://doi.org/10.1038/s41467-022-30206-w

Download citation

Received: 07 September 2021

Accepted: 14 April 2022

Published: 16 June 2022

DOI: https://doi.org/10.1038/s41467-022-30206-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Frontiers of Physics (2023)

Frontiers of Physics (2023)

Nature Italy (2022)

Nature Physics (2022)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.