1. Introduction
Over the last two decades, the experimental observation of many new hadronic states is challenging our current understanding of hadrons as conventional mesons and baryons with valence contents of quark–antiquark and three quarks, respectively, since most of them do not fit in the well-known quark model. This difficulty brought back a long-standing discussion on the exotic hadronic structures, i.e., multiquark configurations that might have quantum numbers beyond those assigned to the conventional mesons and baryons [
1,
2].
Exotic quark configurations, such as tetraquarks [
3,
4], hadron–hadron molecules [
5], glueballs, and hybrids [
6,
7], among others, have been suggested to describe suitably most of the properties of these new states, such as the
quantum numbers, mass, and decay width, especially for those lying in the charmonium and bottomonium spectra.
On the other hand, distinguishing the exotic states from the conventional hadrons is a more complicated task in the light quark sector. Many states have their masses close to each other, and the possibility of mixing brings additional difficulty to the problem. The situation improves as the quantum numbers do not fall into those allowed by the conventional quark model. It seems to be the case of the newly discovered state, dubbed
, by the BESIII Collaboration [
8,
9], observed in the invariant mass distribution of the
meson pair in the
decay channel with a significance of
. Its mass and width reported by BESIII are
MeV and
MeV, respectively, with likely
quantum numbers, which cannot be formed by a pair of quark and antiquark. The
is not the only state experimentally found with that set of quantum numbers. As of today, three other hadronic structures, called
,
, and
, with
, were observed by several collaborations [
7,
10].
From the theoretical point of view, the hybrid model has been used to investigate these exotic meson states, in particular the
ones. Lattice quantum chromodynamics (QCD) calculations have pointed out hybrid supermultiplets with exotic
quantum numbers, including the
one [
11,
12,
13,
14,
15,
16]. In this picture, however, the mass of the lightest
state and decay modes are inconsistent with the corresponding experimental results, while the
and
structures can fit into the nonets predicted by lattice QCD [
7].
The newly observed
state has also been the focus of some studies. In particular, the authors in Ref. [
17] proposed two hybrid nonet schemes in which the
resonance can be either the lower or higher mass state with isospin
. In Ref. [
18], an effective Lagrangian respecting flavor, parity, and charge conjugation symmetries is used to study the hybrid nonet decays into two-body meson states. The authors have fixed the couplings to those two-body meson states by performing a combined fit to the experimental and lattice results available. As a result, the decay width value estimated for the isoscalar member of the hybrid nonet agrees with the one observed for
state. Additionally, addressing the same picture, Ref. [
19] applied the approach of QCD sum rules to describe the
mass. By contrast, within the same approach, the
resonance is described as a tetraquark state in Ref. [
20].
The
resonance also supports a meson–meson molecule interpretation due to its proximity to the
threshold, as put forward by Refs. [
21,
22]. In particular, the authors in Ref. [
21] have investigated the
interaction through the one-boson exchange model. According to their findings, the
system binds for cut-off values above 2 GeV with a monopole form factor. In addition, the comparison between their result for the branching fraction
to the experimental one led them to conclude that the
molecule can explain the
structure.
An important point to be addressed is the meson–meson interaction around the threshold for the quantum numbers. In this sector, many meson–meson pairs may contribute to that interaction, so a coupled-channel treatment seems appropriate to take these contributions into account. In particular, hadron–hadron interactions in coupled channels have been studied in many works to describe the properties of the new hadronic systems experimentally observed. In those cases, these hadronic structures are called dynamically generated states.
Following this approach, in this work, we aim to explore the
,
, and
hadronic systems as dynamically generated states from pseudoscalar-axial vector meson interactions in coupled channels. Specifically, the low-energy interactions are given by the Weinberg–Tomozawa (WT) term from chiral Lagrangians at the leading order of the chiral expansion by treating the axial vector mesons as matter fields and the pseudoscalar mesons as the pseudo-Nambu–Goldstone bosons of the spontaneous breaking of chiral symmetry. Such Lagrangians have been used to study many hadron structures stemming from meson–meson and meson–baryon interactions in coupled channels in light and heavy sectors, see, e.g., Refs. [
23,
24,
25,
26,
27]. In our case, the amplitudes obtained from the WT term are unitarized via the Bethe–Salpeter equation from which bound states/resonances manifest as poles in the physical/unphysical Riemann sheets of the scattering matrices. The existence of a whole family of kaonic bound states has been pointed out in Ref. [
28] based on unitarizing the WT term for the scattering of the kaon off isospin-1/2 matter fields taking heavy mesons and doubly-charmed baryons as examples. As we shall show in this work, the newly observed
structure may correspond to a dynamically generated state from the pseudoscalar-axial vector interaction in the isospin
sector coupling strongly to the
channel. Moreover, the
and
, may be assigned as the
SU(3) partners which are also dynamically generated from the pseudoscalar-axial vector meson interactions in the
sector. The former resonance couples mainly to the
channel, and the latter has the
as its main coupled channel.
In addition, we have also found two poles around 1.7 GeV in the
sector. These poles are particularly interesting as they could be the origin of the
structure observed experimentally [
10], which is the main component of the
contribution to the
mass distribution in the
decays recently measured by LHCb [
29].
This paper is organized as follows. In
Section 2, we discuss the relevant channels contributing to the pseudoscalar-axial vector meson interactions and the use of the chiral unitary approach (ChUA) for the evaluation of the transition amplitudes among those channels. In
Section 3 and
Section 4, we investigate the dynamical generation of poles stemming from those interactions in the
and
sectors and discuss their possible decay channels. Finally, in
Section 5, we also explore the dynamical generation of poles for
and their connection to the vector
structure observed experimentally.
Section 6 gives a summary.
4. The Dynamical Generation
The WT amplitudes for the pseudoscalar-axial vector meson interactions with
are given by Equation (
16), with the corresponding
coefficients listed in
Table 7. In this case, from Equation (
19), we obtain two
poles shown in
Table 11.
Similar to the previous section, we also provide the couplings of these dynamically generated states to the channels listed in
Table 2.
Table 11 shows a broad
pole at
GeV, and a width of about
GeV.
1 This state is above the
and
thresholds. Its large width stems from the large coupling to the
and the fact that this channel is open for decaying. The
channel is also open. However, according to
Table 7, the corresponding WT term in Equation (
16) is zero for the diagonal
transition. On the other hand, the next
pole in
Table 11 has a sizeable dependence on the mixing angles. Using set
A, we find that pole at
GeV. It couples most strongly to the
channel, which is closed for decaying. Nonetheless, the state can decay into
and
, albeit their corresponding couplings are small compared to the
one, but still large enough to provide a sizeable width for the pole. In contrast, when set
B is adopted, the higher
pole is now located at
GeV, above the
threshold, which is now open. One might think that the width should increase since now three channels are open for decaying. However, although the coupling to the
has increased in this case, at the same time the couplings to the other open channels have decreased. Hence, the overall effect leads to a smaller width compared to the previous case.
The lower pole mass is slightly higher than the mass of the
state listed in RPP,
MeV [
10]. Notice that we use the same subtraction constant for all channels. In principle, it can take different values and lead to a shift of the poles. In addition, we did not include in the loops the
width, that is relatively large and whose effects could influence the pole position. However, it is expected to affect more the imaginary part of the pole than the real one (see
Figure 5a below). We can obtain a rough estimate of this change by adding the
width to the previous result for
, with
the lower
pole, i.e.,
which is close to the
width reported in RPP,
MeV [
10]. From these results, we are led to claim that the lower
pole may explain the
resonance; in other words, the
is suitably described in our approach as a dynamically generated state with the
as its main component.
Alternatively, following the prescription used in
Section 3, we can also study the changes in the results caused by the inclusion of the finite widths for the axial-vector mesons by looking at the line shape for the relevant
T-matrix elements. In
Figure 5a, we show the line shapes for the
T-matrix element corresponding to the elastic
transition, which is the one we would expect the lower pole in
Table 11 manifests most due to its large coupling to the
channel. It becomes clear that the bumps become broader when the widths of axial-vector mesons are taken into account. A similar behavior can be seen in
Figure 5b for the
T-matrix element associated with the scattering of
, which is the channel to which the higher
pole couples most strongly. The peak mass and width extracted from the line shape of the diagonal
T-matrix element for the dominantly coupled channel are also listed in
Table 11 when the axial-vector meson widths are considered.
The higher
pole, denoted now by
, has a mass consistent with that of the
, whose pole mass has been reported to be
MeV in Ref. [
54] and
MeV in Ref. [
55]. It can decay into the
and
channels. The corresponding diagrams for both amplitudes are illustrated in
Figure 6, from which we have
and
with
and
the polarization vectors of the
and
mesons. Here,
,
, and
are the effective coupling of the
pole to the corresponding couplings, and
and
are the loops involving the
and
mesons, respectively. Notice that the effective couplings are computed from the residues of the
T matrix elements; thus they contain contributions from all coupled channels.
In order to compare our findings with the experimental information, we evaluate the ratio
where
q and
are the momentum in the c.m. frame of the
and
pairs, respectively. Numerically, Equation (
44) gives
The ratio is slightly bigger for the mixing angles in the set
A. Nevertheless, the result in Equation (
45) is consistent to the corresponding ratio
reported by the E852 Collaboration [
56]. This good agreement with the experimental data supports the molecular picture for the
state.
6. Conclusions
We have studied the interactions between the pseudoscalar and axial-vector mesons in coupled channels with quantum numbers for the isospin 0, 1, and sectors. Using the chiral unitary approach, we describe the interaction with the Weinberg–Tomozawa term derived from chiral Lagrangians. The transition amplitudes among all the relevant channels are unitarized using the Bethe–Salpeter equation from which resonances (bound states) manifest themselves as poles on the (un)physical Riemann sheets of the complex energy plane.
We consider the physical isoscalar axial-vector states as mixtures of the corresponding SU(3) singlets and octets. In addition, the and physical states are also mixtures of the and mesons, which are the strange partners of the and resonances, respectively. We group into two sets, called A and B, the mixing angles accounting for such mechanisms and investigate their influence on the pole positions.
According to our findings, we obtain poles with
quantum numbers in the energy range from
to
GeV, in each isospin sector studied (
). The
quantum numbers are exotic in the sense that they cannot be formed from a pair of quark and antiquark. In particular, we have found an isoscalar state that may correspond to the
state, newly observed by the BESIII Collaboration [
8]. In addition, we have also found two dynamically generated isovector states that we assign to be the
and
resonances. Hence, within our formalism, they are dynamically generated through the pseudoscalar-axial vector meson interactions, with the
state coupling mostly to
channel, while the
couples strongly to the
, and
structure couples most strongly to the
. We also find two
states with a mass around 1.7 GeV. They combined together could be responsible to the observed
structure.
In addition, we also evaluate the decays of the
and the
. We find that the three-body decay channel
has a significantly larger branching fraction than the
, which is the channel where the observation of the
was made. The obtained ratio between the
and
decays, given by Equation (
45), is consistent with the corresponding experimental value.
We suggest searching for two additional exotic mesons with masses of about 1.4 and 1.7 GeV, respectively. In particular, the latter should be relatively narrow with a width around 0.1 GeV and one of its main decay channels is .