Next Article in Journal
A Fuzzy Gain-Based Dynamic Ant Colony Optimization for Path Planning in Dynamic Environments
Next Article in Special Issue
Nonequilibrium Dynamics of the Chiral Quark Condensate under a Strong Magnetic Field
Previous Article in Journal
More Effective Conditions for Oscillatory Properties of Differential Equations
Previous Article in Special Issue
Strong-Interaction Matter under Extreme Conditions from Chiral Quark Models with Nonlocal Separable Interactions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coupling Hadron-Hadron Thresholds within a Chiral Quark Model Approach

by
Pablo G. Ortega
1,† and
David R. Entem
2,*,†
1
Department of Fundamental Physics and IUFFyM, U. Salamanca, E-37008 Salamanca, Spain
2
Grupo de Física Nuclear and IUFFyM, U. Salamanca, E-37008 Salamanca, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2021, 13(2), 279; https://doi.org/10.3390/sym13020279
Submission received: 18 December 2020 / Revised: 23 January 2021 / Accepted: 2 February 2021 / Published: 6 February 2021
(This article belongs to the Special Issue Advances in Chiral Quark Models)

Abstract

:
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 give an overview of a way in which 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.

1. 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 q q ¯ 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] a few years before.
Afterwards, heavy meson spectroscopy developed very fast during the following years. The Υ ( 1 S ) , the first state with bottom quarks discovered, was found at Fermilab [4] in 1977. Already in 1980, there were 11 new mesons included in the Particle Data Group (PDG) 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 evidence of charmed baryons came six months after the discovery of the J / ψ , in 1975 [8], but, in 1980, only the Λ c and Σ c baryons were 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 Υ ( 4 S ) 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 make many important contributions to heavy hadron physics.
From a theoretical point of view, the description of the hadron spectroscopy has mostly relied on potential models since the discovery of the J / ψ meson. Regarding heavy quarkonium spectroscopy, potential models have been very successful at describing the charmonium and bottomonium spectra (see, e.g., [9,10,11,12,13,14,15,16,17]). It is worth mentioning that new theoretical approaches directly derived from QCD, such as Effective Field Theories (EFTs) like non-relativistic QCD (NRQCD) [18,19,20] or lattice gauge theories [21,22,23] have been recently employed to study hadron spectroscopy. However, due to their complexity, many features of the hadron spectroscopy are more easily addressed with potential model approaches. Already in 1978, the Cornell model [24] was developed to understand heavy meson spectroscopy. This model is a non-relativistic approach for heavy quarks with interactions that are governed by S U ( 3 ) color gauge symmetry, with the 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 that 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 [25]. 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 [26] and more recently with NRQCD up to N 3 L O [27].
For the baryon spectra, also very soon, in 1979, quark models developed for the light sector were applied in the heavy quark sector [28]. In addition, Stanley and Robsen [29] 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 [24], the coupling with two-meson thresholds was considered for excited states, although it was found to be 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 [30]. It was very soon confirmed by the CDF [31], D0 [32], and BaBar [33] Collaborations. It has some intriguing properties that are 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 [34,35]. 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 c c ¯ 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 overview 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.

2. The Naive Chiral Quark Model

The quark model we use is a constituent quark model based on spontaneous chiral symmetry breaking [36]. It was first applied to the light-quark sector [37] and then extended to the heavy sector [38].
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 [39], 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 [40] and, in unquenched QCD, this string is broken [41] when there is enough energy to produce a quark–antiquark pair. Finally, we introduce QCD perturbative effects through the one-gluon exchange interaction [42]. The model has been reviewed in many works, and all of the details can be found in Refs. [43,44].
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 Rayleigh–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 whether 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 straightforwardly 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 [45] and has been applied to many few body problems [46,47].
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
φ n l m ( r ) = N n l r l e ν n r 2 Y l m ( r ^ )
N n l = 2 l + 2 ( 2 ν n ) l + 3 / 2 π ( 2 l + 1 ) ! ! 1 / 2 .
The basis function is generated by taking several values for the Gaussian parameter ν n . In the GEM, these parameters are taken in geometrical progression as
ν n = 1 r n 2
r n = r min a n 1
a = r max r min 1 N 1
where r min and r max are the minimum and maximum radius and N the number of Gaussians. r min and r max should be chosen such that the range of the wave functions is well represented. For excited states, usually the range increases, so states with a higher range than r max will not be correctly represented. r min should be below the range of the potential. Numerically, it is important that the parameter a is not too close to 1 so that the problem does not become singular, so N can not be increased arbitrarily for fixed values of r min and r max . Then, one can use the expansion
φ ( r ) = n l c n l φ n l m ( r )
and ends up with the generalized eigenvalue problem
n l ( H n l , n l E N n l , n l ) c n l = 0
with
H n l , n l = φ n l m | H | φ n l m
N n l , n l = φ n l m | φ n l m .
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
φ ( r ) = n l c n l φ n l ( r ) χ S J χ F ξ c
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 is needed, typical accurate three-body Fadeev calculations of Triton binding energy need up to 38 of such partial waves. The GEM approaches the problem in a different way; it also considers 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 lower number of partial waves are 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
R c m = m 1 u 1 + m 2 u 2 + m 3 u 3 M T
r i = u j u k
R i = u i m j u j + m k u k m j + m k
where u i is the position vector of particle i, and ( i j k ) is one of the three even permutations of ( 123 ) . The orbital wave function is taken as
φ = i = 1 3 n i l i N i L i c n i l i N i L i φ n i l i ( r i ) φ N i L i ( R i ) L T
where ( 1 ) l i + L i gives the parity, and L T is the total orbital angular momentum. 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
P i j φ n i l i ( r i ) φ N i L i ( R i ) L T = ( 1 ) l i φ n i l i ( r j ) φ N i L i ( R j ) L T
P i j φ n j l j ( r j ) φ N j L j ( R j ) L T = ( 1 ) l j φ n j l j ( r i ) φ N j L j ( R i ) L T
P i j φ n k l k ( r k ) φ N k L k ( R k ) L T = ( 1 ) l k φ n k l k ( r k ) φ N k L k ( R k ) L T .
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
φ = n 1 l 1 N 1 L 1 c n 1 l 1 N 1 L 1 φ n 1 l 1 ( r 1 ) φ N 1 L 1 ( R 1 ) L T
φ = n 23 l 23 N 23 L 23 c n 23 l 23 N 23 L 23 φ n 23 l 23 ( r 2 ) φ N 23 L 23 ( R 2 ) L T ± φ n 23 l 23 ( r 3 ) φ N 23 L 23 ( R 3 ) L T
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. In principle, one could only work with one basis including all possible angular momentums but numerically is more efficient to use both with the lower angular momentum. 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 Table 1 and Table 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 such as the Coulomb one, the GEM works very well.
If the three particles are identical, then the wave function to be used is
φ = n l N L c n l N L i = 1 3 φ n l ( r i ) φ N L ( R i ) L T .
As explained above, we have to include the spin–flavor–color wave function for baryons in a quark model 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 [49] model, calculated with the GEM and compared with a Fadeev calculation by B. Silvestre-Brac [50]. The GEM calculation only includes wave functions with l i = L i = 0 . In addition, the matter radius square and the charge radius square defined by
R m 2 = Ψ | i = 1 3 m i M ( u i R c m ) 2 | Ψ
R c 2 = Ψ | i = 1 3 e i ( u i R c m ) 2 | Ψ
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 [51], 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.

3. The 3 P 0 Model

The quark–pair creation model or 3 P 0 model is a microscopic model that allows for coupling channels with different number of quarks (Figure 1). 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 [54] and, afterwards, Le Yaouanc et al. applied it to the strong decays of mesons [55] and baryons [56]. These authors also evaluated strong decay partial widths of the three charmonium states ψ ( 3770 ) , ψ ( 4040 ) , and ψ ( 4415 ) within the same model [57,58].
The 3 P 0 model is usually formulated in terms of the Hamiltonian operator
H I = 3 g s d 3 x ψ ¯ ( x ) ψ ( x )
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.
It can be also formulated in terms of a transition operator given by
T = 3 μ , ν d 3 p μ d 3 p ν δ ( 3 ) ( p μ + p ν ) g s 2 m μ 2 5 π × Y 1 p μ p ν 2 1 2 1 2 1 0 a μ ( p μ ) b ν ( p ν )
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 l m ( p ) = p l Y l m ( 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. [59]. The relation is given as
γ = g s / 2 m
m being the mass of the pair created, which is usually a light pair.
We consider a process 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 c m = 0 , with P A the momentum of the initial hadron and P c m the total momentum of the final hadrons. Then, the matrix element of the transition operator is written as
B C | T | A = δ ( 3 ) ( P c m ) M A B C .
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
| M = ϕ M ( 1 , 2 ) b 1 d 2 | 0 d 3 p 1 d 3 p 2
| B = N B Ψ B ( 1 , 2 , 3 ) b 1 b 2 b 3 | 0 d 3 p 1 d 3 p 2 d 3 p 3
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,
| ϕ M ( 1 , 2 ) | 2 d 3 p 1 d 3 p 2 = 1
| Ψ B ( 1 , 2 , 3 ) | 2 d 3 p 1 d 3 p 2 d 3 p 3 = 1
where the sum is for discrete degrees of freedom. With these states, two hadron states with correct quantum numbers are constructed.
With the transition amplitude, the widths for strong decays can be evaluated
Γ A B C = 2 π E B ( k 0 ) E C ( k 0 ) m A k 0 | M A B C ( k 0 ) | 2
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,
A = k c k A k
B = i c i B i
C = j c j C j
where | A k , | B k , and | C k are the different GEM components for the hadrons and the c’s the coefficients for the solution of the GEM, one can evaluate the matrix element in terms of matrix elements in the GEM basis as
B C | T | A = i j k c i * c j * c k B i C j | T | A j
where the matrix elements 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, a scale-dependent parameter was considered in Ref. [60] as
γ ( μ ) = γ 0 log ( μ μ 0 )
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.

4. The Unquenched Quark Model

The unquenched quark model [61] or unitarized quark models [62] include the effects of two-hadron thresholds in one hadron state through hadron loops. They are usually implemented with approximate wave functions, and the interaction between the two-hadrons in the loops is not considered. Here, we present how these effects can be incorporated in an easy way.
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
Ψ = α c α ψ α + β χ β ( P ) ψ H 1 ψ H 2 β
where | ψ α are naive quark model one-hadron states with α quantum numbers and | ψ H 1 ψ H 2 β are two-hadron state H 1 H 2 with β quantum numbers and relative momentum P. We define now the transition amplitude
ψ H 1 ψ H 2 β | T | ψ α = P δ ( 3 ) ( P c m ) h β α ( P ) .
If we impose the Schrödinger equation
H Ψ = E Ψ
and solve for the one-hadron amplitudes, we find
c α = 1 E M α β h α β ( P ) χ β ( P ) P 2 d P
where χ β ( P ) is given by the equation in the two-hadron sector
β H β β H 1 H 2 ( P , P ) + V β β eff ( P , P ) χ β ( P ) P 2 d P = E χ β ( P ) .
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 depend 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
V β β eff ( P , P ) = α h β α ( P ) h α β ( P ) E M α .
It is interesting to notice that this effective potential has special relevance at an energy close to the threshold of a two meson channel β , E E β thr , where one expects the molecules to appear. Furthermore, one should expect attraction for E E β thr < M α since V β β e f f ( P , P ) < 0 and there is repulsion in the other case, it being more intense when the bare state is closer to the threshold M α E β thr . Thus, 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
T β β ( E ; P , P ) = V T β β ( P , P ) + β d P P 2 V T β β ( P , P ) 1 E E β ( P ) T β β ( E ; P , P )
with
V T β β ( P , P ) = V β β ( P , P ) + V eff β β ( P , P )
H β β H 1 H 2 ( P , P ) = δ β β δ ( P P ) P 2 P 2 2 μ β + V β β ( P , P ) .
The solution to this equation is given in Ref. [63]
T β β ( E ; P , P ) = T V β β ( E ; P , P ) + α , α ϕ β α ( E ; P ) Δ α α 1 ( E ) ϕ ¯ α β ( E ; P ) .
The first term on the right-hand side is the non-resonant contribution given by the solution of the equation
T V β β ( E ; P , P ) = V β β ( P , P ) + β d P P 2 V β β ( P , P ) 1 E E β ( P ) T V β β ( E ; P , P ) .
The resonant part include the dressed vertex functions
ϕ α β ( E ; P ) = h α β ( P ) β T V β β ( E ; P , q ) h α β ( q ) q 2 / 2 μ β E q 2 d q ,
ϕ ¯ α β ( E ; P ) = h α β ( P ) β h α β ( q ) T V β β ( E ; q , P ) q 2 / 2 μ β E q 2 d q
and the dressed two hadron propagator defined as the inverse of
Δ α α ( E ) = ( E M α ) δ α α + G α α ( E ) ,
G α α ( E ) = β d q q 2 ϕ α β ( q , E ) h β α ( q ) q 2 / 2 μ β E .
The dressed propagator has singularities at the energies of the resonance states so, to find these energies, we solve the equation
Δ α α ( E ¯ ) = ( E ¯ M α ) δ α α + G α α ( E ¯ ) = 0 .
Once the resonance energies E ¯ are known, we find one-hadron amplitudes by solving
M α δ α α G α α ( E ¯ ) c α ( E ¯ ) = E ¯ c α ( E ¯ )
and the two-hadron wave function is given by
χ β ( P ) = 2 μ β α ϕ β α ( E ; P ) c α P 2 k β 2
with k β 2 = 2 μ β ( E E β thr ) . Notice that the normalization of the wave-function requires
α | c α | 2 + β χ β | χ β = 1 .
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. [64], 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 c c ¯ meson, getting very similar results to the present approximation.

5. Coupled Channel Effects

As mentioned before, since 2003, it has been clear that the naive quark model is not enough to understand the heavy hadron spectra. In some cases, such as the pentaquarks, the energy scale of its mass makes it unavoidable to include higher Fock components. However, as we will see, in other cases, threshold effects can easily explain properties that are 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.

5.1. 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 [65], 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 [66], which is an isospin 0 channel. The ratio between these two decay modes was found to be
B ( X ( 3872 ) π + π π 0 J / ψ ) B ( X ( 3872 ) π + π J / ψ ) = 1.0 ± 0.4 ( stat ) ± 0.3 ( syst ) .
Thus, 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 χ c 1 ( 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 D D ¯ * 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 D D ¯ * channel. In any case, the crucial property of this state is that its mass is very close to the D 0 D ¯ * 0 with a binding energy given by
(57) B m D * 0 + m D ¯ 0 m X ( 3872 ) = 1 . 1 0.4 0.3 + 0.6 + 0.1 MeV (58) = ( 0.00 ± 0.18 ) MeV (59) = ( 0.07 ± 0.12 ) MeV
where the first number is from Ref. [67], the second from Ref. [68], and the third is from Ref. [69].
Since the mass is so close to the D D ¯ * 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
m D * + + m D m D * 0 m D ¯ 0 = 8.2 ± 0.1 MeV B
to be 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 [70] 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. [71] in the framework of an Effective Field Theory, where the coupling of the states to the different final channels could be evaluated. The couplings obtained for D D ¯ * states were g = 2982 MeV and g = 3005 MeV for charged and neutral channels, respectively, showing an isospin breaking of less than 1 % . In fact, although the ratio given in Equation (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 [71]. Excluding phase space effects, the decay in the ρ channel is only around 3.2 % of the ω channel decay. This was clarified in Ref. [72] 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. [73] in the framework of the chiral quark model previously described. Within the model, the naive quark model χ c 1 ( 2 P ) state has a mass of 3947 MeV, which is far above the X ( 3872 ) . Coupling with D D ¯ * states makes the naive quark model masses change slightly. However, the important effect is that a new state appears in the D D ¯ * threshold with properties in overall good agreement with those of the X ( 3872 ) . This is in contrast to other unquenched quark models [74] where the X ( 3872 ) is not an additional state but just a corrected naive quark model state by two meson loops.
The ratio of Equation (56) was analyzed in Ref. [75], 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 D D ¯ * charged and neutral components should have the same probabilities, while we see that, close to the D D ¯ * neutral threshold, this component dominates.
This effect, seen on the X ( 3872 ) state, may appear in any hadron–hadron molecule close to the threshold. Of particular interest are the famous pentaquark states [34,35], which are close to the Σ c D ¯ , Σ c * D ¯ , and Σ c D ¯ * 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 [76], and these effects could be magnified in the P c ( 4457 ) [77]. For this state, the binding energy is 2 . 5 4.2 + 1.8 MeV for the charged channel with lower threshold and 6 . 9 4.1 + 1.8 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
R Δ + / p = | B ( P c ( 4457 ) + J / ψ Δ + ) | | B ( P c ( 4457 ) + J / ψ p ) |
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.

5.2. 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 not exactly being realized, the ratio M D * M D M D * + M D 3.6 % shows that the breaking is, indeed, small. In the hidden charm sector, we have M J / ψ M η c M J / ψ + M η c 1.8 % and M χ c 2 ( 1 P ) M χ c 0 ( 1 P ) M χ c 2 ( 1 P ) + M χ c 0 ( 1 P ) 2.0 % , so even smaller breaking. 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 D D ¯ * molecule, HQSS [78,79,80] leads to unavoidable predictions. The interaction between D D ¯ * ( 1 + + ) and D * D ¯ * ( 2 + + ) channels is the same, so, if the X ( 3872 ) is a D D ¯ * ( 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 ) . In addition, HFS requires the interaction between charmed mesons to be the same as for bottom mesons [81], 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 [82]
2 3 D * D ¯ * ( 0 + + ) | H | D D ¯ ( 0 + + ) = D D ¯ ( 0 + + ) | H | D D ¯ ( 0 + + ) D * D ¯ * ( 0 + + ) | H | D * D ¯ * ( 0 + + )
D D ¯ * ( 1 + + ) | H | D D ¯ * ( 1 + + ) = D * D ¯ * ( 2 + + ) | H | D * D ¯ * ( 2 + + ) ,
D D ¯ * ( 1 + + ) | H | D D ¯ * ( 1 + + ) = 3 2 D D ¯ ( 0 + + ) | H | D D ¯ ( 0 + + ) 1 3 D * D ¯ * ( 0 + + ) | H | D * D ¯ * ( 0 + + )
2 D D ¯ * ( 1 + ) | H | D D ¯ * ( 1 + ) = D D ¯ ( 0 + + ) | H | D D ¯ ( 0 + + ) + D * D ¯ * ( 0 + + ) | H | D * D ¯ * ( 0 + + )
where H represents the interacting Hamiltonian. In the chiral quark model of previous sections, the results for diagonal matrix elements are shown in Figure 3 and Figure 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 2 S + 1 L J , we may have the states 1 P 1 ( J P C = 1 + ), 3 P 0 ( 0 + + ), 3 P 1 ( 1 + + ) and 3 P 2 ( 2 + + ) coupled to 3 F 2 . Under HQSS, all these states are degenerated (within a small deviation due to the breaking). The largest mass splitting is given by the difference of the ground states M χ c 2 ( 1 P ) M χ c 0 ( 1 P ) 141 MeV, the splitting for excited states being smaller in the naive quark model picture. The χ c J ( 2 P ) and h c ( 2 P ) states are in the region of the D ( * ) D ¯ ( * ) thresholds; however, the threshold difference M D * D ¯ * M D D ¯ 282 MeV is larger. Taking only S-wave two-meson states, only D D ¯ and D * D ¯ * can have 0 + + quantum numbers, while only D D ¯ * 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 Figure 5, where we can see that, for the 1 + + channel, the naive quark model state is above the D D ¯ * threshold, giving additional attraction, while, in the 2 + + channel, the P-wave state is below the threshold, giving repulsion (The state above D * D ¯ * threshold is an F state.). This explains why the 1 + + channel has an additional state, while the 2 + + does not, which is against HQSS expectations. A systematic study of this effect was performed at the hadron level in Ref. [83], and a more elaborate study at the quark level was performed in Ref. [84].
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 [82]. A more elaborate calculation including more thresholds is underway (Figure 6).

5.3. 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 [85]. It was first noticed in Ref. [86] and then proposed to be measured by Cabibbo [87]. Due to the very precise experimental data available, a very precise determination of this combination was performed [88] that was in very good agreement with χ PT predictions.
This effect has been also used to explain strong energy dependencies near the threshold of invariant mass distributions, as in the case of Z c and Z b states [89,90]. The effect should be present when the interaction is attractive [91]; however, it has been argued [92] 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 P C = 1 + + structure observed in the ϕ J / ψ invariant mass spectrum by many collaborations, such as CDF [93], D0 [94], CMS [95], Belle [96], BaBar [97], and LHCb [98]. Within the chiral quark model described above, a coupled calculation of the main open-charm channels [99] showed that the structure just above the ϕ J / ψ threshold is not caused by the effect of a nearby c c ¯ pole, but it is associated with the presence of the D s D ¯ s * channel. The residual D s D ¯ 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. For the Zc(3900), there exist analyses of the BESIII data considering two-body coupled channels and triangle singularities [100,101,102]. In the chiral quark model, the Z c states have been studied in Ref. [103]. 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 c c ¯ states, and, at least, four quarks are needed. The nature of these charged states is still an open question.
The Z c ( 3900 ) was discovered by the BESIII [104] and Belle [105] Collaborations in the π J / ψ invariant mass distribution of the reaction e + e π + π J / ψ . It was then seen by the BESIII Collaboration [106] in the D D ¯ * invariant mass distribution of the e + e π ± ( D D ¯ * ) reaction with a lower mass and was referred to as the Z c ( 3885 ) , although now it is 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 [107]. Later on, BESIII also reported about the neutral partner [108], completing the isospin triplet.
In Ref. [103], D * D ¯ ( * ) was analyzed in the sector I G ( J P C ) = 1 + ( 1 + ) , within the formalism previously mentioned. Although here there is no q q ¯ 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 are 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 D D ¯ * and D * D ¯ * thresholds, respectively, and in the second Riemann sheet corresponding to virtual states. Despite these poles already emerging 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 Figure 7 and Figure 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.

6. Conclusions

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 described 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 has been used during the last few years to study the meson spectrum, and we have shown a few examples where 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 such as HFS or HQQS, highlighting the importance of analyzing such effects especially in the heavy meson and baryon sectors.

Author Contributions

Conceptualization and methodology, P.G.O. and D.R.E.; Software, P.G.O. and D.R.E.; Validation, P.G.O. and D.R.E.; Writing—original draft preparation, P.G.O. and D.R.E.; Writing—re-view and editing, P.G.O. and D.R.E.; Funding acquisition, P.G.O. and D.R.E.; All authors have read and agreed to the published version of the manuscript.

Funding

This work has been funded by the Ministerio de Economía, Industria y Competitividad under Contract No. FPA2016-77177-C2-2-P and Ministerio de Ciencia, Innovación y Universidades under Contract No. PID2019-105439GB-C22, and by the EU Horizon 2020 research and innovation program, STRONG-2020 project, under Grant No. 824093.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
HQQSHeavy Quark Spin Symmetry
HFSHeavy Flavor Symmetry
χ PTChiral Perturbation Theory
QCDQuantum Chromodynamics

References

  1. Aubert, J.J.; Becker, U.; Biggs, P.J.; Burger, J.; Chen, M.; Everhart, G.; Goldhagen, P.; Leong, J.; McCorriston, T.; Rhoades, T.G.; et al. Experimental Observation of a Heavy Particle J. Phys. Rev. Lett. 1974, 33, 1404–1406. [Google Scholar] [CrossRef] [Green Version]
  2. Augustin, J.E.; Boyarski, A.M.; Breidenbach, M.; Bulos, F.; Dakin, J.T.; Feldman, G.J.; Fischer, G.E.; Fryberger, D.; Hanson, G.; Jean-Marie, B.; et al. Discovery of a Narrow Resonance in e+e Annihilation. Phys. Rev. Lett. 1974, 33, 1406–1408. [Google Scholar] [CrossRef] [Green Version]
  3. Glashow, S.L.; Iliopoulos, J.; Maiani, L. Weak Interactions with Lepton-Hadron Symmetry. Phys. Rev. D 1970, 2, 1285–1292. [Google Scholar] [CrossRef]
  4. Herb, S.W.; Hom, D.C.; Lederman, L.M.; Sens, J.C.; Snyder, H.D.; Yoh, J.K.; Appel, J.A.; Brown, B.C.; Brown, C.N.; Innes, W.R.; et al. Observation of a Dimuon Resonance at 9.5 GeV in 400-GeV Proton-Nucleus Collisions. Phys. Rev. Lett. 1977, 39, 252–255. [Google Scholar] [CrossRef]
  5. Kelly, R.L.; Horne, C.P.; Losty, M.J.; Rittenberg, A.; Shimada, T.; Trippe, T.G.; Wohl, C.G.; Yost, G.P.; Barash-Schmidt, N.; Bricman, C.; et al. Review of particle properties. Rev. Mod. Phys. 1980, 52, S1–S286. [Google Scholar] [CrossRef] [Green Version]
  6. Hagiwara, K.; Montanet, L.; Barnett, R.M.; Groom, D.E.; Trippe, T.G.; Wohl, C.G.; Armstrong, B.; Wagman, G.S.; Murayama, H.; Stone, J.; et al. Review of Particle Properties. Phys. Rev. D 2002, 66, 010001. [Google Scholar] [CrossRef] [Green Version]
  7. Particle Data Group; Zyla, P.A.; Barnett, R.M.; Beringer, J.; Dahl, O.; Dwyer, D.A.; Groom, D.E.; Lin, C.J.; Lugovsky, K.S.; Pianori, E.; et al. Review of Particle Physics. Prog. Theor. Exp. Phys. 2020, 2020, 083C01. [Google Scholar] [CrossRef]
  8. Cazzoli, E.G.; Cnops, A.M.; Connolly, P.L.; Louttit, R.I.; Murtagh, M.J.; Palmer, R.B.; Samios, N.P.; Tso, T.T.; Williams, H.H. Evidence for ΔS=−ΔQ Currents or Charmed-Baryon Production by Neutrinos. Phys. Rev. Lett. 1975, 34, 1125–1128. [Google Scholar] [CrossRef]
  9. Eichten, E.; Gottfried, K.; Kinoshita, T.; Lane, K.D.; Yan, T.M. Charmonium: The Model. Phys. Rev. D 1978, 17, 3090–3117, Erratum in 1980, 21, 313. [Google Scholar] [CrossRef]
  10. Eichten, E.; Gottfried, K.; Kinoshita, T.; Lane, K.D.; Yan, T.M. Charmonium: Comparison with Experiment. Phys. Rev. D 1980, 21, 203–233. [Google Scholar] [CrossRef]
  11. Gupta, S.N.; Radford, S.F.; Repko, W.W. Semirelativistic Potential Model for Charmonium. Phys. Rev. D 1985, 31, 160. [Google Scholar] [CrossRef] [PubMed]
  12. Barnes, T.; Close, F.E.; Page, P.R.; Swanson, E.S. Higher quarkonia. Phys. Rev. D 1997, 55, 4157–4188. [Google Scholar] [CrossRef] [Green Version]
  13. Ebert, D.; Faustov, R.N.; Galkin, V.O. Properties of heavy quarkonia and Bc mesons in the relativistic quark model. Phys. Rev. D 2003, 67, 014027. [Google Scholar] [CrossRef] [Green Version]
  14. Eichten, E.; Godfrey, S.; Mahlke, H.; Rosner, J.L. Quarkonia and their transitions. Rev. Mod. Phys. 2008, 80, 1161–1193. [Google Scholar] [CrossRef] [Green Version]
  15. Danilkin, I.V.; Simonov, Y.A. Channel coupling in heavy quarkonia: Energy levels, mixing, widths and new states. Phys. Rev. D 2010, 81, 074027. [Google Scholar] [CrossRef] [Green Version]
  16. Ferretti, J.; Santopinto, E. Higher mass bottomonia. Phys. Rev. D 2014, 90, 094022. [Google Scholar] [CrossRef] [Green Version]
  17. Godfrey, S.; Moats, K. Bottomonium Mesons and Strategies for their Observation. Phys. Rev. D 2015, 92, 054034. [Google Scholar] [CrossRef] [Green Version]
  18. Caswell, W.E.; Lepage, G.P. Effective Lagrangians for Bound State Problems in QED, QCD, and Other Field Theories. Phys. Lett. B 1986, 167, 437–442. [Google Scholar] [CrossRef]
  19. Pineda, A.; Soto, J. Effective field theory for ultrasoft momenta in NRQCD and NRQED. Nucl. Phys. B Proc. Suppl. 1998, 64, 428–432. [Google Scholar] [CrossRef] [Green Version]
  20. Brambilla, N.; Pineda, A.; Soto, J.; Vairo, A. Effective Field Theories for Heavy Quarkonium. Rev. Mod. Phys. 2005, 77, 1423. [Google Scholar] [CrossRef] [Green Version]
  21. Dudek, J.J.; Edwards, R.G.; Mathur, N.; Richards, D.G. Charmonium excited state spectrum in lattice QCD. Phys. Rev. D 2008, 77, 034501. [Google Scholar] [CrossRef] [Green Version]
  22. Gray, A.; Allison, I.; Davies, C.T.H.; Dalgic, E.; Lepage, G.P.; Shigemitsu, J.; Wingate, M. The Upsilon spectrum and m(b) from full lattice QCD. Phys. Rev. D 2005, 72, 094507. [Google Scholar] [CrossRef] [Green Version]
  23. Meinel, S. The Bottomonium spectrum from lattice QCD with 2+1 flavors of domain wall fermions. Phys. Rev. D 2009, 79, 094501. [Google Scholar] [CrossRef] [Green Version]
  24. Gottfried, K. The spectroscopy of massive quark-antiquark systems, Progress in Particle and Nuclear Physics. Prog. Part. Nucl. Phys. 1982, 8, 49–71. [Google Scholar] [CrossRef]
  25. Gostev, V.B.; Mineev, V.S.; Frenkin, A.R. The inverse problem of quantum mechanics for a linear potentia. Theor. Math. Phys. 1983, 56, 682–686. [Google Scholar] [CrossRef]
  26. Sumino, Y. QCD potential as a “Coulomb-plus-linear” potential. Phys. Lett. B 2003, 571, 173–183. [Google Scholar] [CrossRef] [Green Version]
  27. Mateu, V.; Ortega, P.G.; Entem, D.R.; Fernández, F. Calibrating the Naïve Cornell Model with NRQCD. Eur. Phys. J. C 2019, 79, 323. [Google Scholar] [CrossRef]
  28. Copley, L.A.; Isgur, N.; Karl, G. Charmed baryons in a quark model with hyperfine interactions. Phys. Rev. D 1979, 20, 768–775. [Google Scholar] [CrossRef]
  29. Stanley, D.P.; Robsen, D. Do Quarks Interact Pairwise and Satisfy the Color Hypothesis? Phys. Rev. Lett. 1980, 45, 235–238. [Google Scholar] [CrossRef]
  30. Choi, S.K.; Olsen, S.L.; Abe, K.; Abe, T.; Adachi, I.; Ahn, B.S.; Aihara, H.; Akai, K.; Akatsu, M.; Akemoto, M.; et al. Observation of a Narrow Charmoniumlike State in Exclusive B±K±π+πJ/ψ Decays. Phys. Rev. Lett. 2003, 91, 262001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Acosta, D.; Affolder, T.; Ahn, M.H.; Akimoto, T.; Albrow, M.G.; Ambrose, D.; Amerio, S.; Amidei, D.; Anastassov, A.; Anikeev, K.; et al. Observation of the Narrow State X(3872)→J/ψπ+π in p ¯ p Collisions at s = 1.96 TeV. Phys. Rev. Lett. 2004, 93, 072001. [Google Scholar] [CrossRef] [Green Version]
  32. Abazov, V.M.; Abbott, B.; Abolins, M.; Acharya, B.S.; Adams, D.L.; Adams, M.; Adams, T.; Agelou, M.; Agram, J.L.; Ahmed, S.N.; et al. Observation and Properties of the X(3872) Decaying to J/ψπ+π in p p ¯ Collisions at s = 1.96 TeV. Phys. Rev. Lett. 2004, 93, 162002. [Google Scholar] [CrossRef] [Green Version]
  33. Aubert, B.; Barate, R.; Boutigny, D.; Couderc, F.; Gaillard, J.M.; Hicheur, A.; Karyotakis, Y.; Lees, J.P.; Tisserand, V.; Zghiche, A.; et al. Study of the BJ/ψKπ+π decay and measurement of the BX(3872)K branching fraction. Phys. Rev. D 2005, 71, 071103. [Google Scholar] [CrossRef] [Green Version]
  34. Aaij, R.; Adeva, B.; Adinolfi, M.; Affolder, A.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; et al. Observation of J/ψp Resonances Consistent with Pentaquark States in Λ b 0 J/ψKp Decays. Phys. Rev. Lett. 2015, 115, 072001. [Google Scholar] [CrossRef] [Green Version]
  35. Aaij, R.; Beteta, C.A.; Adeva, B.; Adinolfi, M.; Aidala, C.A.; Ajaltouni, Z.; Akar, S.; Albicocco, P.; Albrecht, J.; Alessio, F.; et al. Observation of a Narrow Pentaquark State, Pc(4312)+, and of the Two-Peak Structure of the Pc(4450)+. Phys. Rev. Lett. 2019, 122, 222001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Manohar, A.; Georgi, H. Chiral quarks and the non-relativistic quark model. Nucl. Phys. B 1984, 234, 189–212. [Google Scholar] [CrossRef]
  37. Fernandez, F.; Valcarce, A.; Straub, U.; Faessler, A. The nucleon-nucleon interaction in terms of quark degrees of freedom. J. Phys. G Nucl. Part. Phys. 1993, 19, 2013–2026. [Google Scholar] [CrossRef]
  38. Vijande, J.; Fernández, F.; Valcarce, A. Constituent quark model study of the meson spectra. J. Phys. G Nucl. Part. Phys. 2005, 31, 481–506. [Google Scholar] [CrossRef] [Green Version]
  39. Burgio, G.; Schröck, M.; Reinhardt, H.; Quandt, M. Running mass, effective energy, and confinement: The lattice quark propagator in Coulomb gauge. Phys. Rev. D 2012, 86, 014506. [Google Scholar] [CrossRef] [Green Version]
  40. Bali, G.S. QCD forces and heavy quark bound states. Phys. Rep. 2001, 343, 1–136. [Google Scholar] [CrossRef] [Green Version]
  41. Bali, G.S.; Neff, H.; Düssel, T.; Lippert, T.; Schilling, K. Observation of string breaking in QCD. Phys. Rev. D 2005, 71, 114513. [Google Scholar] [CrossRef] [Green Version]
  42. De Rújula, A.; Georgi, H.; Glashow, S.L. Hadron masses in a gauge theory. Phys. Rev. D 1975, 12, 147–162. [Google Scholar] [CrossRef] [Green Version]
  43. Segovia, J.; Yasser, A.M.; Entem, D.R.; Fernández, F. JPC=1−− hidden charm resonances. Phys. Rev. D 2008, 78, 114033. [Google Scholar] [CrossRef]
  44. Segovia, J.; Entem, D.R.; Fernandez, F.; Hernandez, E. Constituent quark model description of charmonium phenomenology. Int. J. Mod. Phys. E 2013, 22, 1330026. [Google Scholar] [CrossRef] [Green Version]
  45. Kamimura, M. Nonadiabatic coupled-rearrangement-channel approach to muonic molecules. Phys. Rev. A 1988, 38, 621–624. [Google Scholar] [CrossRef] [PubMed]
  46. Hiyama, E.; Kino, Y.; Kamimura, M. Gaussian expansion method for few-body systems. Prog. Part. Nucl. Phys. 2003, 51, 223–307. [Google Scholar] [CrossRef]
  47. Hiyama, E. Gaussian expansion method for few-body systems and its applications to atomic and nuclear physics. Prog. Theor. Exp. Phys. 2012, 2012, 01A204. [Google Scholar] [CrossRef] [Green Version]
  48. Morton, D.; Wu, Q.; Drake, G.W. Energy Levels for the Stable Isotopes of Atomic Helium (4He I and 3He I). Can. J. Phys. 2006, 84, 83–105. [Google Scholar] [CrossRef]
  49. Bhaduri, R.K.; Cohler, L.E.; Nogami, Y. A unified potential for mesons and baryons. Il Nuovo Cimento A 1981, 65. [Google Scholar] [CrossRef]
  50. Silvestre-Brac, B. Spectrum and Static Properties of Heavy Baryons. Few-Body Syst. 1996, 20. [Google Scholar] [CrossRef]
  51. Kandula, D.Z.; Gohle, C.; Pinkert, T.J.; Ubachs, W.; Eikema, K.S.E. Extreme Ultraviolet Frequency Comb Metrology. Phys. Rev. Lett. 2010, 105, 063001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Martin, W.C. Energy Levels of Neutral Helium (4He I). J. Phys. Chem. Ref. Data 1973, 2, 257–266. [Google Scholar] [CrossRef] [Green Version]
  53. Tech, J.L.; Ward, J.F. Accurate Wavelength Measurement of the 1s2p3P0 − 2p23P Transition in 4He I. Phys. Rev. Lett. 1971, 27, 367–370. [Google Scholar] [CrossRef]
  54. Micu, L. Decay rates of meson resonances in a quark model. Nucl. Phys. 1969, B10, 521–526. [Google Scholar] [CrossRef]
  55. Le Yaouanc, A.; Oliver, L.; Pène, O.; Raynal, J.C. “Naive” Quark-Pair-Creation Model of Strong-Interaction Vertices. Phys. Rev. D 1973, 8, 2223–2234. [Google Scholar] [CrossRef]
  56. Le Yaouanc, A.; Oliver, L.; Pène, O.; Raynal, J.C. Naive quark-pair—Creation model and baryon decays. Phys. Rev. D 1974, 9, 1415–1419. [Google Scholar] [CrossRef]
  57. Yaouanc, A.L.; Oliver, L.; Pene, O.; Raynal, J.C. Strong decays of ψ(4028) as a radial excitation of charmonium. Phys. Lett. B 1977, 71, 397–399. [Google Scholar] [CrossRef]
  58. Yaouanc, A.L.; Oliver, L.; Pène, O.; Raynal, J. Why is ψ(4414) so narrow? Phys. Lett. B 1977, 72, 57–61. [Google Scholar] [CrossRef]
  59. Ackleh, E.S.; Barnes, T.; Swanson, E.S. On the mechanism of open-flavor strong decays. Phys. Rev. D 1996, 54, 6811–6829. [Google Scholar] [CrossRef] [Green Version]
  60. Segovia, J.; Entem, D.; Fernández, F. Scaling of the P03 strength in heavy meson strong decays. Phys. Lett. B 2012, 715, 322–327. [Google Scholar] [CrossRef] [Green Version]
  61. Bijker, R.; Santopinto, E. Unquenched quark model for baryons: Magnetic moments, spins, and orbital angular momenta. Phys. Rev. C 2009, 80, 065210. [Google Scholar] [CrossRef] [Green Version]
  62. Heikkilä, K.; Törnqvist, N.A.; Ono, S. Heavy c c ¯ and b b ¯ quarkonium states and unitarity effects. Phys. Rev. D 1984, 29, 110–120. [Google Scholar] [CrossRef]
  63. Baru, V.; Hanhart, C.; Kalashnikova, Y.S.; Kudryavtsev, A.E.; Nefediev, A.V. Interplay of quark and meson degrees of freedom in a near-threshold resonance. Eur. Phys. J A 2010, 44, 93. [Google Scholar] [CrossRef] [Green Version]
  64. Ortega, P.G.; Entem, D.R.; Fernández, F. Unquenching the Quark Model in a Nonperturbative Scheme. Adv. High Energy Phys. 2019, 2019, 3465159. [Google Scholar] [CrossRef]
  65. Abulencia, A.; Acosta, D.; Adelman, J.; Affolder, T.; Akimoto, T.; Albrow, M.G.; Ambrose, D.; Amerio, S.; Amidei, D.; Anastassov, A.; et al. Measurement of the Dipion Mass Spectrum in X(3872)→J/ψπ+π Decays. Phys. Rev. Lett. 2006, 96, 102002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Abe, K. Evidence for X(3872) —> γ J / ψ and the sub-threshold decay X(3872) —> ω J / ψLepton and photon interactions at high energies. In Proceedings of the 22nd International Symposium, LP 2005, Uppsala, Sweden, 30 June–5 July 2005. [Google Scholar]
  67. Aushev, T.; Eidelman, S.; Gabyshev, N.; Shwartz, B.; Usov, Y.; Zhulanov, V.; Zyukova, O.; Drutskoy, A.; Goldenzweig, P.; Lange, J.S.; et al. Study of the BX(3872)(→D*0 D ¯ 0 )K decay. Phys. Rev. D 2010, 81, 031103. [Google Scholar] [CrossRef] [Green Version]
  68. Guo, F.K. Novel Method for Precisely Measuring the X(3872) Mass. Phys. Rev. Lett. 2019, 122, 202002. [Google Scholar] [CrossRef] [Green Version]
  69. Aaij, R.; Abellan Beteta, C.; Ackernley, T.; Adeva, B.; Adinolfi, M.; Afsharnia, H.; Aidala, C.A.; Aiola, S.; Ajaltouni, Z.; Akar, S.; et al. Study of the ψ2(3823) and χc1(3872) states in B+Jψπ+πK+ decays. JHEP 2020, 08, 123. [Google Scholar] [CrossRef]
  70. Swanson, E.S. Diagnostic decays of the X(3872). Phys. Lett. B 2004, 598, 197–202. [Google Scholar] [CrossRef] [Green Version]
  71. Gamermann, D.; Oset, E. Isospin breaking effects in the X(3872) resonance. Phys. Rev. D 2009, 80, 014003. [Google Scholar] [CrossRef] [Green Version]
  72. Gamermann, D.; Nieves, J.; Oset, E.; Arriola, E.R. Couplings in coupled channels versus wave functions: Application to the X(3872) resonance. Phys. Rev. D 2010, 81, 014029. [Google Scholar] [CrossRef] [Green Version]
  73. Ortega, P.G.; Segovia, J.; Entem, D.R.; Fernández, F. Coupled channel approach to the structure of the X(3872). Phys. Rev. D 2010, 81, 054023. [Google Scholar] [CrossRef] [Green Version]
  74. Ferretti, J.; Galatà, G.; Santopinto, E. Interpretation of the X(3872) as a charmonium state plus an extra component due to the coupling to the meson-meson continuum. Phys. Rev. C 2013, 88, 015207. [Google Scholar] [CrossRef] [Green Version]
  75. Ortega, P.G.; Entem, D.R.; Fernández, F. Molecular structures in the charmonium spectrum: TheXYZpuzzle. J. Phys. G Nucl. Part. Phys. 2013, 40, 065107. [Google Scholar] [CrossRef] [Green Version]
  76. Burns, T.J. Phenomenology of Pc(4380)+, Pc(4450)+ and related states. Eur. Phys. J. A 2015, 51, 152. [Google Scholar] [CrossRef] [Green Version]
  77. Guo, F.K.; Jing, H.J.; Meißner, U.G.; Sakai, S. Isospin breaking decays as a diagnosis of the hadronic molecular structure of the Pc(4457). Phys. Rev. D 2019, 99, 091501. [Google Scholar] [CrossRef] [Green Version]
  78. Nieves, J.; Pavón Valderrama, M. Heavy quark spin symmetry partners of the X(3872). Phys. Rev. D 2012, 86, 056004. [Google Scholar] [CrossRef] [Green Version]
  79. Hidalgo-Duque, C.; Nieves, J.; Valderrama, M.P. Light flavor and heavy quark spin symmetry in heavy meson molecules. Phys. Rev. D 2013, 87, 076006. [Google Scholar] [CrossRef] [Green Version]
  80. Baru, V.; Epelbaum, E.; Filin, A.A.; Hanhart, C.; Nefediev, A.V. Molecular partners of the X(3872) from heavy-quark spin symmetry: A fresh look. EPJ Web Conf. 2017, 137, 06002. [Google Scholar] [CrossRef] [Green Version]
  81. Guo, F.K.; Hidalgo-Duque, C.; Nieves, J.; Pavón Valderrama, M. Consequences of heavy-quark symmetries for hadronic molecules. Phys. Rev. D 2013, 88, 054007. [Google Scholar] [CrossRef] [Green Version]
  82. Entem, D.R.; Ortega, P.G.; Fernández, F. Partners of the X(3872) and heavy quark spin symmetry breaking. AIP Conf. Proc. 2016, 1735, 060006. [Google Scholar] [CrossRef] [Green Version]
  83. Cincioglu, E.; Nieves, J.; Ozpineci, A.; Yilmazer, A.U. Quarkonium Contribution to Meson Molecules. Eur. Phys. J. C 2016, 76, 576. [Google Scholar] [CrossRef] [Green Version]
  84. Ortega, P.G.; Segovia, J.; Entem, D.R.; Fernández, F. Charmonium resonances in the 3.9 GeV/c2 energy region and the X(3915)/X(3930) puzzle. Phys. Lett. B 2018, 778, 1–5. [Google Scholar] [CrossRef]
  85. Batley, J.; Lazzeroni, C.; Munday, D.J.; Slater, M.W.; Wotton, S.A.; Arcidiacono, R.; Bocquet, G.; Cabibbo, N.; Ceccucci, A.; Cundy, D.; et al. Observation of a cusp-like structure in the π0π0 invariant mass distribution from K±π±π0π0 decay and determination of the ππ scattering lengths. Phys. Lett. B 2006, 633, 173–182. [Google Scholar] [CrossRef]
  86. Budini, P.; Fonda, L. Pion-Pion Interaction from Threshold Anomalies in K+ Decay. Phys. Rev. Lett. 1961, 6, 419–421. [Google Scholar] [CrossRef]
  87. Cabibbo, N. Determination of the a0a2 Pion Scattering Length from K+π+π0π0 Decay. Phys. Rev. Lett. 2004, 93, 121801. [Google Scholar] [CrossRef] [Green Version]
  88. Cabibbo, N.; Isidori, G. Pion-pion scattering and the K→3π decay amplitudes. J. High Energy Phys. 2005, 2005, 021. [Google Scholar] [CrossRef] [Green Version]
  89. Bugg, D.V. An explanation of Belle states Z b (10610) and Z b (10650). EPL (Europhys. Lett.) 2011, 96, 11002. [Google Scholar] [CrossRef] [Green Version]
  90. Swanson, E.S. Zb and Zc exotic states as coupled channel cusps. Phys. Rev. D 2015, 91, 034009. [Google Scholar] [CrossRef] [Green Version]
  91. Dong, X.K.; Guo, F.K.; Zou, B.S. Why there are many threshold structures in hadron spectrum with heavy quarks. arXiv 2020, arXiv:2011.14517. [Google Scholar]
  92. Guo, F.K.; Hanhart, C.; Wang, Q.; Zhao, Q. Could the near-threshold XYZ states be simply kinematic effects? Phys. Rev. D 2015, 91, 051504. [Google Scholar] [CrossRef] [Green Version]
  93. Aaltonen, T. Observation of the Y(4140) structure in the J/ψϕ Mass Spectrum in B±J/ψϕK decays. arXiv 2011, arXiv:1101.6058. [Google Scholar] [CrossRef]
  94. Abazov, V.M. Inclusive Production of the X(4140) State in p p ¯ Collisions at D0. Phys. Rev. Lett. 2015, 115, 232001. [Google Scholar] [CrossRef] [PubMed]
  95. Chatrchyan, S. Observation of a peaking structure in the J/ψϕ mass spectrum from B±J/ψϕK± decays. Phys. Lett. 2014, B734, 261–281. [Google Scholar] [CrossRef]
  96. Shen, C.P.; Yuan, C.Z.; Aihara, H.; Arinstein, K.; Aushev, T.; Bakich, A.M.; Balagura, V.; Barberio, E.; Bay, A.; Belous, K.; et al. Evidence for a new resonance and search for the Y(4140) in the gamma gamma —> phi J/psi process. Phys. Rev. Lett. 2010, 104, 112004. [Google Scholar] [CrossRef] [PubMed]
  97. Lees, J.P. Study of B±,0J/ψK+KK±,0 and search for B0J/ψϕ at BABAR. Phys. Rev. 2015, D91, 012003. [Google Scholar] [CrossRef] [Green Version]
  98. Aaij, R. Observation of J/ψϕ structures consistent with exotic states from amplitude analysis of B+J/ψϕK+ decays. arXiv 2016, arXiv:1606.07895. [Google Scholar]
  99. Ortega, P.G.; Segovia, J.; Entem, D.R.; Fernández, F. Canonical description of the new LHCb resonances. Phys. Rev. D 2016, 94, 114018. [Google Scholar] [CrossRef] [Green Version]
  100. Albaladejo, M.; Guo, F.K.; Hidalgo-Duque, C.; Nieves, J. Zc(3900): What has been really seen? Phys. Lett. B 2016, 755, 337–342. [Google Scholar] [CrossRef] [Green Version]
  101. Pilloni, A.; Fernández-Ramírez, C.; Jackura, A.; Mathieu, V.; Mikhasenko, M.; Nys, J.; Szczepaniak, A. Amplitude analysis and the nature of the Zc(3900). Phys. Lett. B 2017, 772, 200–209. [Google Scholar] [CrossRef] [Green Version]
  102. Ikeda, Y.; Aoki, S.; Doi, T.; Gongyo, S.; Hatsuda, T.; Inoue, T.; Iritani, T.; Ishii, N.; Murano, K.; Sasaki, K. Fate of the Tetraquark Candidate Zc(3900) from Lattice QCD. Phys. Rev. Lett. 2016, 117, 242001. [Google Scholar] [CrossRef] [PubMed]
  103. Ortega, P.G.; Segovia, J.; Entem, D.R.; Fernández, F. The Zc structures in a coupled-channels model. Eur. Phys. J. C 2019, 79, 78. [Google Scholar] [CrossRef]
  104. Ablikim, M.; Achasov, M.N.; Albayrak, O.; Ambrose, D.J.; An, F.F.; An, Q.; Bai, J.Z.; Ferroli, R.B.; Ban, Y.; Becker, J.; et al. Observation of a Charged Charmoniumlike Structure in e+eπ+πJ/ψ at s = 4.26 GeV. Phys. Rev. Lett. 2013, 110, 252001. [Google Scholar] [CrossRef] [Green Version]
  105. Liu, Z.Q.; Shen, C.P.; Yuan, C.Z.; Adachi, I.; Aihara, H.; Asner, D.M.; Aulchenko, V.; Aushev, T.; Aziz, T.; Bakich, A.M.; et al. Study of e+eπ+πJ/ψ and Observation of a Charged Charmoniumlike State at Belle. Phys. Rev. Lett. 2013, 110, 252002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Ablikim, M.; Achasov, M.N.; Albayrak, O.; Ambrose, D.J.; An, F.F.; An, Q.; Bai, J.Z.; Baldini Ferroli, R.; Ban, Y.; Becker, J.; et al. Observation of a Charged ( D D ¯ * )± Mass Peak in e+eπ D D ¯ * at s = 4.26 GeV. Phys. Rev. Lett. 2014, 112, 022001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  107. Ablikim, M.; Achasov, M.N.; Albayrak, O.; Ambrose, D.J.; An, F.F.; An, Q.; Bai, J.Z.; Baldini Ferroli, R.; Ban, Y.; Becker, J.; et al. Observation of a Charged Charmoniumlike Structure Zc(4020) and Search for the Zc(3900) in e+eπ+πhc. Phys. Rev. Lett. 2013, 111, 242001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  108. Ablikim, M.; Achasov, M.N.; Ai, X.C.; Albayrak, O.; Albrecht, M.; Ambrose, D.J.; Amoroso, A.; Lou, X.; BESIII Collaboration. Observation of e+eπ0π0hc and a Neutral Charmoniumlike Structure Zc(4020)0. Phys. Rev. Lett. 2014, 113, 212002. [Google Scholar] [CrossRef] [Green Version]
  109. Ablikim, M.; Achasov, M.N.; Ai, X.C.; Albayrak, O.; Albrecht, M.G.; Ambrose, D.J.; BESIII Collaboration. Confirmation of a charged charmoniumlike state Zc(3885) in e+eπ±( D D ¯ * ) with double D tag. Phys. Rev. 2015, D92, 092006. [Google Scholar] [CrossRef] [Green Version]
  110. Ablikim, M.; Achasov, M.N.; Ai, X.C.; Albayrak, O.; Albrecht, M.; Ambrose, D.J.; Amoroso, A.; An, F.F.; An, Q.; Bai, J.Z.; et al. Determination of the Spin and Parity of the Zc(3900). Phys. Rev. Lett. 2017, 119, 072001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  111. Ablikim, M.; Achasov, M.N.; Albayrak, O.; Ambrose, D.J.; An, F.F.; An, Q.; Bai, J.Z.; Ferroli, R.B.; Ban, Y.; Becker, J.; et al. Observation of a charged charmoniumlike structure in e+e→( D * D ¯ * )±π at s = 4.26 GeV. Phys. Rev. Lett. 2014, 112, 132001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Diagrams of the 3 P 0 model that contributes to the decay of a meson into two mesons (a) and a baryon into a baryon and a meson (b).
Figure 1. Diagrams of the 3 P 0 model that contributes to the decay of a meson into two mesons (a) and a baryon into a baryon and a meson (b).
Symmetry 13 00279 g001
Figure 2. Probabilities of different components in the physical X ( 3872 ) state obtained in the unquenched chiral quark model as a function of the mass. The mass is varied using the 3 P 0 strength parameter. Lines are D D ¯ * , charged neutral in black and charged channel in red, χ c 1 ( 2 P ) in dashed magenta and χ c 1 ( 1 P ) in thin black.
Figure 2. Probabilities of different components in the physical X ( 3872 ) state obtained in the unquenched chiral quark model as a function of the mass. The mass is varied using the 3 P 0 strength parameter. Lines are D D ¯ * , charged neutral in black and charged channel in red, χ c 1 ( 2 P ) in dashed magenta and χ c 1 ( 1 P ) in thin black.
Symmetry 13 00279 g002
Figure 3. Diagonal matrix elements of the two meson interaction in momentum space for the D ( * ) D ¯ ( * ) sector (left panel) and B ( * ) B ¯ ( * ) sector (right panel). For the left panel, the dashed blue line gives the D * D ¯ * ( 0 + + ) matrix element, the dashed red line the D D ¯ ( 0 + + ) , the solid blue line the right-hand side of Equation (61), and the solid red line left-hand side of the same equation. For the right panel, it is the same but for B ( * ) B ¯ ( * ) channels.
Figure 3. Diagonal matrix elements of the two meson interaction in momentum space for the D ( * ) D ¯ ( * ) sector (left panel) and B ( * ) B ¯ ( * ) sector (right panel). For the left panel, the dashed blue line gives the D * D ¯ * ( 0 + + ) matrix element, the dashed red line the D D ¯ ( 0 + + ) , the solid blue line the right-hand side of Equation (61), and the solid red line left-hand side of the same equation. For the right panel, it is the same but for B ( * ) B ¯ ( * ) channels.
Symmetry 13 00279 g003
Figure 4. Diagonal matrix elements of the two meson interaction in momentum space for the D ( * ) D ¯ ( * ) sector (left panel) and B ( * ) B ¯ ( * ) sector (right panel). For the left panel, the dashed blue line gives the D * D ¯ * ( 0 + + ) matrix element, the dashed red line the D D ¯ ( 0 + + ) , the solid blue line the right-hand side of Equation (62), the solid red line the left-hand side of the same equation, and the solid green line the right-hand side of Equation (63). At exact HQSS, the three solid lines should be equal. For the right panel, it is the same but for B ( * ) B ¯ ( * ) channels.
Figure 4. Diagonal matrix elements of the two meson interaction in momentum space for the D ( * ) D ¯ ( * ) sector (left panel) and B ( * ) B ¯ ( * ) sector (right panel). For the left panel, the dashed blue line gives the D * D ¯ * ( 0 + + ) matrix element, the dashed red line the D D ¯ ( 0 + + ) , the solid blue line the right-hand side of Equation (62), the solid red line the left-hand side of the same equation, and the solid green line the right-hand side of Equation (63). At exact HQSS, the three solid lines should be equal. For the right panel, it is the same but for B ( * ) B ¯ ( * ) channels.
Symmetry 13 00279 g004
Figure 5. Charmonium spectrum in the energy region of 1 P and 2 P states. Blue boxes show the states in the Particle Data Group [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 D D ¯ * but not in D D . The states in red are naive Q Q ¯ states predicted by the model.
Figure 5. Charmonium spectrum in the energy region of 1 P and 2 P states. Blue boxes show the states in the Particle Data Group [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 D D ¯ * but not in D D . The states in red are naive Q Q ¯ states predicted by the model.
Symmetry 13 00279 g005
Figure 6. Bottomonium spectrum in the energy region of 2 P and 3 P states. Blue boxes show the states in the Particle Data Group [7]. The states in red are the pure Q Q ¯ states predicted by the model.
Figure 6. Bottomonium spectrum in the energy region of 2 P and 3 P states. Blue boxes show the states in the Particle Data Group [7]. The states in red are the pure Q Q ¯ states predicted by the model.
Symmetry 13 00279 g006
Figure 7. Line shapes for D D ¯ * (left panel) and π J / ψ (right panel) at s = 4.26 GeV [103]. Experimental data are from Refs. [109,110], respectively. The theoretical line shapes have been convoluted with the experimental resolution. The line-shape’s 68 % uncertainty is shown as a shadowed area.
Figure 7. Line shapes for D D ¯ * (left panel) and π J / ψ (right panel) at s = 4.26 GeV [103]. Experimental data are from Refs. [109,110], respectively. The theoretical line shapes have been convoluted with the experimental resolution. The line-shape’s 68 % uncertainty is shown as a shadowed area.
Symmetry 13 00279 g007
Figure 8. Line shapes for D * D ¯ * (left panel) and η c ρ (right panel) at s = 4.26 GeV [103]. Experimental data for D * D ¯ * are from Ref. [111]. The theoretical line shapes have been convoluted with the experimental resolution. The line-shape’s 68 % -uncertainty is shown as a shadowed area.
Figure 8. Line shapes for D * D ¯ * (left panel) and η c ρ (right panel) at s = 4.26 GeV [103]. Experimental data for D * D ¯ * are from Ref. [111]. The theoretical line shapes have been convoluted with the experimental resolution. The line-shape’s 68 % -uncertainty is shown as a shadowed area.
Symmetry 13 00279 g008
Table 1. First levels of parahelium ( S = 0 ) with orbital angular momentum L 2 . The experimental data are taken from NIST, Ref. [48]—results with the GEM method from this work. All energies are in eV.
Table 1. First levels of parahelium ( S = 0 ) with orbital angular momentum L 2 . The experimental data are taken from NIST, Ref. [48]—results with the GEM method from this work. All energies are in eV.
TermJNISTGEM
1 s 2 1 S 00.000000000.00
1 s 2 s 1 S 020.615775133420.61
1 s 2 p 1 P 121.218023021821.21
1 s 3 s 1 S 022.92031768222.91
1 s 3 d 1 D 223.0740751194123.07
1 s 3 p 1 P 123.087018852823.08
1 s 4 s 1 S 023.673570913323.67
1 s 4 d 1 D 223.7363353578623.73
1 s 4 p 1 P 123.742070391823.74
1 s 5 s 1 S 024.011215312924.00
1 s 5 d 1 D 224.04280373493024.04
1 s 5 p 1 P 124.045800729724.04
1 s 6 s 1 S 024.191160598224.18
1 s 6 d 1 D 224.20925011689324.20
Table 2. First levels of orthohelium ( S = 1 ) with orbital angular momentum L 2 . The experimental data are taken from NIST, Ref. [48]—results with the GEM method from this work. All energies are in eV.
Table 2. First levels of orthohelium ( S = 1 ) with orbital angular momentum L 2 . The experimental data are taken from NIST, Ref. [48]—results with the GEM method from this work. All energies are in eV.
TermJNISTGEM
1 s 2 s 3 S 119.8196148420319.82
1 s 2 p 3 P 220.9640872067520.98
120.96409668230
020.96421916817
1 s 3 s 3 S 122.71846674222.71
1 s 3 p 3 P 223.007073467323.01
123.0070761918
023.0071097475
1 s 3 d 3 D 323.0736510299023.07
223.07365134140
123.07365682165
1 s 4 s 3 S 123.59395903623.59
1 s 4 p 3 P 223.707891551123.71
123.7078926664
023.7079063452
1 s 4 d 3 D 323.7360905124723.73
223.73609066143
123.73609295768
1 s 5 s 3 S 123.971971741323.97
1 s 5 p 3 P 224.028225387024.02
124.0282259477
024.0282328220
1 s 5 d 3 D 324.04266256481924.04
224.042662644310
124.042663817021
1 s 6 s 3 S 124.168998546324.16
1 s 6 p 3 P 224.200815777624.20
124.2008160981
024.2008200312
1 s 6 d 3 D 324.20916343333524.22
224.209163480258
124.209164158016
Table 3. Results for the Bhaduri potential obtained with a Fadeev calculation (FD) [50] and the GEM method (this work). Masses are given in MeV, while matter radius square R m 2 and charge radius square R c 2 are in fm 2 .
Table 3. Results for the Bhaduri potential obtained with a Fadeev calculation (FD) [50] and the GEM method (this work). Masses are given in MeV, while matter radius square R m 2 and charge radius square R c 2 are in fm 2 .
StateM(FD)M(GEM) R m 2 (FD) R m 2 (GEM) R c 2 (FD) R c 2 (GEM)
Λ c + ( c u d ) 23002298.50.0970.09840.1170.1180
Σ c 0 ( c u d ) 24732475.00.1110.1116−0.224−0.2247
Σ c + ( c u d ) 0.1340.1347
Σ c + + ( c u d ) 0.4940.4941
Λ b ( b u d ) 56535649.60.0430.04350.1150.1169
Σ b ( b u d ) 58585859.80.0510.0509−0.280−0.2804
Σ b 0 ( b u d ) 0.1380.1383
Σ b + ( b u d ) 0.5550.5571
Ξ c 0 ( n s c ) 24902490.90.0970.0978−0.145−0.1463
Ξ c + ( n s c ) 0.1600.1617
Ω c 0 ( c s s ) 27002701.00.1000.0999−0.111−0.1111
Ξ b ( n s b ) 58265824.80.0450.0459−0.193−0.1951
Ξ b 0 ( n s b ) 0.1510.1517
Ω b ( b s s ) 60466046.70.0500.0505−0.164−0.1642
Ξ c c + ( n c c ) 36313632.20.0760.0766−0.034−0.0330
Ξ c c + + ( n c c ) 0.2850.2852
Ξ b b ( n b b ) 10,19710,197.40.0310.0309−0.128−0.1279
Ξ b b 0 ( n b b ) 0.2150.2140
Ω c c + ( s c c ) 37393738.70.0730.07390.0080.0091
Ω c b + ( s c b ) 70237024.20.0430.0430−0.023−0.0232
Ω b b ( s b b ) 10,27110,271.30.0300.0304−0.083−0.0829
Ω c c c + + ( c c c ) 48064807.20.0620.06190.1240.1239
Ω c c b + ( c c b ) 80328030.90.0380.03780.0890.0891
Ω c b b 0 ( c b b ) 11,22011,218.60.0260.02640.0320.0318
Ω b b b ( b b b ) 14,37014,371.80.0190.0192−0.019−0.0192
Table 4. Levels of orthohelium ( S = 1 ) above the first open threshold. The experimental data are taken from NIST, Refs. [52,53]—results with the GEM method from this work. All energies are in eV.
Table 4. Levels of orthohelium ( S = 1 ) above the first open threshold. The experimental data are taken from NIST, Refs. [52,53]—results with the GEM method from this work. All energies are in eV.
TermNISTGEM
2 s 2 p 3 P 58.311 [52]58.31
2 p 2 3 P 59.67378 [53]59.66
2 p 3 p 3 D 63.120 [52]
2 p 3 d 3 D 63.78658 [52]
2 p 3 d 3 P 64.0719 [52]
Table 5. The S-matrix pole positions, in MeV / c 2 , for different coupled-channels calculations [103]. 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.
Table 5. The S-matrix pole positions, in MeV / c 2 , for different coupled-channels calculations [103]. 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 Z c ( 3900 ) PoleRS Z c ( 4020 ) PoleRS
D D ¯ * 3871.37 2.17 i (S)--
D D ¯ * + D * D ¯ * 3872.27 1.85 i (S,F) 4014.16 0.10 i (S,S)
ρ η c + D D ¯ * 3871.32 0.00 i (S,S)--
ρ η c + D D ¯ * + D * D ¯ * 3872.07 0.00 i (S,S,F) 4013.10 0.00 i (S,S,S)
π J / ψ + ρ η c + D D ¯ * + D * D ¯ * 3871.74 0.00 i (S,S,S,F) 4013.21 0.00 i (S,S,S,S)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ortega, P.G.; Entem, D.R. Coupling Hadron-Hadron Thresholds within a Chiral Quark Model Approach. Symmetry 2021, 13, 279. https://doi.org/10.3390/sym13020279

AMA Style

Ortega PG, Entem DR. Coupling Hadron-Hadron Thresholds within a Chiral Quark Model Approach. Symmetry. 2021; 13(2):279. https://doi.org/10.3390/sym13020279

Chicago/Turabian Style

Ortega, Pablo G., and David R. Entem. 2021. "Coupling Hadron-Hadron Thresholds within a Chiral Quark Model Approach" Symmetry 13, no. 2: 279. https://doi.org/10.3390/sym13020279

APA Style

Ortega, P. G., & Entem, D. R. (2021). Coupling Hadron-Hadron Thresholds within a Chiral Quark Model Approach. Symmetry, 13(2), 279. https://doi.org/10.3390/sym13020279

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop