1. Introduction
The scalar meson
was established experimentally 40 years ago; however, so far, its nature has been a long-standing puzzle. This scalar meson dominantly decays to
final states, but also decays to
final states with a small branching fraction [
1]. Resonant parameters of
were obtained via measuring decay channels, such as
by SND [
2] and KLOE [
3], and
by BESII [
4], or via analyzing scattering processes, such as
by BABAR [
5],
scattering data and
decay data [
6,
7]. Based on those experimental measurements and data analysis, the mass and width of
estimated by the Particle Data Group (PDG) were
MeV/
and
to 100 MeV until 2016 [
8]. It should be noted that not all uncertainties come from experimental data. The shape of the
varies for different decay processes, and that is why the mass and width are quoted from their process-independent pole position at an unphysical sheet of the complex energy plane, defined as
.
In 2012, the BESIII Collaboration first observed anomalously narrow widths of about 10 MeV of
via isospin-violating decays of
[
9]. The isospin-violating ratio between
and
is up to
. Three years later, the BESIII Collaboration once again observed the similarly narrow widths of
via isospin-violating decays of
, in which the isospin-violating ratio is measured to be
[
10]. Since then the width of
in PDG has been updated to be
10 to 100 MeV [
1]. In 2018, the anomalously narrow width of
was further confirmed by the BESIII Collaboration via an isospin-violating decay of
[
11]. It is worth noting that the narrow-width
mesons are only produced in isospin-violating decays, which apparently disagrees with the normal width of
in isospin-conserving decays.
Theoretically, the internal structure of the
is not only considered as the conventional quark-antiquark, but also proposed to be tetraquarks [
12],
molecule [
13], or quark-antiquark gluon hybrid [
14]. However, explanations about the nature of
have been controversial to date. The most famous theoretical study about
is the
-
mixing mechanism, which was first proposed in the late 1970s [
15]. Because
and
both can decay into
, the charged and neutral kaon mass thresholds are different by about 8 MeV due to isospin-violating effects. The mixing amplitude between
and
is dominated by the unitary cuts of the intermediate two-kaon system and proportional to the phase-space difference between charged and neutral kaon systems. As a consequence, a narrow peak of about 8 MeV in width is predicted between the charged and neutral
mass thresholds. The Feynman diagram of
mixing in the decay of
is shown in
Figure 1(a), where
X can be
,
or
. The
-
mixing mechanism has been investigated extensively for a long time, and many decay processes have been discussed [
16,
17,
18,
19,
20,
21,
22,
23,
24,
25,
26,
27]. There was no experimental results until the BESIII Collaboration reported
-
mixing via the decays of
and
[
11,
28]. The mixing intensity of
,
i.e. the isospin-violating ratio, is measure to be
, which is less than
. Obviously, the
-
mixing mechanism can not completely describe the large isospin-violating ratio in the decays of
and
.
Since the anomalously narrow widths of about 10 MeV of
was reported by the BESIII experiment [
9], a novel scenario called triangle singularity mechanism was proposed to explain the internal behavior in isospin-violating decay processes [
29]. Taking the decay
as an example, the intermediate
c.c. pair can exchange an on-shell kaon, the three interaction vertices satisfy the energy-momentum conservation, the physical amplitude has a logarithmic triangle singularity, and the kinematic effects result in a narrow peak in the
invariant mass distributions. The Feynman diagram of the triangle singularity mechanism is shown in
Figure 1(b). The triangle singularity mechanism can well explain the narrow width of the
and the large isospin-violating ratio in the decays of
. The triangle singularity mechanism was continuously discussed more extensively and deeply [
30,
31,
32], and applied to other decay processes [
33,
34,
35,
36,
37].
In this article, firstly, we use the non-relativistic Breit-Wigner formula to simultaneously fit five
invariant mass distributions reported by the BESIII Collaboration, thereby obtaining the accurate mass and width of the
for isospin-violating decays in
Section 2. Secondly, the energy-dependent Flatté formula as the shape of the
is used to simultaneously fit the
invariant mass distributions to determine the coupling constants of
and
,
i.e.,
and
in
Section 3. Thirdly, we obtain the joint confidence region of the two coupling constants between
and
as a way to provide quantitative constraints on different theoretical models of the
meson in
Section 4. Finally, we summarize the whole article in
Section 5.
2. Simultaneous fit to invariant mass distributions in isospin-violating decays
So far, the BESIII Collaboration has reported anomalously narrow width of the
meson in five independent isospin-violating decay processes, they are
[
9],
[
10], and
[
11], respectively. In order to accurately determine the mass and width of the
in isospin-violating decays, we perform a simultaneous fit to the
invariant mass spectra of the above five decay channels. In the simultaneous fit, the signal shape of the narrow
meson is described by the same non-relativistic Breit-Wigner function convolved with a Gaussian mass resolution of each decay channel. The mass resolution has an important influence on determining the width of the
meson, therefore the mass resolution is obtained in advance using Monte Carlo simulations, listed in
Table 1. The backgrounds are represented by first-order or second-order Chebyshev polynomials, where the polynomial order is kept consistent with the fit in original papers published by BESIII.
The results of the simultaneous fit are illustrated in
Figure 2.
Figure 2 (a) and (b) are the fit results of the
invariant mass distributions in the decays of
,
and
;
Figure 2 (c) and (d) are the fit results of the
invariant mass spectra in the decays of
,
and
;
Figure 2 (e) is the fit result of the
invariant mass distribution in the decay of
,
. The mass and width of the
meson obtained from the simultaneous fit, as well as the results reported by the BESIIII Collaboration, are listed in
Table 1. The mass and width from the simultaneous fit are
MeV/
and
MeV, respectively. The fit results are consistent with the mass and width reported by the BESIIII Collaboration for each channel, but the errors of the mass and width are improved remarkably. We can also obtain the joint confidence regions of the mass and width from the simultaneous fit, as shown in
Figure 2 (e), in which the circles from inside to outside represent the confidence level of the two parameters from 1 to 6 standard deviations in order.
3. Determination of coupling constants of and
In general, the
meson can be descibed by the Flatté form of the propagator [
29]:
where,
s is the square of
invariant mass,
is the square of the nominal mass of
in PDG [
1],
and
are energy-dependent partial widths of
and
, respectively. They are defined as:
where,
and
are coupling constants of
and
, and
is the momentum of the particle A or B in the center-of-mass frame of the two-body decay.
In order to obtain the coupling constants of
and
, a similar simultaneous fit to the
invariant mass spectra of the five decay channels is performed. In the new fit, only the shape of the
meson is replaced by the Flatté formula from the original Breit-Wigner function, and the remaining parts keep the same as the previous fit. The two coupling constants of
and
from the simultaneous fit are determined to be
and
, respectively, which are listed in
Table 2.
4. Joint confidence regions of the coupling constants and
Theoretically, the
is not only considered as the conventional quark-antiquark (
) meson, but also the internal structures are proposed to be tetraquarks (
) state [
12],
molecule [
13], or quark-antiquark gluon (
) hybrid [
14]. For those different theoretical models, the correspending coupling constants,
and
, in the decays of
and
have been predicted [
27], which are also listed in
Table 2 for comparison.
In order to distinguish the different theoretical models using the experimental data, we also draw the joint confidence regions of these two coupling constants from the simultaneous fit, shown in
Figure 3. The dashed lines from inside to outside represent the confidence level of the two coupling constants from 1 to 5 standard deviations in order. The predicted positions from different theoretical models are also marked on the same plot, we find that only the prediction of the
molecule model is located within the region of 5 standard deviations, while the predictions from other theoretical models are all outside the region of the five standard deviations, and the position of the quark-antiquark (
) model is very close to the boundary of the five standard deviations. The predicted positions of the tetraquarks (
) and the quark-antiquark gluon (
) hybrid are far from the region of the five standard deviations. Therefore, the experimental data tend to support the
molecule model and the quark-antiquark (
) model, but tend not to support the tetraquarks (
) model and the quark-antiquark gluon (
) hybrid model.