Exact unified tetraquark equations

Recently we formulated covariant equations describing the tetraquark in terms of an admixture of two-body states $D\bar D$ (diquark-antidiquark), $MM$ (meson-meson), and three-body-like states where two of the quarks are spectators while the other two are interacting [Phys. Rev. D 107, 094014 (2023)]. A feature of these equations is that they unify descriptions of seemingly unrelated models of the tetraquark, like, for example, the $D\bar D$ model of the Moscow group [Faustov et al., Universe 7, 94 (2021)] and the coupled channel $D \bar D-MM$ model of the Giessen group [Heupel et al., Phys. Lett. B718, 545 (2012)]. Here we extend these equations to the exact case where $q\bar{q}$ annihilation is incorporated explicitly, and all previously neglected terms (three-body forces, non-pole contributions to two-quark t matrices, etc.) are taken into account through the inclusion of a single $q\bar{q}$ potential $\Delta$.


I. INTRODUCTION
In 1992, Khvedelidze and Kvinikhidze (KK) formulated exact covariant equations for four particles interacting via pairwise interactions [1].A feature of these equations is that pairwise interactions enter the theory not only via single two-body kernels K a describing the scattering of a particle pair a while particle pair a ′ is spectating, but also via subtractions of products of kernels K a K a ′ , which are needed to avoid overcounting, see Eq. (32) below.
However, the equations are rearranged in such a way that all two-body kernels, including these subtraction terms, disappear, with two-body interactions entering the full theory as a sum of single two-body t matrices T a and products of t matrices T a T a ′ .In particular, the final equations use as input the amplitudes T aa ′ , defined as the full 4-body t matrix where all interactions are switched off except those within the pairs labelled by a and a ′ , and which themselves are given in terms of T a and T a ′ as where G 0 a and G 0 a ′ are the Green functions describing the free propagation of particle-pairs a and a ′ , respectively [1,2].
In 2012, members of the Giessen group, Heupel et al. [2], applied KK's formulation to the case of two quarks (2q) and two antiquarks (2q), introducing the further simplifications where the first two terms on the right hand side (rhs) of Eq. ( 1) are neglected, while the two-body t matrices in the last term T a T a ′ are approximated by their meson (M), diquark (D), and antidiquark ( D) pole contributions.The resulting equations were then reduced to a set of coupled two-body equations for the MM-tetraquark and D D-tetraquark amplitudes φ M and φ D , respectively, as illustrated in Fig. 1.A feature of the Giessen group's approach is that it is based on a rigorous field-theoretic derivation for the 2q2q system where all approximations are clearly specified.Numerical solutions of these equations have been pursued by this group in a series of recent publications [2][3][4][5].
A separate approach has been followed for many years by the Moscow group, Faustov et al. [6][7][8][9], who model tetraquarks as a diquark-antidiquark system.The Moscow group's model can be viewed as being based on the solutions of the bound-state equation for the D D-tetraquark amplitude φ D , as illustrated in Fig. 2. As seen from this figure, the kernel of the equation consists of a single term where a q q pair scatters elastically in the presence of spectating q and q quarks.More specifically, the Moscow model corresponds to the case 1. Tetraquark equations of the Giessen group [2][3][4][5].Form factor φ M couples the tetraquark to two mesons (dashed lines), and form factors φ D couples the tetraquark to diquark-antidiquark states (double-lines).Quarks (antiquarks) are represented by left (right) directed lines.
where T q q, the t matrix describing the mentioned q q scattering, is expressed as a sum of two potentials where V gluon is the q q one-gluon-exchange potential and V conf is a local confining potential. 1  Recently we have shown how these two seemingly unrelated approaches to the tetraquark can be unified by providing a common theoretical basis for the two separate models [10].
The key idea is based on the observation that the first two terms of Eq. ( 1), if retained in the application of KK's formalism to the 2q2q system, would generate three-body-like intermediate states of the form qq(T qq ), q q(T qq ), and q q(T q q), where two of the quarks are Diquark-antidiquark bound state equation encompassing the Moscow group's approach [6][7][8][9].The form factor φ D couples the tetraquark to diquark and antidiquark states (both represented by double-lines).Shown is the general form of the kernel where one q q pair interacts (the red circle representing the corresponding t matrix T q q) while the other q q pair is spectating.
1 To be precise, the Moscow group uses quasipotential bound state form factors instead of the D → qq form factor Γ 12 (p, P ) and the D → q q form factor Γ 34 (p, P ), appearing as small blue circles in Fig. 2. Formally, this is equivalent to assuming that Γ 12 (p, P ) and Γ 34 (p, P ) do not depend on the longitudinal projection of the relative 4-momentum p with respect to the total momentum P of the two quarks or two antiquarks.
spectators while the other two are interacting, the last of which, namely q q(T q q), would account for the kernel of the Moscow group's bound state equation.Thus, in Ref. [10], we applied KK's equations to the 2q2q system, similarly to Heupel et al. [2], but we retained the first two terms on the rhs of Eq. (1).In this way we obtained equations for the tetraqaurk which can be expressed in matrix form as where φ is a column of the tetraquark form factors, D is a diagonal matrix consisting of the two-meson propagator MM and the diquark-antidiquark propagator D D, namely, where we have used the same graphical representation of the tetraquark form factors as in Figs. 1 and 2, and where the kernel of the equation (term in square brackets) consists of an infinite series, the first two terms of which, V (0) and V (1) , are defined in terms of Feynman graphs as shown in Fig. 3 and Fig. 4, respectively.The ellipsis in Eq. ( 3) represents terms consisting of Feynman diagrams, similar to those of V (1) , but describing the full multiple scattering of quarks in intermediate state.We refer to Eq. (3) as the "unified tetraquark equations" [10].
That these equations unify the works of the Giessen and Moscow groups can be seen from the fact that the equations of Heupel et al. correspond to keeping just the term V (0) term in the kernel of Eq. ( 3): while the Moscow group model follows by keeping just the term V (1) term in the kernel of Eq. (3): but with all but the D D → D D elements of V (1) set to zero.Our formulation suggests that the Giessen and Moscow groups have been investigating two non-overalpping aspects of the same model of the tetraquark.
3. Feynman diagrams making up the elements of the coupled channel M M − D D kernel matrix V (0) of Eq. (79): Solid lines with leftward (rightward) arrows represent quarks (antiquarks), dashed lines represent mesons, and double-lines represent diquarks and antidiquarks.4. Feynman diagrams making up the elements of the coupled channel M M − D D kernel matrix V (1) of Eq. (81): 12 = 2 Γ1 In the present work, we extend the unified tetraquark equations of Ref. [10], namely Eq. ( 3), to include all contributions that have been neglected in the derivation of these equations.In particular, we explicitly include q q annihilation within the same formalism as used in Eq. ( 3), and at the same time, take into account all possible other contributions (e.g.three-body forces, non-pole contributions to the t matrices within the products T a T a ′ , correction terms to overcounted or undercounted contributions, etc.) through the introduction of a q q amplitude ∆ that is defined to consist of all corrections needed to restore the exactness of the model.That this is possible has been demonstrated by us in Ref. [11] on the example where the first two terms of Eq. ( 1) have been neglected.By following the same approach as [11], but retaining the first two terms of Eq. (1), we obtain equations of the same form as the unified tetraquark equations of Eq. ( 3), but that include q q annihilation in the way prescribed by quantum field theory, and that are otherwise exact due to the inclusion of the ∆ amplitude.Similarly to Eq. ( 3), these equation can be expressed in matrix form as where ϕ is a rank-3 column of tetraquark form factors, D is a diagonal matrix consisting of the two-meson propagator MM, the diquark-antidiquark propagator D D, and the q q propagator Q Q, namely, where Φ M , Φ D , and Γ * are form factors (together with their graphical representations) that couple the tetraquark to MM, D D and q q states, respectively.The kernel of Eq. ( 7) (term in square brackets) consists of two parts: (i) the kernel V, defined formally as the sum of the infinite series which is the kernel of the corresponding equation without q q annihilation, i.e.Eq. ( 3), but expressed in rank-3 form as illustrated for the first two terms V (0) and V (1) in Fig. 5 and Fig. 6, respectively, and (ii) the kernel W, illustrated in Fig. 7, which consists of q qirreducible amplitudes NM , ND , N M , and N D connecting MM and D D states to q q states [as described in the figure caption of Fig. 7], and the aforementioned q q elastic scattering amplitude ∆.
A feature of the tetraquark description as modelled by Eq. ( 7), is that it is formally exact no matter what model or approximations are used in calculating the kernel V, or what model is used for the "N amplitudes" NM , ND , N M , and N D .This is due to the fact that the ∆ appearing in W, by definition, consists of compensating terms that keep the tetraquark description exact.To illustrate this point, let us consider the model where , corresponding to the Giessen model, but where q q annihilation is taken into account through the inclusion of the kernel W where the N amplitudes are modelled as in Fig. 7.We can then formally adjust the ∆ amplitude to the amplitude ∆ 0 , with the resulting kernel W being denoted by W 0 , such that the tetraquark equation is exact; that is, having the same solution ϕ as Eq. ( 7).That this is the case can be seen from the discussion in Sec.2.1 below, where we show that the three equations represented by Eq. ( 7) are equivalent to the single equation for the tetraquark → q q form factor Γ * , Eq. ( 23), whose q q kernel K (2) , given by Eq. ( 22), can be made exact by adjusting ∆, no matter what model is chosen for the N amplitudes or potential V .This means that ∆ 0 contains within it information about V (1) and all the rest of the contributions not explicitly appearing in Eq. ( 10).Similarly, we can consider the model where V = V (0) + V (1) , and then adjust the ∆ amplitude to the amplitude ∆ 1 , with the resulting kernel W being denoted by W 1 , such that the tetraquark equation is exact.The amplitude ∆ 1 then contains within it information about V (2) and all the rest of the contributions not explicitly appearing in Eq. (11).Equation ( 11) is then the exact tetraquark equation that includes q q annihilation into the unified description of the Giessen and Moscow models of the tetraquark as described by Eq. (61) of [10].It is clear from these examples that the tetraquark equations, Eq. ( 7), have the generality to unify a large variety of other possible tetraquark models.Moreover, the same equations can be applied generally to four-body systems that couple to 2-body channels, as for example the two-electron plus two-positron system, the two-nucleon plus two-antinucleon system, etc.

II. DERIVATION
The main results of this paper have been summarised above, and are embodied in the exact tetraquark equations of Eq. (7).In this section, we aim to provide as short a derivation of these equations as possible, but without sacrificing rigour or clarity.For this purpose we consider only the non-identical meson case, noting that the additional step of symmetrising MM states for the case of identical mesons, can be implemented simply by referring to Ref.
[10] where this procedure is described.
To derive Eq. ( 7), we incorporate q q absorption into the unified tetraquark equations of Ref. [10], by exploiting the same procedure as we previously used in Ref. [11] to incorporate q q absorption into the tetraquark model of the Giessen group.This procedure is described in Sect.2.1.In order to derive the precise mathematical expressions corresponding to the diagrams of Figs. 1 -7, it is also necessary to discuss the derivation of the unified tetraquark equations themselves, at least for the case of non-identical mesons.Although this derivation has been previously given in a part of Ref. [10], for completeness we reproduce this derivation in Sect.2.2 below, although using a somewhat simpler notation.
Elements of the kernel matrix W: W 31 = NM and W 32 = ND , are q q-irreducible amplitudes for M M → q q and D D → q q, respectively; W 13 = N M and W 23 = N D , are the corresponding amplitudes for q q → M M and q q → D D. Illustrated are models for these amplitudes as given by Eqs.(21).W 33 = ∆ is the q q amplitude consisting of correction terms needed to make Eq.( 7) exact.
A. Describing the tetraquark in the presence of q q annihilation In the absence of q q annihilation (either because annihilation is being neglected or because quantum numbers forbid it), the tetraquark can be simply defined as the four-body bound state of two quarks and two antiquarks as signalled by a pole in the full 2q2q Green function G (4) .However, in the presence of q q annihilation (when annihilation is physically realisable and taken into account), a pole in G (4) becomes an insufficient criterion for a tetraquark, as can be seen from the exact field-theoretic expression for G (4) : where G ir is the q q-irreducible part of G (4) , G 0 is the disconnected part of the two-body q q Green function G (2) (corresponding to the independent propagation of q and q in the s channel), and G ) is the sum of all q q-irreducible diagrams corresponding to the transition q q ← 2q2q (2q2q ← q q).Equation (12) shows that any pole in the q q Green function G (2) , even if not associated with a tetraquark, is automatically a pole in G (4) .
Here we consider the case of the tetraquark in the presence of annihilation, meaning that the transition Green functions G (4−2) ir and G (2−4) ir in Eq. ( 12) are not zero.For such a case, both G (4) and G (2) will display simultaneous poles corresponding to a tetraquark of mass M, so that as P 2 → M 2 where P is the total four-momentum of each system, In Eq. ( 13), Ψ is the tetraquark 4-body (2q2q) bound state wave function, while Γ * is the form factor for the disintegration of a tetraquark into a q q pair.
As previously introduced in Ref. [11], our key idea for describing a tetraquark in the presence on q q annihilation, is to: (i) express G ir in terms of some model for the 2q2q system in absence of q q annihilation, and (ii) to then express the q q kernel K (2) , defined through the Dyson equation in terms of G ir as where operators A (4−2) (A (2−4) ) are q q-irreducibe amplitudes describing the transitions 2q2q ← q q (q q ← 2q2q) in the chosen model, and ∆ is defined as the q q-irreducible q q four-point function consisting of all contributions not accounted for by the last term of Eq. (15).A feature of this approach is that no matter what model is chosen for G (4) ir and A (4−2) (A (2−4) ), Eq. ( 15) is exact in view of the definition of ∆.Moreover, all terms contained in ∆ are known thanks to its clear definition; for example, it contains sums of multiple-gluon exchanges.This fact enables one to systematically improve the approximations used in modelling the tetraquark, as for example, to take into account one-gluon exchange explicitly.
Furthermore, as the purpose of this work is to extend the unified tetraquark equations to include q q annihilation and other corrections, we assume that G ir (which is the Green function determining the unified tetraquark equations) has a pole at P 2 = M 2 0 , where M 0 is the tetraquark mass described by these equations; thus where Ψ 0 is the corresponding tetraquark wave function, and R is a background term.It is then straightforward to prove 2 2 It follows from the second of Eqs.(13), Eq. ( 14) and Eq. ( 15) that which shows that in the absence of annihilation and other corrections (described by ∆), the tetraquark's mass M 0 plays the role of a "bare" mass, and that the inclusion of annihilation and other corrections, then shifts the bare mass to the physical tetraquark mass M.
Here we apply this approach to the case where is expressed in terms of the tetraquark model as described by the unified tetraquark equations, Eq. ( 3), and as derived specifically in Sect.2.2 below.In particular, we write this term as where G is the Green function in MM − D D space generated by the kernel of Eq. ( 3), thus and N (simlarly N) is a matrix of transition amplitudes NM (for q q ← MM) and ND (for q q ← D D).We follow Ref. [12] where the simplest model for such amplitudes was realised, namely those described graphically as and given analytically in the notation of Sect.2.2.3 below, as where S 23 is the quark propagator connecting quark lines 2 and 3. Expressions for the kernels making up V are also derived in Sect.2.2.3, with V (0) being given by Eq. ( 79), and V (1) by Eq. (81).Thus It is also evident from Eq. ( 14) and the second of the relations in Eq. ( 13), that the tetraquark state will also satisfy the two-body equation Using Eq. ( 16) with P 2 = M 2 in the above expression, and defining X = Ψ0 A (4−2) G (2) 0 Γ * , one obtains which then implies Eq. ( 17).
Using Eq. ( 22) in Eq. ( 23), one obtains where with φ being the column of tetraquark form factors φ M and φ D as in Eq. ( 4).By introducing a column of tetraquark form factors that also includes the form factor Γ * , and defining G (2) 0 ≡ Q Q, Eq. ( 24) and Eq. ( 25) can be expressed succinctly as the rank-3 matrix equation Eq. ( 7).

B. Unified tetraquark equations
The exact unified tratraquark equations, Eq. ( 7), take as an input the kernel of the unified tetraquark equations of Ref. [10], including the separate form factors for the M ↔ q q, D ↔ qq, and D ↔ q q processes that may be used to model the input MM ↔ q q and D D ↔ q q transition amplitudes, as for example in Eqs.(21).For this reason, it is essential to refer to the derivation of the unified tetraquark equations where all the expressions used in the current work as input, are developed.Although this derivation can be found in full in Ref. [10], for the purposes of the current work, it is sufficient to consider this deriviation just for the case of non-identical mesons.For completeness, we present this derivation here, following closely the presentation of Ref. [10], but using a somewhat simpler notation.

Four-body equations for distinguishable quarks
We consider the 4-body system consisting of 2q and 2q treated as distinguishable particles where q q annihilation (or creation) is not allowed.The 4-body Green function describing this system is then G (4) ir , the q q-irreducible part of the full Green function G (4) , as in Eq. ( 12).One can then introduce the 4-body kernel K and corresponding t matrix T (both q q-irreducible) through the equations where G 0 describes the free propagation of all four particles.We assign labels 1,2 to the quarks and 3,4 to the antiquarks and following Ref.[1], introduce an index a ∈ {12, 13, 14, 23, 24, 34} which enumerates the six possible pairs of particles.Similarly, we introduce the double index aa ′ ∈ {(13, 24), (14, 23), (12,34)} which enumerates the three possible two-pairs of particles, and use the Greek index α as an abbreviation for aa ′ such that α = 1 denotes aa ′ = (13, 24), α = 2 denotes aa ′ = (14, 23), and α = 3 denotes aa ′ = (12,34).
Using this labelling scheme, it is useful to consider the Green function G aa ′ defined as the part of G ir where all interactions are switched off except those within the pairs a and a ′ .One can then introduce the corresponding kernel K aa ′ and t matrix T aa ′ through the equations Note that we drop the superscript indicating the number of particles described, when this is otherwise implied by the meaning of the subscripts.It is evident that in the theory where only pairwise interactions are allowed, one has that This is a key expression for K as it is of similar form to that describing the kernel of a 3-body system interacting via pairwise interactions, thereby leading to the formulation of 4-body equations by analogy to those of 3-body equations [1].It is also clear that where G a (similarly G a ′ ) is the full 2-body Green function for the scattering of particles of pair a, with corresponding 2-body kernel K a and t matrix T a defined through where G 0 a describes the free propagation of the particles in pair a.Using the above equations it is then easy to show Eq. ( 1) and similarly that where the presence of a minus sign in the last term of Eq. ( 32) is necessary to avoid overcounting.
The 2q2q bound state form factor for distinguishable quarks is then The four-body kernels K α can be used to define the Faddeev components of Φ as From Eq. (33) follow Faddeev-like equations for the components, where δαβ = 1 − δ αβ .

Four-body equations for indistinguishable quarks
The bound state equation for two identical quarks 1, 2 and two identical antiquarks 3, 4, is modified from Eq. (33) to where now the kernel K is antisymmetric with respect to swapping quark or antiquark quantum numbers either in the initial or in the final state; that is, where the exchange operator P ij swaps the quantum numbers associated with particles i and j in the quantity on which it is operating.The factor 1 4 in Eq. ( 37) is a product of the combinatorial factors 1  2 , one for identical quarks and another for identical antiquarks.
In order to distinguish the quantities Φ, K, and T in the present case of indistinguishable quarks, from their distinguishable-quark counterparts discussed in the previous section, we shall simply use a superscript d to denote the quantities that apply to the case of distinguishable quarks.Thus Eq. (33) for distinguishable quarks shall now be written as Eq. ( 36) as and Eq. ( 1) as The kernel K that is antisymmetric in the way specified by Eq. ( 38), can be represented as where K d is symmetric with respect to swapping either quark or antiquark quantum numbers in the initial and final states simultaneously, P 12 K d P 12 = P 34 K d P 34 = K d .This symmetry property of K d can be written in the form of commutation relations and follows directly from the following relations implied by Eqs.(32): Due to the antisymmetry properties of K as specified in Eq. ( 38), the solution of the identical particle bound state equation, Eq. (37), is correspondingly antisymmetric; namely, P 34 Φ = P 12 Φ = −Φ.However, because K d usually corresponds to a fewer number of diagrams than K, rather than solving Eq. (37), it may be more convenient to determine Φ by antisymmetrising the solution Φ d of the bound state equation for distinguishable quarks, as Then, in view of the commutation relations of Eq. ( 43), if the solution Φ d exists, its antisymmetrized version as given by Eq. ( 45), also satisfies the bound state equation for distinguishable quarks, Eq. (39), as well as the one for indistinguishable ones, Eq. (37): As Φ satisfies the same bound state equation as Φ d , the kernels K d α can again be used to define Faddeev components, but this time for Φ: where In view of Eqs.(44), the Faddeev components Φ α have the following properties: Since Φ satisfies the same bound state equation as Φ d , the components Φ α satisfy the same Faddeev-like equations as for distinguishable quarks, Eq. (40), Without loss of generality, we can assume that the solution of Eq. ( 51) has the symmetry properties of Eq. (50) [10].We also note that the input 2-body t matrices T d 12 and T d 34 can be antisymmetrized by defining so that By defining further that we have made Eq. ( 1) apply, with unchanged notation, also for the case of indistinguishable quarks.Moreover, one can now write Eq. ( 51) simply as where dropping the superscript d from T d α for the case α = 3 can be justified by multiplying Eq. ( 51) for this case by (1 − P 12 ) and using the symmetry properties of Eq. (50): The advantage of using T 3 rather than T d 3 in this equation stems from the fact that the physical (antisymmetric) t matrices for qq and q q scattering are T qq = (1−P 12 )T d 12 = 2T 12 and T q q = (1 − P 34 )T d 34 = 2T 34 , respectively, therefore it is convenient to use the antisymmetric T 12 and T 34 as the input qq and q q t matrices.
Writing out Eqs. (55) in full, but to save notation without explicitly showing G 0 's, one has For physical (antisymmetric) solutions of Eqs.(57), only two of these three equations are independent.For example, Eq. (57b) can be written as where Eq. (50) and T 2 = P 12 T 1 P 12 have been used; then, after a further application of P 12 , one obtains Eq. (57a).Choosing Eq. (57a) and Eq.(57c) as the two independent equations, we can use Φ 2 = −P 12 Φ 1 to obtain closed equations where, necessarily, P 12 Φ 3 = −Φ 3 .Again, without loss of generality, we can assume that the solutions of Eq. ( 59) have the symmetry properties of Eq. (50) [10].Equation (59b) can be further simplified using P 12 Φ 1 = P 34 Φ 1 and the assumption that T 12 and T 34 are antisymmetric in their labels, so that In this way Eqs.(59) take the form Again, without loss of generality, we choose a solution of Eqs. ( 61) which has all the symmetry properties of Eq. (50).
3. Tetraquark equations with exposed q q(T q q ), qq(T q q), and q q(T qq ) channels Choosing Eqs.(61) as the four-body equations describing a tetraquark, they may be expressed in matrix form as where Writing where we have that where Thus and consequently To be close to previous publications we choose a separable approximation for the two-body t matrices in T × 1 and T × 3 (but not necessarily in T + 1 and T + 3 ); namely, for a ∈ {13, 24, 12, 34} we take where D a = D a (P a ) is a propagator whose structure can be chosen to best describe the twobody t matrix T a , and Γ a is a corresponding vertex function.In the simplest case, one can follow previous publications and choose the pole approximation where D a (P a ) = 1/(P 2 a −m 2 a ) is the propagator for the bound particle (diquark, antidiquark, or meson) of mass m a .In view of Eq. (53), note that We can thus write where the 2 × 2 matrix potential (with reinserted G 0 ) is Expanding the term in square brackets in powers of T + (i.e., with respect to the contribution of intermediate states q q(T q q), qq(T q q), and q q(T qq )), it turns out that each of the first two terms of this expansion corresponds to different existing approaches to modelling tetraquarks in terms of MM −D D coupled channels.In particular, the lowest order term where consists of Feynman diagrams illustrated in Fig. 3, and corresponds to the Giessen group (GG) model of Heupel et al. [2] where tetraquarks are modelled by solving Eq. ( 5).
Similarly, the first order correction (without the lowest order term included) is which consists of Feynman diagrams illustrated in Fig. 4, and corresponds to the Moscow group (MG) model of Faustov et al. [7] where they modelled tetraquarks by solving Eq. ( 6) albeit, with only diquark-antidiquark channels retained.It is an essential feature of this formulation, that it is the sum of the potentials V (0) and V (1) , each associated with the separate approaches of the Giessen and Moscow groups, with tetraquarks modelled by the bound state equation φ = [V (0) + V (1) ]Dφ, that results in a complete MM − D D coupled channel description up to first order in T + [i.e., up to first order in intermediate states where one 2q pair (qq, q q, or q q) is mutually interacting while the other 2q pair is spectating].
Finally we need to point out that for the case of identical mesons, the MM states in the above MM − D D coupled channel description need to be symmetrised.This does not change the structure of any of the above equations, or the essential results.The details of this symmetrisation are given in Ref. [10].