Decay widths of the spin-2 partners of the X(3872)

We consider the X(3872) resonance as a JPC=1++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J^\mathrm{{PC}}=1^{++}$$\end{document}DD¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D\bar{D}^*$$\end{document} hadronic molecule. According to heavy quark spin symmetry, there will exist a partner with quantum numbers 2++\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2^{++}$$\end{document}, X2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_{2}$$\end{document}, which would be a D∗D¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D^*\bar{D}^*$$\end{document} loosely bound state. The X2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_{2}$$\end{document} is expected to decay dominantly into DD¯\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D\bar{D}$$\end{document}, DD¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D\bar{D}^*$$\end{document} and D¯D∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{D} D^*$$\end{document} in d-wave. In this work, we calculate the decay widths of the X2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_{2}$$\end{document} resonance into the above channels, as well as those of its bottom partner, Xb2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_{b2}$$\end{document}, the mass of which comes from assuming heavy flavor symmetry for the contact terms. We find partial widths of the X2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_{2}$$\end{document} and Xb2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_{b2}$$\end{document} of the order of a few MeV. Finally, we also study the radiative X2→DD¯∗γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_2\rightarrow D\bar{D}^{*}\gamma $$\end{document} and Xb2→B¯B∗γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_{b2} \rightarrow \bar{B} B^{*}\gamma $$\end{document} decays. These decay modes are more sensitive to the long-distance structure of the resonances and to the DD¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$D\bar{D}^{*}$$\end{document} or BB¯∗\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\bar{B}^{*}$$\end{document} final state interaction.


Introduction
In the infinite quark mass limit, heavy quark spin-flavor symmetry (HQSFS) implies that the dynamics involving heavy quarks is independent of their spin or flavor. In this way, charm and bottom spectra can be related, up to corrections suppressed by 1/m Q with m Q the heavy quark mass. It should also be possible, in a given heavy flavor sector, to relate states with different spins. These relations are very useful in the study of composites that mix heavy and light quarks. In this work, we focus on hadronic molecular states composed by a heavy-light meson and a heavylight antimeson (P ( * )P ( * ) , P = D,B). These molecular states were predicted in the mid 70s [1,2]. So far, the best experimental candidate to fit this molecular description is the ting between the first radially excited charmonia with 2 ++ and 1 ++ in the well-known Godfrey-Isgur quark model is 30 MeV [17], which is much smaller than the value in Eq. (1). In a quark model calculation with screened potential, the 2 ++ − 1 ++ mass splitting for the 2P charmonia is around 40 MeV [18]. As for the tetraquark states, the corresponding mass splitting predicted in Ref. [19] is 70 MeV, which is again much smaller than M D * − M D . Notice that it is generally believed that the χ c2 (2P) has been discovered [20,21], and its mass is much lower than 2M D * . Therefore, we conclude that a possible discovery of a 2 ++ charmonium-like state with a mass around 4013 MeV as a consequence of HQSS [22] would provide a strong support for the interpretation that the X (3872) is dominantly a DD * hadronic molecule. It is thus very important to search for such a tensor resonance, as well as the bottom analogs, in various experiments and in lattice QCD (LQCD) simulations.
Some exotic hidden charm sectors on the lattice have been recently studied [23][24][25][26][27], and evidence for the X (3872) from DD * scattering on the lattice has been found [24], while the quark mass dependence of the X (3872) binding energy was discussed in Refs. [28,29]. The 2 ++ sector has not been exhaustively addressed yet in LQCD, though a state with these quantum numbers and a mass of (m η c + 1041 ± 12) MeV= (4025 ± 12) MeV, close to the value predicted in Refs. [7,8], was reported in Ref. [23]. The simulation used dynamical fermions, novel computational techniques, and the variational method with a large basis of operators. The calculations were performed on two lattice volumes with pion mass 400 MeV. There exists also a feasibility study [30] of future LQCD simulations, where the EFT approach of Refs. [7,8] was formulated in a finite box.
On the other hand, despite the theoretical predictions on their existence, none of these hypothetical particles has been observed so far. Nevertheless, they are being and will be searched for in current and future experiments such as BESIII, LHCb, CMS, Belle-II, and PANDA. It is thus of paramount importance to provide theoretical estimates on their production rates in various experiments, as well as the dominant decay modes and widths. 3 The production of these states in hadron colliders and electron-positron collisions has been studied in Refs. [31,32]. In this work, we will investigate the dominant decay modes of the spin-2 partners of the X (3872), i.e. the X 2 (4013) and X b2 , and provide an estimate of their decay widths.
Besides, we will also discuss the radiative X 2 → DD * γ and X b2 →B B * γ transitions. These decay modes are more sensitive to the long-distance structure of the resonances and might provide valuable details of their wave functions. The situation is similar to that of the X (3872) → D 0D0 π 0 decay studied in Refs. [33,34]. Also here the widths will be affected by the DD * or BB * final state interaction (FSI). FSI effects are expected to be large because they should be enhanced by the presence of the isovector Z c (3900) and Z b (10610) resonances located near the DD * andB B * thresholds, respectively. Besides, FSI corrections will be also sensitive to the negative C-parity isoscalar DD * orB B * interaction. Eventually, precise measurements of these radiative decay widths might provide valuable information on the interaction strength in this sector, which would be important in understanding the P ( * )P ( * ) system and other exotic systems related to it through heavy quark symmetries [15,35].
The structure of the paper is as follows. First in Sect. 2, we briefly discuss the relation of the charm and bottom 2 ++ states with the X (3872) resonance, and in Sect. 3 we present our predictions for the X 2 → DD, DD * hadron decays and the X b2 → BB, BB * ones in the bottom sector. In Sect. 4, the X 2 and X b2 radiative decays are investigated, paying special attention to the loop mechanisms responsible for the FSI contributions. The conclusions of this work are outlined in Sect. 5 and in addition, there are three appendices. In the first one (Appendix A), we collect different heavy meson Lagrangians used through this work, while the validity of the perturbative treatment of the DD for the X 2 is discussed in the second one (Appendix B). Finally, in Appendix C, we give some details of the evaluation of different three-point loop functions that appear in the computation of the hadronic and radiative decays.

X (3872)
As mentioned, we start assuming the X (3872) to be a positive C-parity DD * bound state, with quantum numbers J PC = 1 ++ . At very low energies, the leading order (LO) interaction between pseudoscalar and vector charmed (D 0 , D + , D * 0 , D * + ) and anti-charmed (D 0 , D − ,D * 0 , D * − ) mesons can be described just in terms of a contact-range potential, which is constrained by HQSS [7,8,15]. Pion exchange and particle coupled-channel 4 effects turn out to be subleading [7,36]. For the case of the X (3872), isospin breaking is important [37] as this bound state is especially shallow. The energy gap between the D 0D * 0 and D + D * − channels is around 8 MeV, which is much larger than the X (3872) binding energy with respect to the D 0D * 0 threshold. As a consequence, the neu-tral (D 0D * 0 ) and charged (D + D * − ) channels should be treated independently. The coupled-channel 5 contact potential in the 1 ++ sector is given by [8] (see also Appendix A1) with C 0X and C 1X low-energy constants (LECs) that need to be fixed from some input. This interaction is used as kernel of the Lippmann-Schwinger equation (LSE) in the coupledchannel space in the 1 ++ sector, with M 1 and M 2 the masses of the involved mesons, μ −1 12 = M −1 1 + M −1 2 , E the center of mass (c.m.) energy of the system and p ( p ) the initial (final) relative three-momentum of the DD * pair in the c.m. frame. The used normalization is such that above threshold [E > (M 1 + M 2 )], the single channel elastic unitary condition is Im T −1 (E) = μ 12 The discussion is similar for any other J PC sector. Due to the use of contact interactions, the LSE shows an ill-defined ultraviolet (UV) behavior, and requires a regularization and renormalization procedure. We employ a standard Gaussian regulator (see, e.g. [38]) with C I X any of the LECs of Eq. (3) in the case of the X (3872), or the relevant ones for any other J PC sector. We take cutoff values = 0.5-1 GeV [7,8], where the range is chosen such that will be bigger than the wave number of the states, but at the same time it will be small enough to preserve HQSS and prevent that the theory might become sensitive to the specific details of short-distance dynamics. 6 The dependence of the results on the cutoff, when it varies within this window, provides a rough estimate of the expected size of subleading corrections. Bound states correspond to poles of the T -matrix below threshold on the real axis in the first Riemann sheet (RS) of the complex energy, while the residues at the pole give the s-wave couplings of the state to each channel (D 0D * 0 and D + D * − in the case of the X (3872) resonance 7 ). 5 Actually, positive C-parity combinations in both the neutral D 0D * 0 and the charged D +D * − channels are being considered. 6 However, as will be shown later on, the situation is more complicated in the two-body d-wave hadronic decays. 7 For instance, in the case of the X (3872), we have The LECs C 0X and C 1X can in principle be determined [8] from M X (3872) = (3871.69 ± 0.17) MeV (mass average quoted by the PDG [5]) and the isospin violating ratio of the decay amplitudes for the X (3872) → J/ψπ π and X (3872) → J/ψπππ, R X (3872) = 0.26 ± 0.07 [39]. The ratio is given by (see Eq. (80) of Ref. [8]) whereˆ n,c give the average of the neutral and charged wave function components in the vicinity of the origin, and are related to the LECs introduced in Eq. (3) by [37] where the neutral and charged loop functions, where T i j are the matrix elements of the T -matrix solution of the UV regularized LSE. Fig. 1 X 2 binding energy histograms obtained from the interaction of Eq. (9) using LECs distributions determined from the X (3872) resonance inputs (blue) or using LECs distributions additionally modified to account for the HQSS systematic error (red). Left and right plots correspond to UV cutoffs of 0.5 and 1 GeV, respectively. MC sample distributions are first generated and are used to evaluate the X 2 mass. X (3872) mass trials above threshold are rejected. To evaluate the red shaded histogram, and to account for the HQSS 20 % uncertainty in the 2 ++ interaction, each of the members of any MC sample (C 0X + C 1X , C 0X − C 1X ) pair is multiplied by independent N (μ = 1, σ = 0.2) Gaussian distributed random quantities r ± mass, identical to that in the X (3872) sector (1 ++ ) and given by Eq. (3) [7,8]. Thus, in the 2 ++ sector, the potential in the (D * 0D * 0 ), (D * + D * − ) coupled-channel space reads (see also Appendix A1) with the same structure and involving the same LECs that in the X (3872) channel. Besides, in the above equation m c ∼ 1.5 GeV is the charm quark mass and q ∼ QCD , a scale related to the light degrees of freedom. Taking QCD ∼ 300 MeV [5], corrections of the order of 20 % to the interaction predicted by HQSS cannot be discarded. 8 Nevertheless, it seems natural to expect a 2 ++ D * D * loosely bound state (X 2 ), the HQSS partner of the X (3872), and located in the vicinity of the D * 0D * 0 threshold (∼ 4014 MeV) [7,8,15]. This is illustrated in the X 2 binding energy distributions depicted in Fig. 1. Neglecting the O(q/m c ) corrections to the LECs, and using those obtained from the X (3872) resonance, we find a clear signal (blue histograms) of a weakly bound state. However, the case is less robust when the latter corrections are taken into account. Thus, because of the additional 20 % HQSS uncertainty, the area below the red shaded = 0.5 GeV ( = 1 GeV) histogram is only 0.77 (0.68). This means that approximately a 23 % (32 %) of X (3872) events [(M X (3872) , R X (3872) ) MC samples] do not produce a X 2 pole in the first RS, since the strength of the resulting interaction in the 2 ++ sector would not be attractive enough 8 The two type of uncertainties (HQSS subleading corrections and errors inherited from the inputs) affecting the determination of the LECs in the 2 ++ sector are statistically uncorrelated, and should be accordingly added up. To that end, we use MC techniques as explained in the caption of Fig. 1. to bind the state, though a virtual state in the second RS will be generated instead. Given the existence of the X (3872) as a DD * molecule, if the X 2 resonance exists, we would expect its mass (binding energy) to lie most likely in the interval [2m D * 0 , 4006 MeV] ([0, 8] MeV), as shown in Fig. 1. Note that the discussion in Ref. [15] was simpler, because there we worked in the isospin symmetric limit and used the averaged masses of the heavy mesons, which are larger than those of the physical D 0 and D * 0 mesons.
For later use, we also need the couplings of the X 2 to its neutral (D * 0D * 0 ) and charged (D * + D * − ) components, g X 2 0 and g X 2 c , respectively. They turn out to be slightly different because the X 2 resonance is an admixture of isospin 0 and 1, since its binding energy is much smaller than the energy difference between the two thresholds [8,37]. Considering the HQSS uncertainties, we find Owing to the heavy flavor symmetry, the LO 2 ++ B * B * interaction is given by Eq. (9) as well, and thus we should also expect a 2 ++ B * B * bound state (X b2 ), the HQSFS partner of the X (3872), located close to the B * B * threshold (∼10650 MeV) [15]. The X b2 binding energy distributions are shown in Fig. 2 Fig. 1 but in the bottom sector. To better appreciate the distribution details, the = 0.5 (1) GeV red histogram, which includes the 20 % HQSS error, has been multiplied by a factor of 5 (10) 5325.2 MeV. Note that, according to the PDG [5], MeV, from where we might expect isospin breaking effects for the B * mesons to be significantly smaller than in the charm sector. In Ref. [15] we found that the binding energy of the X b2 state is significantly larger than that of its counterpart in the charm sector (X 2 ), around a few tens of MeV. Thus we do not expect any significant isospin breaking effects and the X b2 resonance would be a pure isoscalar (I = 0) state.
This bound state, being isoscalar, equally couples to the neutral and charged components and, therefore: g X b2 Our predictions in Eq. (12), both for the mass and the B * B * coupling of the resonance show some dependence on the UV cutoff, which is to some extent diminished when HQSS uncertainties are taken into account. Nevertheless, this dependence might hint to non-negligible subleading corrections (among others, pion-exchange and coupled-channel effects [7], which can be larger here than in the charm sector due to the larger binding energy and larger meson masses). We will compute the decay widths for both UV regulators, and the spread of results will account for this source of uncertainty.

The hadronic X 2 and X b2 decays
The quantum numbers, J PC = 2 ++ , of these resonances constrain their possible decay channels. In this work, for hadronic decays we only consider the decays into two heavy hadrons: X 2 → DD and X 2 → DD * (D * D ), and the analogous processes X b2 → BB and X b2 → BB * (B * B ). We expect that these d-wave decay modes should largely saturate the widths of these states. Because the DD couples in a d-wave to the 2 ++ system, its contribution to the mass renormalization of the X 2 is of higher order (see Appendix B). We thus did not include the DD as a coupled channel in the T -matrix, but treat it perturbatively. This means that the transitions are mediated by the exchange of a pion. The relevant π P ( * ) P ( * ) vertices are taken from the LO Lagrangian of heavy meson chiral perturbation theory [40][41][42][43] (see Appendix A2). At LO, besides the pion decay constant, f π = 92.2 MeV, there appears only one additional D * Dπ coupling (g). We take g = 0.570 ± 0.006 as inferred from the new value of = (83.4 ± 1.8) keV for the D * + decay width quoted by the PDG [5]. This is mostly determined by the recent BABAR Collaboration measurement [44] of this width, which is approximately a factor of 12 times more precise than the previous value, = (96 ± 4 ± 22) keV by the CLEO Collaboration [45]. Thus, we end up with an uncertainty of the order 1 % for g. Though the hadronic X 2 and X b2 widths evaluated in this section will be proportional to g 4 , this source of the error (∼4 %) will be much smaller than the others and it will be ignored in the following.

X 2 (4013) → DD
We will first consider the X 2 (4013) → D + D − (D 0D0 ) decay, which proceeds through the Feynman diagrams Fig. 3 Feynman diagrams for the X 2 (4013) → D + D − decay. Diagrams for the X 2 (4013) → D 0D0 transition are similar, with the appropriate changes of the exchanged pion charges depicted in Fig. 3. We treat charm mesons non-relativistically, and neglect p D * ,D * /m D * terms and the temporal components in the D * ,D * propagators. We obtain for the X 2 (4013) → D + D − process, in the resonance rest frame and with q and k the four-momenta of the D andD final mesons ( q = − k, where i j (λ) is the symmetric spin-2 tensor with λ denoting the polarization of the X 2 state 10 and N = 8M X 2 m 2 D * √ m D m D * 2 accounts for the normalization of the heavy meson fields 11 and some additional factors needed when the couplings g X 2 c,0 , as determined from the residues at the pole of the EFT T -matrix, are used for the X 2 D * D * vertex. For the neutral and charged pion masses, we have used the values quoted by the PDG [5] and heavy meson isospin averaged masses to compute N . Besides, I i j is a three-point loop function, the detailed evaluation of which is relegated to Appendix C1. 12 The loop is seemingly logarithmically 10 For instance, it could be given by with i 0 = (0, 0, 1), and i ±1 = ∓(1, ±i, 0)/ √ 2, and ( j 1 , j 2 , j|m 1 , m 2 , m) a Clebsch-Gordan coefficient. 11 We use a non-relativistic normalization for the heavy mesons, which differs from the traditional relativistic one by a factor √ M H . 12 In the computation of I i j , we are consistent with the former approximations, and we use non-relativistic charm meson propagators. divergent. However, since the X 2 polarization is traceless, the divergent part which comes with a Kronecker delta does not contribute. This is because the decay occurs in a dwave, thus the loop momentum is converted to external momenta, and the remaining part of the integral is convergent. Nevertheless, we will include two different form factors in the computation of the three-point loop function. One is inherited from the UV regularization/renormalization procedure sketched in Eq. (5) and employed to make the LSE T -matrix finite. In addition, we will include a second form factor to account for the large virtuality of the pion in the loop. We will discuss this at length below and in Appendix C1.
Analogously, the X 2 (4013) → D 0D0 amplitude is, The two-body decay width in the X 2 rest frame reads [5]: The sum over the X 2 polarizations can easily be done in the c.m. frame, Table 1 X 2 (4013) → DD, DD * decay widths using different UV Gaussian regulators for the D * D * X 2 form factor and with/without including a pion-exchange vertex form factor (FF) in each of the D * Dπ and D * D * π vertices in the three-point loop function. The decay widths of the X 2 (4013) →D D * modes are the same thanks to C-parity. Uncertainties are obtained from a Monte Carlo simulation using the X 2 binding energy histograms displayed in Fig. 1 (red shaded) and the g X2 0 and g X2 c couplings given in Eqs. (10) and (11). Note that the procedure takes into account 20 % HQSS uncertainties and the correlations between X 2 masses (binding energies) and g X2 0 and g X2 c couplings. Errors on the widths provide 68 % CL intervals Without pion-exchange FF With pion-exchange FF As discussed in Appendix C1, the three-point loop function has a tensor structure of the type The I 1 term carries the UV divergence, which, however, does not contribute to the width, because it vanishes after the contraction with the traceless spin-2 polarization tensor, as mentioned above. Therefore, only the I 0 term is relevant, which is free from UV divergences. Moreover, the contrac- (16)] leads to a factor of 2 q 4 /3. Thus, the integration over the solid angle d (q) trivially gives 4π , and the width scales like | q | 5 as expected for a d-wave process.
Our predictions for the X 2 (4013) → D + D − , D 0D0 decays are compiled in Table 1. If we look at the first two columns of results in the table, we find widths of the order of a few MeV, with asymmetric errors that favor larger values. This is mostly due to the similar asymmetry of the uncertainties quoted for the g X 2 0 and g X 2 c couplings in Eqs. (10) and (11).
Our scheme is based on a low-energy EFT, in which the momenta should be smaller than a hard scale which serves as a momentum cutoff [see Eq. (5)]. The high-momentum modes are out of control in the low-energy EFT. Therefore in the computation of the width, we include a Gaussian regulator at the D * D * X 2 vertex, as discussed in Eq. (C6). The cutoff should be the same as the one used in generating the X 2 as it is related to the same unitary cut in the D * D * system. In Fig. 4, we display, as an example, the dependence In the left plot we see that in spite of including the Gaussian X 2 form factor, large momenta above 1 GeV provide a sizable contribution to the integral ( 14, 30 and 45 % for = 0.5, 1 GeV and ∞, respectively), which is an unwanted feature within the low-energy EFT scheme and signals a sizeable short-distance contribution. Indeed, the momentum of the exchanged pion, which peaks at around 750 MeV, is a somehow large value in the sense that the hard scale for the chiral expansion which controls the pionic coupling is χ ∼ 1 GeV. We see that below the peak, the curves for both cutoff values are very close to each other, and they are also close to the curve corresponding to the case without any regulator. This is the region where the low-energy expansion works and thus model-independent conclusions can be made. The curves start deviating from each another after the peak, that is, in the region with a pion momentum χ . Because the loop integrals are not completely dominated by momentum modes well below χ , the widths of interest will bear an appreciable systematic uncertainty. This is reflected in the fact that the widths in the second column in Table 1 are larger than those in the first column by a factor around 2. 13 On the other hand, the fact that the pion could be quite far off-shell should be reflected in the D * Dπ vertex, which should be corrected, similarly as it is done in the case of the N Nπ one. Thus, to give an estimate of the hadronic decay widths, in the spirit of the Bonn potential [46], we have included a monopole pion-exchange vertex form factor, with a hadron scale of the order of 1 GeV, in each of the D * Dπ vertices [Eq. (C7)]. Its effect on the internal pion momentum dependence of the Fig. 4. The large pion momenta contribution (| l | > 1 GeV), which is not reliable in a low-energy EFT calculation, is reduced now to 6.5, 13 and 16 % for = 0.5, 1 GeV and ∞, respec- With pion form-factor (Λπ Diagrams for the X 2 (4013) → D 0D * 0 transition are similar, with the appropriate changes of the exchanged pion charges tively. This makes also more appropriate the non-relativistic treatment of the charmed mesons adopted here. Besides, the dependence of the width on the UV Gaussian regulator is significantly softer, though the widths are further reduced by almost an order of magnitude.
We believe that the most realistic estimates are those obtained with the inclusion of the pion-exchange form factor and the spread of results compiled in Table 1 give a conservative estimate of the systematic uncertainties, beyond the mere existence of the X 2 (4013) state, as discussed in Sect. 2.2. We recall here that because of the additional 20% HQSS uncertainty, approximately a 23 % (32 %) for = 0.5 GeV ( = 1 GeV) of the X (3872) events [(M X (3872) , R X (3872) ) MC samples] do not produce the X 2 as a bound state pole, since the strength of the resulting interaction in the 2 ++ sector is not attractive enough to bind the D * D * .
Nevertheless, assuming the existence of the X 2 state, and in view of the results given in Table 1, we estimate the X 2 → DD partial width (including both the charged and the neutral channels) to be where the first error accounts for the dependence on the UV Gaussian regulator used in the D * D * X 2 vertex, while the second one is obtained from the uncertainties given in Table 1.
This latter error includes both some additional systematic (HQSS violations) and statistical (X (3872) input used to fix the properties of the X 2 resonance) uncertainties. Notice that, as discussed above, the calculation is probably already beyond the valid range of the EFT due to the large contribution of high-momentum modes. We thus have adopted a more phenomenological strategy and used the pion-exchange form factor with a cutoff of 1 GeV to make an estimate of the decay widths. The values presented in Eq. (18) refer only to the last two columns in Table 1 that include the effect of the pion-exchange form factor. Thus, and considering both sets of errors, we cover the range 0.5-2.8 MeV, which accounts for all values including errors quoted in the table.
Here, we will study the D + D * − , and D 0D * 0 channels, which proceeds through the Feynman diagrams depicted in Fig. 5. This is also a d-wave decay so that both angular momentum and parity are conserved. The decay widths are expected to be comparable to those found for the X 2 (4013) → DD decays, despite the phase space is considerably more reduced. This extra enhancement is caused by the extra multiplicity due to the spin of the final D * (D * ) meson.
As commented before, we treat charm mesons nonrelativistically and obtain the decay amplitude for the X 2 (4013) → D + D * − process as in the resonance rest frame. Here, q is the four-momenta of the D + meson, and n (λ * ) is the polarization vector of the final D * − meson with helicity λ * , N * = The two-body decay width in the X 2 rest frame in this case reads [5] The sum over theD * and X 2 polarizations can easily be done, and we get The above tensor structure should be contracted with I im ( q ) I rl ( q ). We see that the sum over polarizations of Eq. (22) is orthogonal to δ im and δ rl , which guarantees also here that only the UV finite I 0 term of the three-point loop function is relevant. The contraction leads to a q 4 factor, 14 and thus the width scales like | q | 5 , as expected for a d-wave decay.
Results for the X 2 → DD * decay widths are also compiled in Table 1. We only show predictions for the 14 In the DD mode, studied in Sect. 3.1.1 a factor of 2 q 4 /3 is obtained instead. Thus, neglecting the D − D * mass difference, because the loop integrals are the same, we would find This extra factor 3/2 due to the spin-1 polarization vector produces an enhancement of the DD * decay mode with respect to the DD one, which partially compensates the smaller available phase space. X 2 (4013) → D + D * − (D 0D * 0 ) decays, because the X 2 being an even C-parity state, its decay modes into the chargeconjugated final states have the same decay widths. In what respects to the effect of the form factors, the discussion runs in parallel to that in the Sect. 3.1.1, though the effect of the pion-exchange form factor is significantly smaller here (a factor 4 or 5 at most). As expected, the widths are comparable to those found for the X 2 → DD decays. Finally, we estimate the partial X 2 → DD * (D * D ) width (including both the charged and the neutral channels as well as the charge-conjugated modes) to be where the errors have been estimated as in Eq. (18). The above result, together with that obtained previously for the DD channel, leads to a total X 2 width of the order of 2-8 MeV, assuming its existence.

Charm decays: further analysis of uncertainties
The uncertainties on the results compiled in Table 1 account both for reasonable estimates of HQSS corrections, as well as for the statistical errors on the inputs used to fix the LEC's that determine the properties (mass and D * D * couplings) of the X 2 resonance. Moreover, we are using an EFT to describe these decays, which means that there is an intrinsic uncertainty that can be determined systematically. We indicated the size of this error in Eqs. (18) and (24). For obtaining the EFT uncertainty we combined the predictions obtained for two different UV Gaussian cutoffs ( = 0.5 and 1 GeV) in the D * D * X 2 vertex, and used the spread of values to quantify this error. The idea is that the residual dependence of the results on the cutoff should provide insight into the size of subleading corrections. Now, we try to further test the robustness of the systematic errors quoted in Eqs. (18) and (24). To that end, we have examined: • The dependence of our results on the functional form of the UV regulators, both in the D * D * X 2 and D * D ( * ) π vertices.
-We have a contact theory with a Gaussian regulator and a cutoff between 0.5 and 1 GeV. This theory, though very simple, is the LO of an EFT expansion describing the low-energy dynamics of heavy hadrons (see Ref. [36] for details). Within the EFT we can expand observable quantities as a power series of the type 2 X 2 (4013) → DD, DD * decay widths (in MeV units) using different UV regularization schemes for the D * D * X 2 vertex and pionexchange form factors. FFX 2 -G and FFX 2 -L stand for the results obtained with Gaussian (Eq. (5)) and Lorentzian (Eq. (26)) regularized local interactions, respectively. On the other hand, the widths in the columns FFπ 1 and FFπ 2 were obtained inserting the pion form factor of Eq. (C7) and F( l 2 , π ) = e − l 2 / 2 π in each of the two D * D ( * ) π vertices, respectively. In the latter case, we take π = 1.2 GeV, as determined in the QCDSR calculation of Ref. [49] for the D * Dπ coupling.
In the DD * decay mode, we also show results (in brackets) obtained when a larger cutoff, π = 1.85 GeV, is used in the D * D * π vertex, as obtained in the QCDSR study carried out in Ref. [50] for this coupling. The results presented in Table 1 correspond to the [FFX 2 -G & FFπ 1] columns and for the rest of the choices we only provide central values where p is the momenta of the hadrons and M the scale at which hadrons stop behaving as pointlike particles (about the ρ mass). A LO calculation only keeps the first term in the series above. Hence it should have a relative error of order ( p/ M ). We stress that this is expected to be so regardless of the exact regulator employed (gaussian, lorentzian, etc.), provided that the cutoff is at least of the order of M . The reason for that is that the calculations we show are renormalizable: once the counter-terms are fixed, 15 they only contain negative powers of the cutoff when we expand on the large cutoff limit. Hence, the uncertainty in the calculation is of order ( p/ ). By taking of the order of M , observables are guaranteed to have a cutoff uncertainty of the order ( p/ M ), equivalent to the EFT uncertainty. There are several methods for making a more precise estimation of the EFT error: the one we use in Table 1 is to vary the cutoff around values of the order of M (hence the choice of the 0.5-1 GeV cut-off window). EFT predictions for two different cutoffs differ by a factor of ( p/ M ) and that is why we chose this particular way of assessing the errors. Alternatively, one could use different regulators or form factors to assess the size of this error. This idea also gives a cross-check of the previous error estimates based on varying the cut-off. We have employed a different regulator to check that our former estimation of the EFT errors is correct and to show that the particular regulator employed is not important. We closely follow here the discussion of 15 See for instance the discussion of Eqs. (23) and (24) of Ref. [37].
Sect. VII of Ref. [37] and consider a Lorentz form for the regulator with two values of the cutoff, namely = 0.4 and 0.8 GeV, which were obtained by multiplying the Gaussian cutoffs by a factor of √ 2/π [37]. The resulting 16 widths are presented in Table 2 and turn out to be rather insensitive to the form of the regulator (this is to be understood within the limits of the expected EFT uncertainty). -Next we have studied the dependence of the widths on the pion form factor that accounts for the offshellness of the pion in the mechanisms depicted in Figs. 3 and 5. To that end we used the results from the QCD sum rule (QCDSR) calculations performed in Refs. [49,50]. The first of these two works considers the D * Dπ vertex: it was found the form factor is harder if the off-shell meson is heavy, implying that the size of the vertex depends on the exchanged meson. This means that a heavy meson will see the vertex as point-like, whereas the pion will see its extension. The authors of Ref. [49] find an on-shell value g = 0.48 ± 0.05 (note the different definition used in this reference), around 1-2 sigma below the value of 0.57 used in this work. In addition, they adjust their results for off-shell pions to a Gaussian form e l 2 / 2 π , with l μ the pion four-momentum, and find π = 1.2 GeV. This form factor, 17 in the region of interest for this work, turns out to be in good agreement with that used to obtain the results presented in Table 1. This can be seen in the new results showed in Table 2 and obtained with this new pion off-shell form factor (FFπ 2). Deviations from our previous estimates of the X 2 (4013) → DD, DD * decay widths are both much smaller than the (HQSS & exp) uncertainties quoted in Table 1, and well comprised within the systematic uncertainty generated when the D * D * X 2 cutoff varies in the 0.5-1 GeV window.
In the X 2 → DD * decay, there also appears the D * D * π coupling in one of the vertices; see Fig. 5.
The off-shell behavior of this vertex might differ from that of the D * Dπ one. This coupling was studied using QCDSR in [50] where, translating the definition used therein to the one used here, it was found an on-shell value g = 0.56 ± 0.07 in good agreement with the HQSS expectations. The off-shell behavior was described by a Gaussian, as in the case of the D * Dπ vertex, but with a significantly larger cutoff, π = 1.85 GeV. This significant difference is somehow surprising from the HQSS point of view, and we should note that the QCDSR actual calculations in [50] were carried out for significantly larger values of l 2 > 4 GeV 2 than in the case of the D * Dπ coupling analyzed in Ref. [49]. Nevertheless, we used this softer dependence for the D * D * π vertex and re-computed the DD * widths. Results are shown in brackets in Table 2. Changes are now larger, and in general are of the order of 50 %, though they could be still accommodated within the HQSS and EFT uncertainties already considered in our original calculations. The large momenta of the external mesons, which can even exceed 0.5 GeV, make it possible that the short-distance details of the decay mechanisms could be relevant. This is the weakest point in our scheme. The reason is that the EFT uncertainty is expected to be 0.5 GeV/ M 1/2, as the variations of the DD * widths in Table 2 seem to confirm, and the calculation is really on the limit of validity of this kind of description and should only be considered as a reasonable estimate. 17 To use the form factor of Ref. [49] in the computation of the widths, we have approximated the pion four-momentum squared by − l 2 , which is sufficiently accurate because the dominant part of the integration comes from regions where the two virtual D * andD * mesons are almost on-shell. In this region, the energy of both heavy-light vector mesons is approximately M X2 /2, which coincides with that of the external heavy mesons, and hence l 0 is much smaller than | l |.
• The contribution of decay mechanisms driven by the exchange of shorter-range mesons (heavier) than the pion.
Since the momenta of the external charmed mesons can exceed 0.5 GeV, one might think that shorter-range contributions such as the ρ and ω exchanges could be sizable, and even comparable to those of the diagrams depicted in Figs. 3 and 5 for the exchange of a pion. We will focus on the X 2 → DD decay mode, where the momenta of the external mesons is the largest and we will begin by studying the effects of the exchange of a ρ meson. If we use the phenomenological D * Dρ Lagrangian given in Eq. (3e) of Ref. [51], we find that the amplitudes of this new contribution can be obtained from those driven by pion exchange, and given in Eqs. (13) and (14), by replacing m π 0 and m π + by m ρ 0 and m ρ + in the loop integrals and where we have neglected | q | 2 /m 2 D ( * ) terms, with q the c.m. three-momentum of the external D andD mesons. The coupling constant g D * Dρ has been computed in various schemes [51][52][53][54][55] (ordinary and light cone QCDSR, vector dominance model, SU(4), etc.) and it varies in the range [2.8 ± 0.1, 4.3 ± 0.9] GeV −1 (see Table 5 of Ref. [55]). Taking g D * Dρ ∼ 5 GeV −1 , in the highest part of the interval of calculated values, we would have g 2 D * Dρ /(g 2 / f 2 π ) ∼ 2/3, while a direct calculation of the loop integrals shows that those calculated using the ρ mass instead of the pion mass are around a factor of 2 smaller. Altogether, this indicates that the absolute values of the ρ-exchange amplitudes are about a factor of 3 smaller than those driven by the pion exchange. If one uses g D * Dρ ∼ 3 GeV −1 , now in the lowest part of the interval of values, the ρ contribution will be, at the level of amplitudes, around eight times smaller than those considered in the present work. In any case, these effects are smaller than the HQSS and EFT uncertainties quoted in Table 1, and therefore it seems justified to neglect them. On the other hand, the ω-exchange contributions are even smaller, around a factor of 3, because this meson is neutral, and it cannot generate mechanisms where a light charged meson is exchanged. To estimate the size of the ρ exchange contribution in the X 2 → DD * decay mode, in addition to the D * Dρ vertex, we have used the phenomenological D * D * ρ Lagrangian given in Eq. (1f) of Ref. [51]. The amplitudes of this new mechanism can be obtained from those driven by pion exchange, and given in Eqs. (19) and (20), by replacing m π 0 and m π + by m ρ 0 and m ρ + , and P 2 (x) → (P 2 (x) + x| q |/l) in the computation of the loop integral I 0 (see Eq. (C5) of Appendix C1) and using the appropriate coupling constants where we have neglected as before the | q | 2 /m 2 D ( * ) terms. (At this point, the phase in Eq. (28) is arbitrary, because the phase convention used in [51] and in this work for the heavy meson fields is different.) The g D * D * ρ coupling is set to 2.52 in [51], while it is estimated to be 4.7 ± 0.2 in [53,56] (QCDSR). On the other hand, there exist some partial cancellations in the loop integral in this case, and it turns out to be at least a factor of 5 smaller than that calculated for the pion driven contribution. Thus, and having in mind the previous discussion for the expected range of values of the g D * Dρ coupling, we conclude that the ρ-exchange contribution can be safely neglected in the X 2 → DD * decay mode.
In view of the results discussed in this subsection, we should acknowledge that as a result of the contribution from highly virtual pions, which is certainly in the limit of applicability of the low-energy EFT, the hadronic decay widths bear a large systematic uncertainty. Nevertheless, we have given arguments to be reasonably convinced that the results quoted in Eqs. (18) and (24) provide sensible estimates for the X 2 (4013) → DD, DD * widths.
In the next subsection, we will study these hadronic decays in the bottom sector. There, the considerations are parallel to those discussed here for the charm sector.

Bottom decays
Thanks to the heavy flavor symmetry, the results of the previous subsection can be trivially extended to the bottom sector. There, we have a robust prediction, even when HQSS uncertainties (20 %) are taken into account, for the X b2 resonance (see Fig. 2). Moreover, all sort of non-relativistic approximations adopted in the current scheme are now more suited, since the range of variation of the internal pion momentum in the loops is similar to that shown in Fig. 4 for the charm sector.
On the other hand, as discussed in Sect. 2.3, we do not expect any significant isospin breaking effects and the X b2 resonance would be a pure I = 0 state, with equal coupling to its neutral and charged components. For simplicity, we will also neglect the tiny difference between B 0 and B ± masses, and we will use a common mass m B = (m B 0 + m B ± )/2 = 5279. 42 MeV. Yet, for the pion mass that appears in the loop integral, we take the isospin averaged value m π = (2m π ± + m π 0 )/3. Note that the relevant internal pion momentum is around 750 MeV. With all these approximations, we find for any charge channel (B + B − , B 0B0 , B + B * − , B 0B * 0 ) or charge-conjugation mode (B * B ). Our results for these decay widths are presented in Table 3. We notice in passing that following heavy flavor symmetry we use the same value of g = 0.570 ± 0.006 in the charm and bottom decays. It agrees very well with a recent lattice calculation with relativistic bottom quarks which gives g b = 0.569 ± 0.076 [47] (we have added the systematic and statistic errors in quadrature). Yet, lattice calculations with static heavy quarks tend to give smaller values; see Ref. [48] and references therein. For instance, the ALPHA Collaboration presented a precise computation with the result of g ∞ = 0.492 ± 0.029 [48]. Thus we expect that the decay widths of Table 3 slightly overestimate the real ones.
For the BB mode we find a pronounced dependence on the Gaussian cutoff employed in the dynamical generation of the resonance. This is inherited from the strong dependence of the X b2 mass on this UV cutoff, as discussed in Eq. (12), which affects the available phase space for the decay. With all these shortcomings, we expect a partial width in the 1-10 MeV range, when both charge modes are considered.
In the BB * decay mode, the impact of the Gaussian regulator is even larger, because it turns out that for = 1 GeV, the central value of the resonance mass M X b2 = 10594 +22 −26 MeV is located below the threshold m B + m B * ∼ 10604 MeV. Thus, in that case, the decay will be forbidden. For = 0.5 GeV, we estimate a width also in the 4-12 MeV range, when the four possible decay modes (B + B * − , B 0B * 0 , B − B * + , B 0 B * 0 ) are considered.

The X 2 and X b2 radiative decays
In this section, we will study the X 2 → DD * γ and X b2 → B B * γ decays. The interaction of the photon with the s-wave heavy mesons contains two contributions that correspond to the magnetic couplings to the light and heavy quarks [57] (see also Appendix A3). Both terms are needed to understand the observed electromagnetic branching fractions of the D * + and D * 0 because a cancellation between the two terms accounts Table 3 X b2 → BB, BB * decay widths (here BB ( * ) refers to any of the neutral or charged modes, but it is not the sum of both). The results are given for different treatments of the three-point loop function. The errors have been obtained using the Monte Carlo analysis explained in Table 1, but now considering the X b2 mass histograms displayed in Fig. 2 and the coupling given in Eq. (12). The decay width of the X b2 →B B * mode is the same because of charge-conjugation symmetry for the very small width of the D * + relative to the D * 0 [5]. Actually, one finds [57,58]: where E γ is the photon energy, m c the charm quark mass and α ∼ 1/137.036 the fine-structure constant. In the nonrelativistic constituent quark model β = 1/m q ∼ 1/330 MeV −1 , where m q is the light constituent quark mass. Heavy meson chiral perturbation theory allows one to improve upon this approximation by including corrections from loops with light Goldstone bosons, which give O( √ m q ) corrections [57].
Using isospin symmetry to relate (D * 0 → D 0 π 0 ) and (D * + → D 0 π + ), correcting by the slightly different available p-wave phase space in each of the two decays, and taking into account the experimental D * 0 and D * + widths and radiative branching fractions quoted by the PDG [5], we find (D * 0 → D 0 γ ) = (22.7 ± 2.6) keV, These values differ from those used in Ref. [58] because of the recent accurate BABAR measurement of the D * + decay width, mentioned in Sect. 3, which is around 10 % smaller than the previous CLEO one used in Ref. [58]. Fixing the charm quark mass to m c = 1.5 GeV, we fit the parameter β to the above experimental values and find β −1 = (293 ± 11) MeV.
In what follows, we will study decays of the type X 2 → PP * γ , P and P * being pseudoscalar and vector heavy-light mesons, respectively. Let us define p Since, as we will see, the Feynman amplitudes depend only on the invariant masses m 2 12 and m 2 23 of the final γ P and PP * pairs, respectively, we can use the standard form for the Dalitz plot [5] with |T | 2 the absolute value squared of the decay amplitude with the initial and final polarizations being averaged and summed, respectively. Thus, we readily obtain where for a given value of m 2 23 , the range of m 2 12 is determined by its values when p P is parallel or anti-parallel to p γ [5]: with E * P = (m 2 23 + m 2 P − m 2 P * )/2m 23 and E * γ = (M 2 X − m 2 23 )/2m 23 the energies of the P meson and photon in the m 23 c.m. frame, respectively, and p * P,γ the moduli of the corresponding three-momenta.
Because of parity conservation, this is a p-wave decay and hence the photon momentum appears always in the amplitudes. In the X 2 rest frame, it is given by We will first consider the X 2 (4013) → D 0D * 0 γ decay, which proceeds according to the Feynman diagrams depicted in Fig. 6. This decay can take place directly through the radiative transition of the constituent D * 0 as shown in Fig. 6a, which is the tree-level approximation. However, there are other mechanisms driven by the DD * FSI. After emitting (a) Fig. 6 Feynman diagrams for the X 2 (4013) → D 0D * 0 γ decay. Diagrams for the D + D * − γ transition are similar the photon, the vector meson D * 0 transits into the D 0 , and it can interact with the other constituent in the X 2 as shown in Fig. 6b. There is a third (c) mechanism in which the photon is emitted from theD * 0 meson, and the virtual D * 0D0 rescatter into D 0D * 0 . Finally, Fig. 6d, e present another possibility, namely the decay can also occur through the charged D * + D * − component of the X 2 resonance, and the virtual charged D + D * − and D * + D − pairs then rescatter into D 0D * 0 . Because the X 2 has a well-defined charge parity (+), the decay width into the charge-conjugated modē D 0 D * 0 γ is the same. The Feynman diagrams contributing to the D + D * − γ decay mode are similar with obvious replacements.
One could also construct other FSI mechanisms by replacing the D * Dγ vertices in Fig. 6 with the D * D * γ ones. These contributions are small and will be discussed below.

Tree level approximation
The Feynman amplitude of the mechanism depicted in Fig. 6a reads (as in the previous sections, we treat the charm mesons non-relativistically) with m 12 the invariant mass of the final γ D 0 pair. Besides, i (λ γ ) is the polarization vector of the final photon with helic-ity λ γ , p γ is its three-momentum and N γ = 8M X 2 m 2 D * √ m D m D * accounts for the normalization of the heavy meson fields and the X 2 D * 0D * 0 coupling. Finally, with p 2 12 m D * 0 (M X 2 − m D * 0 − m 12 ) and γ 2 = m D * 0 (M X 2 −2m D * 0 ) < 0. The Gaussian form factor is inherited from the D ( * )D( * ) EFT UV renormalization scheme.
We have neglected the D * 0 width in the above propagator because, since it is quite small, its inclusion only leads to small numerical variations in the decay rate which are certainly smaller than uncertainties induced by the errors in the coupling g X 2 0 and the mass of the X 2 (4013) resonance. Similarly, the use of the non-relativistic D * 0 propagator instead of m 2 12 − m 2 D * 0 + iε −1 leads also to numerically negligible differences, as compared to the HQSS corrections. The sum over theD * 0 , γ , and X 2 polarizations can easily be done, and we get The amplitude for D +D * − γ decay is readily obtained from Eq. (38) making the obvious replacements: g X 2 0 → g X 2 c , β 1 → β 2 , and (m D 0 , m D * 0 ) → (m D + , m D * + ). Performing the phase-space integration, we find at tree level (assuming the existence of the X 2 state) where the values outside and inside the parentheses are obtained with = 0.5 and 1 GeV, respectively. The errors account for the uncertainty, both in the mass of the X 2 state and in its couplings g X 2 0,c , derived from the X (3872) input (M X (3872) and the ratio R X (3872) of the decay amplitudes for the X (3872) → J/ψρ and X (3872) → J/ψω decays) and the HQSS corrections, as explained in the caption of Table 1. We have neglected additional uncertainties stemming from the error on β ( 3 %), because it is totally uncorrelated to those discussed above, and it is much smaller than those affecting for instance the g X 2 0,c couplings. The neutral mode width is much larger than the charged one, thanks to the bigger magnetic D * Dγ coupling and a larger available phase space.
In analogy with the discussion of Eqs. (20) and (21) in Ref. [33] for the X (3872) → D 0D0 π 0 decay, in the X 2 radiative processes the relative distance of the D * D * pair can be as large as allowed by the size of the X 2 (4013) resonance, since the final state is produced by the one body decay of thē D * meson instead of by a strong two-body transition. Thus, the radiative DD * γ decays might provide details on the longdistance part of the resonance wave function. For instance, the d /d| pD * 0 | [d /d| p D * − |] distribution is related to the X 2 (4013) wave function ( pD * 0 ) [ ( p D * − )] [33]. This is in sharp contrast to the DD and DD * decay modes studied in the previous sections, which turned out to be strongly sensitive to short-distance dynamics of the resonance, as revealed by the notorious dependence on the UV form factors.
However, all these considerations are affected by the DD * FSI effects to be discussed next.

DD * FSI effects
To account for the FSI effects, we include in the analysis the DD * → DD * and D * D → DD * T -matrices, which are obtained by solving the LSE [Eq. (A10)] in coupled channels with the V D ( * )D( * ) potential given in Eq. (A8). Some isospin symmetry breaking effects are taken into account because the physical masses of the neutral (DD * ) and charged (D + D * − ) channels are used in Eq. (A10). The X 2 → D 0D * 0 γ decay amplitude for the mechanisms depicted in Fig. 6b, c reads −iT (λ, λ * , λ γ ) (43) where m 23 is the invariant mass of the final D 0D * 0 pair, and the three-point loop integral function J M 1 , M 2 , M 3 , p γ is given in the Appendix C2. The integral is convergent, however, for consistency it is evaluated using the same UV regularization scheme as that employed in the D ( * )D( * ) EFT. In sharp contrast with the hadronic decays studied above, the momenta involved in these integrals are rather low.
On the other hand, we see that in the (b)+(c) contribution there appears the combination T D 0D * 0 →D 0D * 0 (m 23 ) + T D * 0D0 →D 0D * 0 (m 23 ). In the isospin limit, when the mass differences between the neutral (D 0D * 0 ) and charged (D + D * − ) channels are neglected, we will find T 00→00 = (T I =0 C=−1 + T I =1 C=−1 )/2. From Eqs. (A8) and (A10), we find the C-parity odd isospin amplitudes, 18 with G DD * G D 0D * 0 G D +D * − a common loop function. Note that the kernel of this LSE is fixed by the isoscalar (C 0Z ) and isovector (C 1Z ) C(charge conjugation) = −1 terms of V D ( * )D( * ) . This is a trivial consequence of the conservation of this symmetry, taking into account that the X 2 and the photon are even and odd C-parity states, respectively. The (d) and (e) FSI contributions of Fig. 6 are similar, with obvious replacements. We find −iT (λ, λ * , λ γ ) where (47) and in the isospin limit, we would have T +−→00 = with and T +−→+− (m 23 ) = T 00→00 (m 23 ) in the isospin limit. Taking into account isospin corrections, induced by the meson mass differences, all the needed T -matrices ( T 00→00 , T +−→00 , and T +−→+− ) can be calculated by solving the coupled channel LSE, Eq. (A10), with the V D ( * )D( * ) potentials of Eq. (A8), as mentioned above. Thanks to the conservation of C-parity, the FSI corrections will depend only on C 0Z and C 1Z . Indeed, one finds with F (E) = Diag f neu (E), f ch (E) , where the Gaussian form factors are defined after Eq. (A13). The Z b (10610) observed in Ref. [59] carries electric charge, and its neutral partner was also reported by the Belle Collaboration [60]. It lies within a few MeV of the BB * threshold and it is tempting to speculate about it as a hadronic molecule. Belle also reported the discovery of a second exotic electrically charged bottomonium state [59], Z b (10650) in the vicinity of the B * B * threshold. For both the Z b (10610) and the Z b (10650) states, J P = 1 + are favored from angular analyses.
Within our scheme, we assume that the Z b (10610) resonance is an isovector BB * + B * B / √ 2 s-wave bound state with J PC = 1 +− [15]. Note that HQSS predicts the interaction of the B * B * system with I = 1, J PC = 1 + quantum numbers to be identical to that of the BB * pair in the Z b (10610) sector. Thus, HQSS naturally explains [61] the approximate degeneracy of the Z b (10610) and Z b (10650). Taking for the Z b (10610) binding energy (2.0 ± 2.0) MeV [62], we could fix a third linear combination of the LECs that appear in the LO Lagrangian of Eq. (A6) 10 −0.17 −0.30 +0.02 −0.04 fm 2 (52) for = 0.5(1.0) GeV. Errors have been obtained with a procedure similar to that described in the discussion of Eq. (8) and used in the case of the X (3872). Assuming heavy quark flavor symmetry, the above value of C 1Z could be also used in the charm sector, subject to corrections of the order of O( QCD /m c ) 20 %, which we take into account. Therefore, we predict the existence of the isovector charmonium partners of the Z b (10610) and the Z b (10650), though as virtual states in the second Riemann sheet [15], which probably correspond to the recently observed charged charmoniumlike 19 Z c (3900) and Z c (4025) states [63][64][65][66][67]. These resonances lie close to the DD * and D * D * thresholds, respectively, while J P = 1 + quantum numbers are favored from some angular analyses.
Due to the presence of the Z c (3900) close to threshold, one should expect the loop (FSI) mechanisms depicted in Fig. 6 to be important since T I =1 C=−1 must have a pole. However, the value of C 0Z = (C 0 A − C 0B ) is still unknown. It is not determined by the inputs deduced from the X (3872) and Z b (10610) states used in the present analysis. Depending on the value of C 0Z , there can be an isoscalar J PC = 1 +− DD * s-wave bound state or not. For instance, considering the case for = 0.5 GeV and taking the central value for C 0Z ∼ −2.5 fm 2 one finds a bound state pole in the DD * system bound by around 10 MeV; if C 0Z ∼ −1.5 fm 2 , there will be a DD * bound state almost at threshold; if the value of C 0Z is larger, there will be no bound state pole any more. Therefore, the information of C 0Z will be crucial in understanding the DD * system and other exotic systems related to it through heavy quark symmetries [15,35]. Conversely, as we will see, the X 2 radiative decay width could be used to extract information on the fourth LEC, C 0Z , thanks to the FSI effects.
To investigate the impact of the FSI, in Fig. 7 we show the dependence of the partial X 2 (4013) → DD * γ decay widths on C 0Z . For comparison, the tree-level results are also shown 19 The simple EFT scheme employed in this work does not allow for finite width resonances, only for virtual or bound states. The merit of that EFT is actually making a connection between the Z b (10610/10650) and Z c (3900/4025) resonances on the basis of heavy flavor symmetry, not so much predicting the exact location. In that simple theory, we treat the Z b (10610) as a stable bound state and that is the reason why the Z c are predicted to be virtual states. Had we used a more complete EFT that takes into account the finite width of the Z b (10610) and includes next-to-leading corrections, then we would have predicted the Z c as resonances and located them more accurately in the complex plane.   1.8 Fig. 7 Dependence of the X 2 (4013) → D 0D * 0 γ and the X 2 (4013) → D + D * − γ partial decay widths on the low-energy constant C 0Z . The UV cutoff is set to = 0.5 GeV (1 GeV) in the top (bottom) panels.
The red error bands contain the DD * FSI effects, while the three horizontal blue lines stand for the tree-level predictions of Eqs. (41) and (42). Besides the uncertainties on the mass and the couplings of the X 2 resonance, the errors on C 1Z quoted in Eq. (52), together with the expected 20 % heavy quark flavor symmetry corrections when it is used in the charm sector, are also accounted in the 68 % CL bands displayed in the panels. The red dash-dotted (full calculation) and solid blue (treelevel) lines stand for the results obtained with the central values of the parameters in the same plots. As expected, for the decay into the D 0D * 0 γ channel, the FSI effects turn out to be important, and for some values of C 0Z , they dominate the decay width. The maximum effects of the FSI mechanisms approximately occur for values of C 0Z , which give rise to an isoscalar 1 +− DD * bound or virtual state close to threshold. One can see an apparent deviation from the tree-level results in this region. When C 0Z takes smaller values, the binding energy of the bound state increases and moves apart from threshold; when C 0Z takes larger values, the pole moves deeper into a non-physical RS and becomes a virtual state further from the threshold. In both situations, the FSI corrections turn out to be less important.
On the other hand, the FSI corrections are always important in the D + D * − γ channel. This is because the tree-level amplitude involves only the D * ± D ± γ magnetic coupling, while FSI brings in the neutral magnetic coupling, which is much larger than the former one. This is also the reason, besides phase space, why the tree-level width is much larger in the neutral mode than in the charged one, as commented above.
Notice that in the above calculations, we did not include the contribution from the coupled-channel FSI D * D * → DD * which can come from replacing the D * Dγ vertices in Fig. 6 by the D * D * γ ones. 20 We have checked that this is a good approximation when the re-summation of the charmed meson scattering is switched off. This is partly because the loop integral defined in Eq. (C8) takes a much larger absolute value for the considered diagrams than for those with the D * D * γ vertices in most regions of the phase space. The only exception is when the photon has a very low energy so that the D * D * are almost on-shell. However, after integrating over the phase space, this region provides a small contribution to the partial decay width because of the p 3 γ suppression due to relative p-wave between the photon and the 1 +− DD * system. The re-summation would introduce a 20 The electric part of the D * D * γ vertex does not contribute to the FSI X 2 decay width amplitude when the quantum numbers of the final D * D * pair are 1 +− , with the two heavy mesons in relative s-wave. Thus, we are left with the contribution from the magnetic coupling, as in the D * Dγ case.

BE (MeV)
, It contributes mainly to the region close to the pole position of the Z c (4025). Again, this corresponds to the region with the very low-energy photon and thus it is suppressed due to phase space. Therefore, we expect that the neglected contribution discussed here has little impact, and its numerical effects should be safely covered by the sizable HQSS uncertainty exhibited in the results.
To better understand the dependence of the X 2 (4013) → DD * γ decay widths on C 0Z in the bottom panel of Fig. 8 we show the pole positions of the odd C-parity DD * → DD * T -matrix as functions of this LEC and for an UV cutoff of 0.5 GeV. There are two coupled channels, the neutral one which has the lowest threshold, and the charged channel. As a consequence, there are three relevant RS's (among four). The first of them [(++)] is the physical one, while the two non-physical RS's are reached by changing the sign of the imaginary part of the momentum inside of the loop functions G in Eq. (51) of the first channel [(−+)] or the momenta of both channels [(−−)]. Each of these non-physical RS's are continuously connected with the physical one on the real axis above the relevant thresholds. The solid red and dashed blue curves stand for poles in the (++) and (−+) RS's, respectively. For sufficiently large negative values of C 0Z , there is a bound state (a pole below threshold and located in the physical RS), which becomes less bound when |C 0Z | decreases. For a value of this LEC around −1.5 fm 2 , marked in Fig. 8 with a vertical black line, this state reaches the threshold. When C 0Z is further increased, the pole jumps into the (−+) RS, becoming thus a virtual state (a pole below threshold and located in an non-physical RS) and it moves away from the threshold. In the medium panel we show the central values of the X 2 (4013) radiative widths for both decay modes as a function of C 0Z , which were already presented in Fig. 7. As can be seen, it is around this critical value C 0Z = −1.5 fm 2 , when the FSI effects are larger for both decays, due to the vicinity of the pole to the threshold. This happens regardless of whether it is a bound or a virtual state, since the presence of the pole in both situations greatly enhances the odd C-parity DD * → DD * T -matrix near threshold. This can be appreciated in the top panel of Fig. 8, where the dependence of the T 00→00 and T +−→+− scattering lengths on C 0Z is shown.
Thus we have understood why in the region of values of C 0Z around −1.5 fm 2 , FSI corrections strongly affect the D 0D * 0 γ decay width: this channel has the lowest threshold and the bound or virtual state is located on or nearby it. For the charged decay mode, the width exhibits a maximum for values of C 0Z also in this region, followed by a clear minimum placed now in the vicinity of C 0Z ∼ C 1Z . Notice that when C 0Z = C 1Z , T 00→+− vanishes [see Eq. (51)] and therefore the contribution due to the neutral mesons, driven by the largest magnetic coupling (β 1 ), in the FSI loops disappears. The exact position of the minimum is modulated by the further interference between the tree-level and the FSI charged loops, which are comparable.
In the bottom panel, we also observe a virtual state pole, the position of which is rather insensitive 21 to C 0Z . It orig-inates with the interaction in the I = 1 sector, C 1Z , and it should be related to the Z c (3900) exotic charmonium-like state reported by the BESIII and Belle Collaborations. Within our LO EFT scheme, we do not find a DD * bound state but, instead, a pole located near threshold in a non-physical RS [15].
The expressions found in the charm sector can be readily used here, having in mind the following correspondence: D 0 ↔ B − , D + ↔B 0 ,D 0 ↔B + , and D − ↔ B 0 . Since the heavy quark flavor symmetry ensures that g and β 1,2 are the same in the b and c systems (up to corrections of the order QCD /m c ), the expressions of Eqs. (31) and (32) can be used to predict the widths for the B * radiative decays [57], where we have taken the value m b = 4.8 GeV for the bottom quark mass. As in the study of its hadronic decays, we assume the X b2 be a pure I = 0 state, with equal coupling to its neutral and charged components, g X b2 The isospin breaking effects for the B * mesons are expected to be small and the tiny difference between the B 0 and B ± masses can safely be neglected as well. In this limit we find at tree level where the values have been obtained with = 0.5 GeV. We recall that for = 1 GeV, the central value of the resonance mass M X b2 is located below the threshold (m B + m B * ) ∼ 10604 MeV and the decay is forbidden. The errors reflect the uncertainty in the inputs from the X (3872) and the HQSS breaking corrections, as outlined in the caption of Table 3; they are quite large and are dominated by those quoted for g X b2 in Eq. (12). The widths are of the order of a few eV, significantly smaller than those of the radiative decays of the B * 0 and B * − meson because of the quite reduced phase space available (∼20 MeV) for this p-wave decay. They are also orders of magnitude smaller than (X 2 → D 0D * 0 γ ) as a result of (D * 0 → D 0 γ ) (B * → Bγ ). The amplitude for the FSI mechanisms is readily evaluated and we find where the +− (−+) combination stands for the B − B * + γ (B 0 B * 0 γ ) decay mode and m 23 is now the invariant mass of theB B * pair. The C-parity odd isospin amplitudes are obtained by solving Eq. (45) using the bottom sector loop function GB B * . FSI corrections turn out to be important, as can be appreciated in Fig. 9. This is because we are generating in the T I =1 C=−1 amplitude a bound state [Z b (10610)], almost at threshold (binding energy (2.0 ± 2.0) MeV [62]), that enhances the loop mechanisms, as we discussed in the charm sector. If we pay attention for instance to the charged B − B * + γ mode, we could appreciate a distinctive feature: there appears a destructive interference pattern between the tree-level and the FSI amplitudes. Thanks to our MC procedure where correlations are consistently propagated, we also observe a reduction of the size in the uncertainties. Besides the uncertainties on the mass and the couplings of the X b2 resonance, the errors on  Fig. 9 Dependence of the X b2 → B − B * + γ and X b2 →B 0 B * 0 γ partial decay widths on the low-energy constant C 0Z . The red error bands contain theB B * FSI effects, while the three horizontal blue lines stand for the tree-level predictions of Eq. (56) C 1Z quoted in Eq. (52) are also accounted for in the 68 % CL bands displayed in the panels. Actually, these latter uncertainties should have also an important impact on the total CL bands. This is because variations of C 1Z allow for situations where the pole is located precisely at threshold (zero binding energy) or bound by about 4 MeV. In the first case the FSI contribution should be larger than that obtained with the central value of C 1Z , which correspond to a binding energy of 2 MeV. These big 68 % CL bands makes hard to disentangle any further dependence on C 0Z , which in this case turns out to be quite mild.

Conclusions
In this work we have studied the hadronic and radiative decays of a molecular P * P * state with quantum numbers J PC = 2 ++ in the charm (X 2 ) and bottom (X b2 ) sectors using an EFT approach. We have considered the X (3872) resonance as a J PC = 1 ++ DD * hadronic molecule. The X 2 and the X b2 states will be HQSFS partners of the X (3872) with masses and couplings to the P * P * heavy meson pair determined by the properties of the X (3872) resonance. The hadronic d-wave X 2 → DD and X 2 → DD * twobody decays are driven by one-pion exchange. We observed that as a result of the contribution from highly virtual pions, which is out of control in the low-energy EFT, these hadronic decay widths (hence the total width of the X 2 as well) bear a large systematic uncertainty. Even though the momenta involved in these decays probably lie outside the range of applicability of EFT the calculations are still valuable as a way to find reasonable estimates of these partial decay widths, which we expect to almost saturate the X 2 decay width. To this end and in analogy to the Bonn potential, we have included a monopole pion-exchange form factor, with a cutoff around 1 GeV, in each of the D * Dπ and D * D * π vertices to suppress the contribution of large momenta. We finally estimate the partial widths of both processes to be of the order of a few MeV. The analysis runs in parallel in the bottom sector with the assumption that the bare contact terms in the Lagrangian are independent of the heavy flavor. In this sector, we also find widths of the order of a few MeV.
We discussed the radiative X 2 → DD * γ and X b2 → B B * γ decays as well. The widths are small, of the order of keV's (eV's) in the charm (bottom) sectors. Furthermore, they are affected by the DD * or BB * FSI mechanisms. FSI effects are large because they are enhanced by the presence of the isovector Z c (3900) and Z b (10610) resonances located near the D 0D * 0 andB B * thresholds, respectively. In the charm sector, FSI corrections turn out also to be sensitive to the negative C-parity isoscalar DD * interaction (C 0Z ). Thus, future precise measurements of these radiative decay widths might provide valuable information on this LEC, which cannot be in principle determined from the properties of the X (3872), Z b (10610), and Z b (10650) resonances. Constraints on the latter LEC are important in order to understand the dynamics of the P ( * )P ( * ) system.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided 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. Funded by SCOAP 3 .

Appendix A: Heavy meson Lagrangians: s-wave interactions and pionic and electromagnetic decays
We collect in this appendix the Lagrangians used in this work. We use the matrix field H (Q) [H (Q) ] to describe the combined isospin doublet of pseudoscalar heavy-meson P ] (see for example [40]),  [40], this is to say: Their hermitian conjugate fields are defined bȳ and transform as [40]: The definition for H (Q) a also specifies our convention for charge conjugation, which is C P (Q) a C −1 = P (Q)a and C P * (Q) aμ C −1 = −P * (Q)a μ , and thus it follows that with c the Dirac space charge-conjugation matrix satisfying cγ μ c −1 = −γ T μ .

Appendix A1: Quadruple-heavy-meson contact interaction
At very low energies, the interaction between a heavy and anti-heavy meson can be accurately described just in terms of a contact-range potential. The LO Lagrangian respecting HQSS reads [68] with τ the Pauli matrices in isospin space, and C

(τ )
A,B light flavor independent LECs, which are also assumed to be heavy flavor independent. Note that in our normalization the heavy or anti-heavy meson fields, H (Q) or H (Q) , have dimensions of E 3/2 (see [69] for details). This is because we use a nonrelativistic normalization for the heavy mesons, which differs from the traditional relativistic one by a factor √ M H . For later use, the four LECs that appear above are rewritten C 0 A , C 0B and C 1A , C 1B , which stand for the LECs in the isospin I = 0 and I = 1 channels, respectively. The relations read The LO Lagrangian determines the contact interaction potential V = −L/4, which is then used as kernel of the two-body elastic LSE (see Eq. (4) and the related discussion).
The LECs that appear in the J PC = 1 ++ and 2 ++ sectors [Eqs. (3), (9)] turn out to be C 0X ≡ C 0 A + C 0B and C 1X ≡ C 1A + C 1B . The contact interaction in the Z b (10610) sector On the other hand, the interaction in the D 0D * 0 , D * 0D0 , 22 In the bottom sector, the corresponding basis is B − B * + , B * − B + ,B 0 B * 0 ,B * 0 B 0 .
with C 0Z = C 0 A − C 0B and the orthogonal matrix A given by (A9) Equation (A8) trivially follows from the fact that the L 4H interaction of Eq. (A6) is diagonal in the isospin basis and the charge conjugation is well defined. 23 The interaction given in Eq. (A8) can be used as the kernel of an UV finite LSE to obtain the T -matrix that we use to account for the FSI in the radiative decays studied in Sect. 4, with the two particle regularized matrix propagator defined as where trivially G D 0D * 0 = G D * 0D0 and G D + D * − = G D * + D − . In addition, the on-shell UV Gaussian form factor matrix reads

Appendix A2: P ( * )P ( * ) π interactions
The relevant term in the LO Lagrangian of the heavy meson chiral perturbation theory [40][41][42][43] that provides the D * Dπ and D * D * π couplings is 23 For instance in the charm sector, the C-parity states are [DD * ] 1,2 = DD * ±D * D √ 2 (1 ↔ +, 2 ↔ −). In our convention, the C-parity of these states is independent of the isospin and it is equal to ∓1.
with φ a relativistic field that describes the pion, 24 g is the heavy flavor independent P P * π coupling and f π = 92.2 MeV the pion decay constant. Note that in our normalization the pion field has the dimension of energy, while the heavy meson or antimeson fields H (Q) or H (Q) have dimensions of E 3/2 , as we already mentioned.

Appendix A3: H Hγ interactions
The magnetic coupling of the photon to the s-wave heavy mesons is described by the Lagrangian [57,58] where F μν = ∂ μ A ν − ∂ ν A μ , with A μ the photon field (C A μ (x)C −1 = −A μ (x)), Q = Diag(2/3, −1/3) is the light quark charge matrix, and Q is the heavy quark electric charge (in units of the proton charge e = √ 4πα). For charm (bottom) Q = 2/3 (Q = −1/3). Besides, m Q is the heavy quark mass and β is the parameter introduced in Ref. [57]. These two terms describe the magnetic coupling due to the light (preserves HQSS) and heavy quarks (suppressed by 1/m Q ), respectively. Both terms are needed to understand the observed electromagnetic branching fractions of the D * + and D * 0 because a cancellation between the two terms accounts for the very small width of the D * + relative to the D * 0 [5].
In the non-relativistic constituent quark model β = 1/m q ∼ 1/330 MeV −1 , where m q is the light constituent quark mass. Heavy meson chiral perturbation theory provides contributions from Goldstone boson loops, which give O( √ m q ) corrections to the decay rates [57]. If these loop corrections are evaluated in an approximation where heavy hadron mass differences are neglected, the correction to the above formulas can be incorporated by making the following replacements [57]: with f K ∼ 1.2 f π . 24 We use a convention such that φ = φx −iφy √ 2 creates a π − from the vacuum or annihilates a π + , and the φ z field creates or annihilates a π 0 . We adopt the usual convention C( τ · φ)C −1 = ( τ · φ) T .
Appendix B: Validity of the perturbative treatment of the DD for the X 2 In this appendix, we will argue that the d-wave DD may be treated perturbatively in the 2 ++ system. Even though this was already discussed in Ref. [7], we have included here a new argument grounded on a different EFT to make a more compelling case on the smallness of this contribution to the X 2 mass. We will compare the power counting of the selfenergy diagrams of the X 2 from the d-wave DD and the swave D * D * two-point loops; see Fig. 10. If the DD loop is suppressed in comparison with the D * D * one, it will validate the perturbative treatment of the DD. Because in our case the heavy mesons are non-relativistic, we can apply a velocity counting for the loops analogous to the power counting of the heavy meson loops in heavy quarkonium transitions [70,71]. For the D * D * loop, the velocity counting of the selfenergy reads where g S denotes the value of the s-wave coupling of the X 2 to the D * D * , v denotes the velocity of the D * meson, v 5 is for the loop integral measure, since the non-relativistic energy is counted as O(v 2 ), and 1/(v 2 ) 2 accounts for the two non-relativistic propagators. Similarly, for the DD loop, denoting the velocity of the D meson by w, the velocity counting is given by where g D is the d-wave coupling constant normalized to have the same dimension as g S , and the factor of w 4 in the denominator comes from the two d-wave vertices.
Therefore, we obtain the ratio The question is now how g D compares with g S . We can estimate g D by considering the one-pion exchange diagram considered in this work as illustrated in Fig. 11. Because the X 2 is very close to the D * D * threshold, we should count each of the D * (D * ) propagators as 1/v 2 . This is equivalent to affirm that the cut due to the D * D * in the triangle diagram in Fig. 11 is the same as that in the D * D * bubble diagram of the X 2 self-energy. Thus, we can count the D * in both diagrams in the same way. But the pion propagator should be counted differently. The reason is that because the X 2 couples to the DD in a d-wave, the momenta in the D * Dπ vertices of the one-pion exchange diagram should become the momenta of the D andD, q D = m D w, and the pion momentum is of the same order as we discussed in Sect. 3.1.1. This is to say that the pion propagator should be counted as 1/w 2 rather than 1/v 2 . Thus, expressing the content of Fig. 11 in terms of power counting gives where g is the axial coupling constant in Eq. (A14), and χ = 4π f π is the hard scale for the chiral expansion. Numerically, for the case of the X 2 , we have w (M X 2 −2m D )/m D 0.38, and v (M X 2 −2m D * )/m D * ∼ 0.06 if we take 7 MeV as the binding energy (recall that we have the charged D * D * channel explicitly whose threshold is around 7 MeV above the neutral one). With these values, we use Eqs. (B3) and (B4), which leads to g D ∼ 0.6 g S , to The above value suggests a high suppression of the d-wave DD in comparison with the s-wave D * D * . We notice that the power counting of Ref. [7] indicates that the size of the DD loop is N 4 LO (next-to-next-to-next-to-next-to-leading order), in line with the velocity power counting arguments.