The interplay between compact and molecular structures in tetraquarks

Due to the cluster reducibility of multiquark operators, a strong interplay exists in tetraquarks between the compact structure, resulting from the direct confining forces acting on quarks and gluons, and the molecular structure, dominated by the mesonic clusters. This issue is studied within an effective field theory approach, where the compact tetraquark is treated as an elementary particle. The key ingredient of the analysis is provided by the primary coupling constant of the compact tetraquark to the two mesonic clusters. Under the influence of this coupling, an initially formed compact tetraquark bound state evolves towards a new structure, where a molecular configuration is also present. In the strong-coupling limit, the evolution may end with a shallow bound state of the molecular type. The strong-coupling regime is also favored by the large-$N_c$ properties of QCD. The interplay between compact and molecular structures may provide a natural explanation of the existence of many shallow bound states.

The reason these two competing alternatives are arising is related to the fact that the multiquark operators that generate multiquark states are not color-irreducible, in contrast to the ordinary hadron case, in the sense that they are decomposable along combinations of clusters of ordinary hadron operators [16,46]. There are, therefore, two different ways of considering the construction of a multiquark state, as depicted above. The main issue is which one reproduces the most faithful description of reality. The existence or emergence of hadronic clusters inside a multiquark state might be an indication of a kind of instability when the state is built out of confining forces. It might, at some stage, dislocate into the clusters, or, in the case of a bound state, evolve towards another state, dominated by the clusters, that is, towards a molecular-type state.
A theoretical hint to analyze the problem is provided by the study of the energy balance of the two types of configuration [29]. This is most easily done for heavy or static quarks, for which lattice calculations are available [47][48][49][50][51][52][53]. In the strong coupling limit of lattice theory, analytic expressions are obtained by means of Wilson-loop expectation values, which satisfy the area law, or, more generally, are saturated by minimal surfaces; the corresponding predictions have been verified by direct numerical calculations on the lattice (cf. previous references). The qualitative result that emerges from the latter calculations is the following. When the quarks and antiquarks of the multiquark state are gathered into a small volume, it is the compact multiquark configuration that is energetically favored, while in situations where the quarks and antiquarks are seperated from each other at larger distances, it is the cluster-type configurations that are energetically favored. Therefore, there are nonzero probabilities for each type of configuration to occur for the description of the multiquark state. However, since quarks and antiquarks are moving objects and generally reaching, even with small probabilities, large distances, whose integrated volume may be much larger than the small volume of the compact configuration, one expects that an initially formed compact state would gradually evolve towards a cluster-type configuration, typical of a molecular state. In coordinate space, the core of the multiquark state would be better described by the compact representation, whereas the outer layer would be better described by the molecular representation. The relative weight of each representation would, of course, depend on specific parameters, such as the quark masses and the quantum numbers that are involved 1 . The above results have led, in spectroscopic calculations, to the introduction of the concept of "configurationspace-partitioning" (or for short "geometric partitioning") which is realized in the so-called "flipflop" potential model [55][56][57][58][59][60][61][62][63][64][65][66][67], which takes into account more faithfully the role played by each configuration in the formation of multiquark states. However, because of the complicated nature of the constraints, which are coordinate dependent, this model, apart from simplified cases, has not yet led to full spectroscopic results, to be compared, on quantitative grounds, with experimental data.
In principle, if the tetraquark bound state problem could have been solved with high precision, taking into account all interactions that act between quarks and gluons, one would obtain the exact knowledge about its structure. Unfortunately, this is not currently the case; one is obliged to adopt approximations and proceed step by step by including additional inputs to improve the predictions. As mentioned above, the simplest approximations are either the compact scheme or the molecular scheme. While the latter scheme, based on hadron-hadron interactions, has a sufficiently developed theoretical background, the former one needs further analysis. In the diquark model, the diquark being considered in particular in its color-antisymmetric representation (ignoring here spin degrees of freedom) within a very small volume (pointlike or almost pointlike approximation), one always has tetraquark (or multiquark) bound states [68][69][70][71][72][73][74]; this is due to the fact that, in that approximation, all forces acting on the various small volumes (or points) are of the attractive confining types. This is not the case of the molecular scheme, where the occurrence of a bound state depends on the strength of the attractive forces. Therefore, in the compact scheme, one is entitled to start with a tetraquark (or multiquark) candidate, with all its accompanying multiplicities. The main problem that is encountered here is the evaluation of the effect the mesonic (or hadronic) clusters could have on that bound state. The key ingredient which enters in the description of that effect is the effective coupling constant of the compact tetraquark to the meson clusters. One easily guesses that stronger the latter quantity is, more important is the transformation of the compact state into a cluster-like state, which ultimately might take the appearance of a molecular-type state.
It is the main objective of the present paper to evaluate the interplay between the compact and molecular structures of possibly existing tetraquark states. For this, we shall adopt methods of effective field theories, remaining at the same time at the level of simple qualitative features.
In effective field theories of mesons, interactions are described by meson exchanges and by heavy quarks and on the basis of the properties of the confining interactions at short distances, the existence of a tetraquark bound state. They, however, recognized that the large-distance dynamics should be better described by meson-meson interactions and switched for the description of that domain to chiral perturbation theory.
contacts. At lower energies, the exchanged meson fields can be integrated out and one remains only with a theory with contact-type interactions [75,76], which we call here the "lower-energy" theory. The correspondence between the parameters of the two types of theory is not, however, simple and a physical understanding of the results necessitates a more detailed investigation. We devote Sec. II to a presentation of this aspect of the problem. Taking into account the various physical conditions and known results, we propose, in Sec. III, an empirical formula, which relates, in an explicit way, the coupling constant of the lower-energy theory to that of the Yukawa-type theory, and allows an easy understanding of the conditions in which a bound state may emerge.
In order to emphasize the qualitative features of the approach, we neglect spin effects and limit ourselves to scalar interactions with scalar particles, ultimately considered in the nonrelativistic limit. The resulting effective theory is then used to study the meson-meson interaction through the scattering amplitude and the determination of the possibly existing bound state properties. The corresponding scattering length and effective range are evaluated in Sec. IV. Some of the results of Secs. III and IV are well known in the literature and are presented here for the purpose of introducing the method of approach that is applied for more general cases.
The case of compact tetraquarks is studied in Sec. V. With respect to the meson clusters, the compact tetraquark can be represented, in first approximation, as an elementary particle, whose internal structure would be relevant only at short-distance scales. It is then essentially characterized by its mass and quantum numbers and described by means of its propagator. The tetraquark, because of its internal structure, has necessarily interactions with meson pairs, and in particular with those lying closest to its mass. In the simplest case of one meson-pair, one may introduce a bare coupling constant for the interaction tetraquark-two-mesons and analyze its influence on the properties of the tetraquark through the radiative corrections it induces. For the bound state case, it is assumed that the bare compact tetraquark mass lies below the two-meson threshold. It turns out that, in general, the compact-tetraquark-two-meson interaction shifts the binding energy of the tetraquark to lower values. In the strong coupling limit, the shift may even transform the compact tetraquark into a shallow bound state, typical of loosely bound hadronic molecules. This phenomenon is best represented by means of the "elementariness" parameter Z, introduced by Weinberg [36], which measures the probability of a bound state to be considered as elementary, the complementary quantity, (1 − Z), representing its "compositeness". In the strong coupling limit, described above, Z takes small values, approaching zero. In parallel, the physical coupling constant of the tetraquark to two mesons tends also to zero in the same limit. The value of Z is also measured by means of the scattering length and the effective range parameter, the latter taking negative values when Z = 0.
These results, which are the main outcome of the present paper, provide a more refined understanding of the structure of observed tetraquark candidates. Many shallow bound states, which are typical of molecular states, might have a compact origin, provided Z = 0. They would be the result of the deformation, under the influence of the mesonic clusters, of the initially formed compact tetraquark. This is another illustration of the dual representation of the tetraquark found in lattice theory on the basis of the energy balance analysis [47][48][49][50][51][52][53].
On more general grounds, shallow bound states are usually considered as belonging to universality classes, whose binding energy values are not naturally explained by means of the interaction scales of the system [77]. The previous results bring a new lighting to that problem in the case of tetraquarks. The coupling constant compact-tetraquark-meson-clusters introduces an additional scale parameter, whose strong-coupling limit naturally explains the origin of the shallowness.
The case of resonances occurs when the bare mass of the compact tetraquark lies above the twomeson threshold. Here, however, contrary to the bound state case, additional constraints appear for the existence of a physical resonance. The latter may exist only in the weak-coupling regime.
Large values of the coupling constant resend the state to the bound state domain, while for a finite interval of the coupling constant, occurring prior to the strong-coupling regime, the tetraquark state may disappear from the spectrum.
Section VI brings complementary informations with respect to Sec. V, by also considering, in addition to the tetraquark-two-meson interaction effect, the influence of meson-meson interactions, which now renormalize the primary (bare) coupling constant and introduce a competing effect coming from direct molecular-type forces. The qualitative conclusions drawn in Sec. V remain, however, valid. Section VII is devoted to an analysis of the problem in the large-N c limit of QCD.
In that limit, the theory provides additional support to the dominance of the strong-coupling regime in the effective interaction compact-tetraquark-meson-clusters. Conclusion follows in Sec. VIII. A few detailed analytic expressions, approximating energy eigenvalues, are gathered in the Appendix.

II. REDUCTION TO CONTACT-TYPE INTERACTIONS
Effective field theories, which result from the integration of fields operating mainly at high energies, are generally characterized by the presence of contact-type interactions. Often, depending on the energy scale that is considered, these coexist with ordinary-type interactions, whose prototype is the Yukawa interaction, responsible of meson exchanges between interacting particles. At lower energies, the exchange-meson fields themselves are integrated out and one remains only with contact-type interactions. A representative example of such a theory is chiral perturbation theory [38][39][40], in which the interacting particles are the pseudo-Goldstone bosons and where all other massive particle fields have been integrated out. However, explorations of more refined properties, related to spectroscopic problems and to the physical interpretation of numerical values of parameters, may require a more detailed knowledge of the connection between the lower-energy theory and its higher-energy generator.
A. Spectroscopic properties of the higher-energy theory To illustrate the above aspect of the question, we shall consider, following Reference [76], a higher-energy theory, where two massive scalar particles, with masses m 1 and m 2 , interact by means of the exchange of a scalar particle with mass µ, described by the interaction Lagrangian density where, for simplicity, after factorizing, for dimensionality reasons, the mass terms 2m i , we have chosen the same (dimensionless) coupling constant g for the two particles. The fields φ i (i = 1, 2) correspond to those of the external massive particles, while ϕ corresponds to the exchanged particle field.
The scattering amplitude of the process {1(p ′ 1 ) + 2(p ′ 2 ) → 1(p 1 ) + 2(p 2 )} is given, in lowest order in g, by the Born term, or the ladder diagram: where q is the momentum transfer: This term iteratively generates higher-order diagrams, representing the series of ladder-type diagrams, which plays a basic role in a possible production of bound states. They are represented in Fig. 1. A complete study of the theory requires the evaluation of the effects of self-energy and vertex corrections. These, however, do not play a fundamental role in the extraction of the main qualitative features that are of interest for our purpose and will be ignored.
The bound state problem in the nonrelativistic limit can be studied by means of the Schrödinger equation. The latter is obtained from the Bethe-Salpeter equation, with the kernel considered in the ladder approximation, represented by the right-hand side of Eq. (2). One has first to use the instantaneous approximation, where one neglects, in the cm frame, the temporal component of the momentum transfer q, and then to take the nonrelativistic limit. The Schrödinger equation in x-space is: where V (r) is the well-known Yukawa potential: Making in (3) the changes of variable and parameter [78] x = x ′ /µ, one reduces the Schrödinger equation to a dimensionless equation with a single parameter, g 2 : Compared to the Coulomb potential, the Yukawa potential (4) is of the short-range type and hence not all values of g 2 may produce bound states. There exists a critical value of it, g 2 cr = 1.68 × (4π) [76,78], below which bound states do not exist. Their existence is ensured by the inequality ( At the critical value of the coupling constant, a zero-energy bound state appears. When g 2 increases, gradually new bound states appear, while the ground state binding energy itself increases.
In more general cases, one may have the sum of several Yukawa potentials with different meson exchanges. In such cases, one loses the notion of a universal coupling constant squared g 2 . To continue exploring qualitative aspects of the problem, it would be advantageous to approximate the sum by a single Yukawa potential with a mean exchanged mass and a mean coupling constant.
General qualitative properties of the Yukawa potential, concerning the related spectroscopy and the poles of the corresponding S-matrix, could be obtained by considering the soluble model of the spherically symmectric rectangular potential well, which is a prototype of short-range potentials [79,80]. If V 0 is the depth of the potential and R its width, then there is a correspondence between g 2 /(4π) and the product 2m r V 0 R 2 ≡ A 2 . The critical value of A 2 is equal to π 2 /4 = 2.46, which is of the same order of magnitude as g 2 cr /(4π) [Eq. (7)].
B. The lower-energy effective field theory A lower-energy effective field theory can be obtained by integrating out the mediator field ϕ.
One then obtains a theory where the only remaining fields are the φ i s. Their interactions are represented by an infinite series of contact terms with increasing powers of products of φ i s, also possibly containing derivative couplings. They are classified according to their dimensionality and a systematic expansion is organized according to definite power counting rules. The coupling constants of the lower-energy theory are usually determined from matching conditions of the scattering amplitudes calculated in the two theories. At low energies, it is the lowest-dimension operator that is expected to provide the leading contribution. We stick in the following to that term. The corresponding Lagrangian density, for the mutual interaction of particles 1 and 2, is The problem that interests us is to what extent the lower-energy theory in its leading-order approximation, represented by the Lagrangian density (8), provides a faithful description of the higher-energy theory, considered in the ladder-approximation, in the bound state domain.
Similar searches as above have been undertaken long ago in the past, although with a different viewpoint. In this respect, Reference [81], provides a rather detailed account of the corresponding approach. The objective has been to establish an equivalence theorem between the two theories described by Lagrangian densities of the types of (1) and (8), respectively, the external particles being here fermions. Putting aside the question of the nonrenormalizability of the four-Fermi interaction theory, the approach has consisted in considering the first ladder diagram of Fig. 1, together with the radiative corrections of the exchanged meson propagator, and searching for conditions of equivalence of the chain of diagrams of Fig. 2, considered in the t-channel. Our line of investigation is rather different, being based on the effective field theory approach, searching for matching conditions for the leading terms of the scattering amplitude in the s-channel.
Coming back to the present approach with its matching conditions, it turns out that all diagrams of Fig. 1, besides contributing to the coupling constants of higher-dimensional operators, also give contributions to the coupling constant h of (8). Therefore, h, considered as a function of g 2 (or g 2 ), has an expression in the form of an infinite series of it: where h (n) represents the contribution coming from the n-loop diagram (proportional to ( g 2 4π ) n h (0) ). Notice that h (0) is positive. While the first few terms of the series of h are explicitly calculable [76], the complete knowledge of the expression of h in terms of g 2 is not known. It is evident, from the lower-bound (7), that, for the bound state problem, one cannot be satisfied by a perturbative expansion of h, but rather a full nonperturbative expression of it will be needed in order to understand and interpret the bound state properties of the lower-energy theory in terms of those of the higher-energy theory.

LOW-ENERGY THEORY
The nonperturbative properties of the coupling constant h can be deduced from the study of the bound state problem of the low-energy theory. The result of the latter problem is well known in the literature [75,76] and we briefly sketch the corresponding procedure. The scattering amplitude resulting from the chain of diagrams of Fig. 2 is a geometric series expressible in terms of the two-point loop integral J(s): J(s) is divergent, but its divergence can be absorbed in h by renormalization. Defining J r as the regularized (finite) part of J and J div as its diverging part, such that one obtains where h r is the renormalized finite part of h. Notice that this renormalization implies that the unrenormalized h is a vanishing quantity: In the limit where, for finite h r , iJ div → ∞, one has h → 0.
Regularizing J by dimensional regularization and keeping in the diverging part only massindependent terms, J r takes the following form: µ is the mass scale introduced by the dimensional regularization, whose value can be chosen at will; Q(s) is defined in the complex s-plane, cut on the real axis from s = (m 1 + m 2 ) 2 to +∞ and from s = (m 1 − m 2 ) 2 to −∞, and has the expression where (The square-roots are defined for complex s as √ s − m 2 = e iα/2 |s − m 2 |, where α is the angle made by (s − m 2 ) with the positive real axis starting from m 2 .) On the real axis, Q(s) takes the following forms: (The square-roots in (18) are defined with positive values, their arguments representing moduli.
J r (s) is finite at s = 0.) Physical quantities should be independent of µ. Since T eff is such a quantity, this implies that h r itself should be µ dependent and should cancel the µ-dependence of J r . The absorption of the term ln into a redefined 1/h r , which amounts to also absorbing the same quantity into J div , is the simplest way of ensuring this. We shall adopt henceforth that procedure. It is advantageous to define h r as having a simple physical interpretation. A natural choice is the value of T eff at the two-particle threshold [75]. We shall continue using the same notation for the redefined h r and J r ; then the redefined J r obtains the form The scattering amplitude (10) then takes the form, according to (13), Bound states of T will be identified as tetraquark states of the molecular type, whose parameters and ingredients will be labeled by the indices tm; they correspond to the solutions of the equation (We shall often omit, for the simplicity of notation and when no ambiguity is present, the index r from h r and J r .) We restrict ourselves to the case of possible nonrelativistic solutions, located near the two-particle threshold. We introduce the nonrelativistic energy E tm through the definition and, to simplify notations, we shall use henceforth the reduced dimensionless energy variable e and bound state energy e tm , respectively: where m r is the reduced mass [Eq. (3)]. Retaining, in the second of Eqs. (18), the leading term in √ −e tm , which is contained in its first piece, one finds a unique possible solution: .
The existence of the solution is conditioned by a negative value of h. This means that h must have changed sign with respect to the perturbative expression h (0) [Eq. (9)]. To ensure the nonrelativistic interpretation of the solution, one must have a large value of |h|, such that When |h| → ∞, the bound state approaches the two-particle threshold. On the other hand, when h → 0 with negative values, the binding energy increases and tends to ∞.
The stability of the solution (24) in the presence of higher-order terms in √ −e tm can be studied by also considering in J r (s) terms of order (−e), which are contained in the second term of Q(s) [Eq. (18)] and the remaining terms of J r (s) [Eq. (17)]. Accepting condition (25), the resulting equation continues reproducing the solution (24), but yields also a second solution, which lies far from the two-particle threshold and does not have the required nonrelativistic limit. Solution (24) remains therefore the only stable nonrelativistic solution of Eq. To interpret correctly its role, one should try to connect it more explicitly to the original coupling constant g of the higher-energy theory. In the latter theory, the lowest bound state approaches the two-particle threshold when g approaches its critical value g cr from above. Therefore, negative values of h correspond to values of g greater than g cr . Large and negative values of h would correspond to the approach g → g cr , at which value h would have a singularity. Vanishing of h with negative values would correspond to the increasing of g up to infinity. When g is lower than g cr , h should change sign and the system would enter in a phase characterized by the absence of bound states. Finally, the vanishing of h with positive values would correspond to the approach to the perturbative regime. A schematic representation of these correspondences is shown in Fig. 3 [29].
The exact relationship between h and g not being available at present, it would be useful to have an approximate or an empirical relation which qualitatively reproduces its main properties and allows for an easy understanding of the physical situation that is considered. For this, we propose the following formula between h and g 2 [Eq. (5)], the latter having a more universal meaning than where b is an empirical parameter, whose approximate value has been determined by numerical tests. By expanding in (26) g 2 around g 2 cr , one obtains from (24) the further approximation of e tm near the two-meson threshold: The quadratic behavior of e tm with respect to the departure of g 2 from g 2 cr is typical of short-range potentials and can be verified on soluble models.
As a side remark, let us notice that the extension of the formula (26) to the relativistic domain necessitates more elaborate comparisons. The reason is that in the latter domain one has no longer a single energy unit and the corresponding generalization is not straightforward. Furthermore, in the higher-energy theory, the Bethe-Salpeter equation itself has difficulties to describe correctly the first relativistic corrections in the ladder approximation with covariant propagators [82]; one is obliged to use either the instantaneous approximation, or other quasipotential-type approaches, which reduce the Bethe-Salpeter equation to a three-dimensional equation. We shall be content, in the present work, to stick to the nonrelativistic domain, where many experimental data still require a detailed investigation.
One can also obtain, from the expression (20) of the scattering amplitude, the coupling constant of the bound state to the constituent mesons, appearing in the residue of the bound state pole.
Designating by M 1 and M 2 the two constituent mesons, the scattering amplitude has the following behavior near the pole position of the bound state: where g T M 1 M 2 represents the dimensionless coupling constant of the tetraquark to the two mesons, defined with the accompanying mass factor (m 1 + m 2 ). Expanding in (20) J r (s) around s tm and using (21), one obtains giving One notices that the coupling constant decreases when the bound state approaches the two-particle threshold.
Let us finally emphasize, as is evident from the previous results and, in particular, from the structure of T [Eq. (10)], that the low-energy theory can reproduce, in the present scalar theory, only the ground state of the higher-energy theory.

IV. SCATTERING LENGTH AND EFFECTIVE RANGE
The expression of the scattering amplitude (10) where we have defined On the other hand, T can be expressed in terms of the S-wave phase shift δ 0 as The factor cot δ 0 (k) is itself expressed through a low-energy expansion in terms of the scattering length a and the effective range r e [83]: yielding the identifications (The relativistic correction coming from the expansion of √ s in the numerator of T , Eq. (33), has been neglected.) One finds that a is proportional to −h and, therefore, has the same type of behavior as −h in terms of the coupling constant squared g 2 (Fig. 3). (This has also been shown in [78].) According to whether g 2 is greater or smaller than g 2 cr , a is positive or negative. On the other hand, the parameter d [Eq. (32)] is positive and in general, for physical applications, smaller than 1; the effective range r e is then predicted positive and small.
Of particular interest is the case of resonances, which appear as bumps in the cross section above the two-particle threshold. They correspond to complex poles of the scattering amplitude, lying below the cut of the real axis. To check the possible presence of complex poles, we go back to the expression (16) of Q(s) in the complex plane, which can also be rewritten in the following form: Using a definition of the type of (22), one has in approximate form, neglecting quadratic terms in E, The first equation shows that the complex variable E, considered as a vector, is parallel to (s − (m 1 + m 2 ) 2 ) in the s-plane and therefore one can transpose to E the complex-plane analysis, with a right-hand cut starting at E = 0. Making expansions in √ e and retaining terms up to order e, one obtains (α defined in Eq. (24).) J(s) [Eq. (19)] then takes the form where d has been defined in (32). In the nonrelativistic domain, the second term is negligible in front of the first and the resonance equation takes the form which yields a purely imaginary solution for √ e R and gives back the bound state solution (24).
Formal resonance solutions can be obtained outside the nonrelativistic domain for small negative , by including also the second term of the right-hand-side of (39); however, such values of h correspond, according to the correspondence (26), to the strong-coupling limit of g 2 , which goes beyond the validity of the present nonrelativistic approximation.
The previous results mean that the present model does not produce resonances in the vicinity of the two-particle threshold. Resonances can be produced when there are derivative-type couplings [84][85][86], which we have discarded in the present approach. One can also refer to the rectangular well model of [79], where resonances in the S-wave are generally produced far from the real axis.

V. COMPACT TETRAQUARKS
We have considered in the previous sections, in an effective theory approach, the bound state formation problem of a molecular state, or a hadronic molecule, which we also called a moleculartype tetraquark, from two mesons, interacting by short-range Yukawa-type forces, approximated in the effective theory by a contact-type interaction. We have noticed that the smallness of the binding energy is sharpened when the coupling constant of the lower-energy theory takes large negative values, corresponding in the higher-energy theory to the proximity of the coupling constant to the critical value, below which no bound state exists.
The latter mechanism is not the only one that may produce tetraquark states. Another mechanism, based on the direct internal interaction of four-quark systems (more precisely, made of two quarks and two antiquarks) by means of the confining forces, might also produce bound states, in analogy to what happens with the formation of ordinary hadrons [34,35]. Because of the strong nature of the confining forces, one expects that such bound states would have more compact sizes than the molecular-type bound states and are distinguished from the latter in the literature under the terms of "compact tetraquarks".
However, the formation of compact tetraquarks as definite stable bound states (with respect to the strong interactions) remains a matter of debate. This is related to the "cluster reducibility" problem, in the sense that the multiquark operators that create tetraquarks are reducible to a combination of mesonic clusters, and hence the compact tetraquark state would rapidly dislocate into them and would be transformed into a molecular-type object [16,46,87,88].
Another argument which is advocated in favor of the molecular scheme is the proximity of many of the observed tetraquark candidate states to two-meson thresholds. In the molecolar scheme, the two-meson threshold is a natural reference of energy levels. In the compact tetraquark scheme, the elementary confining forces do not refer to meson states and hence, at first sight, no natural justification is proposed for the appearance of tetraquark states near two-meson thresholds.
We shall analyze, in the present section, these questions with the aid of the effective field theory approach, adopted in Secs. III and IV.

A. Compositeness
The comparison of the molecular and compact schemes is reminiscent of a general problem, already raised in the past in the case of the deuteron state, denoted under the term of "composite-ness" [36]. The binding energy of the deuteron, referred to the proton-neutron threshold, is very small as compared to the mass scale involved in the strong interaction dynamics of the nucleons.
One is inclined to consider the deuteron as a loosely bound composite object, or a molecule, made of a neutron and a proton. On the other hand, there might still exist some probability, that should be quantified, that it might be an elementary particle, or a compact object. Weinberg has shown that this question can receive, in the nonrelativistic limit, a precise and model-independent answer, by relating the latter probability to observable quantities, represented by the scattering length and the effective range of the neutron-proton S-wave isospin-0 scattering amplitude [36]. Designating by Z the probability of finding the deuteron in an elementary, or compact, state, Weinberg has found the following relations for the scattering length a and the effective range r e : where R is the deuteron radius, m r the reduced mass of the proton-neutron system [Eq. (3)], E d (≡ 2m r e d ) the deuteron nonrelativistic energy (opposite of its binding energy) and m π the pion mass; O(m −1 π ) represents the scale of the hadronic corrections that are negligible in front of R. On the other hand, introducing the dimensionless deuteron-neutron-proton coupling constant g dnp , accompanied by the factor (m n + m p ), one has the relationship One notices that the effective range is the most sensitive quantity to Z, which, in case Z = 0, is manifested by a sizeable negative value. Using Eqs. (41), one can also express the compositeness factor in a combined form with respect to a and r e : In the case of the deuteron, the experimental data about a and r e rule out a nonzero value of Z and confirm its composite nature [36].
Equations (41) and (42) can also be used, with appropriate relabeling of the parameters, to check the consistency of the results obtained in Secs. III and IV. Equations (35) show that r e has a small negligible value (as compared to R = 1/(2m r √ −e tm )), which could be interpreted as representing the higher-order hadronic corrections. Thus, with respect to the second of Eqs. (41), the main value of r e is 0, which entails that Z = 0, in accordance with the molecular nature of the bound state. Also, the comparison of Eq. (30) with (42) confirms the latter conclusion.
For later reference, using notations adapted to the tetraquark problem, we display here the expression of the scattering amplitude obtained in Weinberg's analysis: where g T M 1 M 2 has been defined in (28) and e and e t in (23). The main assumption underlying this result concerns the absence of zeros in T , at least in the vicinity of the bound state [89].

B. Compact bound states
We consider, in this subsection, the case of possibly existing compact tetraquarks. We shall not enter, for the analysis of the problem, into the details of the mechanism that produces such states, but merely shall assume their existence. If experimental data were sufficiently precise concerning an observed tetraquark candidate, providing us with its coupling amplitude to the nearby twomeson states, as well as the scattering length and the effective range of the related two-meson elastic scattering amplitude, then, for a nonrelativistic state, Eqs. (41) and (42) would allow us to reach a conclusion about the internal structure of the tetraquark. In the absence of high precision data, we proceed by successive steps. In first approximation, we assume that the compact tetraquark is a pointlike object in comparison to a loosely bound molecular state. At this stage, we assume that the mass of the tetraquark has been evaluated by the sole mechanism of the confining forces, from which clustering forces or effects have been removed. (The small volume or pointlike approximations of the diquark system satisfy this requirement.) Furthermore, we assume that the mass of the bound state under study, designated by m tc1 , where the labels tc refer to the compact tetraquark, is rather close to the nearest two-meson threshold mass (m 1 + m 2 ) and, therefore, a nonrelativistic energy of the bound state, E tc1 , can be defined by means of the equation E tc1 remaining a small quantity with respect to the two-meson reduced mass. However, we do not assume that the bound state is shallow. Considering the case of the deuteron as an example of a shallow bound state, whose binding energy is of the order of 2 MeV, E tc1 might have values of the order of 20-30 MeV.
Another point worth emphasizing is that, in general, when one has many different quark flavors inside the tetraquark state, the latter has two different two-meson clusters [29,90,91] and one should, in that case, use a coupled-channel formalism. However, in order to display in a clearer way the main qualitative aspects of the problem, we stick here to a single-channel formalism (which describes an exact situation when two quarks, or two antiquarks, have the same flavor).
Because of the existence of internal two-meson clusters inside the tetraquarks, the compact tetraquark has necessarily a coupling to the two mesons M 1 and M 2 . The corresponding (dimensionless) coupling constant is designated by g ′ , factored by the mass term (m 1 + m 2 ). However, the latter coupling generates, through meson loops, radiative corrections inside the tetraquark propagator, thus modifying the parameters of the bare propagator. They are graphically represented in Designating by m tc0 the bare tetraquark mass, the full tetraquark propagator becomes where s stands for p 2 and J is the same loop function as the one met in Eqs. (11)- (19). The divergence of J is now absorbed by the bare mass term, yielding the renormalized mass m tc1 : The renormalized tetraquark propagator is now Notice that g ′ does not undergo any renormalization. The mass term m tc1 does not yet represent the physical mass of the tetraquark. The latter is determined from the pole position of the propagator.
Sticking to the nonrelativistic limit, one can use for s and m tc1 expansions of the types of (22) and (45). Retaining in J r (s) the dominant contribution, one finds for the nonrelativistic energy of the tetraquark, designated by E tc , the equation −e tc + e tc1 + g ′2 16π −e tc = 0, whose solution is (≡ −E tc1 /(2m r )) has been fixed at 0.01.
(This expression could also be obtained directly from (49) by neglecting in it e tc in front of e tc1 .) Because of the nonlinear relationship between e tc and e tc1 , there appears a strong decrease in e tc .
One has Considering for example −E tc1 ≃ 20 MeV and 2m r ≃ 2 GeV, one has −E tc1 /(2m r ) ≃ 0.01; values of g ′2 /(16π) of the order of or greater than 0.5 produce ratios E tc /E tc1 ≤ 0.04, or equivalently −E tc ≤ 0.8 MeV. These values of g ′2 /(16π) are not exceptional and we may conclude that we are not in the presence of a fine tuning effect. There is a nonnegligible probability, in many physical cases, to meet such a situation.
The contribution of the tetraquark state, in the s-channel, to the two-meson elastic scattering amplitude is obtained by inserting the tetraquark propagator between two tetraquark-two-meson couplings, as shown in Fig. 7. One finds Proceeding as in the molecular case [Eq. (28)], one can obtain the physical coupling constant g T M 1 M 2 of the tetraquark to the two mesons: Comparing this expression with (42), one obtains Z: which, after taking into account Eq. (49), can also be expressed as One notices, in particular from (55), that Z decreases in the strong-coupling limit and takes small values. With the numerical example considered above, one has Z ≤ 0.075. On the other hand, in the same limit, the physical coupling constant g T M 1 M 2 also decreases, like √ −e tc . The variation of Z as a function of g ′2 /(16π), taking into account (50), is represented in Fig. 8. One notices the rapid decrease of Z.
It is also possible to calculate from (53) , (d is defined in (32).) Eliminating g ′2 /(16π) and e tc1 in favor of Z and e tc , one obtains the expected (The second term in r e , which is small, has been neglected.) It is of importance to have a precise interpretation of the values of Z. When g ′ = 0, Z = 1 and the tetraquark is completely decoupled from the two mesons. This would mean that either its mass scale or quark content are different from those of the two mesons. In the opposite case, corresponding to the strong-coupling limit, Z decreases and approaches the value 0. This does not mean, however, that the tetraquark's nature becomes molecular. Nowhere, in the present model, did we consider direct interactions between mesons; the existence of the bound state is entirely due to the confining forces that are responsible of its compact nature. The value of (1 − Z) simply reflects the strength of the (bare, but finite) coupling of the tetraquark to the two meson clusters.
As g ′ grows, the latter, through the radiative corrections, plays, in the internal structure of the tetraquark, an increasingly determinant role, leading to a strong decrease of the binding energy and to a corresponding increase of the radius R [Eq. (58)]. The tetraquark, though of compact nature, is gradually deformed into a molecular-type state.

C. Resonances
A resonance may occur when the renormalized (real) mass of the compact tetraquark [Eq. (47)] lies above the two-meson threshold. Its nonrelativistic energy, defined in (45), is now positive.
The position of the complex mass of the possibly existing resonance can be searched for with the same method and the same approximations as in the bound state case. Designating by E T R the nonrelativistic (complex) energy of the tetraquark resonance, the equivalent of Eq. (49) is Its solutions are The imaginary part of √ e T R comes out negative, which means that the two solutions lie in the second Riemann sheet. The expression of e T R is The two solutions are complex conjugate to each other, the resonance corresponding to the negative imaginary part. One verifies that the modulus of e T R is equal to e tc1 : The condition of the positivity of the real part of e T R (resonance above the threshold) requires This means that one is rather in the weak-coupling regime. As the coupling constant increases, the real part of the resonance energy approaches the threshold, the upper bound of (63) corresponding to its merging with the threshold. On the other hand, the imaginary part remains always different from zero; at threshold, it is only the latter that survives.
The scattering amplitude, due to the resonance, is We have represented, in Fig. 9, the shape of |T | 2 in the vicinity of the threshold, for real E and for several values of g ′2 /(16π) in its allowed domain.
Expanding T around E T R [Eq. (61)] and using for the resonance the plus sign in (60), one with Comparing (66) with (42), the former being considered as a formal extension of the latter to the resonance region, one obtains the compositeness coefficient We notice that, due to the inequality (63), which guarantees the occurrence of the resonance above the threshold, (1 − Z) is real, positive and smaller than 1. This allows us to continue giving it a probabilistic interpretation.
Writing the (reduced) energy of the resonance in the form where e T Rr is the real part of the (reduced) energy and γ the reduced dimensionless width of the resonance, defined as (see also (23)) Γ being the dimensionful width, one has from (61) which allows us to express (1 − Z) in terms of observable quantities 2 : Of particular interest is the case of narrow resonances, characterized by the inequality Γ ≪ 2E T Rr , which entails It is evident from this result that the narrow resonance case favors values of Z close to 1, that is, dominance of the compact nature of the tetraquark. In the opposite case, when the inequality (63) is saturated, one has vanishing of e T Rr and merging of the resonance with the threshold, with Z = 0; the trace of the compact origin of the resonance is then completely lost.
The expressions of the scattering length and the effective range are the same as in Eqs. (57), except that e tc1 is now positive: (The second term in r e , which is small, has been neglected.) One notices that now, as compared to the bound state case (57), both the scattering length and the effective range are negative. In terms of a and r e , the compositeness coefficient (67) takes the form a relation also obtained in [84]. It can be compared with Eq. (43), valid in the bound state case.
In the narrow resonance case, one has the simplified expressions Let us also comment on the case of the strong-coupling limit, where, in particular, the inequality (63) is no longer satisfied. We have to distinguish here two cases. The first corresponds to the domain 2e tc1 < (g ′2 /(16π)) 2 < 4e tc1 , in which case (1 − Z) > 1 and the real part of the energy becomes negative [Eq. (70)]. This is a sign of the instability of the initial system under the influence of the two-meson clusters and might signify the disappearance of the compact tetraquark from the spectrum. The second corresponds to the domain 4e tc1 < (g ′2 /(16π)) 2 . In that case, √ e T R [Eq. (60)] becomes imaginary and e T R becomes real and negative, falling back in the boundstate regime. Considering then the new value of e T R as a starting value (equal to a new e tc1 ), one continues remaining in the bound-state regime for any value of the coupling constant, as we have seen in Eq. (50) and Fig. 6. Therefore, genuine resonances may occur, in the present model, only in the weak-coupling regime, satisfying the inequality (63).
Let us notice, as a final remark, that due to the fact that the internal structure of the compact elementary particle, assumed here as being a tetraquark, was not specified, one is also entitled to apply the previous approach, when the flavor and other quantum numbers are compatible, to ordinary mesons. In that case, its field theoretic basis is even more robust. This is supported by the large-N c limit of QCD [29,[92][93][94]. In that limit, the spectrum of the theory is composed of free mesons, which are made of one quark-antiquark pair. They do not have any internal other meson clusters. Therefore, their elementary nature with respect to the other mesons is well justified.
The couplings to the other mesons appear only at nonleading order in N c , putting the clustering phenomenon at a perturbative level. This is in contrast to the multiquark case, where the clustering occurs already at leading order in N c and becomes even stronger in the large-N c limit [29]. (Cf. also Sec. VII.) Detailed investigations about the compositeness criterion and its applicability to various tetraquark candidates, as well as to ordinary hadrons, can be found in References [84,[95][96][97][98][99][100][101][102][103][104][105][106][107][108][109].

VI. PRESENCE OF MESON-MESON INTERACTIONS
In the model considered in Sec. V, where the influence of the coupling of a compact tetraquark to two mesons was studied, the presence of meson-meson interactions as a background effect was not taken into account. This had the advantage of exhibiting in a clearer way the role of the aforementioned coupling on the observable properties of the tetraquark state, in particular its gradual deformation, in the strong coupling limit, towards a molecular-type object. To complete the previous study, we include, in this section, the effect of the meson-meson interaction into the dynamical process.

A. The meson-meson scattering amplitude
The meson-meson interaction, in the present effective field theory description, was considered in Secs. III and IV. When a compact tetraquark state is present, its effect is first manifested through a vertex renormalization related to the coupling constant g ′ . This is graphically represented in Fig. 10. The coupling constants at the vertices are indicated. g ′ is accompanied by the mass factor (m 1 + m 2 ).
The corresponding vertex function, denoted Γ  (14)) We have seen that the divergence contained in J(s) [Eq. (12)] is absorbed by a renormalization of the coupling constant h, according to (13) or (14). Here, the occurrence of the same divergence necessitates a similar renormalization of g ′ , involving, however, also h: Equivalently, one has the following renormalizations: These ensure the renormalization of Γ A second effect of the meson-meson interactions is manifested through the radiative corrections of meson-one-loop diagrams, occurring in the tetraquark propagator (cf. Fig. 5). Each loop receives radiative corrections, as represented in Fig. 11.
The chain of such diagrams can be summed and yields the full one-loop contribution, which we designate by Γ T T (s). One obtains FIG. 11. Chain of meson-one-loop radiative corrections to the one-loop diagram of the tetraquark propagator.
The coupling constants at the vertices are indicated. One finds for the full propagator Using the renormalizations given in (78) and (79), and after taking the limit iJ div → ∞, one finds One notices that after the renormalizations of the coupling constants h and g ′ have been realized, the radiative corrections of the tetraquark propagator are now finite. This is in contrast to the case where the meson-meson interactions had been ignored (h = 0) and the divergence of the radiative corrections had been absorbed by the mass renormalization, while the coupling constant g ′ had remained finite (cf. (46) and (47)). Nevertheless, the radiative corrections in (83) contain a singularity in h r . When h r → 0, one recovers the divergence that exists in the aforementioned case.
To remedy this defect, one has to subtract that singularity from the global radiative corrections and associate it with a renormalization of the bare mass m tc0 . Designating by m tc1 the renormalized bare mass, one has As long as h r is nonzero, this mass renormalization is finite. When h r → 0, one recovers the situation of Eq. (47).
The full tetraquark propagator takes now the following form: .
The meson-meson scattering amplitude is obtained by inserting the tetraquark propagator inside two vertices of the type of (77), which are finite [Eq. (80)], and by adding the contribution generated by the contact interaction (8) [Fig. 2 and Eqs. (10) and (13)]. This is represented in Fig. 13. One obtains The above expression could also have been obtained by starting from the integral equation iT = K + iKJT , with the kernel K given by After the renormalizations of the coupling constants and of the mass m tc0 are done, one finds (86).

B. Bound states
The singularities of the scattering amplitude (86) are the same as those of the tetraquark propagator (85). The separate molecular-type singularity, present in T (h) , has been cancelled by a similar singularity resulting from the radiative corrections in T (g ′ ,h) .
We first focus on the bound state problem. The tetraquark mass, m t , is given by the equation where s t = m 2 t . Sticking to the nonrelativistic limit and using definitions similar to (22) and (23), and Eq. (88) reduces to (omitting henceforth the renormalization label r from the coupling constants) where α is defined in (24) 3 . We notice that in case g ′ = 0 (absence of tetraquark-two-meson coupling), the equation splits into two independent equations yielding the molecular-type solution (24) (for h < 0), on the one hand, and the bare compact tetraquark mass (84), on the other. In case h = 0 (absence of molecular-type forces), the equation reduces to that of the compact tetraquark case (49). Equation (90) cannot be solved analytically, but accurate analytic approximate solutions can be found for it. One has to distinguish two cases, according to the sign of h.
We first consider the case h > 0. According to the empirical relationship (26) and Fig. 3, this case corresponds to subcritical values of three-meson coupling constants, for which no genuine molecular-type bound states can exist. Even though h may take large values, approaching +∞, one may characterize this domain as globally representing a weak-coupling regime. One therefore expects that the tetraquark bound state originates entirely from a compact configuration, the molecular forces mainly introducing deformations.
A leading approximate solution to (90) is which generalizes (50). This expression, together with its next-to-leading term, which is presented Its next-to-leading term is given in the Appendix . The analytic approximation (A.2) reproduces the behavior of the exact solution with an error of less than one per mil. The behavior of −e t, (1) with respect to variations of g ′2 /(16π), for fixed h, is again similar to that of Fig. 6; however, now, it starts, when g ′2 /(16π) = 0, from √ −e tm [Eq. (24)], instead of √ −e tc1 (Fig. 15). In the limit Expanding the scattering amplitude (86) around the bound state pole (cf. (28)), one obtains the physical coupling constant g T M 1 M 2 , expressed in two equivalent ways, using (90): Comparison of these expressions with (42) yields Z: which is manifestly a positive quantity bounded by 1.
The value of Z and its behavior under variations of the coupling constants depend on the specific bound states that we have met above. When h > 0, we have one bound state with the squareroot of the binding energy having the approximate expressions (90) and (A.1). For this case, the qualitative features of Z are more transparent in the second expression of (95). When g ′ → 0, for fixed h, √ −e t → √ −e tc1 (cf. Fig. 14); then Z → 1. When g ′ increases, √ −e t decreases rapidly and tends to zero. The behavior of Z remains very similar to that of Fig. 8, with the difference that the presence of h slows down the decrease of Z. We display, in Table I configuration in all the interval of variation of g ′ . The main influence of the latter is reflected in the continuous increase of the binding energy (Fig. 15). We display, in Table III,   To determine the scattering length and the effective range, one expresses the scattering amplitude T in the scattering region E > 0: from which one deduces, with the aid of (33) and (34), where d is defined in (32) and α in (24). The relationship of a and r e with Z [Eq. (95)] is not, however, as simple as in the previous cases [Eqs. (41) and (73) if they are measured experimentally or on the lattice, allow the calculation of Z. In particular, the presence of the coupling g ′ between the tetraquark and the two mesons provides r e in general with a negative value, whatever the sign of h is. We notice that for h < 0, a is positive, as it receives contributions from two independent bound states. For h > 0, a is positive for αh < −g ′2 /e tc1 and negative for αh > −g ′2 /e tc1 . For αh = −g ′2 /e tc1 , a vanishes and simultaneously r e tends to −∞.

C. Resonances
Resonances may occur when the renormalized (real) compact tetraquark mass, given by (84), lies above the two-meson threshold. In that case, the nonrelativistic energy E tc1 of (89) is positive.
The poles of the scattering amplitude (86) When h = 0, this equation has to be solved numerically. Introducing the definitions where e T Rr and e T Ri are the real and imaginary parts of e T R , respectively, one can study the behavior of e T R under variations of g ′2 and h. The physical conditions to be imposed on the solutions are v < 0, e T Ri < 0 and e T Rr > 0 (the resonance lies above the two-meson threshold). In the case h = 0, the upper bound (63) had been found for g ′2 /(16π).
When h = 0, one has to distinguish two main domains of h (cf. Fig. 3): (i) h > 0 and small, corresponding to the weak-coupling regime of the meson-meson interaction; (ii) |h| large, corresponding to the strong-coupling regime, with the possibility of existence of a genuine molecular-type bound state (for h < 0). In the first case, the qualitative behavior of the solutions is naturally close to that of the case with h = 0, studied in Sec. V, the only changes being small quantitative ones. Thus, for h/(16π) ≤ 0.2/ √ −e tc1 , a conservative upper bound for When g ′2 increases from 0 to its upper bound, the resonance approaches the two-meson threshold with its imaginary part remaining finite, but relatively decreasing. In the second case, for the same variation of g ′2 , the resonance approches the threshold for h < 0 and moves away from it for h > 0.
There are also regions of h for which large values of g ′2 are possible, however they are not obtained by continuous variations of all interactions.
For completeness, we display here the equation satisfied by the imaginary part of √ e T R : |e T R | being the modulus of e T R , which generally favors negative values of v (when h < 0 or when h > 0 but small).
Expanding T [Eq. (86)] around E T R , as in (65), one obtains the expression of the tetraquarktwo-meson coupling constant squared: Generally, the multiplicative factor accompanying 32π √ e T R is identified with the compositeness (42), (67), (95)]. Here, however, the corresponding coefficient is complex when h = 0 and, therefore, a probabilistic interpretation of it is no longer possible (cf. also [100,101] situation has also been found in Sec. V C and has been interpreted as a sign of the instability of the initial system, leading to the disappearance of the compact tetraquark from the spectrum. On the other hand, large values of g ′2 generally send back the resonance to the bound state domain. We therefore adopt the modulus prescription, with the restriction that g ′2 respects its upper bound: The general properties of Z are then similar to those met in the case h = 0. When the resonance approaches the two-meson threshold, Z tends to 0, while when the resonance stays in the vicinity of its primary position, this mainly corresponding to small values of g ′2 , Z remains close to 1. On the other hand, large and negative values of h have the tendency to push the resonance towards the two-meson threshold, while large and positive values of h repel the resonance from the threshold. It is worthwile recalling that in the case h < 0, the spectrum also contains a molecular-type bound state, whose binding energy has been approximately estimated, in the two-bound-state case, by means of Eq. (92), where e tc1 < 0. The same formula could also be used for the evaluation of the binding energy in the resonance case, where now e tc1 > 0: We find that the bound state exists only when −αh/(16π) is greater than g ′2 /(16πe tc1 ); the condition of the vicinity of the bound state to the two-meson threshold actually requires much larger values. For instance, with e tc1 = 0.01 and g ′2 /(16π) = 0.1, one would need −αh/16π > 10. Values of −αh/16π of the order of 30, which were frequently considered throughout the present work, would then produce a binding energy −e t0,(1) of the order of 0.0025.
The expressions of the scattering length and the effective range are the same as in Eqs. (97), with the only difference that e tc1 is now positive. The contribution of the term proportional to g ′2 in a is now negative and in the case h < 0 the same competition as in the case of the bound-state binding energy (103) arises between the contributions of h and g ′2 . In case −αh/(16π) − g ′2 /(16πe tc1 ) is negative, the bound state disappears and the scattering length becomes negative. While the effective range is insensitive to the sign of that combination, it nevertheless has a singularity when the latter vanishes. Therefore, the latter can act on them only as perturbations, not affecting their compact structure.
This is not the case of tetraquarks, because of the existence of internal mesonic clusters. The compact structure of primarily existing tetraquarks is not protected by the large-N c limit. This is due to the fact that the interaction forces acting for the formation of mesonic clusters are (N c − 1) times larger than the forces forming diquark compact objects [20,29] (assuming that the confining forces have the same color structure as one-gluon exchange terms). This means that the primary coupling constant g ′ that connects the compact tetraquark to the two mesons should behave at large N c like N 1/2 c , assuming that the compact tetraquark state and energy are of order N 0 c . For the squared quantity, one has Concerning the four-meson contact-type coupling constant h, we have emphasized in Sec. III that it is not an elementary coupling constant and should rather be related to the three-meson coupling constants of the higher-energy theory. The latter coupling constants generically behave, at large N c , like N −1/2 c [29,93,94,112] and vanish in that limit. Compared to the critical coupling constant, introduced in (7), they lie, in the subcritical domain. Meson-meson interactions cannot, therefore, produce on their own bound states in the large-N c limit. Going back to the empirical formula (26), we deduce that h is positive and lies in its perturbative domain, behaving like g 2 : The meson-meson interaction has, therefore, only a subleading effect with respect to the direct compact-tetraquark-two-meson interaction. From Eqs. (91) and (52) From (95) and (56), one also finds the behavior of the elementariness coefficient: At large N c , the compact tetraquark is thus transformed into a shallow, molecular-type, bound state.
From Eqs. (94) and (54), one deduces the behavior of the physical coupling constant g T M 1 M 2 : This result is in accordance with other general estimates done in the large-N c limit, which predict the vanishing of the coupling constant in that limit [29,90,91,[114][115][116][117][118]. The behavior (108) should only be considered as a generic one. The power of the decrease may slightly change according to the detailed flavor content of the tetraquark state or other more refined analyses.
Concerning the resonance states, we had found that they show up only in the weak-coupling regime. However, the large-N c limit imposes the strong-coupling regime. Therefore, in that limit, resonances should disappear from the spectrum, at least from the neighberhood of the two-meson threshold.
In conclusion, in the large-N c landscape, tetraquarks may survive only in the form of shallow, molecular-type, bound states, which are relics of primarily created compact states.

VIII. CONCLUSION
The compact tetraquark scheme, considered in its simplest version, where diquarks and antidiquarks are separately gathered in very small volumes, can be considered as a starting point for the analysis of the tetraquark properties. In this situation, because of the dominance of the attractive confining forces, one usually finds confined bound states, in parallel to the case of ordinary hadrons. However, the very existence of underlying meson-clustering interactions in the general system forces the initial compact state to evolve towards a more complicated structure, where now molecular-type configurations are also present. In an effective field theory approach, where the compact tetraquark is represented as an elementary particle, this evolution was studied, within a scalar interaction framework, by means of the primary compact-tetraquark-two-meson coupling constant. Another quantity, more related to physical observables, is the elementariness coefficient Z, which varies between 1, corresponding to the pure elementary case, and 0, corresponding to the completely composite case. Stronger is the primary coupling constant, smaller is the value of Z. In the strong-coupling limit, the system tends to a dominant molecular configuration, characterized by a shallow structure. In the case of resonances, only the weak-coupling regime provides a stable framework for their existence. For higher values of the coupling constant, either the resonance disappears from the spectrum, or is resent to the bound state domain.
The consideration of the large-N c limit of QCD provides an additional support to the dominance of the strong-coupling regime of the coupling constant, with all its consequences.
Many of the tetraquark candidates, whose elementariness has been evaluated in the literature from experimental data, have led to values of Z which are neither 1, nor 0. This is a clear indication of the mixture of configurations that results from the evolution described above. Nonzero values of Z, even small ones, are indicative of the existence of a primary compact state. Shallowness of bound states, with small values of Z, may receive a natural explanation as resulting from the strong-coupling limit of the interaction compact-tetraquark-meson-clusters, also supported by the large-N c limit of QCD.
The analysis undertaken in the present work was based on the simplest qualitative approach, considering scalar interactions, ignoring spin degrees of freedom and details of quark flavors, and using nonrelativistic limit and single-channel formalism for the meson clusters. A more general and refined quantitative analysis, concerning definite candidates, should include the missing ingredients.