Coupling hadron-hadron thresholds within a chiral quark model approach

Heavy hadron spectroscopy was well understood within the naive quark model until the end of the past century. However, in 2003, the $X(3872)$ was discovered, with puzzling properties difficult to understand in the simple naive quark model picture. This state made clear that excited states of heavy mesons should be coupled to two-meson states in order to understand not only the masses but, in some cases, unexpected decay properties. In this work we will review how the naive quark model can be complemented with the coupling to two hadron thresholds. This program has been already applied to the heavy meson spectrum with the chiral quark model and we show some examples where thresholds are of special relevance.


Introduction
Heavy hadron spectroscopy started in November 1974, when Brookhaven National Laboratory announced the discovery of a new particle called J [1] and, at the same time, the Stanford Linear Accelerator reported the existence of another new particle, called ψ [2]. Very soon both particles were seen as the same state, which we know now as the J/ψ state of the charmonium spectrum. This state was understood as a qq bound state in the naive quark model and its discovery was the confirmation of the existence of the charm quark, that was predicted by the GIM mechanism [3] few years before.
Afterwards, heavy meson spectroscopy developed very fast during the following years. The Υ(1S), the first state with bottom quarks discovered, was found at Fermilab [4] in 1977. Already at 1980 there were 11 new mesons included in the Particle Data Group table [5] on these energy ranges, but that rate decreased to only 15 new states added during the period 1980 to 2003 [6]. However, during these last 17 years 35 new states have been added [7], considering only unflavored mesons.
In the case of the baryon spectra, the first evidences of charmed baryons came six months after the discovery of the J/ψ, in 1975 [8], but in 1980 only the Λ c and Σ c baryons where included in the heavy baryon spectrum of the PDG [5]. In 2003, only 14 states were identified [6], while again during the last 17 years 33 new states have been included [7].
This impressive development of the heavy hadron spectrum has been possible thanks in part to the so called B-factories, like Belle and BaBar, which are electron-positron colliders tuned to the center of mass energy of the Υ(4S) that decays into two B mesons. Other facilities like BESIII, with a lower energy electron-positron collider, have contributed. Many impressive results have also been obtained and are underway at LHC by the LHCb, CMS and ATLAS Collaborations. The next generation SuperB-factory Belle II is running from 2018, and it is expected to give many important contribution to heavy hadron physics.
From a theoretical point of view, already in 1978, the Cornell model [9] was developed to understand heavy meson spectroscopy. This model is a non-relativistic approach for heavy quarks with interactions that are governed by SU(3) color gauge symmetry, with flavor only broken by the quark masses. The main pieces of the model are a Coulomb-like interaction, inspired by the one-gluon exchange, and a linear term, which describes the confining effect. It also took into account two important features than one expects from QCD, Heavy Flavor Symmetry (HFS) and Heavy Quark Spin Symmetry (HQSS), considering terms that are flavor and spin independent. The naive quark model from Cornell was fitted to the only 11 states that were known in 1978 and gave a quite good description of the charmonium and bottomonium spectrum at that time [10]. Not only that, the predictions were in quite good agreement with the experiments up to 2003, giving a prediction for 15 new states in the correct energy range. The Cornell potential has been related with the QCD static potential by Sumino [11] and more recently with NRQCD up to N 3 LO [12].
For the baryon spectra, also very soon, in 1979, quark models developed for the light sector were applied in the heavy quark sector [13]. Besides, Stanley and Robsen [14] extended the Cornell model to study heavy baryons. Many new states were predicted, but it took a long time to be seen on experiments.
Until 2003, the simple naive quark model picture was in quite good agreement with experiments. However, already in the original Cornell model [9], the coupling with two-meson thresholds was considered for excited states, although was found of no relevance for the states considered at that time.
The key event in 2003 for heavy hadron spectroscopy was the discovery of the X(3872) by the Belle Collaboration [15]. It was very soon confirmed by the CDF [16], D0 [17] and BaBar [18] Collaborations. It has some intriguing properties difficult to understand in the naive quark model picture but easily explained when coupled channel effects are included.
Nevertheless, one of the clearest indications that coupled-channel effects have to be considered are the famous pentaquarks measured by LHCb [19,20]. These states are unflavored baryons in the region of 4.5 GeV which rules out a three-light quark baryon interpretation. The only possible explanation is a pentaquark with three light quarks and a cc pair. Whether these states are compact pentaquarks or baryon-meson molecules is a matter of intense debate, although the closeness of these states to different meson-baryon thresholds is seen as a clear indication of the second possibility.
In this work we will make a brief review of coupled-channel effects in the framework of the quark model, in the same spirit as the original Cornell model. Thus, we will use naive quark-model states coupled to two-hadron channels, also built from naive quark-model states. The coupling between these two different sectors will be obtained using the microscopic 3 P 0 creation model. The paper is organized as follows: In Section 2 we will give a brief introduction to the naive chiral quark model and how to calculate the spectrum in this picture. In Section 3 we will give the basis of the 3 P 0 model and how to evaluate the transition amplitude. Section 4 will be devoted to present the formalism to couple one and two-hadron states. In Section 5 we will show a few examples in the meson spectrum where such effects are relevant and show some results in the quark model picture. We will end with some conclusions.

The Naive Chiral Quark Model
The quark model we use is a constituent quark model based on spontaneous chiral symmetry breaking [21]. It was first applied to the light-quark sector [22] and then extended to the heavy sector [23].
The main ingredients of the model are the following. The spontaneous chiral symmetry breaking generates two important effects in the light-quark sector. On one side, the light quarks acquire a dynamical mass that, at zero momentum, is of the order of 330 MeV for the u and d quarks, and 550 MeV for the s quark. This dynamical effect has been seen on the lattice [24] and it is the idea introduced phenomenologically in the constituent quark model. On the other side, it introduces the interaction of light quarks through the exchange of pseudo-Goldstone bosons. Another important non-perturbative effect is confinement, seen as the fact that hadrons are only seen in color singlets. We include it phenomenologically using a linear screened confinement interaction. This effect has also been observed in lattice QCD where, in quenched QCD, a linear rising of the energy of two static sources with increasing distance is clearly seen [25] and, in unquenched QCD, this string is broken [26] when there is enough energy to produce a quark-antiquark pair. Finally we introduce QCD perturbative effects through the one-gluon exchange interaction [27]. The model has been review in many works and all the details can be found in Refs. [28,29].
Once the model interaction is settled, in order to obtain the hadron wave functions one has to solve a non-relativistic bound equation for the two-body problem, in the case of mesons, as quark-antiquark pairs or the three body problem, for baryons, as three-quark states. There are many different approaches to solve these systems. The problem can be solved in coordinate space, in momentum space or using the Raileigh-Ritz variational principle using a certain basis function. Sometimes the potential involved poses special problems for some technique. For example, the use of a Coulomb-like potential makes the diagonal part of the potential in momentum space logarithmically divergent, which introduces a numerical problem. On the other hand, if we use a non-local interaction in coordinate space, the Schrödinger equation ends up being an integro-differential equation, which is also numerically more demanding. Another important fact is if coupled channels are considered or not. Momentum space calculations or variational calculations are very easily extended to such a case, while coordinate space calculations are not straight-forwardly extended.
Nonetheless, the method that can be usually used in any case is the variational method. This method is specially interesting for our purposes since we will be able to calculate transition amplitudes based on the 3 P 0 model with matrix elements of the basis functions used.
The main problem of the method is to find the appropriate basis functions. There are many different options for many different systems. However, the Gaussian Expansion Method (GEM) has been shown to be a very good approach in almost any case. The method was firstly proposed by Kamimura [30] and has been applied to many few body problems [31,32].
In the case of the two-body problem, the employed basis functions are a set of Gaussians multiplied by a solid-spherical harmonic in the relative coordinate, to take into account the correct behavior of the wave function at the origin ϕ nlm ( r) = N nl r l e −ν n r 2 Y lm (r) (1) The basis function is generated by taking several values for the Gaussian parameter ν n . In the GEM, this parameters are taken in geometrical progression as r n = r min a n−1 (4) where r min and r max are the minimum and maximum radius and N the number of Gaussians. Numerically, it is important that the parameter a is not too close to 1 so that the problem does not become singular. Then, one can use the expansion and ends up with the generalized eigenvalue problem ∑ nl (H n l ,nl − EN n l ,nl )c nl = 0 with H n l ,nl = ϕ n l m |H|ϕ nlm (8) N n l ,nl = ϕ n l m |ϕ nlm (9) In the case of the naive quark model one has to include the spin-flavor-color degrees of freedom so the total wave function is where χ S is the spin wave function, χ F is the flavor wave function and ξ c is the color singlet quark-antiquark wave function and where the orbital angular momentum l is coupled with the spin S to total angular momentum J.
As a matter of fact, the GEM is more interesting when we have more than two interacting particles. Usually, one chooses a set of coordinates that includes the center of mass, so one solves for the relative motion of the interacting particles. Then, one builds the most general wave function with the desired total quantum numbers. This is usually done considering an expansion in angular momentum of the different coordinates and one assumes that, for short-range interactions, only the lowest partial waves will be needed. To have a feeling of what it is needed, typical accurate three-body Fadeev calculations of Triton binding energy needs up to 38 of such partial waves. The GEM approaches the problem in a different way, it considers also the lowest partial waves although not only in one set of possible Jacobi coordinates, but in different sets. This has been shown to have a much faster convergence than the previous approach, which numerically is less demanding because radial wave functions in a less number of partial waves is needed. The drawback is that now different partial waves are not orthogonal, so they cannot be considered separately.
In the case of the three-body problem, there are three different sets of Jacobi coordinates given by where u i is the position vector of particle i and (ijk) is one of the 3 even permutation of (123). The orbital wave function is taken as where (−1) l i +L i gives the parity. If there are identical particles, some relations between different modes i may be needed. These relations can be easily obtained from the action of the permutation operator If we consider for example the Helium atom, with particle 1 being the nuclei and particles 2 and 3, the electrons then we can consider two different wave functions with definite symmetry against the P 23 operator where the symmetry to the exchange of the electrons is (−1) l 1 in the first case and ±(−1) l 23 for the second case. Then, the spin wave functions of the electrons have to be considered to solve the wave functions for parahelium (S = 0) and orthohelium (S = 1). In Tables 1 and 2 we show the ground state and first excited states for parahelium and orthohelium, respectively, compared to experimental data from the NIST database. Here we only include Coulomb interactions so we should expect deviations of the order of α 2 ∼ 10 −4 . As we can see, even with a long-range interaction as the Coulomb one, the GEM works very well.
If the three particles are identical then the wave function to be used is As explained above, we have to include the spin-flavor-color wave function where, again, the spin is coupled with the total orbital angular momentum L T to give a total angular momentum J and the color wave function corresponds to the color singlet.
The GEM is again very accurate and in Table 3 we give the result of ground state heavy baryons in the Bhaduri [34] model, calculated with the GEM and compared with a Fadeev calculation by B. Silvestre-Brac [35]. The GEM calculation only includes wave functions with l i = L i = 0. Besides, the matter radius square and the charge radius square defined by are given. From the results, one expects to have a good approximation to the solution of the three-body problem using the GEM. However, not all the states are below the lowest open threshold. If we again consider the Helium atom, the first ionization occurs at an energy of 24.587389011 eV [36], when the continuum of a ground state of a He + atom and a free electron starts. As the GEM takes boundary conditions for bound states, one can still find these bound states as shown in Table 4, although it is more difficult. These states can decay into a He + + e − so they are resonances and can be seen on scattering processes. In the case of hadrons there is a similar situation, however quarks can not abandon a hadron and a quark-antiquark pair is produced to generate two hadron states.

The 3 P 0 Model
The quark-pair creation model or 3 P 0 model is a microscopic model that allows to couple channels with different number of quarks. The name comes from the fact that a quark-antiquark pair is created with quantum numbers of the vacuum. It was first proposed by Micu [39] and, afterwards, Le Yaouanc et al applied it to the strong decays of mesons [40] and baryons [41]. These authors also evaluated strong decay partial widths of the three charmonium states ψ(3770), ψ(4040) and ψ(4415) within the same model [42,43].
The 3 P 0 model is usually formulated in terms of the Hamiltonian operator where the only parameter of the model is g s . The factor √ 3 is usually not included but in our case cancels the color factor in the meson sector, so g s has the usual definition. Table 2. First levels of orthohelium (S = 1) with orbital angular momentum L ≤ 2. The experimental data are taken from NIST, Ref. [33]. Results with the GEM method from this work.  Table 3. Results for the Bhaduri potential obtained with a Fadeev calculation (FD) [35] and the GEM method (this work). Masses are given in MeV, while matter radius square R 2 m and charge radius square R 2 c are in fm 2 . 1617 1517  It can be also formulated in terms of a transition operator given by where µ (ν) are the spin, flavor and color quantum numbers of the created quark (antiquark). The spin of the quark and antiquark is coupled to one. The Y lm ( p ) = p l Y lm (p) is the solid harmonic defined as a function of the spherical harmonic. It is common to give the transition operator in terms of the strength of the quark-antiquark pair creation from the vacuum γ as in Ref. [44]. The relation is given as being m the mass of the pair created, which is usually a light pair. We consider a processes where an initial hadron A decays into two final hadrons B and C. When we work in the center of mass system of the initial hadron we have P A = P cm = 0, with P A the momentum of the initial hadron and P cm the total momentum of the final hadrons. Then, the matrix element of the transition operator is written as The matrix element is taken between hadron states written in terms of quark degrees of freedom in second quantization. Meson and baryon states are written as where b † is a quark creation operator and d † an antiquark creation operator, N B is a factor in terms of the number of identical quarks in the baryon, and with the normalization convention for the wave functions, where the sum is for discrete degrees of freedom. With this states, two hadron states with correct quantum numbers are constructed.
With the transition amplitude, the widths for strong decays can be evaluated with k 0 the relative momentum of the two final hadrons. The transition amplitude is basically given in terms of the wave function of the naive quark model considered. There are many factors that are given in terms of the quark model symmetries, but also an important form factor is given in terms of the overlap of the initial and final hadron wave functions with the transition operators. The orbital part of the matrix element can be difficult to compute and this is the reason why the use of the GEM to solve the internal wave function of mesons and baryons is of special interest. The linearity of the operator allows that, using the expansion of the wave function in the GEM basis one can evaluate the matrix element in terms of matrix elements in the GEM basis as where the matrix element in the GEM basis B i C j |T|A j are evaluated analytically.
Once the naive quark model is fixed, the only unknown to determine the transition amplitude is the 3 P 0 strength parameter γ. In the case of meson decays, an scaled dependent parameter was considered in Ref. [45] as with γ 0 = 0.821 ± 0.02 and µ 0 = (49.84 ± 2.58) MeV. The scale is taken as the reduced mass of the quarks on the initial meson. The parameters were fixed to the strong width of a few open-charm, charmonium and bottomonium states and, then, applied to many different states on these sectors. However it is interesting to notice that this scale dependence was able to predict the strong decay widths of open-bottom mesons correctly without including this sector on the fit.

The unquenched quark model
In the previous section we showed how one and two-hadron states are connected and give rise to the strong decay widths. However, the same transition amplitude has as a consequence: the one-hadron and two-hadron states connected gets mixed. As the origin is the strong force, this mixing can be sizable.
For this reason, in some cases, it is important to consider the effect and, for this purpose, we consider the physical state as where |ψ α are naive quark model one-hadron states with α quantum numbers and |ψ H1 ψ H2 β are two-hadron state H 1 H 2 with β quantum numbers and relative momentum P. We define now the transition amplitude If we impose the Schrödinger equation and solve for the one-hadron amplitudes we find where χ β (P) is given by the equation in the two-hadron sector Here H H 1 H 2 β β (P , P) is the Hamiltonian generated by the kinetic energy of all the quarks and interaction between pairs of quarks that depends on the relative momentum of the hadrons, since the other degrees of freedom are fixed by the hadron states. There is also a part of the interaction which is generated by the coupling with one-hadron states and is given by It is interesting to notice that this effective potential has special relevance at E ∼ M α . Furthermore, one should expect attraction for E < M α since V e f f ββ (P, P) < 0 and repulsion in the other case. This is important when one considers a certain threshold, since states above threshold will help to bind a molecule while states below threshold will help to unbind it. This analysis helps to know when a dynamically generated state can appear.
This formalism is suitable for bound states. However, in order to solve the scattering or consider resonances, it is convenient to work with the equivalent Lippmann-Schwinger equation written as The solution to this equations is given in Ref. [46] T β β (E; P , P) = T β β The first term in the right hand side is the non-resonant contribution given by the solution of the equation The resonant part include the dressed vertex functions and the dressed two hadron propagator defined as the inverse of The dressed propagator has singularities at the energies of the resonance states so, to find these energies, we solve the equation Once the resonance energiesĒ are known, we find one-hadron amplitudes by solving and the two-hadron wave function is given by Notice that the normalization of the state requires The solution with naive quark model states |ψ α is only exact if all the state α are included. However, including only those states close to the energy range under consideration has been shown to be a good approximation. In Ref. [47], the unquenched quark model for charmonium mesons was considered but solving not only for the relative two-hadron wave function, but for the wave function of the ψ α meson as a cc meson, getting very similar results to the present approximation.

Coupled channel effects.
As mentioned before, since 2003 it is clear that the naive quark model is not enough to understand the heavy hadron spectra. In some cases, as the pentaquarks, the energy scale of its mass make unavoidable to include higher Fock components. However, as we will see, in other cases threshold effects can easily explain properties very difficult to understand in the naive quark model. In this section we give a few examples of such cases and we will summarize results obtained using the model previously introduced.

Isospin breaking effects.
The X(3872) was discovered in the J/ψππ invariant mass distribution of the B + → K + π + π − J/ψ decay. The two-pions in this decay came from a ρ 0 meson [48], which is an isospin 1 final state. The ratio of the decay into three pions was also measured and the three-pions came from the decay of an ω meson [49] which is an isospin 0 channel. The ratio between these two decay modes was found to be So this state can decay into final states with two different values of isospin, which implies that either isospin is violated in the decay process or the isospin of the X(3872) is not well defined. The X(3872) is now included in the PDG as the χ c1 (3872). Quark models usually predict this state at higher energies, although the deviation can be explained if one considers that this theoretical state is close and above the DD * threshold. One important question is whether the X(3872) is the state expected in this energy region by quark models, or if it is an additional state dynamically generated in the DD * channel. In any case, the crucial property of this state is that its mass is very close to the D 0D * 0 with a binding energy given by where the first number is from Ref. [50], the second from Ref. [51] and the third is from Ref. [52]. Since the mass is so close to the DD * threshold, one would expect a molecule or a mixing with some charmonium state. Considering the large isospin breaking, the most promising source is the mass splitting between charge and neutral states of D and D * mesons, finding a larger scale than the binding energy, which suggests a big effect. Notice that the masses by themselves do not suggest it, since the breaking is only of 0.13 % and 0.08 % for the D and D * mesons, respectively. This effect was introduced by Swanson [53] in a coupled channel calculation in which an isospin 1 channel J/ψρ was introduced. The important point to notice is that the binding energy for the charged channel is around 8 MeV, so the size of this component is around 1.5 fm, while for the neutral the small binding energy B makes the size of the order of 4 fm or bigger. This effect generates a big isospin breaking on the wave function out of the interaction region. This assertion generated some confusion since the isospin breaking effect is small in the interaction region. The isospin breaking was further analyzed in Ref. [54] in the framework of an Effective Field Theory, were the coupling of the states to the different final channels could be evaluated. The couplings obtained for DD * states were g = 2982 and g = 3005 for charged and neutral channels respectively, showing an isospin breaking of less than 1 %. In fact, although the ratio given in Eq. (56) suggests a big isospin breaking, this is only due to the big phase space effects that enhances the isospin 1 channel against the isospin 0 [54]. Excluding phase space effects, the decay in the ρ channel is only around 3.2 % of the ω channel decay. This was clarified in Ref. [55] relating the couplings with the probability of the wave function in the interacting region.
The microscopic calculation at the quark level was performed in Ref. [56] in the framework of the chiral quark model previously described. Within the model, the naive quark model χ c1 (2P) state has a mass of 3947 MeV, which is far above the X(3872). Coupling with DD * states makes the naive quark model masses change slightly. However, the important effect is that a new state appears in the DD * threshold with properties in overall good agreement with those of the X(3872). The ratio of Eq. (56) was analyzed in Ref. [57], finding a value close to the experimental result. The big isospin breaking on the wave function is represented in Figure 2. For exact isospin symmetry, the DD * charged and neutral components should have the same probabilities, while we see that, close to the DD * neutral threshold, this component dominates.
This effect, seen on the X(3872) state, may appear in any hadron-hadron molecule close to threshold. Of particular interest are the famous pentaquark states [19,20], which are close to the Σ cD , Σ * cD and Σ cD * thresholds. This is the reason why it is widely accepted that the nature of these states is more likely to be a hadron-hadron molecule than a compact pentaquark state. Being close to the threshold, isospin breaking effects were studied [58] and these effects could be magnified in the P c (4457) [59]. For this state, the binding energy is 2.5 +1.8 −4.2 MeV for the charged channel with lower threshold and 6.9 +1.8 −4.1 MeV for the higher threshold. Analogously to the decays into J/ψρ and J/ψω, the pentaquark could decay into an isospin 3/2 channel J/ψ∆ + or an isospin 1/2 channel J/ψp. Within an EFT framework, the ratio was evaluated and showed to be up to 30 %. The measurement of this isospin violating decay of the pentaquark could be the best indication of its molecular nature.

HQSS and HFS breaking.
Heavy Quark Spin Symmetry (HQSS) and Heavy Flavor Symmetry (HFS) are good approximate symmetries of QCD, so one would expect them to be realized in the heavy hadron spectrum.
If we look into the heavy-light sector, under exact HQSS, the D and D * mesons should have the same mass. Despite it is not exactly realized, the ratio shows that the breaking is, indeed, small. In the hidden charm sector, we have ∼ 2.0%, so even smaller breakings. HFS implies that interactions do not depend on the heavy quark mass, so when we find a state in the charm sector, there must exist an analog in the bottom sector.
If we consider the X(3872) to be a DD * molecule, HQSS [60][61][62] leads to unavoidable predictions. The interaction between DD * (1 ++ ) and D * D * (2 ++ ) channels is the same, so if the X(3872) is a DD * (1 ++ ) molecule it implies that there should be a bound state with very similar binding energy in the D * D * (2 ++ ), which was dubbed X(4012). Besides, HFS requires the interaction between charmed mesons to be the same as for bottom mesons [63], so the same molecules observed in the hidden-charm sector should appear in the hidden-bottom sector.
HQSS is usually fulfilled by heavy quark models, since the heavy quark mass only appears in fine structure terms that are suppressed as 1/m Q corrections. Within HQSS, one finds for S partial-waves [ where H represents the interacting Hamiltonian. In the chiral quark model of previous sections, the results for diagonal matrix elements are shown in Figures 3 and 4, showing that HQSS is approximately fulfilled. Additionally, comparing the matrix elements of the interactions in the charmed and bottom sectors, we can see that HFS is also fulfilled. Nevertheless, open thresholds appear in energy regions where we can also find naive quark model states, so one-hadron and two-hadron states can be mixed as shown in the previous sections. This effect can have important consequences for states that are mainly dynamically generated molecules. Let's see an example and consider the P-wave charmonium states. Considering the allowed spins, using the spectroscopic notation 2S+1 L J , we may have the states 1 P 1 (J PC = 1 +− ), 3 P 0 (0 ++ ), 3 P 1 (1 ++ ) and 3 P 2 (2 ++ ) coupled to 3 [7]. The X(3915) has been included in the 0 ++ and in the 2 ++ channels since the J quantum number is not known. The quantum numbers of the X(3940) are also not known but it has been included in the 1 ++ channel since it has been seen in DD * but not in DD. The states in red are naive QQ states predicted by the model.
states are in the region of the D ( * )D( * ) thresholds, however the threshold difference M D * D * − M DD ∼ 282 is larger. Taking only S-wave two-meson states, only DD and D * D * can have 0 ++ quantum numbers, while only DD * can have 1 ++ and D * D * can have 2 ++ . This means that the relevant threshold in each channel will have a different relative position with respect to naive quark model states. This is shown in Fig. 6, where we can see that, for the 1 ++ channel, the naive quark model state is above the DD * threshold, giving additional attraction, while in the 2 ++ channel the P-wave state is below the threshold, giving repulsion 1 . This explains why the 1 ++ channel has an additional state, while the 2 ++ has not, which is against HQSS expectations. A systematic study of this effect was performed at hadron level in Ref. [65] and a more elaborate study at the quark level was performed in Ref. [66]. In bottomonium we have a different situation. In the 1 ++ channel the naive quark model state generates repulsion and no additional state appears, while in the 2 ++ channel there is attraction from the state above threshold and repulsion from the state below and the final result is that an additional state appears. This result is, again, against HFS expectations [64]. A more elaborate calculation including more thresholds in underway.

Threshold cusps
These are enhancements of the cross sections near the opening of a threshold. One famous example is the measurement of the scattering lengths difference a 0 − a 2 in ππ scattering using a cusp-like enhancement in the π 0 π 0 invariant mass distribution in the K ± → π ± π 0 π 0 decay [67]. It was first noticed in Ref. [68] and then proposed to be measured by Cabibbo [69]. Due to the very precise experimental data available, a very precised determination of this combination was performed [70] that was in very good agreement with χPT predictions. 1 The state above D * D * threshold is an F state. This effect has been also used to explain strong energy dependencies near threshold of invariant mass distributions, as in the case of Z c and Z b states [71]. However it has been argued [72] that such big effects may not appear without the existence of a nearby pole (bound, virtual or resonance state).
One example of a threshold cusp effect in the hidden-charm sector is the X(4140) resonance, a J PC = 1 ++ structure observed in the φJ/ψ invariant mass spectrum by many collaborations, such as CDF [73], D0 [74], CMS [75], Belle [76], BaBar [77] and LHCb [78]. Within the chiral quark model described above, a coupled calculation of the main open-charm channels [79] showed that the structure just above the φJ/ψ threshold is not caused by the effect of a nearby cc pole, but it is associated to the presence of the D sD * s channel. The residual D sD * s interaction is strong enough to show a rapid increase of the experimental counts, but too weak to develop a bound or virtual state.
Even more interesting is the case of Z c and Z b states. In the chiral quark model the Z c states have been studied in Ref. [80]. The Z c (3900) ± /Z c (3885) ± and Z c (4020) ± are meson states in the charmonium energy range. The fact that they are charged rules out the possibility of being cc states, and, at least, four quarks are needed. Which is the nature of these charged states is still an open question.
The Z c (3900) was discovered by the BESIII [81] and Belle [82] Collaborations in the π J/ψ invariant mass distribution of the reaction e + e − → π + π − J/ψ. It was then seen by the BESIII Collaboration [83] in the DD * invariant mas distribution of the e + e − → π ± (DD * ) ∓ reaction with a lower mass and was referred to as the Z c (3885), although now are seen as the same state. Soon after this discovery, the BESIII Collaboration reported the discovery of another charged state, the Z c (4020), in the reaction e + e − → π + π − h c [84]. Later on, also BESIII, reported about the neutral partner [85], completing the isospin triplet.
In Ref. [80] D * D( * ) was analyzed in the sector I G (J PC ) = 1 + (1 +− ), within the formalism previously mentioned. Although here there is no qq state coupled to two-meson components, this system is interesting for another reason. There are two close-by channels that one would not expect to have an important effect, the π J/ψ and ρη c , since the interactions between these mesons is expected to be small. However, the non-diagonal interaction D * D( * ) − π J/ψ and D * D( * ) − ρη c are dominant, and they do not generate bound states but virtual states.
In Table 5 we give the pole position for the states corresponding to the Z c (3900) ± and Z c (4020) ± . The poles are below the DD * and D * D * thresholds, respectively, and in the second Riemann sheet Table 5. The S-matrix pole positions, in MeV/c 2 , for different coupled-channels calculations [80]. The included channels for each case are shown in the first column. Poles are given in the second and fourth columns by the value of the complex energy in a specific Riemann sheet (RS). The RS columns indicate whether the pole has been found in the first (F) or second (S) Riemann sheet of a given channel. Each channel in the coupled-channels calculation is represented as an array's element, ordered with increasing energy.
Calculation  [80]. Experimental data are from Ref. [86,87], respectively. The theoretical line shapes have been convoluted with the experimental resolution. The line-shape's 68% uncertainty is shown as a shadowed area. corresponding to virtual states. Despite these poles already emerge when the main open-charm channels are included, other channels are important in order to describe the experimental lineshapes. Lineshapes for different reactions are shown in Figures 7 and 8. Although we find poles on the S-matrix that produce structures in the lineshapes, in some cases they don't seem enough, which supports the idea that a threshold cusp effect might not be sufficient to describe the experimental data without the existence of some associated pole.

Conclusion
The naive quark model has been very successful in describing heavy hadron phenomenology for a very long time. However, since 2003, it seems clear that the coupling with two hadron states are of relevance to describe the phenomenology of new discovered states.
In this work we have describe how a microscopic quark model can be used, and in particular the chiral quark model, to describe systems in which conventional quark model states can couple to two-hadron states in a consistent framework. Although these effects are of no relevance in many low-lying states, keeping the validity of the naive quark model, in some cases the effects can be very important.
This framework have been used during the last years to study the meson spectrum and we have shown a few examples were deviations from naive quark model expectations are of special relevance. In fact, threshold effects can generate deviations from expected results predicted by well-known symmetries  [80]. Experimental data for D * D * are from Refs. [88]. The theoretical line shapes have been convoluted with the experimental resolution. The line-shape's 68%-uncertainty is shown as a shadowed area.
such as HFS or HQQS, remarking the importance of analyzing such effects specially in the heavy meson and baryon sectors.