The isospin and compositeness of the T cc (3875) state

We perform a ﬁt to the LHCb data on the T cc (3875) state in order to determine its nature. We use a general framework that allows to have the D 0 D ∗ + , D + D ∗ 0 components forming a molecular state, as well as a possible nonmolecular state or contributions from missing coupled channels. From the ﬁts to the data we conclude that the state observed is clearly of molecular nature from the D 0 D ∗ + , D + D ∗ 0 components and the possible contribution of a nonmolecular state or missing channels is smaller than 3%, compatible with zero. We also determine that the state has isospin I = 0 with a minor isospin breaking from the diﬀerent masses of the channels involved, and the probabilities of the D 0 D ∗ + , D + D ∗ 0 channels are of the order of 69% and 29% with uncertainties of 1%. The diﬀerences between these probabilities should not be interpreted as a measure of the isospin violation. Due to the short range of the strong interaction where the isospin is manifested, the isospin nature is provided by the couplings of the state found to the D 0 D ∗ + , D + D ∗ 0 components, and our results for these couplings indicate that we have an I = 0 state with a very small isospin breaking. We also ﬁnd that the potential obtained provides a repulsive interaction in I = 1, preventing the formation of an I = 1 state, in agreement with what is observed in the experiment.


I. INTRODUCTION
The discovery of the T cc (3875) reported in [1,2] was a turning point in hadron physics, showing the first evidence of a meson state clearly exotic with two open charm quarks.The state, showing as a very narrow peak in the D 0 D 0 π + mass distribution very close to the D 0 D * + and D + D * 0 thresholds, has spurred lots of work in the hadron community, trying to find out its origin and nature.It should be quoted that prior to this discovery there were many works looking for tetraquark structures which predicted the existence of a state of this nature with two heavy quarks [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17], yet with a large span in the mass from about 250 MeV below to about 250 MeV above the actual T cc (3875) mass.There were also works predicting such an object as a molecular state stemming from the interaction of the D 0 D * + and D + D * 0 pairs in [18][19][20].Actually, the proximity of the mass of the T cc (3875) to the D 0 D * + and D + D * 0 thresholds will guarantee that a molecule is made provided a reasonable attractive force between the pseudoscalar and vector-meson components appears [21].The existence of such a state could also be concluded from the results of the D * D * interaction in [22], where an extrapolation of the local hidden gauge approach for meson-meson interaction [23][24][25][26] was used and sufficient attraction was found to form a bound state.The extreme proximity of the T cc (3875) state to the D 0 D * + threshold, of about 360 keV, has motivated many works to explain this state as a molecular state of the D 0 D * + and D + D * 0 channels .The compact tetraquark structure has also had its supporters in [51][52][53][54].Some works suggest that there should be a molecular state with a small admixture of a compact tetraquark [55].A comparative study of the tetraquark and molecular structures is done in [56] showing that the compact tetraquarks are about 100 MeV more bound than the molecular ones.Other works suggest to measure magnetic moments to elucidate the nature of the states [49,57].As mentioned, the proximity of the mass of the T cc (3875) state to the D 0 D * + , D + D * 0 thresholds, makes the molecular picture very natural and this can explain the large amount of works supporting this picture, mostly by using effective meson-meson potentials, or quark model interaction [50], QCD sum rules [38,49] or lattice QCD simulations [40,44].Further discussions concerning the T cc (3875) and related states can be seen in [54,[58][59][60][61].
While the proximity of the T cc (3875) to the D 0 D * + and D + D * 0 thresholds, together with a moderate attractive interaction, will lead to a molecular state of these components [21], one cannot rule out that there could be a mixture, even if small, of a genuine compact tetraquark structure, as suggested in [55].The issue of the compositeness and elementariness of states has drawn much attention, starting by the well known work of Weinberg [62] supporting the deuteron as a bound state of a proton and a neutron.Improvements on the original idea of [62] have been done in [41,.In [84] the effect of the finite range of the interaction, implicitly assumed to be of zero range in coordinate space in [62], is investigated and found to be relevant in the determination of the compositeness.
For bound states the probability to have a certain component in the wave function is given by −g2 ∂G/∂s, where g is the coupling of the state to the given component and G the loop function of the propagators of the two hadrons of the component, and for a potential that is energy independent the sum rule −g 2 i ∂G i /∂s = 1 holds, indicating that the state is fully molecular and the sum of probabilities of the different components is unity [68,78,83].However, if the potential is energy dependent, the former sum rule can be smaller than 1, accounting for the possibility of having other components in the wave function.In fact, it is possible to eliminate some channel resulting in an energy dependent effective potential for the remaining channels, and the sum rule of these latter channels gives unity minus the probability of the eliminated channel [69,74,83,85].Similarly, the existence of a contribution from a genuine, or compact tetraquark in this case, can also be taken into account eliminating this channel in favor of an energy effective potential for the remaining channels [86].
The works assuming a molecular picture for the T cc (3875) ordinarily assume energy independent interaction potentials between the components of the molecule, and hence, the sum rule is satisfied and the sum of probabilities of the D 0 D * + and D + D * 0 channels is unity.In the case of [34], the probabilities of these states are 71% and 29% respectively1 , or the case of [42], where the probabilities are also of the order of 72% and 27% respectively, or in [87] with 69% and 31% 2 , or in [33] where the probabilities are 80% and 20%.There is a similarity in all these results, and the probability of the first channel is bigger than the second channel because its threshold is closer to the energy of the state.This difference has not to be interpreted as a large violation of isospin.Indeed the T cc (3875) is assumed to have isospin I = 0 and no trace of a peak is seen in the partner isospin I = 1 channel D + D * + .Isospin is a symmetry of strong interactions, which are of short range, and what matters in reactions where isospin is evidenced is the wave function at the origin (for s-wave), g G according to [68].With values of G similar for the two channels here, the couplings are what determine the isospin nature of the states.From [87] we obtain values of g i for each channel which are of opposite sign and equal in size within 2-3%, as it corresponds to an I = 0 state.Since the probability is given by −g 2 ∂G/∂s and ∂G/∂s goes to infinity as we approach the threshold of the channel, it comes natural that, even having good isospin, the probability of the channel whose threshold is closer to the mass of the state will be larger than the one of the other channel.This is also what happens with the D 0 D * 0 and D + D * − components in the X(3872) [88].The other thing to consider is that what determines the molecular nature of the state is the sum of the probabilities of the two channels that form the state.
Very recently, two papers have addressed the issue of the probability of the D 0 D * + and D + D * 0 components.One of them is [89], where only the probability of the first channel is evaluated with a result of about 50-80%.Note that in this case the Weinberg prescription is used, which according to many posterior works, has been shown to be inaccurate in many cases 3 .The other work is [90], where a general discussion is done based on the value of the binding energy alone, together with several other assumptions made to validate the results.As we shall see, we will base our conclusions on two types of fits to the data, one fitting the binding, scattering lengths and effective ranges, seven data in our case, and another one making a fit to the D 0 D 0 π + distribution itself.This information obviously contains more dynamical information to make conclusions than the binding energy alone, which is clearly insufficient, since one can easily have a genuine state with a given binding energy, provided there is no interaction between the meson-meson components, or a pure molecular state from the interaction of the meson-meson components, provided there is no contribution from a genuine state, or such a state is very far away from the mass of the observed state (see further discussions along this line in [91].We must also stress that in [90] the effect of the interaction range is also investigated and it is found that, in addition to the binding energy of the T cc , if the range of the interaction in momentum space is assumed to be of the order of 770 MeV, as it would correspond to vector exchange, the compositeness of the T cc is much bigger than assuming a range due to pion exchange, and becomes of the order of 94%.
The former discussion shows that the subject is ripe to make a thorough discussion of the issue, which is our purpose here.For this purpose we shall include a general potential for the interaction of the D 0 D * + and D + D * 0 channels, including the necessary terms with energy dependence to account for possible nonmolecular components, and carry fits to data, from where we will evaluate the probabilities of the two channels and Z, the probability of other possible nonmolecular components.We shall see that the results strongly support the molecular picture, with a negligible probability of a nonmolecular component.

II. FORMALISM
The T cc (3875) is found according to Ref. [2] with mass The state is in J P = 1 + and is assumed to be in isospin I = 0, since no signal is seen in the I = 1, D + D * + channel.The two channels close to the T cc mass are D 0 D * + and D + D * 0 , which we consider explicitly in our study.We have the isospin doublets (D + , −D 0 ),(D * + , −D * 0 ) and hence In our study we shall not assume a priori that the T cc is a state of I = 0, it will come out from the analysis of the experimental data.We shall only assume that isospin is a good symmetry for the interaction such that the potential can be diagonalized with the isospin basis.The different masses of D 0 D * + and D + D * 0 , that go in the kinetic energy of the Hamiltonian, can, and will, produce a small violation of isospin in the final wave function.This said, we take a general potential for the interaction of the two channels given by sign to the experimental one (see details in [84]).However, as shown in [65], the prescription works better when the effective range is negative.
from where the scattering matrix is with , where G i are the loop functions for the two channels which we regularize in the cutoff method, with where ω i = q 2 + m 2 i , m 1 is the mass of the D and m 2 that of D * .The value of q max reflects the range of the interaction in momentum space [68,84] and will be obtained from the fits to the data.Since we assume isospin symmetry for the potential we shall have We also find that According to the findings of [2] we should expect that V 11 − V 12 is negative, to allow for a bound state, and V 11 + V 12 should be positive to prevent a bound state, from where we should have However, we could relax these conditions, since, as we shall see, the state that we obtain has I = 0 and all we need is that V 11 − V 12 is negative and V 11 + V 12 relatively smaller in size than V 11 − V 12 , such that it cannot lead to any bound state.
Then we have three parameters in the formalism, V 11 , V 12 and q max .But to take into account the possibility of a genuine component (compact tetraquark) or contribution from other meson-meson channels different than the two explicitly considered, we will add terms in the potential which are energy dependent.According to the discussion in [84,86] it is a good approximation to retain just the terms linear in s.Thus we write where α, β are free parameters, m V = 800 MeV is a typical vector mass, chosen to have α, β dimensionless, and s 0 is the mass squared of the T cc .Our formalism has now five free parameters, but demanding that we have a pole at s 0 , they will be reduced to four.For the experimental data we will follow two different strategies: (a) a fit to the a and r 0 parameters of the effective range expansion for the two channels, D 0 D * + and D + D * 0 , obtained from the experimental analysis of the data in Ref. [2]; (b) a fit to the unitary amplitudes extracted in the analysis of the experimental data in Ref. [2].
The scattering length a and effective range r 0 correspond to the effective range expansion of the scattering matrix Let us also mention here that in our formalism with Eqs. ( 5), (6) for the scattering amplitude we have the relationship which one can easily induce by looking at the imaginary parts of Im T −1 (s) = −Im G (for a real potential) and f −1 (s).
To avoid unnecessary complications due to the small width of the D * , in the strategy (a) we shall disregard the D * width4 .In the direct fit to the unitary amplitude of [2] we shall take into account the D * width.The two fits are complementary and for fit (a) we rely on the fact that the very small width of the D * will not change the character of the T cc state.Actually, we are used to call bound state to atoms and nuclei with unstable particles.For strategy (a) we have a 1 , r 0,1 real for the D 0 D * + channel and a 2 , r 0,2 complex for the D + D * 0 channel.Hence, we have 6 data to be fitted with 4 free parameters, for which we minimize the χ 2 in the fit.The data of [2] neglecting the D * width are a 1 = 6.134 ± 0.51 fm , r 0,1 = −3.516± 0.50 fm , a 2 = (1.707± 0.30) − i (1.07 ± 0.30) fm , r 0,2 = (0.259 ± 0.30) − i (3.769 ± 0.30) fm . ( The errors for a 1 , r 0,1 are taken from [2], and for a 2 , r 0,2 are taken relatively of the same order of magnitude.The errors influence the value of χ 2 but not the parameters of the minimum of χ 2 when this reaches values close to zero, as is the case here. It is interesting to compare the results of Eq. ( 13) with those induced from [87] a 1 = 6.42 fm , r 0,1 = −3.49fm , The results are very similar although not equal.The real parts of r 0,2 are different and even with opposite sign, but let us note that they are one order of magnitude smaller than their corresponding imaginary parts, and not incompatible within the assumed uncertainties.Note that a 2 , r 0,2 are complex, even neglecting the D * width because the decay of D + D * 0 to D 0 D * + is open at the D + D * 0 threshold.The magnitudes of Eq. ( 14) are obtained with a potential where M, M ′ are the initial, final vector masses and m, m ′ the initial, final pseudoscalar masses, with s calculated at the threshold of D 0 D * + and average masses of D 0 , D + and D * 0 , D * + .The actual masses are used in the evaluation of the scattering matrices of Eqs. ( 5), (6).When we fit the magnitudes of Eq. ( 14) with the general potential of Eqs. ( 4), (10) and the T matrix of Eq. ( 5), we obtain a best fit with χ 2 = 0 with α = β = 0 , certainly the potential of Eq. ( 15), and the sum rule is exactly fulfilled indicating that we have a pure molecular state.This is actually a tautology, since we start with a potential that produces a molecule and we get the answer that it is indeed a molecule.However, the test tells us that the inverse problem of getting the potential parameters from the magnitudes a 1 , r 0,1 , a 2 , r 0,2 is feasible and meaningful.The interesting test is to obtain the potential of Eqs. ( 4), (10) from the experimental data of Eq. (13) where no assumptions concerning the interaction have been done, and then evaluate the probabilities P 1 , P 2 for the D 0 D * + and D + D * 0 channels.

III. EVALUATION OF THE SCATTERING MATRIX, THE COUPLINGS AND THE PROBABILITIES OF THE CHANNELS
The evaluation of the T matrix of Eq. ( 4) with two channels is easy and we show the formula below with the determinant of (1 − V G), DET, given by The bound state appears when DET = 0 at s 0 , hence from where we can eliminate V 12 in terms of V 11 , as (note that the α, β terms of Eq. ( 10) vanish at s = s 0 ) and we take the positive root of V 12 according to Eq. ( 9).The couplings are defined from the residues of the T matrix.Since then where in the second step in the former equations we have used L'Hôpital rule.Note that the last equation allows us to get the relative sign of g 1 and g 2 .Then, according to [68,83] we have the probabilities for the D 0 D * + and D + D * 0 channels, respectively, as and consequently, the nonmolecular part, that could also account for possible missing coupled channels, is With the explicit formulae for the couplings, Eq. ( 22), it is now easy to prove explicitly the sum rule of Eq. ( 16) when α = β = 0, although a general proof, extended to any number of channels, can be found in [68,83], identifying each term of the sum rule with the probability of the corresponding channel in the wave function of the state.

IV. EVALUATION OF THE SCATTERING LENGTHS AND EFFECTIVE RANGES
From Eqs. (11), (12) we find Hence, for each of the two channels D 0 D * + and D + D * 0 , we have with and then we can write with s i the square of the threshold energies for the two channels and where to calculate r 0,1 we have used that V. RESULTS

A. Fit to the data of a and r0
We first conduct the fit (a) to the data of a, r 0 of Mikhasenko [89] given in Eq. ( 13).We define where y exp i are a 1 , r 0,1 which are real, and a 2 , r 0,2 which are complex.In the end we have six data (the binding energy, was used in Eq. ( 20) to eliminate V ′ 12 ), and four parameters V ′ 11 , α, β, q max .The reduced χ 2 r is χ 2 /2.We minimize χ 2 r to obtain these parameters and, once obtained, we calculate P 1 , P 2 , Z. Given the fact that there are correlations between the parameters, as we shall see, we find useful in order to find the uncertainties in the output to use a method often used in these works [93][94][95][96][97][98][99][100].This consist in calculating χ 2 of Eq. ( 28), using instead of y exp i a random number within the interval [y exp i − ∆y exp i , y exp i + ∆y exp i ].We conduct several runs in this way, selecting the good ones with χ 2 r < 0.4, and for each run we determine the magnitudes P 1 , P 2 , Z (and also any other magnitude that one wished to obtain).The errors obtained for the different magnitudes in this way do not change if we take χ 2 r < 0.4, 0.3, 0.2, 0.1 and the average values are always within these error bars.Then, we evaluate the average value of these magnitudes and the dispersion in the usual way, for instance The results that we obtain are shown in Table I.The agreement of the values obtained for the scattering lengths and effective range with the experimental ones of Eq. ( 13) is remarkable.
As we see in Table II, we obtain P 1 , P 2 of the order of 69%, 29% with uncertainties of the order of 1%, and Z = 0.03 ± 0.02.The numbers obtained qualify the T cc as a clear molecular state made of the D 0 D * + and D + D * 0 components.We also find there that the couplings g 1 , g 2 are very close to each other and of opposite sign, indicating, according to Eq. ( 3) that we have basically a state of I = 0.It is illustrating to take one typical fit of those obtained.We find This contrasts with the fit obtained in [87] using the potential of Eq. ( 15) and α = β = 0, We should first realize that in a range of about 5 MeV from the energy of the T cc the terms α(s − s 0 ), β(s − s 0 ) are of the order of 16% of V 12 of Eq. ( 30), or 5% compared to V 12 of Eq. ( 31), hence, a relatively small magnitude.Yet, it looks like the solutions of Eqs.(30), (31) are very different.Indeed, the values of q max , although qualitatively similar, are still different, and the strength of the potential is also rather different.However, this only reflects the strong correlations between the parameters.Indeed, if we look only to the binding energy (but we have more information) and assuming we just have the I = 0 state that we have induced from the data, we would have in one channel where, according to Eq. ( 8)  FIG.1: Fit to the data of Fig. 8 of supplementary material in Ref. [2] with the size of the errors taken from the raw data in the way explained in the text.
One should not worry too much about these values and their errors, since we know that there are also correlations between V ′ 11 and V ′ 12 and α and β (the differences between these parameters is what matters if one has an I = 0 state).The relevant thing is what we get for the observables, and from the different fits in the resampling method we obtain the average values of the observables and their dispersion.The results are summarized in Tables III and IV.These results are very similar to those displayed in Table I.Yet, they should be better compared to the results of Ref. [2] when the D * width is explicitly taken into account, which are The agreement is perfect within errors if we assume the relative errors in a exp 2 to be similar to those of a exp 1 .We see that a 1 is now complex, yet with the imaginary part reasonably smaller than the real one.The real parts obtained for a 1 or a exp 1 in this case are also very similar to a 1 of Eqs. ( 13) and ( 14) obtained with fit (a) to the data neglecting the width of the D * .
We have taken advantage of the weak dependence of the obtained magnitudes with the D * decay width and evaluate P 1 , P 2 , P 1 + P 2 (P 1 + P 2 = 0.998 ± 0.025(±0.0004))5and r 0,1 , r 0,2 , g 1 , g 2 in the limit of Γ D * → 0, using the parameters obtained in Eq. ( 37) and the formulas described in section III.
We see that the systematic errors in the effective ranges are bigger than for the scattering lengths, as one might expect, given the dependence of these variables on q max .The values of the couplings obtained are also very similar to those shown in Table II with the fit (a), compatible within errors.Coming now to the compositeness P 1 , P 2 , the values obtained are also remarkably similar to those in Table II, indicating again the molecular nature of the T cc state with P 1 + P 2 essentially 1 in the present fit with an uncertainty of the order of 0.025.

VI. LIMITING CASE OF A NONMOLECULAR STATE
We shall assume that the origin of the T cc state corresponds to a state, which has a very small overlap with the D 0 D * + , D + D * 0 components (a minimum overlap is needed such that the state is observed in the D 0 D * + and D + D * 0 channels).We shall also assume, to be consistent with experiment, that the state has I = 0, although the conclusions are not tied to this fact.The amplitude for the DD * (I = 0) state will be given by depicted in Fig. 2, and g 2 should be small to prevent the overlap of the genuine state with the DD * components.Consistently with the nonmolecular assumption, we will also assume that the D 0 D * + , D + D * 0 components have no interaction stemming from a potential, or t-channel exchange of other hadrons.Yet, the picture for the T matrix is not complete since, for consistency with the small coupling of the genuine state to the DD * components, we have to consider the selfenergy of the state due to its coupling to the DD * components, as depicted in Fig. 3.
By taking into account the I = 0 structure of the DD * state, we find from Fig. 3 T Once again, by considering Eq. ( 25) and that Im T −1 = 1 8π √ s k, we realize that the amplitude of Eq. ( 39) satisfies unitarity in the two channels D 0 D * + , D + D * 0 , reproducing exactly the imaginary part of f −1 (s).We now derive the results for a, r 0 in the limit of g → 0, such that Eq. ( 39) has a pole at s 0 , where it should be.Let us call D the denominator of Eq. ( 39), In order to have a pole of T D 0 D * + ,D 0 D * + at s 0 , we need Let us note that if g ̸ = 0, then s 0 ̸ = s 0 and we can make an expansion of T in powers of s − s 0 , and we would then go back to the previous analysis accounting for the genuine state in terms of the energy dependent terms, the constant parts reabsorbed in the V ij coefficients.Thus, the novel thing should be found in the g 2 → 0 limit.
In this case it is easy to find a and r 0 from Eq. ( 39) Thus, we immediately see that − 1 a → ∞ and hence a → 0; a 1 → 0, a 2 → 0.
For the radius we find 1 2 r 0,1 = lim We see that in the strict limit of the nonmolecular state, the scattering length goes to 0 and the radius to −∞.This contrast with the experimental results of Eq. ( 13) where a is relatively large for hadronic reactions, and r 0 is smaller than a.A discussion on the effects of the genuine state along the line described here, but not addressing the scattering length and effective range can be seen in [33] (see also section 2 of Ref. [86]).

VII. CONCLUSIONS
We have performed fits to the data of the LHCb collaboration on the T cc (3875) state with a general framework that contains explicitly the D 0 D * + and D + D * 0 components, which are allowed to interact with a potential containing energy dependent terms by means of which one can account for contributions of missing coupled channel, as well as effects from some genuine state, like a compact tetraquark state.We have conducted two types of fits to the data.One of them is a fit to data from LHCb for the scattering lengths and effective ranges of the channels D 0 D * + and D + D * 0 , evaluated ignoring the D * width.The other fit is done directly to the data of the mass distribution of D 0 D 0 π + corrected by the experimental resolution and parametrized in terms of a unitary amplitude in [2], where the D * width is explicitly considered.In all fits we have calculated the parameters of the theoretical framework and from them we have looked at the bound state produced and evaluated the probabilities P 1 , P 2 of having D 0 D * + and D + D * 0 in the wave function, as well as the couplings g 1 , g 2 of the state obtained to these two components.Without any assumption of the nature of the state observed, we find that the couplings obtained are very close to each other and of opposite sign, indicating the isospin I = 0 nature of the state observed 6 .The potential obtained produces a repulsive interaction in I = 1, justifying the non observance of any structure with I = 1 in the LHCb experiment.On the other hand, from the couplings obtained, we could calculate the probabilities P 1 , P 2 of D 0 D * + , D + D * 0 and Z = 1−P 1 −P 2 , the non compositeness of the state.We found values for Z of the order of 3% with an uncertainty that makes it compatible with zero or, conversely, in the case of fit (b) P 1 + P 2 = 1 with uncertainties of 2.5%, indicating that we have a very clear case of a molecular state made of the D 0 D * + , D + D * 0 components.We found that the probabilities of the two channels are of the order of 69% and 29% respectively, very similar to what has been obtained in other approaches which, however, assume a priori energy independent potentials which do not allow nonmolecular components to appear.We also clarified that the isospin nature should not be induced by the probabilities P 1 , P 2 , but from the couplings g 1 , g 2 , or related to them, the wave functions at the origin of these two channels, since in strong interactions, where the isospin symmetry holds, the wave functions at short distances is what matters.We also shortly discussed that using all the available experimental information, apart from the mass of the state, was essential to reach the present conclusions, with the information on the binding being clearly insufficient to reach these conclusions.

TABLE I :
The obtained scattering lengths and effective ranges.

TABLE II :
The obtained coupling constants and probabilities.

TABLE III :
The obtained scattering lengths and effective ranges.The values between brackets correspond to the systematic error.

TABLE IV :
The obtained binding energy, width, coupling constants and probabilities.The values between brackets correspond to the systematic error.